0

我想在 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。

Treat 是分类 rootpath 是分类 LWP_meas2 是连续的

在此处输入图像描述

4

2 回答 2

1

正如本所说，您可能想要这样的东西（由于没有提供数据，因此必须未经测试：

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"))
      )
于 2013-02-03T01:19:46.413 回答
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)
于 2013-02-03T03:44:37.490 回答