r - DCA:使用自动绘图或 ggplot2 标记点
问题描述
我发现很难在 autoplot 或 ggplot 中为带有 DCA 的站点添加标签。
我还想根据它们的组来区分 autoplot/ggplot 上的点。
这是我使用的数据和代码,它运行良好,直到自动绘图/ggplot 的命令:
library(vegan)
data(dune)
d <- vegdist(dune)
csin <- hclust(d, method = "single")
cl <- cutree(csin, 3)
dune.dca <- decorana(dune)
autoplot(dune.dca)
这是获得的自动绘图:
我正在使用简单的编码,我尝试了这些代码,但它们并没有引导我到任何地方:
autoplot(dune.dca, label.size = 3, data = dune, colour = cl)
ggplot(dune.dca(x=DCA1, y=DCA2,colour=cl))
ggplot(dune.dca, display = ‘site’, pch = 16, col = cl)
ggrepel::geom_text_repel(aes(dune.dca))
如果有人有一个简单的建议,那就太好了。
解决方案
通过添加的信息(包),我能够深入挖掘。
问题是(简而言之)autoplot.decorana
将数据添加到特定层(geom_point
或geom_text
)。这不会继承到其他层,因此添加其他层会导致空白页。
基本上注意到下面的 2 个代码字符串之一会导致错误,并注意data
参数的位置:
# Error:
ggplot() +
geom_point(data = mtcars, mapping = aes_string(x = 'hp', y = 'mpg')) +
geom_label(aes(x = hp, y = mpg, label = cyl))
# Work:
ggplot(data = mtcars) +
geom_point(mapping = aes_string(x = 'hp', y = 'mpg')) +
geom_label(aes(x = hp, y = mpg, label = cyl))
ggvegan:::autoplot.decorana
像示例中一样放置数据,返回错误。
我看到了解决这个问题的两种方法:
ggplot_build
使用或提取图层数据layer_data
并创建整体或单层映射。- 提取用于生成数据的代码,并手动创建我们的绘图(不使用
autoplot
)。
老实说,我认为第二个更简单,因为我们可能必须提取更多信息才能使我们的数据变得合理。通过查看源代码ggvegan:::autoplot.decorana
(只需通过省略括号将其打印到控制台),我们可以提取以下代码,该代码生成与图中使用的数据相同的数据
ggvegan_data <- function(object, axes = c(1, 2), layers = c("species", "sites"), ...){
obj <- fortify(object, axes = axes, ...)
obj <- obj[obj$Score %in% layers, , drop = FALSE]
want <- obj$Score %in% c("species", "sites")
obj[want, , drop = FALSE]
}
有了这个,我们就可以生成我们想要的任何图,使用适当的映射而不是层单独的映射
dune.plot.data <- ggvegan_data(dune.dca)
p <- ggplot(data = dune.dca, aes(x = DCA1, DCA2, colour = Score)) +
geom_point() +
geom_text(aes(label = Label), nudge_y = 0.3)
p
推荐阅读
- android - 插页式广告的点击率非常高
- php - 删除字符串php中的数字
- arrays - 如何获得两个数组的大小?
- jquery - 清空之前的数据后无法显示预测数据
- javascript - How to check and replace undefined and zero inside Javascript array map?
- xamarin - Xamarin iOS 对 Xcode 11.2.1 GM 的支持
- java - 如何将 xml 视图关联到 java 类?
- python - 除 IDLE 外,其他 IDE 不适用于 openCV
- amazon-web-services - 如何为弹性豆茎容器创建 CNAME?
- java - 有没有办法在假代理中将 HttpServletResponse 作为 RequestParameter 转发?