使用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