我正在尝试使用lme包中的重复测量方差分析nlme来复制结果。aov我已经为一个单因素重复测量实验和一个具有一个主体间因素和一个主体内因素的双因素实验做了这个，但是对于一个有两个内部因素的双因素实验，我很难做到这一点——主体因素。
一个例子如下所示。是固定效应因子，A是随机效应因子。Bsubject
set.seed(1)
d <- data.frame(
Y = rnorm(48),
subject = factor(rep(1:12, 4)),
A = factor(rep(1:2, each=24)),
B = factor(rep(rep(1:2, each=12), 2)))
summary(aov(Y ~ A*B + Error(subject/(A*B)), data=d)) # Standard repeated measures ANOVA
library(nlme)
# Attempts:
anova(lme(Y ~ A*B, data=d, random = ~ 1 | subject)) # not same as above
anova(lme(Y ~ A*B, data=d, random = ~ 1 | subject/(A+B))) # gives error
我在皮涅罗和贝茨的书中没有看到对此的解释，但我可能忽略了它。