blocks|key|431164|text|我已经开始研究一个解决方案，但我有点卡住了。我张贴这篇文章是希望一些核心的数据可视化人员可以为一个真正的解决方案做出贡献(如果这不是解决方案的方法，请让我知道，我会删除这个半生不熟的项目)。|type|unstyled|depth|inlineStyleRanges|entityRanges|data|431165|#+data
require(vegan)
data(dune)
data(dune.env)

#+RDA
RDA+<-+rda(dune+~+A1,+data+=+dune.env)

#+extract+species+scores
df1+<-+data.frame(RDA.scores+=+scores(RDA)$sites[+,+1],+Management+=+dune.env$Management)
df2+<-+data.frame(spec_scores+=+scores(RDA)$species[+,+1],+taxa+=+rownames(scores(RDA)$species))
df2+<-+df2[abs(df2$spec_scores)+>+0.05,+]
df2+<-+cbind(df2,+x+=+rep(1:4,+6))

#+plot
P+<-+ggplot(data+=+df1,+aes(+x+=+Management,+y+=+RDA.scores)+)+%2B+geom_boxplot()+%2B+theme_bw()++%2Bgeom_hline(y=0,+lty=2)+%2B+geom_rug(data=+df2,+mapping=aes(x+=+x+,+y+=+spec_scores),+sides="r")+%2B+geom_text(data=+df2,+mapping=aes(x=5,+y=spec_scores,+label=taxa),+size=2)+%2B+xlim(c("BF",+"HF",+"NM",+"SF",+""))
P|code-block|syntax|javascript|431166|​|431167|📷|atomic|offset|length|431168|431169|entityMap|0|IMAGE|mutability|IMMUTABLE|imageUrl|https://ask.qcloudimg.com/http-save/yehe-900000/6d079b04b097f9c52a536166be5d960a.png|imageAlt^0|0|0|0|0|1|0|0|0^^$0|@$1|2|3|4|5|6|7|X|8|@]|9|@]|A|$]]|$1|B|3|C|5|D|7|Y|8|@]|9|@]|A|$E|F]]|$1|G|3|H|5|6|7|Z|8|@]|9|@]|A|$]]|$1|I|3|J|5|K|7|10|8|@]|9|@$L|11|M|12|1|13]]|A|$]]|$1|N|3|H|5|6|7|14|8|@]|9|@]|A|$]]|$1|O|3|-4|5|6|7|15|8|@]|9|@]|A|$]]]|P|$Q|$5|R|S|T|A|$U|V|W|-4]]]]

I have started to work on a solution, but I'm kinda stuck. I'm posting it in the hope that someone of the hardcore data visualization people can contribute to a real solution (if this is not the way to do it please let me know and I'll delete this half-baked project).

<pre><code># data
require(vegan)
data(dune)
data(dune.env)

# RDA
RDA &lt;- rda(dune ~ A1, data = dune.env)

# extract species scores
df1 &lt;- data.frame(RDA.scores = scores(RDA)$sites[ , 1], Management = dune.env$Management)
df2 &lt;- data.frame(spec_scores = scores(RDA)$species[ , 1], taxa = rownames(scores(RDA)$species))
df2 &lt;- df2[abs(df2$spec_scores) &gt; 0.05, ]
df2 &lt;- cbind(df2, x = rep(1:4, 6))

# plot
P &lt;- ggplot(data = df1, aes( x = Management, y = RDA.scores) ) + geom_boxplot() + theme_bw() +geom_hline(y=0, lty=2) + geom_rug(data= df2, mapping=aes(x = x , y = spec_scores), sides="r") + geom_text(data= df2, mapping=aes(x=5, y=spec_scores, label=taxa), size=2) + xlim(c("BF", "HF", "NM", "SF", ""))
P
</code></pre>

<img src="https://i.stack.imgur.com/WL9nF.png" alt="half-baked project plot">

blocks|key|212923|text|这是仅使用ggplot2的解决方案，因此它不会是最优雅的解决方案。|type|unstyled|depth|inlineStyleRanges|offset|length|style|CODE|entityRanges|data|212924|我从@Eric+Fail建议的第一步开始。|212925|require(vegan)
data(dune)
data(dune.env)

#+RDA
RDA+<-+rda(dune+~+A1,+data+=+dune.env)

#+extract+species+scores
df1+<-+data.frame(RDA.scores+=+scores(RDA)$sites[+,+1],+Management+=+dune.env$Management)
df2+<-+data.frame(y+=+scores(RDA)$species[+,+1],+taxa+=+rownames(scores(RDA)$species))
#+Order+data+according+to+species+scores
df2<-df2[order(df2$y),]

#+define+values+for+rugs+(segments).+4.9+and+5.0+used+because+data+has+4+levels+(%2B1)
df2$x=4.9
df2$xend=5

#+define+coordinates+for+names+(30+is+number+of+species)
df2$yend2<-seq(min(df1$RDA.scores),max(df1$RDA.scores),length.out=30)
df2$xend2=5.3

#+plot
P+<-+ggplot(data+=+df1,+aes(+x+=+Management,+y+=+RDA.scores)+)+%2B+
+++++geom_segment(y=0,yend=0,x=0,xend=5,+lty=2,+size+=+I(0.3))+%2B+
+++++geom_boxplot()+%2B+
+++++#add+extra+levels+to+get+space
+++++scale_x_discrete(limits=c("BF",+"HF",+"NM",+"SF",+"",""))+%2B+
+++++#set+y+scale
+++++scale_y_continuous(limits=c(-3,5),expand=c(0,0))+%2B
+++++#add+rugs+as+segments+and+add+segments+connecting+rugs+and+texts
+++++geom_segment(data=+df2,+mapping=aes(x=x,xend=xend,y=y,yend=y),+size+=+I(0.3))+%2B
+++++geom_segment(data=+df2,+mapping=aes(x=xend,xend=xend2,y=y,yend=yend2),+size+=+I(0.4))+%2B
+++++#add+texts
+++++geom_text(data=df2,mapping=aes(x=xend2%2B0.1,y=yend2,label=taxa),hjust=0,+size=3)+%2B
+++++#add+rectangular+to+imitate+box+around+plot
+++++geom_rect(xmax=5,xmin=0.4,ymax=5,ymin=-3,colour="black",fill=NA)

#+Final+adjustments+of+plot
P%2Btheme(axis.line=element_blank(),
++++++++panel.grid=element_blank(),
++++++++panel.background=element_blank())|code-block|syntax|javascript|212926|​|212927|📷|atomic|212928|212929|entityMap|0|IMAGE|mutability|IMMUTABLE|imageUrl|https://ask.qcloudimg.com/http-save/yehe-900000/dff7a21e14c0dbcaf98ef66577bd600a.png|imageAlt^0|5|7|0|0|0|0|0|1|0|0|0^^$0|@$1|2|3|4|5|6|7|11|8|@$9|12|A|13|B|C]]|D|@]|E|$]]|$1|F|3|G|5|6|7|14|8|@]|D|@]|E|$]]|$1|H|3|I|5|J|7|15|8|@]|D|@]|E|$K|L]]|$1|M|3|N|5|6|7|16|8|@]|D|@]|E|$]]|$1|O|3|P|5|Q|7|17|8|@]|D|@$9|18|A|19|1|1A]]|E|$]]|$1|R|3|N|5|6|7|1B|8|@]|D|@]|E|$]]|$1|S|3|-4|5|6|7|1C|8|@]|D|@]|E|$]]]|T|$U|$5|V|W|X|E|$Y|Z|10|-4]]]]

This is solution only using <code>ggplot2</code>, so it won't be most elegant one.

I started with first steps as suggested by @Eric Fail.

<pre><code>require(vegan)
data(dune)
data(dune.env)

# RDA
RDA &lt;- rda(dune ~ A1, data = dune.env)

# extract species scores
df1 &lt;- data.frame(RDA.scores = scores(RDA)$sites[ , 1], Management = dune.env$Management)
df2 &lt;- data.frame(y = scores(RDA)$species[ , 1], taxa = rownames(scores(RDA)$species))
# Order data according to species scores
df2&lt;-df2[order(df2$y),]

# define values for rugs (segments). 4.9 and 5.0 used because data has 4 levels (+1)
df2$x=4.9
df2$xend=5

# define coordinates for names (30 is number of species)
df2$yend2&lt;-seq(min(df1$RDA.scores),max(df1$RDA.scores),length.out=30)
df2$xend2=5.3

# plot
P &lt;- ggplot(data = df1, aes( x = Management, y = RDA.scores) ) + 
 geom_segment(y=0,yend=0,x=0,xend=5, lty=2, size = I(0.3)) + 
 geom_boxplot() + 
 #add extra levels to get space
 scale_x_discrete(limits=c("BF", "HF", "NM", "SF", "","")) + 
 #set y scale
 scale_y_continuous(limits=c(-3,5),expand=c(0,0)) +
 #add rugs as segments and add segments connecting rugs and texts
 geom_segment(data= df2, mapping=aes(x=x,xend=xend,y=y,yend=y), size = I(0.3)) +
 geom_segment(data= df2, mapping=aes(x=xend,xend=xend2,y=y,yend=yend2), size = I(0.4)) +
 #add texts
 geom_text(data=df2,mapping=aes(x=xend2+0.1,y=yend2,label=taxa),hjust=0, size=3) +
 #add rectangular to imitate box around plot
 geom_rect(xmax=5,xmin=0.4,ymax=5,ymin=-3,colour="black",fill=NA)

# Final adjustments of plot
P+theme(axis.line=element_blank(),
 panel.grid=element_blank(),
 panel.background=element_blank())
</code></pre>

<img src="https://i.stack.imgur.com/8OumT.png" alt="ggplot2 boxplot">

Can I make such a plot with ggplot2?

<pre><code># data
require(vegan)
data(dune)
data(dune.env)

# RDA
RDA &lt;- rda(dune ~ A1, data = dune.env)

# extract species scores
df &lt;- data.frame(spec_scores = scores(RDA)$species[ , 1], 
 taxa = rownames(scores(RDA)$species))
df &lt;- df[abs(df$spec_scores) &gt; 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)
</code></pre>

<img src="https://i.stack.imgur.com/XQflp.jpg" alt="enter image description here">

Basicly I am looking for a way how to plot a labelled rug beneath the boxplot.

Any hints or suggestions?

ggplot2 boxplot with labelled rug

翻译质量差，导致语言生硬或混乱。

没有提供实际的解决方法或示例。

解答不清晰，无法理解或解决问题。

页面排版不美观，阅读体验差。

文章

问答

视频

学习中心

腾讯云实验室

直播

竞赛

腾讯云代码分析专区

腾讯iOA零信任安全管理系统专区

腾讯云架构师技术同盟交流圈

腾讯云数据库专区

腾讯云顾问专区

腾讯云原生专区

腾讯混元专区

腾讯云TCE专区

腾讯云Lighthouse专区

腾讯云HAI专区

腾讯云Edgeone专区

腾讯云存储专区

腾讯云智能专区

腾讯轻联专区 

腾讯云开发专区

TAPD专区

腾讯轻量云游戏服专区

腾讯云最具价值专家

腾讯云架构师技术同盟

腾讯云创作之星

腾讯云开发者先锋

腾讯云代码助手

云原生构建

TAPD 敏捷项目管理

Cloud Studio

SDK中心

API中心

命令行工具

涵盖代码开发、场景应用、自动测试全流程，助你从零构建专属AI助手

一站式MCP教程库，解锁AI应用新玩法

我可以用ggplot2绘制这样的图吗？# datarequire(vegan)data(dune)data(dune.env)# RDARDA <- rda(dune ~ A1, data = dune.env)# extract species scoresdf <- data.frame(spec_scores =...

问带有标签地毯的ggplot2盒子图
EN

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问带有标签地毯的ggplot2盒子图EN

社区

活动

圈层

关于

腾讯云开发者

热门产品

热门推荐

更多推荐

问带有标签地毯的ggplot2盒子图
EN