0

我想在 r 上制作一张地图，并通过颜色的强度显示价格弹性。这是我使用线性回归得到的数据：

在此处输入图像描述

    IN  log_PRICE   -1.1059770  0.05943414  -18.60845   7.186435e-72
    KY  log_PRICE   -1.0459502  0.03410250  -30.67078   2.739798e-196
    OH  log_PRICE   -0.9076732  0.01259117  -72.08806   0.000000e+00
    TX  log_PRICE   -0.2252409  0.01053847  -21.37321   4.117490e-101
    IN  log_PRICE   -1.1059770  0.05943414  -18.60845   7.186435e-72
    KY  log_PRICE   -1.0459502  0.03410250  -30.67078   2.739798e-196
    OH  log_PRICE   -0.9076732  0.01259117  -72.08806   0.000000e+00
    TX  log_PRICE   -0.2252409  0.01053847  -21.37321   4.117490e-101

library(usmap)
library(ggplot2)
library(tmpa)
library(sf)
library(leaflet)

plot_usmap(regions="state", data = CC1, values = "estimate", 
    include = c("IN", "KY", "OH", "TX"), color = "orange") + 
  scale_fill_continuous(low = "white", high = "orange", 
     name = "Price elasticity", label = scales::comma) + 
  labs(title = "Price elasticity for Cold Cereal", 
       subtitle = "States include: IN, KY, OH, TX") +
  theme(legend.position = "right") 

#That is the code I wrote but I keep getting this error
Error in match.arg(regions) : 'arg' must be NULL or a character vector
4

1 回答 1

2

我不确定是什么导致了错误。尽管该函数可能需要字符向量，但似乎使用了非字符向量。您想在您的代码和我的代码之间进行比较。也许，您也想检查数据。以下对我有用。

plot_usmap(data = mydata, values = "estimate",
           regions = "state", include = c("IN", "KY", "OH", "TX")) +
scale_fill_continuous(low = "white", high = "orange", 
                      name = "Price elasticity", label = scales::comma) + 
labs(title = "Price elasticity for Cold Cereal", 
     subtitle = "States include: IN, KY, OH, TX") +
theme(legend.position = "right")

在此处输入图像描述

为了好玩，以下是我自己的带有 albersusa 包的版本。

library(tydyverse)
library(sf)
library(albersusa)
library(viridis)

mysf <- left_join(usa_sf("laea"), mydata, by = c("iso_3166_2" = "state")) 

ggplot() +
geom_sf(data = mysf, aes(fill = estimate)) +
scale_fill_viridis(discrete = FALSE, option = "plasma")

在此处输入图像描述

数据

mydata <- structure(list(state = c("IN", "KY", "OH", "TX"), log = c("log_PRICE", 
"log_PRICE", "log_PRICE", "log_PRICE"), estimate = c(-1.105977, 
-1.0459502, -0.9076732, -0.2252409), se = c(0.05943414, 0.0341025, 
0.01259117, 0.01053847), statistic = c(-18.60845, -30.67078, 
-72.08806, -21.37321), pvalue = c(7.186435e-72, 2.739798e-196, 
0, 4.11749e-101)), class = "data.frame", row.names = c(NA, -4L
))
于 2020-02-21T08:59:50.907 回答