2013-02-02 50 views
0

我想在R中繪製一個情節,用LWP繪製兩個分類變量(treat)。使用ggplot在R中添加錯誤欄

我想改變點的顏色(填充)以與變量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是連續

enter image description here

+1

'SE'是一個函數,而不是一個數值;大概你想'SE(某事)'... –

+0

如何在每組點之後添加文本「ns」和「*」來表示它們是否顯着不同? – user1977802

回答

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) 
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")) 
    )