[R]

2017-12-01 5 views
0

データは508校から33988回の観測を持っている5.[R]

を私はhttps://www.cmm.bris.ac.uk/lemma/login/index.phpから混合モデルを学んでいますし、この質問は、モジュール内の運動に関するれます。学生データは第1レベルであり、学校は第2レベルです。私はヌルモデル

nullmodel <- lmer(score ~ (1|schoolid), data = mydata, REML = FALSE) 

(私はデータセットへのオープンアクセスを見つけることができませんでした。質問は、より概念的であるので、私は推測する、それは問題ではありません)

私たちが実行してもらう切片推定値は万一の実行しますヌルモデルのサマリーは「スコア」の平均と等しくないでしょうか?私はインターセプトの推定値を30.60、平均値[平均値(mydata $ score)]を31.09とします。値は近いが同じではない。

また、学校レベルの残差u0を得る直接的な方法はありますか?

答えて

0

ランダムインターセプト効果u0は、平均0と分散σ_u0^ 2で正規分布します。

σ_u0^ 2は学校間の差です。

σ_u0^ 2は直接的な方法で計算できます。

  1. 分散がσ_u0^ 2に近いこれらの値の
  2. 計算分散

各学校の平均スコアを計算します。

+0

REMLの推定値は、ランダムな影響に対して不偏であるため、最終的なモデル適合のためにREMLメソッドを推奨します。 –