0

正定値でない共分散行列によって特異でない正規ランダム ベクトルを作成できるのはなぜですか? たとえば、実行できませんchol(V)。 

pigeon<-data.frame(response=c(10,19,27,28,9,13,25,29,4,10,20,18,5,6,12,17),
               treatment=factor(rep(1:4,4)),
               subject=factor(rep(1:4,each=4))
               )
m<-dcast(pigeon,subject~treatment, value.var="response")
fit<-lm(as.matrix(m[,-1])~1)
V<-cov(residuals(fit))
eigen(V)$values
rmvnorm(mean=rep(0,4),sigma=V,n=subject,method="chol")
4

2 に答える 2

1

MASS::rmvnorm を次のように実行する場合

rmvnorm(mu=rep(0,4),V=V,nsim=4,method="chol")、

私は得る

chol.default(V) のエラー: 次数 4 の先行マイナーが正定値ではありません

、chol(V) を実行するのと同じです。R 2.15.3 を使用しています。質問を明確にしていただけますか？

于 2013-06-27T15:57:21.347 に答える