我遇到了像这样的SPSS语法
MIXED value BY factor1
/CRITERIA=CIN(95) MXITER(100) MXSTEP(10) SCORING(1) SINGULAR(0.000000000001)
HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE) PCONVERGE(0.000001,
ABSOLUTE)
/FIXED=factor1 | SSTYPE(3)
/METHOD=REML
/REPEATED=factor1 | SUBJECT(participant) COVTYPE(UN).
和斗争找到一个等效lmer/ nlme（或R一般）制剂对于这种模型。
有人知道如何将REPEATED子命令转换为R代码吗？
这概括了答案我在r-SIG-混合模型邮件列表得到：
该REPEATED命令指定要在剩余方差-协方差矩阵（结构- [R矩阵），即所谓的R侧的结构，的模型。对于lme4::lmer()这种结构，固定为单位矩阵的倍数。但是，可以使用和参数指定R端结构，如下所示：weightscorrelationnlme::gls()
gls(value ~ factor1,
correlation = corSymm(form = ~ 1|participant),
weights = varIdent(form = ~1|factor1),
method = "REML",
data = data)
如果除了R侧结构还需要G侧效果，请nlme::lme()提供适当的扩展。
本文收集自互联网，转载请注明来源。
如有侵权，请联系 [email protected] 删除。
我来说两句