r - 在 R : Shapiro Wilk, Jarque Bera, Lilliefors 中重复 1000 次正态性测试

Question

学生和专业人士，

我目前正在尝试为随机样本大小 (T=10,30,50,100,500) 编写正态性测试。

我用于正态性测试的函数如下:

sim1 <- rnorm(10)

sw10 <- shapiro.test(sim1)

这对于每个样本量

这会生成一个包含测试信息的列表，必须以 90%、95% 和 99% 的置信度对其进行解释。

我面临的问题是我需要重复这个过程 1000 次。但是在计算相同 p 值的情况下，使用相同的样本 sim1 在这种情况下没有帮助。

那么我使用以下内容吗？

rsw10 <- replicate(shapiro.test(rnorm(10))

另外我必须计算相对拒绝频率，我如何提取该信息？

此致

Answer 1

如果我猜对了，它会是这样，首先是代表的数量，然后是函数:

sim = replicate(1000,shapiro.test(rnorm(10)))

拒绝是这样的，假设 alpha 为 0.05 :

table(sim["p.value",]<0.05)

FALSE  TRUE 
  961    39