可以在 ggplot stat_smooth 调用后提取模型拟合参数吗？

Question

使用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 创建的对象的层次结构，也许有一些方法可以通过它获得系数？

Answer 1

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

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