我试图用R中的随机块和重复测量的线性混合模型复制SAS模型 . 我使用SAS的教授在SAS中为此测试编写了语法 . 我对SAS的了解很差，因此我在R中的尝试 . 我在将重复测量语法合并到我的R模型中时遇到了麻烦 . 请指教 .
我正在研究动物丰度（模式物种）与三年内在27个网格中收集的许多栖息地变量（例如灌木密度，水分，木材）的关联 . 在这三年中，夏季和 Spring 季（季节）收集了丰富的数据 .
随机效应=网格
固定效果=木材，水分，灌木密度
重复措施=年份和季节
这是SAS代码：
proc mixed data=SHREWDATA ic;
class grid year season treatment;
model species= year season shrubdensity moisture &wood &wood*year &wood*season/htype = 3 solution DDFM=SATTERTHWAITE ;
random grid/solution;
repeated /subject=grid(treatment);
lsmeans year season; run;
每年和季节重复的“单位”（主题）是网格内的治疗 .
这是我使用lme尝试的R代码
modelspecies2 <- lme(modelspecies ~ year + season + shrubdensity + moisture + wood + wood * year + wood * season, random = list(~1|grid), data = SHREWDATA, method = "ML")
我没有看到任何有关重复测量和随机混合块的混合模型的讨论 . 我试图合并一个重复的度量语法，但是我得到一个错误，或者没有得到与SAS中相同数据集相同的结果 .