使用stat_smooth,我可以将模型拟合成数据。例如。ggplot 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)显示了ggplot创建的对象的层次结构,也许有一些方法可以通过它获得系数?
P.S.除了上面的as.numeric等方法之外,还有更好的方法来改变性别吗? – 2010-04-19 11:03:41