Dear R experts,
I am trying to use John Fox's effects package to conduct some
post-estimation analysis.? I think it is a great package for graphing
predictions after GLM type models.? However, I am having problems
figuring out how to produce an adjusted mean.? For example, using the
following two lines, I can get predictions when phd varies from 1 to
5, female set to 0 and enrol set to its mean, and then use the plot to
get predicted probability plot.
logitmod <- glm(hijob~female+enrol+phd, data=mydta,
family=binomial(lin="logit"))
plot(effect("phd", logitmod, xlevels=list(phd=1:5),
given.values=c(female=0)))
But I couldn't figure out how to get just one predicted probability
with a given set of values for independent variables (for example,
female=1, phd=3 and enrol = 7), and it appears that using the effect
command, I have to get multiple predictions?? Thanks a lot!
Jun Xu, PhD
Assistant Professor
Department of Sociology
Ball State University
Muncie, IN 47306
