我在R中使用nor.test函数作为单向测试时遇到问题。我的数据包含按处理(Traitement(分组的产量值(Rdt_pied(。在每次治疗中,我的值在60到90之间。
> describe(Rdt_pied ~ Traitement, data = dataMax)
n Mean Std.Dev Median Min Max 25th 75th Skewness Kurtosis NA
G(1) 90 565.0222 282.1874 535.0 91 1440 379.00 751.25 0.7364071 3.727566 0
G(2) 90 703.1444 366.1114 632.5 126 1628 431.50 1007.75 0.4606251 2.392356 0
G(3) 90 723.9667 523.5872 650.5 64 2882 293.50 1028.50 1.2606231 5.365014 0
G(4) 90 954.1000 537.0138 834.5 83 2792 565.25 1143.75 1.1695460 4.672321 0
G(A) 60 368.0667 218.1940 326.0 99 1240 243.00 420.00 2.2207612 9.234473 0
G(H) 60 265.4667 148.0383 223.5 107 866 148.00 357.25 1.3759925 5.685456 0
G(S) 60 498.8000 280.1277 401.0 170 1700 292.75 617.50 1.6792061 7.125804 0
G(T) 60 521.7167 374.7822 448.5 74 1560 214.00 733.25 1.1367209 3.737134 0
>
为什么nor.test会给我这个答案?
> nor.test(Rdt_pied ~ Traitement, data = dataMax)
Error in shapiro.test(y[which(group == (levels(group)[i]))]) :
sample size must be between 3 and 5000
谢谢你的帮助!
尚未使用该包,但根据文档(以及您的错误(,nor.test默认执行Shapiro-Wilk正态性测试,该测试需要一个数字向量作为输入(至少3个值(。我的猜测是,有一个基于Traitement
的组,其值小于3,或大于5000。试着用之类的东西检查一下
table(dataMax$Traitement)