2010-04-19 27 views
2

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

+0

P.S.除了上面的as.numeric等方法之外,还有更好的方法来改变性别吗? – 2010-04-19 11:03:41

回答

2

不,因为模型只在绘制情节时创建。但是,使用plyr自己做它通常很容易。

为什么你想要将性转换为数字? as.nummeric本身就够了,但如果你要在模型中进行减法,则需要用I()来包围它。

+2

有没有机会看到plyr如何完成这项工作?我有同样的问题,但不知道从哪里开始。谢谢! – Atticus29 2012-07-23 03:20:02

相关问题