Добавление пользовательских строк в ggplot barplot

Question

Я создал этот участок

с этим кодом:

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

Теперь мне нужно добавить пользовательские линии, которые показывают различия между парами значений. Вот пример того, что я хочу сделать (см первые два бара с р = 0,46):

Я понятия не имею, о решении, и вещи, и более сложным для меня, так как я использовал facet_grid. Может ли кто-нибудь мне помочь?

Answer 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 

Теперь вам нужно добавить два новых значения для позиций строк. 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 

Второй кадр данных производится для этикеток. Здесь в каждом положении фаска у снова оригинальный 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 уже суммированы значение geom_bar(stat="identity") вместо stat_summary(). Дополнительные строки добавляются с двумя вызовами geom_segment() - сначала один показывает вертикальные линии, а второй добавляет горизонтальную линию. 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)