[求助] R语言中 shapiro.test() 的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

各位高人，小弟在作多元回归分析，x里有20个变量和28个样本来做回归分析，最优的模型也出来了，但是最后要求用这个shapiro.test(x) 对x做个test，但是在R中却总是报错了

> shapiro.test(x.ini)
错误于`[.data.frame`(x, complete.cases(x)) :
  undefined columns selected

忘各位高人指点一下小弟，万分感谢

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：shapiro test Est API R语言 test 语言 shapiro

shapiro.test是做Shapiro-Wilk正态性检验的，x要求是数值向量，你的X是数据框格式。所以出错！

shapiro是检验正态性的，你要选x里的某个变量来做

但是老师要求对这个28个样本和20个变量整体作 shapiro检验，我应该怎么做这个检验呢？谢谢

shapiro.test只能做单变量正态性检验！

明白了，谢谢高人解答

不客气，互相学习，一起提高！！

这个问题很重要，R中有命令可以对多变量数据进行正态性检验，例子如下:

> library(mvnormtest)
> data(EuStockMarkets)
> EuStockMarkets[1:5,]
         DAX    SMI    CAC   FTSE
[1,] 1628.75 1678.1 1772.8 2443.6
[2,] 1613.63 1688.5 1750.5 2460.2
[3,] 1606.51 1678.6 1718.0 2448.2
[4,] 1621.04 1684.1 1708.1 2470.4
[5,] 1618.16 1686.6 1723.1 2484.7
> C <- t(EuStockMarkets[15:29,1:4])
> mshapiro.test(C)

        Shapiro-Wilk normality test

data:  Z
W = 0.8161, p-value = 0.005955

> R <- t(diff(t(log(C))))
> mshapiro.test(R)

        Shapiro-Wilk normality test

data:  Z
W = 0.8841, p-value = 0.06641

>

[此贴子已经被作者于2009-2-15 21:10:31编辑过]

嗯，对!mvnormtest中的mshapiro.test 是做多变量正态性检验，这个方法好像是Patrick Royston 1982年提出来的。

Stats中的shapiro.test只能做单变量正态性检验。

不过建议这两种方法都试一下，看看最终的检验结果有什么差异！

你的数据不对。shapiro.test(x)
x 为你要检验的数据列。如果x来自数据框，则要用子集提取的形式，把要要研究的数据提取出来  再检验。