利用R画置信椭圆

利用R绘制置信椭圆(confidence ellipse)

有时候我们需要查看两组数据的分布情况，这时候通过可以绘制置信椭圆达到目的。

library(car)
dataEllipse(Duncan$income, Duncan$education, levels=0.95,xlim=c(-30,120),ylim=c(-30,160))

上述代码可以很快速地画出income和education两列数据的95%置信椭圆。结果如下：

把最后一行代码赋值给一个变量，就可以很方便地获得置信椭圆的边界：

boundary <- dataEllipse(Duncan$income, Duncan$education, levels=0.95,xlim=c(-30,120),ylim=c(-30,160))
head(boundary)

更多功能自己通过“?dataEllipse”去探索吧。