如何解决将定量添加到每行ggplot2
我试图在具有三行的这张图上阴影0.025和0.975的分位数。我已经尝试过geom_area，geom_ribbon，但无法突出显示每一行中的每个分位数。
请注意，此密度图中忽略了“ y”。
example <-data.frame(source=c("Leaflitter","Leaflitter","Biofilm","Algae","Algae"),n=c(1,2,3,4,5,6,7,8,9,10,1,10),density=c(0.554786934,0.650578421,0.039317168,0.53537613,0.435081982,0.904056941,0.556284164,0.855319434,0.399169622,0.570246304,0.076722032,0.257427999,0.172736928,0.447424473,0.520976948,0.011720494,0.311348655,0.120698996,0.016336661,0.331741377,0.368491034,0.09199358,0.787945904,0.017199397,0.04394107,0.084222564,0.132367181,0.023981569,0.584493716,0.098012319))
example
一个亚组和分位数
L <- filter(QPA_G_Feb17,source == "Leaflitter")
L <-as.data.frame(L)
Lq025 <- quantile(L$density,.025)
Lq975 <- quantile(L$density,.975)
ggplot(QPA_G_Feb17,aes(x=density,color=source)) +
labs(y="Density",x="Sorce contribution") +
geom_density(aes(linetype = source),size=1.2) +
scale_color_manual(values=c("#31a354","#2c7fb8","#d95f0e")) +
scale_linetype_manual(values = c("solid","dotted","longdash")) +
theme_classic()+
ylim(0,5)+
theme(axis.text.y=element_text(angle=0,size=12,vjust=0.5,color="black")) +
theme(axis.text.x =element_text(angle=0,color="black")) +
theme(axis.title.x = element_text(color="black",size=14))+
theme(axis.title.y = element_text(color="black",size=14))
解决方法
我认为该数据更能代表您的绘图中显示的数据：
ggplot我发现，当您尝试在dens <- lapply(split(QPA_G_Feb17,QPA_G_Feb17$source),function(x) density(x$density,from = 0,to = 1))
df <- do.call(rbind,mapply(function(x,y) {
data.frame(x = x$x,y = x$y,source = y)
},dens,names(dens),SIMPLIFY = FALSE))
df <- df %>%
group_by(source) %>%
mutate(cdf = cumsum(y * mean(diff(x))),lower = cdf < 0.025,upper = cdf > 0.975)
中进行复杂或非标准的操作时，最好的方法是提前计算要绘制的数据。在这种情况下，我们可以计算密度曲线和累积密度（包括它们的0.025和0.975分位数），并将它们全部放入这样的数据框中：
geom_area现在很容易使用ggplot(df,aes(x,y,color = source)) +
geom_area(data = df[df$lower,],aes(fill = source),alpha = 0.5,position = "identity") +
geom_area(data = df[df$upper,position = "identity") +
labs(y = "Density",x = "Source contribution") +
geom_line(aes(linetype = source),size = 1.2) +
scale_fill_manual(values = c("#31a354","#2c7fb8","#d95f0e")) +
scale_color_manual(values = c("#31a354","#d95f0e")) +
scale_linetype_manual(values = c("solid","dotted","longdash")) +
theme_classic() +
ylim(0,5) +
xlim(0,1) +
theme(axis.text.y = element_text(size = 12,vjust = 0.5),axis.text.x = element_text(size = 12,axis.title.x = element_text(size = 14),axis.title.y = element_text(size = 14))
进行绘制：
arrowsize= 
在这里，每条密度曲线的2.5％和97.5％的极限值在每条线下方均以阴影显示。唯一的例外是“ Leaflitter”行，该行明显超出了示例中绘制的0-1范围。
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。
