r - 在 R 中，ggplot geom_point() 具有基于特定离散值的颜色 - 第 2 部分

Question

我的问题类似于this one ，除了我的数据不同。就我而言，我无法使用给出的解决方案。我希望点会根据 cut() 值显示在我的 map 上。有人能指出我正确的方向吗？

> test
# A tibble: 10 × 5
        TC1    TC2       Lat      Long       Country
      <dbl>  <dbl>     <dbl>     <dbl>        <fctr>
1       2.9 2678.0  50.62980 -95.60953        Canada
2    1775.7 5639.9 -31.81889 123.19389     Australia
3       4.4 5685.6 -10.10449  38.54364      Tanzania
4       7.9     NA  54.81822 -99.91685        Canada
5      11.2 2443.0   7.71667  -7.91667 Cote d'Ivoire
6     112.1 4233.4 -17.35093 128.02609     Australia
7       4.4  114.6  45.21361 -67.31583        Canada
8    8303.5 4499.9  46.63626 -81.39866        Canada
9  100334.8 2404.5  46.67291 -93.11937           USA
10       NA 1422.9 -17.32921  31.28224      Zimbabwe

ggplot(data = test, aes(x= Long, y= Lat)) + 
  borders("world", fill="gray75", colour="gray75", ylim = c(-60, 60)) +
  geom_point(aes(size=TC2, col=cut(TC1, c(-Inf, 1000, 5000, 50000, Inf)))) +
  # scale_colour_gradient(limits=c(100, 1000000), low="yellow", high="red") +
  scale_color_manual(name = "TC1",
                     values = c("(-Inf,1000]" = "green",
                                "(1000,5000]" = "yellow",
                                "(5000,50000]" = "orange",
                                "(50000, Inf]" = "red"),
                     labels = c("up to 1", "1 to 5", "5 to 50", "greater than 50")) +
   theme(legend.position = "right") +  
   coord_quickmap()
Warning message:
Removed 10 rows containing missing values (geom_point).

Answer 1

你快到了!这只是不正确的“削减”因素的名称。如果您尝试:

cut(test$TC1, c(-Inf, 1000, 5000, 50000, Inf))
#  [1] (-Inf,1e+03]  (1e+03,5e+03] (-Inf,1e+03]  (-Inf,1e+03]  (-Inf,1e+03] 
#  [6] (-Inf,1e+03]  (-Inf,1e+03]  (5e+03,5e+04] (5e+04, Inf]  <NA>         
# Levels: (-Inf,1e+03] (1e+03,5e+03] (5e+03,5e+04] (5e+04, Inf]

如您所见，关卡名称与您输入的名称略有不同。

library(ggplot2)

ggplot(data = test, aes(x = Long, y = Lat)) +
    borders("world", fill="gray75", colour="gray75", ylim = c(-60, 60)) +
    geom_point(aes(size=TC2, color = cut(TC1, c(-Inf, 1000, 5000, 50000, Inf)))) +
    scale_color_manual(name = "TC1",
    values = c("(-Inf,1e+03]" = "green",
               "(1e+03,5e+03]" = "yellow",
               "(5e+03,5e+04]" = "orange",
               "(5e+04, Inf]" = "red"),
    labels = c("up to 1", "1 to 5", "5 to 50", "greater than 50")) +
    theme(legend.position = "right") +  
    coord_quickmap()
#> Warning: Removed 2 rows containing missing values (geom_point).

数据:

test <- read.table(text = 'TC1    TC2 Lat Long    Country
1       2.9 2678.0  50.62980 -95.60953        Canada
2    1775.7 5639.9 -31.81889 123.19389     Australia
3       4.4 5685.6 -10.10449  38.54364      Tanzania
4       7.9     NA  54.81822 -99.91685        Canada
5      11.2 2443.0   7.71667  -7.91667 "Cote d\'Ivoire"
6     112.1 4233.4 -17.35093 128.02609     Australia
7       4.4  114.6  45.21361 -67.31583        Canada
8    8303.5 4499.9  46.63626 -81.39866        Canada
9  100334.8 2404.5  46.67291 -93.11937           USA
10       NA 1422.9 -17.32921  31.28224      Zimbabwe', header = T)