以下のいくつかの係数のサイズが大きいことに疑問を持っている場合に備えて、私は金融/経済データを扱っています... 私の一般的な質問は、R の線形ランダム効果モデルから出力されるパラメーター係数のシミュレーションに関するものです。 R の同じモデルからモデル係数と分散共分散 (VCOV) 行列を使用して、ベータ係数のランダム サンプルを生成しようとしています。 mvtnorm{} パッケージの rmvnorm() 関数? この警告/問題に対処するにはどうすればよいですか?
#Example call: lmer model with random effects by YEAR
#mlm<-lmer(DV~V1+V2+V3+V2*V3+V4+V5+V6+V7+V8+V9+V10+V11+(1|YEAR), data=dat)
#Note: 5 years (5 random effects total)
#LMER call yields the following information:
coef<-as.matrix(c(-28037800,0.8368619,2816347,8681918,-414002.6,371010.7,-26580.84,80.17909,271.417,-239.1172,3.463785,-828326))
sigma<-as.matrix(rbind(c(1834279134971.21,-415.95,-114036304870.57,-162630699769.14,-23984428143.44,-94539802675.96,
-4666823087.67,-93751.98,1735816.34,-1592542.75,3618.67,14526547722.87),
c(-415.95,0.00,41.69,94.17,-8.94,-22.11,-0.55,0.00,0.00,0.00,0.00,-7.97),
c(-114036304870.57,41.69,12186704885.94,12656728536.44,-227877587.40,-2267464778.61,
-4318868.82,8909.65,-355608.46,338303.72,-321.78,-1393244913.64),
c(-162630699769.14,94.17,12656728536.44,33599776473.37,542843422.84,4678344700.91,-27441015.29,
12106.86,-225140.89,246828.39,-593.79,-2445378925.66),
c(-23984428143.44,-8.94,-227877587.40,542843422.84,32114305557.09,-624207176.98,-23072090.09,
2051.16,51800.37,-49815.41,-163.76,2452174.23),
c(-94539802675.96,-22.11,-2267464778.61,4678344700.91,-624207176.98,603769409172.72,90275299.55,
9267.90,208538.76,-209180.69,-304.18,-7519167.05),
c(-4666823087.67,-0.55,-4318868.82,-27441015.29,-23072090.09,90275299.55,82486186.42,-100.73,
15112.56,-15119.40,-1.34,-2476672.62),
c(-93751.98,0.00,8909.65,12106.86,2051.16,9267.90,-100.73,2.54,8.73,-10.15,-0.01,-1507.62),
c(1735816.34,0.00,-355608.46,-225140.89,51800.37,208538.76,15112.56,8.73,527.85,-535.53,-0.01,21968.29),
c(-1592542.75,0.00,338303.72,246828.39,-49815.41,-209180.69,-15119.40,-10.15,-535.53,545.26,0.01,-23262.72),
c(3618.67,0.00,-321.78,-593.79,-163.76,-304.18,-1.34,-0.01,-0.01,0.01,0.01,42.90),
c(14526547722.87,-7.97,-1393244913.64,-2445378925.66,2452174.23,-7519167.05,-2476672.62,-1507.62,21968.29,
-23262.72,42.90,229188496.83)))
#Error begins here:
betas<-rmvnorm(n=1000, mean=coef, sigma=sigma)
#rmvnorm breaks, Error returned:
警告メッセージ: In sqrt(ev$values) : NaN が生成されました
次の検索文字列を Google で検索すると、「rmvnorm、「警告メッセージ: sqrt(ev$values) 内: NaNs が生成されました」と表示されました: http://www.nickfieller.staff.shef.ac.uk/sheff-only /mvatasksols6-9.pdf 4 ページで、このエラーは「負の固有値」を示していますが、概念的または実際的に、負の固有値が何であるか、またはこの場合に生成される理由がわかりません。
2 番目の検索結果: [ http://www.r-tutor.com/r-introduction/basic-data-types/complex 2 -1 の平方根を取得しようとしたためにこのエラーが発生したことを示します。 「複雑な値ではありません」(-1 の平方根を取ることはできません)。
ここで、ベータのランダム生成で何が起こっているのか、そしてこれをどのように修正できるのかという疑問が残ります。
sessionInfo() R バージョン 3.0.2 (2013-09-25) プラットフォーム: x86_64-apple-darwin10.8.0 (64 ビット)
以下のパッケージ/バージョンを使用 mvtnorm_0.9-9994、lme4_1.1-5、Rcpp_0.10.3、Matrix_1.1-2-2、lattice_0.20-23