嵌套方差分析与多级线性模型

机器算法验证 方差分析 多层次分析 lme4-nlme 嵌套数据
2022-04-04 11:42:53

我正在研究一个在两种治疗条件下围绕单个轴突/神经细胞的问题。所有相关问题都明确在轴突水平上,并且取决于单个轴突的特性(大小、髓鞘形成等)。然而,我认为将每个轴突视为一个独立的观察结果而不考虑它们属于哪些动物会作弊。因此,我想使用嵌套在动物中的轴突的两级方法。

阅读此内容,在我看来,我可以使用嵌套方差分析(aov在 R 中,使用 animal_id 作为随机效应)或多级线性模型(nlme/lme4具有随机截距)。但是我很难在它们之间做出决定(实际上,两者实际上都产生了相当相似的结果)。

这是我到目前为止所学到的:

  • 我理解nlme/lme4使用最大似然方法,这应该使它们更适合不平衡设计,但我认为这对我来说不是一个特别大的问题。

  • 我还听说过嵌套方差分析可能更适合较小的样本量的建议。这是真的?

  • 最后,开发人员lme4似乎强烈反对将 p 值分配给他们的结果的想法。如果我理解正确,这是因为“自由度”在混合线性模型中不是一个非常直接的概念。但是每当我看到人们提出这个论点时,他们似乎也反对一般的 p 值和零假设检验。aov与 相比,信任嵌套中的 p 值是否有更好的情况lme4

我非常感谢任何回答或评论。

谢谢

1个回答

在图书馆呆了几天之后,我仍然觉得离理解这两个选项的所有含义还很遥远。但是,我学到了一些东西,如果其他人也有同样的想法,我会非常推荐 Snijders 和 Bosker 的《多层次分析》一书来指导你的决定。这是迄今为止我能找到的关于这个主题的最佳文本,涵盖了从嵌套方差分析到非常复杂的多级线性模型。

对我来说,它为我上面的三点提供了一些澄清:

1)即使我有相当平衡的群体,我的群体规模也足够小,以至于在微观层面上改变一些值会对群体方差产生显着影响。在这种情况下,多级线性模型应该更稳定,因为一个群体的群体效应会影响其他群体的群体效应。S&B 强调,这种影响对于人数少于 50 人的群体可能特别明显。[+1 传销]

2) S&B 建议在 MLM 分析中至少有 10 个,最好是 20 多个组。我绝对处于这个范围的低端。[+1 方差分析?]

3)是的,嵌套方差分析的 F 检验肯定要简单得多。S&B 似乎认为传销中也存在 p 值的情况,但阅读这本书可以很清楚地知道为什么它要复杂得多。[+1 方差分析]

我希望这可以帮助别人。我仍然非常不确定如何继续,所以如果有人有评论,我会很高兴阅读它。

干杯