模型規格與glmer

Question

我正在嘗試建立一個GLMM，但我有問題表達R與glmer分層數據結構。我的數據結構如下：

y (dependent variable); 
visit (L1: nested independent variable); 
year (L2: nested independent variable); 
site (L3: independent variable); 

其中visit嵌套在year，這是嵌套在site（即有可能達到每年5人次，其中可能有長達10年的數據，每個站點）。

我感興趣的year在y的效果，同時也希望把site作爲隨機效應（即攔截和斜坡[同期]通過網站而異）。

我用下面的命令：

model <- glmer(y ~ year + (year | site) + (visit|year), 
       family=poisson, data=data) 

有我正確表達的模式？

Answer 1

我想我會建議

model <- glmer(y ~ year + 
        (1+year | site) + 
        (1|visit:site:year), 
       family=poisson, data=data)  

我建議你把一年的數字（連續）變量;否則(year|site)術語將嘗試估計每年的站點間變化以及逐年協同性（可能會失敗）。在此模型中，

• year指定相對於年的整體固定線性趨勢;
• 1+year|site指定站點間截距的變化（即總數），以及站點間坡度相對於年份的變化以及斜率截距相關性;
• 1|visit:site:year給出觀察水平的隨機效應，其允許過度分散。

還有其他可能的方式來切這個（你可以考慮一個(1|year)長期佔整體，整個現場周圍的線性趨勢varation相對於年）。

一般來說，我會小心(var|grp)方面，其中var是絕對的，因爲它們會導致估計相當複雜的模型，如果var有超過2或3級...