我想在R中复制SPSS的线性判别分析输出,而且我很难找到执行m-box测试的方法。R Box'M测试homocedasticity
我发现的唯一的事情是在论坛上发布了一些代码来手动实现这个过程,但我想知道是否没有为此目的已经包含在语言本身中的任何东西。
我想在R中复制SPSS的线性判别分析输出,而且我很难找到执行m-box测试的方法。R Box'M测试homocedasticity
我发现的唯一的事情是在论坛上发布了一些代码来手动实现这个过程,但我想知道是否没有为此目的已经包含在语言本身中的任何东西。
有code,可以发现一个简单的rseek搜索。这通常不会完成,因为它非常高的灵敏度会导致显着的p值,这可能没有多大意义。
编辑: 那旧的链接不工作了,但事实证明,试验是在Biotools公司包的功能boxM实现。这仍然是一个相对容易的搜索。而且,你可能不应该像所有这些测试一样使用它。你应该仔细检查你的协方差矩阵和你的假设。
在包Biotools公司你可以找到函数boxM(数据分组)。它根据一个分类因子对从多元正态数据获得的协方差矩阵的同质性进行Box的M检验。该测试基于卡方近似。
我不熟悉“m-box”测试。 SPSS专用术语?我听说过Box-Pierce,Box-Jenkins和Box-Q测试,但这些都与自相关有关。 – 2010-08-03 23:22:30
我拼错了名字,它被称为Box'm测试。无论如何,它不像我那样常用(我刚开始学习这个主题),可能是因为它的高敏感性,正如约翰在下面评论的那样。 它用于测试MANOVA假设之一,看起来大多数人要求它是SPSS用户,因为该软件将其作为标准输出之一。 只是为了防止别人感兴趣,这里是SPSS算法的规范: ftp://ftp.spss.com/pub/spss/statistics/spss/algorithms/ap14boxm.pdf (维基百科有一个定义,但只有在德语版本中) – Sam 2010-08-04 14:15:54