Я не уверен, что вызвало ошибку. Кажется, что использовался не символьный вектор, хотя функция, возможно, требует символьный вектор. Вы хотите сравнить свой код с моим. Возможно, вы также хотите проверить данные. Следующее работает для меня.
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
))