Question

您好我在R中运行混合模型相对较新。但我有在SAS Proc Mixed和SPSS Mixed中使用它们的经验。我有以下代码，我想帮助翻译在R中用于lmer：

PROC MIXED COVTEST; 
CLASS ACTOR PARTNER DYAD GROUP; 
MODEL LEAD = /S DDFM=SATTERTH NOTEST; 
RANDOM INTERCEPT /TYPE=VC SUB=ACTOR; 
RANDOM INTERCEPT /TYPE=VC SUB=PARTNER; 
RANDOM INTERCEPT / TYPE=VC SUB=GROUP; 
REPEATED /TYPE=CS SUB=DYAD;

或者，SPSS中的代码是：

MIXED
LEAD BY GROUP
/FIXED =
/PRINT = SOLUTION TESTCOV
/RANDOM INTERCEPT | SUBJECT(GROUP) COVTYPE(VC) .
/RANDOM INTERCEPT | SUBJECT(ACTOR) COVTYPE(VC)
/RANDOM INTERCEPT | SUBJECT(PARTNER) COVTYPE(VC)
/RANDOM INTERCEPT | SUBJECT(DYAD) COVTYPE(VC) .

Answer 1

我认为这更像是Stack Overflow的主题，但无论如何我都会回答它。看起来你正试图适应Kenny's Social Relations Model，我有一点经验。

简短的回答是，这目前无法在lme4中完成。原因是SAS脚本涉及＆＃34; R侧协方差结构＆＃34; （即，它调用REPEATED语句），但lme4目前不支持这些语句。作为替代方案，您可以查看TripleR包以确定它是否符合您的需求。

作为最终评论，我不认为您的SAS和SPSS模型是等效的。 SAS模型（如果我记得所有语法都正确）使用REPEATED语句来允许残差方差在二元组之间变化。但SPSS模型反而具有随机二元效应（通过RANDOM声明），我认为这在统计上是等价的。我不是100％肯定，因为有时可以使用RANDOM或REPEATED语句获得相同的结果。

在R中将SAS语法转换为lmer

1 个答案: