我正在尝试在R中复制具有随机块和重复测量的线性混合模型的SAS模型。我使用SAS的教授在SAS中为此测试编写了语法。我对SAS的了解很差，因此我在R中的尝试。我在将重复测量语法合并到我的R模型中时遇到了麻烦。请指教。
我正在研究动物丰度（模式物种）与三年内在27个网格中收集的许多栖息地变量（如灌木密度，水分，木材）之间的关系。在这三年中，夏季和春季（季节）收集了丰富的数据。
随机效果=网格
固定效果=木材，水分，灌木密度
重复措施=年份和季节
这是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中相同数据集相同的结果。