読者です 読者をやめる 読者になる 読者になる

多変量指数乱数 in R

rdrr.io

無事、この関数でできました…。ほっ。

とりあえず、ちゃんと動いていそうです。

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)

 

f:id:senedesmus:20170217161649p:plain

ちなみに、lcmixのインストールはこちら。CRANになくて、win var.もないみたい。

install.packages("lcmix", repos="http://R-Forge.R-project.org")