PROSAGA码农传奇-模型可视化-物种丰度图与ggmap

<div class =“post-text”itemprop =“text”>
  <P>
    欢迎来到StackOverflow！为了帮助人们帮助你，你应该努力提供一个
    <a href="https://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example">
      可重复的例子
    </A>
    。例如，这里将是您的数据的一个小代码摘录：
  </p>
   <pre>
    <code>
      library(tidyverse)
Anatomy <- tribble(~Species,  ~Lat, ~Long, ~Station,
                "A",  50,   60,       "I",
                "A",  50,   60,       "I",
                "A",  40,   30,      "II",
                "B", 50,   60,       "I",
                "B",  40,   30,      "II")

</code>
  </pre>
  <P>
    您的数据/代码存在一些问题：
  </p>
  <UL>
    <LI>
      投影：您很可能需要重新投影数据。只需查看坐标，您的点数就是50,60，而地图显示的是-50，-60。找出投影，并使用例如
       <code>
        st_transform
      </code>
       从包
       <code>
        sf
      </code>
    </LI>
    <LI>
      引用变量：你不需要再次调用数据帧，如
       <code>
        Anatomy$Latitude
      </code>
      。只是用
       <code>
        Latitude
      </code>
      。加
       <code>
        latitude
      </code>
       实际上是
       <code>
        lat
      </code>
       在你的数据！？
    </LI>
    <LI>
      聚合：我建议只使用
       <code>
        count()
      </code>
       用于查看每个站点的观测数量。
    </LI>
  </UL>
  <P>
    这是一段代码。注意我只是反转（60,50）到（-60，-50）这显然是错误的！
  </p>
   <pre class="lang-r prettyprint-override">
    <code>
      library(ggmap)
#> Loading required package: ggplot2
#> Google's Terms of Service: https://cloud.google.com/maps-platform/terms/.
#> Please cite ggmap if you use it! See citation("ggmap") for details.
library(tidyverse)

B<-get_map(location = c(lon = -56.731405, lat =-61.4831206),
           zoom = 6, 
           scale = "auto",
           maptype = c("terrain"),
           source = c("google"))

library(tidyverse)
Anatomy <- tribble(~Species,  ~Lat, ~Long, ~Station,
                "A",  50,   60,       "I",
                "A",  50,   60,       "I",
                "A",  40,   30,      "II",
                "B", 50,   60,       "I",
                "B",  40,   30,      "II")

Anatomy_clean <- Anatomy %>% 
  mutate_at(c("Lat", "Long"), funs(-1*.)) %>% 
  count(Species, Lat, Long, Station)
#> Warning: funs() is soft deprecated as of dplyr 0.8.0
#> please use list() instead
#> 
#> # Before:
#> funs(name = f(.)
#> 
#> # After: 
#> list(name = ~f(.))
#> This warning is displayed once per session.

ggmap(B) + 
  geom_point(data=Anatomy_clean,
             aes(x= Lat,y=Long, color=Species, size= n)) +
  scale_size_continuous(range=c(1,12))
#> Warning: Removed 2 rows containing missing values (geom_point).

</code>
  </pre>
  <P>
    <img src =“https://i.imgur.com/YtwXdIo.png”alt =“”/>
  </p>
</DIV>

物种丰度图与ggmap - >如何调整点的大小足够