# NOT RUN {
data(france)
left.mod <- glm(voteleft ~ male + age + retnat +
poly(lrself, 2, raw=TRUE), data=france, family=binomial)
data(france)
left.mod <- glm(voteleft ~ male + age + retnat +
poly(lrself, 2, raw=TRUE), data=france, family=binomial)
out <- probci(left.mod, france, numQuantVals=3,
changeX=c("retnat", "lrself"), calcPW = TRUE)
print(out, filter=list(retnat=c("Better", "Worse")))
print(out, type="pw",
filter=list(retnat=c("Better", "Worse")),
const="lrself")
# }
Run the code above in your browser using DataLab