多変量指数乱数 in R
無事、この関数でできました…。ほっ。
とりあえず、ちゃんと動いていそうです。
1番と4番目の変数の間に0.7の相関を仮定しました。
平均も大体1になってる。
test.r<-rmvexp(1000,c(1,1,1,1),matrix(c(1,0,0,0.7, 0,1,0,0, 0,0,1,0, 0.7,0,0,1),ncol=4))
apply(test.r,2,mean)
pairs.panels(test.r)
ちなみに、lcmixのインストールはこちら。CRANになくて、win var.もないみたい。
install.packages("lcmix", repos="http://R-Forge.R-project.org")