state_choropleth不适用于>1个州。

huangapple
go评论96阅读模式
英文:

state_choropleth not working for >1 state

问题

尝试一个非常基本的州级分布地图：

  1. library(choroplethr)
  2. library(choroplethrMaps)
  3. state_choropleth(df_pop_state)

出现以下错误：

  1. Error in private$zoom == "alaska" || private$zoom == "hawaii" : 
  2.   'length = 51' 在强制转换为 'logical(1)' 时出错

这个错误与这里找到的错误相同，尽管在那种情况下，地图似乎已经生成。

如果我缩放到一个州，它可以正常工作：

  1. state_choropleth(df_pop_state, zoom=c("washington"))

不知道出了什么问题。我相当确定一切都是最新的。几个月前我进行了R升级（也许是从版本3到4？），所以也许与此有关。

sessionInfo()的结果如下：

  1. R 版本 4.3.1 (2023-06-16)
  2. 平台：aarch64-apple-darwin20 (64 位)
  3. 运行在：macOS Monterey 12.4
  5. Matrix products: default
  6. BLAS：   /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib 
  7. LAPACK： /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/lib/libRlapack.dylib;  LAPACK 版本 3.11.0
  9. locale:
  10. [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
  12. 时区：Europe/Rome
  13. tzcode 源：internal
  15. attached base packages:
  16. [1] datasets  utils     stats     graphics  grDevices methods   base     
  18. other attached packages:
  19.  [1] dplyr_1.1.2           zoo_1.8-12            viridis_0.6.3         viridisLite_0.4.2    
  20.  [5] scales_1.2.1          readr_2.1.4           RColorBrewer_1.1-3    plyr_1.8.8           
  21.  [9] plotrix_3.8-2         PBSmapping_2.73.2     maptools_1.1-7        sp_2.0-0             
  23. ...

以上是您提供的代码和信息的翻译。

英文:

Trying a very basic state choropleth:

<!-- language: R -->
library(choroplethr)
library(choroplethrMaps)
state_choropleth(df_pop_state)

It exits with this error:

<!-- language: R -->

  1. Error in private$zoom == &quot;alaska&quot; || private$zoom == &quot;hawaii&quot; : 
  2.   &#39;length = 51&#39; in coercion to &#39;logical(1)&#39;

which is the same as the one found here, though in that case, the map seems to get produced.

If I zoom to one state, it works fine:
<!-- language: R -->
state_choropleth(df_pop_state, zoom=c("washington"))

No idea what's going on. I'm pretty sure everything is up to date. I did have an R upgrade (maybe from version 3 to 4?) a few months ago, so maybe that's related.

TIA.

Results of sessionInfo()

  1. R version 4.3.1 (2023-06-16)
  2. Platform: aarch64-apple-darwin20 (64-bit)
  3. Running under: macOS Monterey 12.4
  5. Matrix products: default
  6. BLAS:   /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib 
  7. LAPACK: /Library/Frameworks/R.framework/Versions/4.3-arm64/Resources/lib/libRlapack.dylib;  LAPACK version 3.11.0
  9. locale:
  10. [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
  12. time zone: Europe/Rome
  13. tzcode source: internal
  15. attached base packages:
  16. [1] datasets  utils     stats     graphics  grDevices methods   base     
  18. other attached packages:
  19.  [1] dplyr_1.1.2           zoo_1.8-12            viridis_0.6.3         viridisLite_0.4.2    
  20.  [5] scales_1.2.1          readr_2.1.4           RColorBrewer_1.1-3    plyr_1.8.8           
  21.  [9] plotrix_3.8-2         PBSmapping_2.73.2     maptools_1.1-7        sp_2.0-0             
  22. [13] lubridate_1.9.2       jsonlite_1.8.5        ggmap_3.0.2           ggplot2_3.4.2        
  23. [17] devtools_2.4.5        usethis_2.2.1         deldir_1.0-9          clipr_0.8.0          
  24. [21] choroplethrZip_1.5.0  choroplethrMaps_1.0.1 choroplethr_3.7.1     acs_2.1.4            
  25. [25] XML_3.99-0.14         stringr_1.5.0         car_3.1-2             carData_3.0-5        
  26. [29] calibrate_1.7.7       MASS_7.3-60           alphahull_2.5         datapasta_3.1.1      
  28. loaded via a namespace (and not attached):
  29.   [1] rstudioapi_0.14       magrittr_2.0.3        spatstat.utils_3.0-3  rmarkdown_2.22       
  30.   [5] fs_1.6.2              vctrs_0.6.3           memoise_2.0.1         base64enc_0.1-3      
  31.   [9] htmltools_0.5.5       Formula_1.2-5         KernSmooth_2.23-21    htmlwidgets_1.6.2    
  32.  [13] cachem_1.0.8          uuid_1.1-0            mime_0.12             lifecycle_1.0.3      
  33.  [17] pkgconfig_2.0.3       Matrix_1.5-4.1        R6_2.5.1              fastmap_1.1.1        
  34.  [21] tidycensus_1.4.1      shiny_1.7.4           digest_0.6.31         colorspace_2.1-0     
  35.  [25] ps_1.7.5              pkgload_1.3.2         Hmisc_5.1-0           fansi_1.0.4          
  36.  [29] timechange_0.2.0      httr_1.4.6            polyclip_1.10-4       abind_1.4-5          
  37.  [33] compiler_4.3.1        proxy_0.4-27          splancs_2.01-43       remotes_2.4.2        
  38.  [37] withr_2.5.0           htmlTable_2.4.1       backports_1.4.1       DBI_1.1.3            
  39.  [41] RgoogleMaps_1.4.5.3   pkgbuild_1.4.1        R.utils_2.12.2        rappdirs_0.3.3       
  40.  [45] sessioninfo_1.2.2     classInt_0.4-9        tools_4.3.1           units_0.8-2          
  41.  [49] foreign_0.8-84        httpuv_1.6.11         nnet_7.3-19           R.oo_1.25.0          
  42.  [53] glue_1.6.2            callr_3.7.3           promises_1.2.0.1      grid_4.3.1           
  43.  [57] sf_1.0-13             checkmate_2.2.0       cluster_2.1.4         generics_0.1.3       
  44.  [61] gtable_0.3.3          spatstat.data_3.0-1   tzdb_0.4.0            R.methodsS3_1.8.2    
  45.  [65] class_7.3-22          tidyr_1.3.0           data.table_1.14.8     hms_1.1.3            
  46.  [69] xml2_1.3.4            utf8_1.2.3            spatstat.geom_3.2-1   pillar_1.9.0         
  47.  [73] later_1.3.1           lattice_0.21-8        tidyselect_1.2.0      miniUI_0.1.1.1       
  48.  [77] knitr_1.43            gridExtra_2.3         xfun_0.39             stringi_1.7.12       
  49.  [81] evaluate_0.21         interp_1.1-4          sgeostat_1.0-27       tibble_3.2.1         
  50.  [85] cli_3.6.1             rpart_4.1.19          xtable_1.8-4          munsell_0.5.0        
  51.  [89] processx_3.8.1        Rcpp_1.0.10           spatstat.random_3.1-5 WDI_2.7.8            
  52.  [93] png_0.1-8             ellipsis_0.3.2        prettyunits_1.1.1     jpeg_0.1-10          
  53.  [97] profvis_0.3.8         urlchecker_1.0.1      bitops_1.0-7          tigris_2.0.3         
  54. [101] e1071_1.7-13          purrr_1.0.1           crayon_1.5.2          rlang_1.1.1          
  55. [105] rvest_1.0.3

答案1

得分: 1

更新：此错误是由 R 语言本身的更改引入的（版本为 4.3.0）。我已经更新了 choroplethr 以解决这个问题。已经修复此问题的 choroplethr 版本是 v3.7.2，您可以在控制台输入 install.packages("choroplethr") 来获取它。我将保留我的原始答案，以便其他人想更好地理解底层问题。

我是 choroplethr 的作者。感谢您报告此错误。我可以在本地重现它，使用最新版本的 R（v4.3.1）。我已经追踪到这个问题是由于 R 版本 4.3.0 中引入的这个更改引起的：

> R 4.3.0 中的更改：
>
> 重要的用户可见更改：
>
> * 使用左操作数（LHS）或（如果已评估）右操作数（RHS）的长度大于一现在总是会报错，报错信息的形式如下：
>
> 'length = 4' in coercion to 'logical(1)'
>
> 环境变量 R_CHECK_LENGTH_1_LOGIC2 不再起作用。
>
（您可以在此处看到整个更改列表。）

这个更改导致这个代码行失败：

  1.       if (private$zoom == "alaska" || private$zoom == "hawaii") {
  2.         choro = self$render_helper(self$choropleth.df, self$scale_name, self$theme_clean()) + ggtitle(self$title)
  3.         if (self$add_state_outline)
  4.         {
  5.           choro + self$render_state_outline(private$zoom)
  6.         }        
  7.       } else {
  8.         # remove AK and HI from the "real" df
  9.         continental.df = self$choropleth.df[!self$choropleth.df$state %in% c("alaska", "hawaii"), ]
  10.         continental.ggplot = self$render_helper(continental.df, self$scale_name, self$theme_clean()) + ggtitle(self$title)
  11.         if (self$add_state_outline)
  12.         {
  13.           continental.regions = subset(private$zoom, private$zoom!="alaska" & private$zoom!="hawaii")
  14.           continental.ggplot = continental.ggplot + self$render_state_outline(continental.regions)
  15.         }

简而言之，当您创建国家地图时，zoom 参数的长度为 51（50 个州加华盛顿特区），这导致了错误。

我的第一反应是我不知道他们为什么要引入这个更改。这不是一个难以解决的错误。但在提交修复之前，我希望花一些时间更好地理解他们为什么要做出这个更改，以及他们是否有一种首选的实现此逻辑的方式。

英文:

Update: This bug was introduced by a change in the R language itself (v4.3.0). I have updated choroplethr to work around the issue. The version of choroplethr that has the fix is v3.7.2 and you can get it by typing install.packages(&quot;choroplethr&quot;) at the console. I am leaving my original answer here in case others want to better understand the underlying issue.

I am the author of choroplethr. Thank you for reporting this bug. I can reproduce it locally, on the latest version of R (v4.3.1). I have tracked the issue down to this change which was introduce in R version 4.3.0:

> CHANGES IN R 4.3.0:
>
> SIGNIFICANT USER-VISIBLE CHANGES:
>
> * Calling && or || with LHS or (if evaluated) RHS of length greater
> than one is now always an error, with a report of the form
>
> 'length = 4' in coercion to 'logical(1)'
>
> Environment variable R_CHECK_LENGTH_1_LOGIC2 no longer has any
> effect.
>
(You can see the entire changelist here.)

That change causes this line of code to fail:

  1.       if (private$zoom == &quot;alaska&quot; || private$zoom == &quot;hawaii&quot;) {
  2.         choro = self$render_helper(self$choropleth.df, self$scale_name, self$theme_clean()) + ggtitle(self$title)
  3.         if (self$add_state_outline)
  4.         {
  5.           choro + self$render_state_outline(private$zoom)
  6.         }        
  7.       } else {
  8.         # remove AK and HI from the &quot;real&quot; df
  9.         continental.df = self$choropleth.df[!self$choropleth.df$state %in% c(&quot;alaska&quot;, &quot;hawaii&quot;), ]
  10.         continental.ggplot = self$render_helper(continental.df, self$scale_name, self$theme_clean()) + ggtitle(self$title)
  11.         if (self$add_state_outline)
  12.         {
  13.           continental.regions = subset(private$zoom, private$zoom!=&quot;alaska&quot; &amp; private$zoom!=&quot;hawaii&quot;)
  14.           continental.ggplot = continental.ggplot + self$render_state_outline(continental.regions)
  15.         }

In short, when you create a national map, the zoom parameter has length 51 (50 states plus DC), which is causing the error.

My first reaction is that I have no idea why they introduced this change. It is not a hard error to work around. But I would like to spend some time better understanding why they made this change and if they have a preferred way to implement this logic before submitting a fix.

huangapple
  • 本文由 发表于 2023年6月26日 00:00:51
  • 转载请务必保留本文链接：https://go.coder-hub.com/76551285.html
  • choroplethr
  • r
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定