使用回归模型(LM,LMER)重复测量方差分析

Posted

技术标签:

【中文标题】使用回归模型(LM,LMER)重复测量方差分析【英文标题】:repeated measure anova using regression models (LM, LMER) 【发布时间】:2013-02-14 10:24:51 【问题描述】:

我想使用回归模型而不是“方差分析”(AOV) 函数在 R 中运行重复测量方差分析。

这是我的 3 个主体内因素的 AOV 代码示例:

m.aov<-aov(measure~(task*region*actiontype) + Error(subject/(task*region*actiontype)),data)  

谁能告诉我使用回归模型运行相同分析的确切语法?我想确保尊重残差的独立性,即使用与 AOV 一样的特定错误术语。

在上一篇文章中,我读到了一个类型的答案:

lmer(DV ~ 1 + IV1*IV2*IV3 + (IV1*IV2*IV3|Subject), dataset))

我真的不确定这个解决方案,因为它仍然将变量视为对象之间的变量,而且我不明白添加随机因素会如何改变这一点。

有人知道如何在考虑剩余独立性的情况下使用 lm/lmer 运行重复测量方差分析吗?

非常感谢, 索伦

【问题讨论】:

【参考方案1】:

如果你的 aov 例子是正确的(也许你不想嵌套东西)你想要这个:

lmer(measure~(task*region*actiontype) + 1(1|subject/(task:region:actiontype))

如果残差独立意味着截距和斜率独立计算,您需要分别指定它们:

+(1|yourfactors)+(0+variable|yourfactors)

或使用符号:

+(1||yourfactors)

无论如何,如果您阅读帮助文件,您会发现 lme4 无法处理最常见的问题。

【讨论】:

以上是关于使用回归模型(LM,LMER)重复测量方差分析的主要内容,如果未能解决你的问题,请参考以下文章

R语言使用回归方法解决方差分析问题

多元线性回归模型用r语言怎么来实现

在线性回归分析中,若检验的结果为不显著,可能原因是啥

生存/回归分析结果的最佳/有效绘图

R语言使用lmPerm包应用于线性模型的置换方法(置换检验permutation tests)使用lm模型构建简单线性回归模型使用lmp函数生成置换检验回归分析模型

R使用lm构建单变量线性回归模型