1

confint.default漸近正規性に基づくlmerモデルが機能しないのはなぜですか?

fit <- lmer(y~(1|operator)+(1|part),data=dat)

Linear mixed model fit by REML ['lmerMod']
Formula: y ~ (1 | operator) + (1 | part)
 Data: dat
REML criterion at convergence: 409.3913
Random effects:
Groups   Name        Std.Dev.
part     (Intercept) 3.2018  
operator (Intercept) 0.1031  
Residual             0.9398  
Number of obs: 120, groups:  part, 20; operator, 3
Fixed Effects:
(Intercept)  
     22.39  


 confint.default(fit)
 Error in as.integer(x) : 
   cannot coerce type 'S4' to vector of type 'integer'

エラーは何を言っていますか? lmer モデルの漸近正規性に基づいて信頼区間を取得するにはどうすればよいですか?

4

2 に答える 2

4

を使用しないconfint.default()で、そのまま使用してくださいconfint()。信頼区間を計算する方法は、モデル タイプによって異なります。でさまざまな方法を確認できますmethods(confint)。関数の「正しい」バージョンは、関数に渡す最初のオブジェクトのクラスに基づいて呼び出されます。通常、メソッドの 1 つを直接呼び出すことはお勧めできません。

オブジェクトの境界を計算する方法にはオプションがありlmerます。?confint.merModそのモデル タイプに固有のオプションについては、 のヘルプ ページを参照してください。

于 2015-08-01T01:12:56.300 に答える
2

confint.merMod()@MrFlick は正しいですが、デフォルトで尤度プロファイル CIを提供することを追加する価値があるかもしれません。confint(.,method="Wald")漸近正規性に基づく信頼区間が得られます。

「Wald」: 推定された尤度面の局所曲率に基づいて、信頼区間を近似します (固定効果パラメーターのみ。すべての分散共分散パラメーター CI は「NA」として返されます)。

(これはヘルプ ページから明らかですが、おそらくここでもう一度説明する価値があります)。

于 2015-08-01T01:44:28.257 に答える