我想在R中用LWP绘制两个分类变量(治疗)。
我想改变点的颜色(填充)以对应变量rootpatH(每个处理中一个点)。
如果差异很大,我还想添加一个表示法。
Plants.means<-aggregate(Plants[10:44],Plants[c(6,9)],mean,na.rm=TRUE)
SE<-function(x)(sd(x,na.rm=TRUE)/ (length(x)-sum(is.na(x)))^0.5)
ggplot(Plants.means, aes(x=Treat, y=Yield, group=1)) +
geom_errorbar(width=.1, aes(ymin=Yield-SE(Plants.means$Yield), ymax=Yield+SE(Plants.means$Yield))) +
geom_point(shape=21, size=3, fill="white") +
ylim(0,1750)
我还想更改误差线以代表95%CI而不是SE。
善待是绝对的 rootpatH是绝对的 LWP_meas2是连续的
答案 0 :(得分:1)
正如Ben所说,你可能想要这样的东西(由于没有提供数据,因此没有必要进行测试:
with( Plants.means,
errbar(x =Treat, y = LWP_Meas2,
yplus = LWP_Meas2 + SE( LWP_Meas2),
yminus = LWP_Meas2 - SE( LWP_Meas2),
add=T,na.rm=TRUE, col=c("blue", "hotpink"))
)
答案 1 :(得分:0)
error <- function(x)(qt(0.975,df=length(x)-1)*sd(x)/sqrt(length(x)))
ggplot(Plants.means,aes(x=Treat, y=Yield, colour=rootpatH, group=rootpatH)) +
geom_errorbar(aes(ymin=Yield-error(Plants.means$Yield),ymax=Yield+error(Plants.means$Yield)), colour="black",width=.1, position=pd) +
geom_point(position=pd, size=4)