GLMM をセットアップしようとしていますが、R で階層データ構造を表現するのに問題がありますglmer。私のデータは次の構造を持っています：
y (dependent variable);
visit (L1: nested independent variable);
year (L2: nested independent variable);
site (L3: independent variable);
wherevisitは 内にネストされyear、これは 内にネストされますsite(つまり、1 年間に最大 5 回の訪問があり、サイトごとに最大 10 年間のデータがある可能性があります)。
私は on の効果に興味がありますがyear、ランダムな効果としてy扱いたいと思ってsiteいます (つまり、切片と勾配 [年間] はサイトによって異なります)。
次のコマンドを使用しました。
model <- glmer(y ~ year + (year | site) + (visit|year),
family=poisson, data=data)
モデルを正しく表現しましたか?