如何在状态模型MANOVA，Python中求出组内方差协方差矩阵？

Question

我试图提取组内方差协方差矩阵来估计状态模型中MANOVA模型中的平面度，使用以下公式：

​

我已经阅读了statsmodels用户指南，但还没有看到它。到目前为止，我所做的是：

resid = manova.mv_test().results['x1']['E']
n * (xbar_tot).T @ np.linalg.solve(?, resid) @ xbar_tot

在R中，也可以通过调用：

n * t(xbar.tot) %*% solve(summary(mod, test="Wilks")$SS$Residuals) %*% xbar.tot

我所缺少的是逆的第一部分，代码中有?的位置。

Answer 1

结果，解决方案只是简单的E matrix：因此，对于我提出的statsmodels问题，可以这样做：

#1. run the statsmodel manova model
manova = MANOVA.from_formula('dv1 + dv2 ~ y', data=dy)

#2. extract the E matrix from the mv_test().results
resid = manova.mv_test().results['y']['E']

#3. obtain Hotelling's T^2
n * (xbar_tot).T @ np.linalg.solve(resid) @ xbar_tot