我正在尝试在 R 中复制 SPSS 的线性判别分析输出,但我很难找到执行 m-box 测试的方法。
我发现的唯一一件事是在论坛上发布的一些代码,用于手动实现该过程,但我想知道是否没有任何用于此目的的内容已经包含在语言本身中。
我正在尝试在 R 中复制 SPSS 的线性判别分析输出,但我很难找到执行 m-box 测试的方法。
我发现的唯一一件事是在论坛上发布的一些代码,用于手动实现该过程,但我想知道是否没有任何用于此目的的内容已经包含在语言本身中。
在包biotools 中,您可以找到函数boxM(data, grouping)。它根据一个分类因子对从多元正态数据获得的协方差矩阵的同质性执行 Box 的 M 检验。该测试基于卡方近似。
可以通过简单的 rseek 搜索找到代码。通常不会这样做,因为它非常高的灵敏度会导致显着的 p 值,这可能并不意味着什么。
编辑: 该旧链接不再起作用,但事实证明该测试是在 biotools 包中使用功能 boxM 实现的。这仍然是一个相对容易的搜索。确实,您可能不应该像所有此类测试一样费心使用它。您可能应该仔细检查您的协方差矩阵和您的假设。