这是你脑子里有什么?
centroids <- aggregate(cbind(x,y)~class,df,mean)
ggplot(df,aes(x,y,color=factor(class))) +
geom_point(size=3)+ geom_point(data=centroids,size=5)
这就产生了一个单独的数据帧,centroids
,与y
列x
,和class
其中x
和y
是由类的平均值。然后我们使用centroid
作为数据集添加第二个点几何图层。
这是一个稍微有趣的版本,可用于聚类分析。
gg <- merge(df,aggregate(cbind(mean.x=x,mean.y=y)~class,df,mean),by="class")
ggplot(gg, aes(x,y,color=factor(class)))+geom_point(size=3)+
geom_point(aes(x=mean.x,y=mean.y),size=5)+
geom_segment(aes(x=mean.x, y=mean.y, xend=x, yend=y))
编辑回应OP的评论。
使用geom_errorbar(...)
和geom_errorbarh(...)
可以添加垂直和水平误差条。
centroids <- aggregate(cbind(x,y)~class,df,mean)
f <- function(z)sd(z)/sqrt(length(z)) # function to calculate std.err
se <- aggregate(cbind(se.x=x,se.y=y)~class,df,f)
centroids <- merge(centroids,se, by="class") # add std.err column to centroids
ggplot(gg, aes(x,y,color=factor(class)))+
geom_point(size=3)+
geom_point(data=centroids, size=5)+
geom_errorbar(data=centroids,aes(ymin=y-se.y,ymax=y+se.y),width=0.1)+
geom_errorbarh(data=centroids,aes(xmin=x-se.x,xmax=x+se.x),height=0.1)
如果你想计算,比方说,95%的置信度,而不是性病。错误,与
f <- function(z) qt(0.025,df=length(z)-1, lower.tail=F)* sd(z)/sqrt(length(z))
这是伟大的替代
,谢谢。是否也可以将水平和垂直条添加到表示x和y值的标准误差的质心? – cyril
请参阅最后的编辑。 – jlhoward
谢谢!太棒了。 – cyril