[英]Position dodge on a flipped axis with error bars in ggplot2
我正在使用以下數據框在ggplot2中創建深度圖:
TTTTOM<- data.frame(Treatment=c("Mineral", "Mineral", "Mineral", "OM", "OM", "OM", "Veg", "Veg", "Veg"), Depth2=c(-7.5, -4.5, -1.5, -7.5, -4.5, -1.5, -7.5, -4.5, -1.5), mean=c(2.83, 3.33, 3.16, 9.16, 11.17, 11.67, 4.83, 5, 5.17), se=c(0.7, 1.12, 0.65, 2.41, 3.28, 3.12, 1.83, 1.71, 0.95))
我已經將此代碼與ggplot2一起使用來創建繪圖,但是位置閃避似乎並沒有像處理線和點一樣使用誤差線。 我需要我的錯誤欄來顯示我的觀點(即,讓它們以相同的數量躲避)。 謝謝你的幫助。
pd <- position_dodge(0.3)
ggplot(TTTTOM, aes(x=Depth2, y=mean, linetype=Treatment)) +
geom_line(aes(linetype=Treatment), colour="black", size=0.75, position=pd) +
geom_errorbar(aes(ymin=mean-se, ymax=mean+se), linetype=1, colour="black", width=.1, position=pd) +
geom_point(colour="black",shape=16, size=2.0, position=pd) +
scale_linetype_manual(values=c(1, 2, 4), labels=c("Organic", "Mineral", "Vegetated"))+
geom_vline(aes(xintercept=0), linetype="dashed")+
theme(legend.position="none",
axis.title.x=element_text(size=18, margin=margin(10,0,0,0)),
legend.title=element_blank(),
legend.text=element_text(size=18),
legend.key.width=unit(2,"line"),
axis.text=element_text(size=18),
axis.title.y=element_text(size=16, margin=margin(0,10,0,0)),
plot.margin = unit(c(.25, .25, .25, .25), "in"))+
labs(x=paste("Distance from \n sediment surface (cm)"), y="% Total OM") +
expand_limits(x=c(-10,0))+
coord_flip()
您對geom_errorbar()的調用需要一個組參數; 在這種情況下, group = Treatment :
pd <- position_dodge(.3)
ggplot(TTTTOM, aes(x=Depth2, y=mean, linetype=Treatment)) +
geom_line(aes(linetype=Treatment), colour="black", size=0.75, position=pd) +
geom_errorbar(aes(ymin=mean-se, ymax=mean+se, group=Treatment),
linetype=1, colour="black", width=.1, position=pd) +
geom_point(colour="black",shape=16, size=2.0, position=pd) +
scale_linetype_manual(values=c(1, 2, 4), labels=c("Organic", "Mineral", "Vegetated"))+
geom_vline(aes(xintercept=0), linetype="dashed")+
theme(legend.position="none",
axis.title.x=element_text(size=18, margin=margin(10,0,0,0)),
legend.title=element_blank(),
legend.text=element_text(size=18),
legend.key.width=unit(2,"line"),
axis.text=element_text(size=18),
axis.title.y=element_text(size=16, margin=margin(0,10,0,0)),
plot.margin = unit(c(.25, .25, .25, .25), "in"))+
labs(x=paste("Distance from \n sediment surface (cm)"), y="% Total OM") +
expand_limits(x=c(-10,0))+
coord_flip()
輸出:
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.