ベクトルがあるとしましょう:
Q<-rnorm(50,mean=0,sd=30)
Q+error term
10 個の変数のそれぞれが相関誤差項を持ち、10 個の変数間の相関が 0.8 などになるような 10 個の変数を作成したいと思います。
通常、次を使用して相関変数を作成します。
vars<-rmvnorm(n=10, mean=c(rep(0,50)),cov.mat, method="chol")
を追加するとどうなりますかQ+vars
。探しているものを手に入れることができますか?