0
我試圖實現家庭數據的EM算法，我假設我的觀察具有多變量t分佈。我每個家庭只有兩個兄弟姐妹，所以所有的家庭組都只有兩個觀察。基本上我試圖按照這篇文章中的E（C）M步驟： https://pdfs.semanticscholar.org/9445/ef865c4eb1431f9cb2abdb5efc1c361172cc.pdf多變量t混合模型的EM算法
但是，現在我不確定EM是否適合這種數據，家庭的對角線。
因此，這裏是我的家庭結構如何
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
...
5日迭代說，相關矩陣的PSI後不EM算法扼流圈的R例如：
在solve.default錯誤（psi_hat ）： 系統在計算上是單數
如果任何人都可以對此有所瞭解，我會非常高興！