Построение прогнозируемых эффектов: извлечение верхних и нижних границ из команды пакета {Effects}

Question

У меня есть модель логистической регрессии, и я пытался получить график предсказанных вероятностей для дихотомического предиктора x1 (что означает "пол"). Теперь есть простой способ, используя пакет effects:

dat<-data.frame(
  y=round(runif(100,0,1)),
  x1=as.factor(round(runif(100,0,1))),
  x2=round(runif(100,20,80))
)

model<-glm(data=dat,
           formula=y~x1+x2, 
           family=binomial(link="logit"))

require(effects)
plot(Effect("x1",model))

Возвращает:

Но так как это не очень причудливо, и было бы неплохо сделать это в ggplot: есть ли способ извлечь верхнюю и нижнюю границы из Effect("x1",model)? Команда сама по себе только возвращает предсказанный эффект человека, являющегося мужчиной или женщиной, но не дает верхних и нижних границ на 95%. И нет никакого упоминания в ?effects.

Мне нужна такая рамка, чтобы я мог позволить ggplot сделать работу:

x1   effect   lower   upper
1    0.39     .       .
0    0.56     .       .

Answer 1

Хорошо, я не ожидал, что это будет так легко. Очевидно, достаточно написать summary() вокруг команды эффектов.

summary(Effect("x1",model))


#x1 effect
x1
        0         1 
0.5572423 0.3984501 

#Lower 95 Percent Confidence Limits
x1
        0         1 
0.4111167 0.2782799 

#Upper 95 Percent Confidence Limits
x1
        0         1 
0.6940889 0.5322441