请问怎么做多元正态检验

winfight2012

7352

收藏 2015-07-28

我想用mshapiro.test这个函数，但结果如下，这是什么情况啊

Error in solve.default(R %*% t(R), tol = 1e-18) :
Lapack例行程序dgesv: 系统正好是奇异的: U[77,77] = 0

如果不用这个函数，用什么检验多元正态好呢
另外请问一下矩阵是怎么放的，我有321个变量，每个变量32个样本，是应该321行32列还是32行321列

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

victorchan0633

2015-8-1 14:49:35

应该是321行32列，按照你的数据结构，出现奇异是必然的，321行32列，矩阵的秩最大为32，再乘以自身的转置，必然要小于等于32，因此自然不会满秩（321*321），一般来说检验多元正态就用mshapiro.test了，不过你这里出现就矩阵可逆，建议先主成分分析，提起主成分，然后做一元正态性检验。