針對模型估算值繪製固定效果？

所以我看了一些類似的主題帖子，但沒有一個看起來完全是我需要的，或者我根本不明白他們提供的解決方案......所以這裏就去...針對模型估算值繪製固定效果？

我用lme4運行混合效果模型來查看一些黑猩猩數據。我有兩個因素（攻擊率;交配率）影響我的依賴（餵養時間）。

我想生成兩個散點圖，顯示每個預測變量和結果變量之間的關係，但我想繪製一條線，該線來自模型估計值（而不是（lm （y〜x））類型，它只給出一條簡單的迴歸線，而不是基於完整LMM的迴歸線）。

我有一種感覺，這是唯一可能的ggplot2，但我一直沒有能夠真正弄清楚如何做到這一點。花了大部分時間閱讀書籍和論壇，我希望這是一個可以有一個相當直接的答案，如果知道他們在做什麼。

感謝您提前提出任何提示！

亞歷

2014-03-31 Alex Georgiev

如果您[讓您的問題可重現]，您將獲得更多幫助（http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example） – alexwhan

'？預測.merMod'可能會幫助你。 –

非常感謝，@BenBolker - 這工作！亞歷克斯 –

要與I開始具有下列模式：

M3reml 
Linear mixed model fit by REML ['lmerMod'] 
Formula: z.feeding_time ~ z.copul_rate + z.agro_given + z.agro_recd + (1 | Male) + ac_term 
    Data: N85

其中變量的z變換的值：雄性黑猩猩餵食時間（z.feeding_time）;與女性的每日交配率（act/hr; z.copul_rate）;每日侵略率（z.agro_given）;和每日侵略率（z.agro_recd）。隨機效應 - 我研究的12名男性的男性ID;和時間自相關術語（ac_term）。

我想根據男性餵養時間的模型估計值產生迴歸線。

獲取估計：

p1<-predict(M3reml)

密謀反對侵略的雄性率（Z變換值）的估計：

plot(p1~z.agro_given, data=N85)

添加回歸線：

abline(lm(p1~z.agro_given, data=N85))

我會在這裏發佈劇情的圖片，但顯然我還沒有被允許。

來源

2014-04-01 18:42:24

不知道爲什麼文本的格式看起來非常殘暴......雖然之前沒有明確這麼做過。 –

啊。感謝編輯這個......我會試着弄清楚下一次這是如何工作的！ –

針對模型估算值繪製固定效果？

回答

相關問題