0

我在使用 nor.test 函数作为 R 中的单向测试时遇到问题。我的数据包含按处理 (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

谢谢您的帮助！

4

1 回答 1

0

没有使用那个包，但是根据文档（和你的错误），nor.test 默认执行 Shapiro-Wilk 正态性测试，它需要一个数字向量作为输入（至少 3 个值）。我的猜测是，有一个基于 的组，Traitement其值少于 3 个或超过 5000 个。尝试使用类似

table(dataMax$Traitement)
于 2020-04-10T09:23:32.643 回答