我想創(chuàng)建一個(gè)具有95％“精確”置信橢圓的二元正態(tài)分布的散點(diǎn)圖.

library(mvtnorm)

library(ggplot2)

set.seed(1)

n

c95

rho

Sigma

我從雙變量法線生成了1000個(gè)觀測(cè)值,平均值為零,方差=西格瑪

x

z

for(i in 1:n){

z[i]

p95[i]

}

我們可以使用stat_ellipse輕松地在生成數(shù)據(jù)的散點(diǎn)圖頂部繪制95％置信度橢圓.得到的數(shù)字是完全令人滿意的,直到你注意到幾個(gè)紅點(diǎn)位于置信橢圓內(nèi).我猜這種差異來(lái)自某些參數(shù)的估計(jì),并隨著樣本量的增大而消失.

data

p

p + stat_ellipse(type = "norm")

有沒(méi)有辦法微調(diào)stat_ellipse(),以便它描繪“精確”置信橢圓,如下圖所示,是使用“手工制作”橢圓函數(shù)創(chuàng)建的？

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')

總結(jié)

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

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。