我能用ggplot2制作这样的图吗?
# data
require(vegan)
data(dune)
data(dune.env)
# RDA
RDA <- rda(dune ~ A1, data = dune.env)
# extract species scores
df <- data.frame(spec_scores = scores(RDA)$species[ , 1],
taxa = rownames(scores(RDA)$species))
df <- df[abs(df$spec_scores) > 0.05, ]
# plot
par(mar = c(5,4,4,8))
# boxplot of sites-scores along A1-axis
boxplot(scores(RDA)$sites[ , 1] ~ dune.env$Management)
abline(h = 0, lty = "dotted")
# add species scores to plot
rug(df$spec_scores, side=4)
linestack(df$spec_scores, labels=df$taxa, at = par("usr")[2], add = TRUE, hoff = 1)
我基本上正在寻找一种在箱线图下方绘制标记地毯图的方法。
有什么提示或建议吗?
autoplot()
。还需要一些头脑风暴 - 只是生成图形还是返回一个有用的对象?等等。 - Gavin Simpson