我想創建一個具有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')

總結

以上是生活随笔為你收集整理的r语言dataellipse_r – 在ggplot2中微调stat_ellipse()的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。