2010-04-19 77 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

相關問題