3

使用stat_smooth，我可以将模型拟合到数据中。例如

g=ggplot(tips,aes(x=tip,y=as.numeric(unclass(factor(tips$sex))-1))) +facet_grid(time~.) 
g=g+ stat_summary(fun.y=mean,geom="point") 
g=g+ stat_smooth(method="glm", family="binomial")

我想知道 glm 二项式拟合的系数。我可以用 重新拟合并dlply得到系数ldply，但我想避免这种重复。

调用str(g)揭示了ggplot2创建对象的层次结构，也许有一些方法可以通过它来获得系数？

4

1 回答 1

3

不，因为模型仅在绘制绘图时创建。但是，通常很容易自己使用plyr.

为什么要将性别转换为数字？usingas.numeric本身就足够了，但如果你要在模型中做减法，你需要用I().

于 2010-04-19T12:30:03.303 回答