2013-09-25 16 views
5

でカスタムの行を追加します。私はこのコードでこのプロット enter image description hereggplotのbarplot

を作成しました:

ggplot(data_long, aes(Trial,value,fill=factor(Trial))) +  
    stat_summary(fun.y=mean,geom="bar") + facet_grid(Task~Gruppo) + labs(x="Trial 
    type",y="Accuracy %") + theme(legend.position="none") 

は今、私は値の夫婦間の違いを示すカスタムの行を追加する必要があります。ここで(P = 0.46で最初の2小節を参照してください)私が何をしたいの例です。

enter image description here

私は解決策については考えている、と私は使用するので、物事が私のためにも、より複雑ですfacet_grid。 誰でも助けてくれますか?

答えて

6

まず、サンプルデータが提供されないため、私自身のサンプルデータを作成しました。これらのデータは、すでにレベルの各組み合わせごとに1つだけ値がある(要約されている。

set.seed(1) 
df<-data.frame(expand.grid(c("Control","Effect"),c("Self","Other"),c("Type1","Type2")), 
    runif(8,0,1)) 
colnames(df)<-c("Treatment","Group","Type","value") 
df 
    Treatment Group Type  value 
1 Control Self Type1 0.2655087 
2 Effect Self Type1 0.3721239 
3 Control Other Type1 0.5728534 
4 Effect Other Type1 0.9082078 
5 Control Self Type2 0.2016819 
6 Effect Self Type2 0.8983897 
7 Control Other Type2 0.9446753 
8 Effect Other Type2 0.6607978 

今、あなたは線の位置のための2つの新しい値を追加する必要があります。ymin値が元の値を加えた小さな定数である。ymax値(Treatment及びグルーピングとしてTypeを使用して)各ファセットについて計算し、それが最大のファセットの値に加え、いくつかの定数である。

library(plyr) 
df<-ddply(df,.(Treatment,Type),transform,ymax=max(value)+0.2) 
df$ymin<-df$value+0.05 
df 
    Treatment Group Type  value  ymax  ymin 
1 Control Self Type1 0.2655087 0.7728534 0.3155087 
2 Control Self Type2 0.2016819 1.1446753 0.2516819 
3 Control Other Type1 0.5728534 0.7728534 0.6228534 
4 Control Other Type2 0.9446753 1.1446753 0.9946753 
5 Effect Self Type1 0.3721239 1.1082078 0.4221239 
6 Effect Self Type2 0.8983897 1.0983897 0.9483897 
7 Effect Other Type1 0.9082078 1.1082078 0.9582078 
8 Effect Other Type2 0.6607978 1.0983897 0.7107978 

第二のデータフレームは、ラベルのために作られている。ここで、各ファセットy位置にすると、再び元であるymax値にある定数とlabには、表示する必要のあるラベルが含まれています。

df.names<-ddply(df,.(Treatment,Type),summarise,ymax=ymax[1]+0.1) 
df.names$lab<-c("p=0.46","**","***","*") 
df.names 
    Treatment Type  ymax lab 
1 Control Type1 0.8728534 p=0.46 
2 Control Type2 1.2446753  ** 
3 Effect Type1 1.2082078 *** 
4 Effect Type2 1.1983897  * 

ようになりましたdfが既に要約されている値の代わりに使用stat_summary()geom_bar(stat="identity")。追加の行には2つのコールが追加されます。最初は縦線をプロットし、2番目の行は水平線を追加します。 geom_text()は、行の上にラベルを追加します。

ggplot(df, aes(Group,value,fill=Group)) +  
    geom_bar(stat="identity") + facet_grid(Type~Treatment) + 
    theme(legend.position="none")+ 
    geom_segment(aes(x=Group,xend=Group,y=ymin,yend=ymax))+ 
    geom_segment(aes(x="Self",xend="Other",y=ymax,yend=ymax))+ 
    geom_text(data=df.names,aes(x=1.5,y=ymax,label=lab),inherit.aes=FALSE) 

enter image description here

+0

あなたの忍耐、Didzisのために非常に感謝します。それは私にとって非常に明確で完全です。 –