受启发与某群的讨论,今天来玩一下用ggplot画个小椭圆
ggplot画置信椭圆
我们之前介绍了置信椭圆的画法,参见置信椭圆与R画法,我们介绍了利用ggplot来画椭圆,
我们如果对于二维的聚类的散点图添加置信椭圆往往用的是stat_ellipse()这个参数,
我们看看它的帮助文档:
我们尝试这对一幅图添加95%置信区间:
ggplot(faithful, aes(waiting, eruptions, color = eruptions > 3)) +
geom_point() +
stat_ellipse(type = "norm",level = 0.95)
## level = 0.95为95%置信区间
## type代表符合的分布
我们还可以将点图换成椭圆的:
library(ggplot2)
ggplot(faithful, aes(waiting, eruptions, color = eruptions > 3)) +
geom_point(shape = 48,size = 5) +
stat_ellipse(type = "norm", linetype = 2)
ggforce画椭圆
还有一个画椭圆的包ggforce,里面的geom_ellipse()可以画,但是不太好表示置信度:
library(ggforce)
ggplot() +
geom_ellipse(aes(x0 = 0, y0 = 0, a = 10, b = 3, angle = 0)) +
coord_fixed()