我正在制作一个箱线图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)
但它将所有点绘制在同一条“垂直线”中。我尝试了此处指出的解决方案,但也失败了:
p + stat_summary(fun.y=mean, colour="red", geom="point", position=position_dodge(width=0.75), shape=2)
那么,我错过了什么?