R：（予測）ゼロ膨張したモデルには、関数zeroinfl（）</p> <p>を使用して設定されている<strong>ゼロ膨張したモデル</strong>のために、私は<strong>信頼区間</strong>を計算しようとしていますse.fit

Question

を返しません。私は機能

predict(glm, newdata, type = "response", se.fit = TRUE) 

を用いた線形モデルまたはGLMからそれらを計算する場合、それは正常に動作し、$フィット、$ se.fit、$ DFと$ residual.scaleの値を返します（私がわかりません最後の2つの出力があります - それについての説明も非常に高く評価されています）。

しかし、glmをゼロ膨張モデルで置き換えると、$ fitの値しか得られません。ここで

例： （私はset.seedを使用する方法がわからない - そのために私を責めないで）誰が間違ってここに何が起こるのか知っていますし、私はどのように計算することができ

set.seed(123) 
a <- data.frame(participant = c(1:10), 
       activity = c(round(abs(rnorm(10)))), 
       METmin = c(round(abs(rnorm(10)))), 
       Var4 = c(round(abs(rnorm(10)))), 
       Var5 = c(round(abs(rnorm(10))))) 

b <- data.frame(participant = c(1:10), 
        activity = c(round(abs(rnorm(10)))), 
        METmin = c(round(abs(rnorm(10)))), 
        Var4 = c(round(abs(rnorm(10)))), 
        Var5 = c(round(abs(rnorm(10))))) 
identical(a,b) #FALSE 
model <- zeroinfl(METmin~activity + Var4 + Var5, data = a) 
predict(model, newdata = b, type = "response", se.fit = T) 

これらの標準エラー？

ありがとうございました。

Answer 1

残念ながら、zeroinfl（およびhurdle）オブジェクトのpredict()メソッドは、現時点ではse.fit引数を持っていません。これを調べたい場合は、モンテカルロ法（係数見積りの分布から描画）を使用するか、lsmeansのzeroinflメソッドを目的に合わせて活用できるかどうかを確認する必要があります。