With multiple predictors, I feel that predict ends up being less elegant.
Plus, I've found a lot of the
simple log-odds plots available in the effect package to be simple and easy
to use as compared to
base predict.
On Mon, Aug 13, 2012 at 3:01 PM, Bert Gunter <gunter.ber...@gene.com> wrote:
> Why not use
>
> ?predict.glm ## with type = "response" ?
>
> -- Bert
>
> On Mon, Aug 13, 2012 at 12:39 PM, Abraham Mathew <abmathe...@gmail.com>
> wrote:
> > I'm trying to run a logit model and plot the probability curve for a
> number
> > of the important predictors. I'm trying to do this
> > with the Effects package.
> >
> >
> > df=data.frame(income=c(5,5,3,3,6,5),
> > won=c(0,0,1,1,1,0),
> > age=c(18,18,23,50,19,39),
> > home=c(0,0,1,0,0,1))
> > str(df)
> >
> > md1 = glm(factor(won) ~ income + age + home,
> > data=df, family=binomial(link="logit"))
> >
> > summary(md1)
> >
> > plot(effect("income", md1), grid=TRUE)
> >
> >
> > But I want to know how to plot a graph so that it shows the probability
> of
> > won (response) based on income (or any of the other predictors).
> >
> > However, what I want to do is generate the same plot, with won don y axis
> > and income on x axis, but the curves showing the probabilities for age
> and
> > home.
> >
> >
> > Not seeing how to do this in the effects documentation. Help!
> >
> >
> > Thanks.
> >
> >
> >
> > --
> > *Abraham Mathew
> > Statistical Analyst
> > www.amathew.com
> > 720-648-0108
> > @abmathewks*
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>
>
>
> --
>
> Bert Gunter
> Genentech Nonclinical Biostatistics
>
> Internal Contact Info:
> Phone: 467-7374
> Website:
>
> http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm
>
--
*Abraham Mathew
Statistical Analyst
www.amathew.com
720-648-0108
@abmathewks*
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
