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)

但它将所有点绘制在同一条“垂直线”中。我尝试了此处指出的解决方案,但也失败了:

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

那么,我错过了什么?

4

1 回答 1

4

我假设你想要每个组合Método的平均值。Estimativa0在这种情况下,您还需要告诉stat_summary分组依据Estimativa0;现在,映射仅在 中完成,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)

产生这个:

在此处输入图像描述

于 2013-09-24T21:07:25.897 回答