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

	[[alternative HTML version deleted]]

