混合效应模型中的“方差成分参数”是什么?


9

贝茨关于混合效果模型的书的第12页,他对模型进行了如下描述:

贝茨的混合效应模型

在屏幕快照的结尾处,他提到了

相对协方差因子 ,这取决于方差分量参数θΛθθ

没有解释到底是什么关系。假设我们给出,我们如何获得Λ θ从它?θΛθ

与此相关的是,这是我发现贝茨的论述缺乏细节的众多例子之一。是否有更好的文字实际经过参数估计的优化过程和测试统计量分布的证明?


1
我认为只是一种手段方差分量的,你将承担什么,如AR(1)或联合国等θ
深北

@DeepNorth我一直在更仔细地阅读本文,在某些时候,作者谈论了关于优化的可能性。因此我认为θ必须是一个实际参数。(第108页,第5.4.2秒)θθ
Heisenberg

您设法弄清楚了吗?,我在理解协方差矩阵和theta之间的关系时遇到了同样的困难。

你放弃了这个问题吗?到目前为止,已经提供了两个答案,但没有对它们进行任何评论。请考虑对答案给出建设性的反馈意见,这样,如果他们没有提供(令人满意的)解决方案,则至少可以展开讨论来缩小问题的范围并寻求解决方案。如果您对问题的答案没有反应,则会阻止进一步的答案。
跳过

Answers:


3

这是层次推理。线性模型中有一堆参数,即b的成分。在一个纯粹的固定效应模型中,您将仅获得这些估计值。取而代之的是,您想象b中的值本身是从一个多元正态分布中提取的,该多元分布具有由theta参数化的协方差矩阵。这是一个简单的例子。假设我们在10个不同位置的五个不同时间段查看动物计数。我们将得到一个线性模型(我在这里使用R talk),它看起来像count〜time + factor(location),这样,在这种情况下,您将对所有回归具有共同的斜率(每个斜率一个)位置),但每个位置的截距不同。我们可以将其称为固定效果模型,然后估计所有截距。然而,如果希望从大量可能的位置中选择10个位置,我们可能不关心这些位置。因此,我们在截距上放置了协方差模型。例如,我们将截距声明为多变量正态且独立,且具有共同方差sigma2。然后sigma2是“ theta”参数,因为它描述了每个位置处的截距的总体特征(因此是随机效应)。


2

θd

Λθq×qθq×q

Λθ=θ×一世q

fm01ML

Λθθ

fm03

与两个嵌套的随机效应项相同(第43页,图2.10,此处未显示)。

Λθ

在此处输入图片说明

Λθ

在此处输入图片说明


补充笔记:

θ一世=σ一世σσ一世σ

lme4merModΛθgetME

image(getME(fm01ML, "Lambda"))
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.