わからない素朴な疑問。相互作用を含むモデルから条件付き予測確率を生成しようとしています。たとえば、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 の天才にとっては簡単なことだと思います。ご協力ありがとうございました!