我想用mshapiro.test这个函数,但结果如下,这是什么情况啊
Error in solve.default(R %*% t(R), tol = 1e-18) :
Lapack例行程序dgesv: 系统正好是奇异的: U[77,77] = 0
如果不用这个函数,用什么检验多元正态好呢
另外请问一下矩阵是怎么放的,我有321个变量,每个变量32个样本,是应该321行32列还是32行321列
|
楼主: winfight2012
|
7410
2
[问答] 请问怎么做多元正态检验 |
|
初中生 61%
-
|
回帖推荐victorchan0633 发表于2楼 查看完整内容 应该是321行32列,按照你的数据结构,出现奇异是必然的,321行32列,矩阵的秩最大为32,再乘以自身的转置,必然要小于等于32,因此自然不会满秩(321*321),一般来说检验多元正态就用mshapiro.test了,不过你这里出现就矩阵可逆,建议先主成分分析,提起主成分,然后做一元正态性检验。
| ||
|
|
加好友,备注cda京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明


