我试图实现一个em算法家庭数据,我假设我的观察有一个多元t分布。我每个家庭只有两个兄弟姐妹,所以所有的家庭成员只有两个观察结果。基本上,我试图遵循本文中的E(C)m步骤:
https://pdfs.semanticscholar.org/9445/ef865c4eb1431f9cb2abdb5efc1c361172cc.pdf
但是,现在我不确定em是否适用于这类数据,因为我的相关矩阵psi对于族应该是块对角的。
这里有一个关于我的家庭结构的例子

fam_id = sort(rep(1:5, 2))

Z= matrix(0, nrow = length(fam_id), ncol = length(unique(fam_id)))
colnames(Z) = unique(fam_id)

k = 1
i = 1
# Random effects dummy matrix
while (k <= ncol(Z)) {
    Z[i:(i+1), k] = c(1, 1)
    k = k +1
    i = i+2
}

> Z
  1 2 3 4 5
[1,] 1 0 0 0 0
[2,] 1 0 0 0 0
[3,] 0 1 0 0 0
[4,] 0 1 0 0 0
...

EM算法在第5次迭代后窒息,表示相关矩阵Psi不是:
解算错误。默认值(磅/平方英寸):
系统在计算上是奇异的
如果有人能解释一下,我会很高兴的!

最佳答案

请在Statschange网站上查看此答案
https://stats.stackexchange.com/questions/76488/error-system-is-computationally-singular-when-running-a-glm
在第五次迭代中,你可能会得到一个不可逆矩阵

关于r - 多元t混合模型的EM算法。,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/43367852/

10-10 15:31