我有一个看起来像这样的数据框。
> head(df)
DGene JGene cdr3_len Sum
1 IGHD1 IGHJ1 0 22
2 IGHD1 IGHJ1 1 11
3 IGHD1 IGHJ1 2 16
4 IGHD1 IGHJ1 3 40
5 IGHD1 IGHJ1 4 18
6 IGHD1 IGHJ1 5 30
...
facet_grid非常简单。
ggplot(df,aes(x=cdr3_len,y=Sum)) + geom_line() + xlim(c(1,42)) + facet_grid(JGene~DGene,scales="free_y")
并得到类似的东西。
我想知道是否有人可以帮助我在每个网格的均值上添加一条线。或者可能如何在右上角打印每个网格的均值。
谢谢,
编辑-完整链接到数据框
这是一种cdr3_len
通过预先计算所需值(根据@ jwillis0720的注释)为平均值添加文本和垂直线的方法:
首先,计算cdr3_len
每个面板的平均值,然后left_join
将该数据框计算为第二个数据框,该数据框计算用于将文本放置在每个面板上的适当y值(因为适当的y值仅根据的水平而变化JGene
)。
library(dplyr)
meanData = df %>% group_by(JGene, DGene) %>%
summarise(meanCDR = sum(Sum*cdr3_len)/sum(Sum)) %>%
left_join(df %>% group_by(JGene) %>%
summarise(ypos = 0.9*max(Sum)))
现在为情节:
ggplot(df,aes(x=cdr3_len, y=Sum)) +
geom_vline(data=meanData, aes(xintercept=meanCDR), colour="red", lty=3) +
geom_line() +
geom_text(data=meanData,
aes(label=round(meanCDR,1), x=40, y=ypos), colour="red",
hjust=1) +
xlim(c(1,42)) +
facet_grid(JGene~DGene,scales="free_y")
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句