r - 在 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() 将其包围。