重新縮放斜率
上一個圖顯示各郡的犯罪率變化幅度不同。這表示你應在模型中同時把 Year 納入為隨機效應與固定效應。以這種方式加入 Year,會同時估計跨所有郡的整體斜率,以及各郡各自的斜率。固定效應的斜率估計所有馬里蘭州郡的重大犯罪變化;隨機效應的斜率則讓模型允許各郡有不同的變化趨勢。
但是,擬合這個模型時會出現警告訊息!為了處理這個警告,請將 Year 從以 2006 起算改為以 0 起算。我們已提供這個新變數 Year2(例如,Year 中的 2006 在 Year2 中是 0)。有時在擬合迴歸時,你需要將截距縮放或中心化到從 0 開始。這能提升模型的數值穩定性。
本練習屬於課程
R 的階層式與混合效果模型
練習說明
- 建立一個
lmer(),以Crime為應變數,並將Year同時作為固定效應與隨機效應的斜率,County作為隨機效應的截距。 - 再建立第二個
lmer(),以Crime為應變數,並將Year2同時作為固定效應與隨機效應的斜率,County作為隨機效應的截距。
動手互動練習
試著完成這個範例程式碼,體驗一下這個練習。
# Fit the model with Year as both a fixed and random-effect
lmer(___ ~ Year + (1 + Year | ___) , data = md_crime)
# Fit the model with Year2 rather than Year
lmer(___ ~ Year2 + (1 + Year2 | ___) , data = md_crime)