lme（）的多个嵌套结构

Question

我想在R中构建一个具有多个嵌套结构的混合效果模型。例如，月份嵌套在一年中，一天嵌套在年份和年份中。响应变量Y是温度，预测变量X是时间。以下SAS代码可以完全满足我的需求。

proc mixed data=dataset;
 class year month day;
 model Y=X;
 random day(year) month(year) year;
 run;

但是，当我使用lme（）时，我不知道如何设置这个结构。它应该像

fitted.model <- lme(Y~X, random=list(~1|year/day, ~1|year/month, ~1|year), data=dataset)

以上命令不起作用。如果我写的是

fitted.model <- lme(Y~X, random=list(~1|year, ~1|month, ~1|day), data=dataset)

我将获得嵌套的随机效果年/月/日，这是我不想要的。