r语言dataellipse_r – 在ggplot2中微调stat_ellipse()

我想创建一个具有95％“精确”置信椭圆的二元正态分布的散点图.

library(mvtnorm)

library(ggplot2)

set.seed(1)

n

c95

rho

Sigma

我从双变量法线生成了1000个观测值,平均值为零,方差=西格玛

x

z

for(i in 1:n){

z[i]

p95[i]

}

我们可以使用stat_ellipse轻松地在生成数据的散点图顶部绘制95％置信度椭圆.得到的数字是完全令人满意的,直到你注意到几个红点位于置信椭圆内.我猜这种差异来自某些参数的估计,并随着样本量的增大而消失.

data

p

p + stat_ellipse(type = "norm")

有没有办法微调stat_ellipse(),以便它描绘“精确”置信椭圆,如下图所示,是使用“手工制作”椭圆函数创建的？

ellips

t

Sigma

a

b

x

y

X

R

data.frame(X%*%R)

}

dat

p + geom_path(data=dat, aes(x=X1, y=X2), colour='blue')