我需要比较许多多元方法。通常,我会使用 Hotelling 的 T 方检验统计量来做到这一点。
原来的霍特林方程是:T^2 = (nxny/nx+ny) (XY)' S^-1 (XY)
其中 X 和 Y 是向量均值,S 是池化协方差矩阵,nx/y 是样本大小。
但是,正常 Hotelling 检验的假设是样本协方差矩阵是相等/齐次的。我从 Box 的测试中知道这对我的数据不正确。这些网站提供了 Hotelling 的 T 方检验的修改版本,它不假定协方差矩阵相等:
修改后的方程为:T^2 = (XY)' ((Sx/nx) + (Sy/ny))^-1 (XY)
其中 X 和 Y 是向量均值,Sx/y 是相应的协方差矩阵,nx/y 是样本大小。
我已经搜索了 R 包,试图找到一个没有运气的等式的修改版本。有谁知道可以在 R 中执行此操作的包?