[R] Using the effects package to plot logit probabilities

[Bert Gunter]

Why not use

?predict.glm ## with type = "response"  ?

-- Bert

On Mon, Aug 13, 2012 at 12:39 PM, Abraham Mathew <abmathewks at 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 at 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