利用R绘制置信椭圆(confidence ellipse)
有时候我们需要查看两组数据的分布情况,这时候通过可以绘制置信椭圆达到目的。
library(car)
dataEllipse(Duncan$income, Duncan$education, levels=0.95,xlim=c(-30,120),ylim=c(-30,160))
上述代码可以很快速地画出income和education两列数据的95%置信椭圆。结果如下:
把最后一行代码赋值给一个变量,就可以很方便地获得置信椭圆的边界: