开发者

Can one extract model fit parameters after a ggplot stat_smooth call?

开发者 https://www.devze.com 2022-12-27 21:46 出处:网络
Using stat_smooth, I can fit models to data. E.g. g=ggplot(tips,aes(x=tip,y=as.numeric(unclass(factor(tips$sex))-1))) +facet_grid(time~.)

Using stat_smooth, I can fit models to data. E.g.

g=ggplot(tips,aes(x=tip,y=as.numeric(unclass(factor(tips$sex))-1))) +facet_grid(time~.) 
g=g+ stat_s开发者_如何学JAVAummary(fun.y=mean,geom="point") 
g=g+ stat_smooth(method="glm", family="binomial")

I would like to know the coefficients of the glm binomial fits. I could re-do the fit with dlply and get the coefficients with ldply, but I'd like to avoid such duplication.

Calling str(g) reveals the hierarchy of objects that ggplot2 creates, perhaps there's some way to get to the coefficients through that?


No, because the models are only created when the plot is rendered. However, it's usually pretty easy to do it yourself with plyr.

Why do you want to convert sex to a number? Using as.numeric should be enough by itself, but if you're going to do the subtraction in the model you'll need to surround it with I().

0

精彩评论

暂无评论...
验证码 换一张
取 消