0

わからない素朴な疑問。相互作用を含むモデルから条件付き予測確率を生成しようとしています。たとえば、x2==1 および x3==0 の場合の予測確率と、x2==0 および x3==1 の場合の予測確率を比較できるようにしたいと考えました。

私は次のようにこれをやろうとしています:

model <- glm(y~x1 + x2 * x3, family=binomial(link="logit"), data=data)
predprob1 <- predict(model, type="response", newdata=(x1=mean(x1) & x2==1 & x3==0))
predprob2 <- predict(model, type="response", newdata=(x1=mean(x1) & x2==0 & x3==1))
probdiff<-predprob1-predprob2

その後、probdiff の 95CI を計算する必要があります。これは、R の天才にとっては簡単なことだと思います。ご協力ありがとうございました!

4

1 に答える 1

0

おそらく(実際の例がないためテストされていません):

model <- glm(y~x1 + x2 * x3, family=binomial(link="logit"), data=data)
predprob1 <- predict(model, type="response", 
                     newdata=with(data, data.frame(x1=mean(x1) & x2=1 & x3=0))
predprob2 <- predict(model, type="response", 
                     newdata=with(data, data.frame(x1=mean(x1) & x2=0 & x3=1))
probdiff <- predprob1-predprob2

もう 1 つの R エラーは、代入に「==」を使用することでした。論理テスト用です。

于 2012-08-06T18:33:15.083 に答える