使用stat_smooth,我可以将模型拟合到数据中。 E.g。
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)揭示了ggplot创建的对象的层次结构,也许有一些方法可以通过它来获取系数?
答案 0 :(得分:3)
不,因为模型仅在渲染绘图时创建。然而,通常很容易用plyr自己做。
为什么要将性别转换为数字? As.numeric本身应该足够了,但是如果你要在模型中进行减法,你需要用I()包围它。