2013-09-24 76 views
0

我正在制作一个带有ggplot2的箱线图,并尝试(没有任何成功)将每个变量的平均值放在它上面。有人可以指出什么是错的?带平均分的箱线图

重复的例子:

data <- data.frame(id=1:5000, Estimativa0=gl(5, 1000), NM.Intra1=rnorm(5000, rep(rnorm(5, 0, 4), each=1000)), NM.Entre2=rnorm(5000, rep(rnorm(5, 0, 4), each=1000)), SPADE3=rnorm(5000, rep(rnorm(5, 0, 4), each=1000)), NCI4=rnorm(5000, rep(rnorm(5, 0, 4), each=1000))) 
data <- reshape(data, direction='long', varying=list(c("NM.Intra1", "NM.Entre2", "SPADE3", "NCI4")), idvar=c('id', 'Estimativa0'), times=c("NM.Intra1", "NM.Entre2", "SPADE3", "NCI4"), timevar='Método') 
levels(data$Estimativa0) <- c('Est 1', 'Est 2', 'Est 3', 'Est 4', 'Est 5', 'Est 6') 
p <- ggplot(data, aes(Método, NM.Intra1)) 
p <- p + geom_boxplot(aes(fill=factor(Estimativa0))) + ylab('Vício') + xlab('Método') + scale_fill_discrete("Estatística") + ggtitle('Gráfico maluco da Greice') + geom_hline(yintercept=0) 
p 

的第一次尝试是:

medias <- with(data, aggregate(data.frame(Media=NM.Intra1), by=list(Estimativa0=Estimativa0, Método=Método), FUN=mean)) 
p + geom_point(aes(x=Método, y=Media), data=medias, shape=2) 

但在相同的 '竖线' ploted所有点。我试图解决方案指出here,但不能太:

p + stat_summary(fun.y=mean, colour="red", geom="point", position=position_dodge(width=0.75), shape=2) 

所以,我错过了什么?

回答

4

我假设你想要每个Método - Estimativa0组合的平均值。在这种情况下,您需要通过Estimativa0告知stat_summary进行分组;现在,映射只能在geom_boxplot中完成,因此不会被stat_summary继承。

因此,删除:

aes(fill=factor(Estimativa0)) 

而在初始映射把fill=factor(Estimativa0)

aes(Método, NM.Intra1, fill=factor(Estimativa0)) 

你会这样结束了:

p <- ggplot(data, aes(Método, NM.Intra1, fill=factor(Estimativa0))) 
p <- p + geom_boxplot() + ylab('Vício') + xlab('Método') + scale_fill_discrete("Estatística") + ggtitle('Gráfico maluco da Greice') + geom_hline(yintercept=0) 
p + stat_summary(fun.y=mean, colour="red", geom="point", position=position_dodge(width=0.75), shape=2) 

产生这样的:

enter image description here