线性混合效应模型的报告结果


12

线性混合效应模型在我的生物学研究中并不常用,我需要报告我试图写的论文中使用的统计检验。我知道生物科学的某些领域(对依赖性的解决方案:使用多层次分析来容纳嵌套数据)已经开始出现对多层次建模的意识,但是我仍在尝试学习如何报告结果!

我的实验设计简而言之:
*将受试者分配到四个治疗组中的一个
*在开始治疗后的不同天对因变量进行测量
*设计不平衡(治疗组中受试者人数不等,并且缺少某些天的测量值)
*治疗A是参考类别
*我将数据集中在治疗的最后一天

我想知道治疗A(参考类别)是否比其他治疗(治疗结束)产生明显更好的结果。

我使用nlme在R中进行了分析:

mymodel <- lme(dv ~ Treatment*Day, random = ~1|Subject, data = mydf, na.action = na.omit, 
+ correlation = corAR1(form = ~1 |Subject), method = "REML")

输出(部分内容;为了简洁起见被截断)为:

>anova(mymodel)
              numDF denDF  F-value p-value
(Intercept)      1   222 36173.09  <.0001
Treat            3    35    16.61  <.0001
Day              7   222     3.43  0.0016
Treat:Day       21   222     3.62  <.0001

>summary(mymodel)
Fixed effects: dv ~ Treatment * Day 
                       Value Std.Error  DF  t-value p-value
(Intercept)         7.038028 0.1245901 222 56.48945  0.0000
TreatmentB          0.440560 0.1608452  35  2.73903  0.0096
TreatmentC          0.510214 0.1761970  35  2.89570  0.0065
TreatmentD          0.106202 0.1637436  35  0.64859  0.5208

因此,我知道Day的效果因治疗而异,并且在治疗的最后一天(数据居中),治疗A中的dv与治疗B或C显着不同。

我想说的是:“正如预测的那样,我们发现接受治疗A的受试者(平均+/- SE)的因变量显着低于接受治疗B(平均+/- SE,p = 0.0096)或治疗的受试者的因变量。 C(平均+/- SE,p = 0.0065),在治疗的最后一天测得。”

但是,我必须指出进行了哪些统计检验。这是描述分析的可接受方法吗?“ [测量方法]在所示的日期进行,并确定了因变量(​​单位);我们使用了以[治疗的最后一天]为中心的线性混合效应模型来分析对数转换后的数据。符号代表平均dv;误差线是标准误差。在治疗的最后一天,治疗A的dv(平均+/- SE)明显低于治疗B的dv(平均+/- SE,p = 0.0096)...”

具体来说,
*这足以说明所使用的统计检验吗?(读者习惯于看到更多类似“平均值+/- SE,p = 0.0096,学生t检验”的方法,但似乎很奇怪地写出“ p = 0.0096,线性混合效应对治疗B与治疗A的系数” [治疗的最后一天]建立模型。“)
*是否有更好的方法来表达这一点?

(“方法”部分将提供有关统计信息的更多信息:“ [测量方法]数据已使用R和R包进行了分析...我们使用线性混合效应模型(以“主题”作为随机效应,并使用了线性混合效应模型来分析了对数转换后的因变量数据)。阶1(AR1)的自相关结构作为固定效应,包括“治疗”和“天”,以及“治疗”和“天”的相互作用,通过对残差图与拟合值进行目视检查来检查正态性和同质性,以评估在进行混合效应分析后,我们进行了似然比测试,将具有固定效应的模型与仅具有随机效应的零模型进行了比较。”)

对于如何向通常不喜欢统计数据的受众报告线性混合效应模型结果的任何建议(由相对统计数据的新手撰写),将不胜感激!


您能否进一步解释研究最后一天“居中”的含义?为什么不将研究的最后一天的天数从0抽到天数呢?另外-在摘要表中-交互作用在哪里?
colin 2014年

Answers:


1

这可能无法帮助回答您的问题,但是我注意到您的实验中有一个重复的度量(天),但是您没有表明这是模型中的一个重复的度量。我会认为您模型中的随机项是这样的:

mymodel <- lme(dv ~ Treatment*Day, random = ~1|Subject/Day, 
               data = mydf, na.action = na.omit,
               correlation = corAR1(form = ~1 |Subject/Day), method = "REML")

至于报告结果,您是否打算在开始看到治疗之间有显着差异的那一天进行报告?如果是这样,那么我认为您也需要查看/报告交互作用术语上的对比。我本人是统计新手,基本上和您有相同的问题:-)

Andy Field的“使用R发现统计信息”介绍了如何报告Ch14中的线性混合效应模型的结果。我手头没有这本书,但是一旦我再次使用它,就可以编辑该帖子。


这取决于,如果所有科目的日期都相同,我想那是一个交叉的因素,而不是嵌套的随机因素,对吗?然后它是因为day是连续的,因此它位于垂直线的前面。评论赞赏!random= Day|Subject
Aura Borealis
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.