R中的正态性测试重复1000次：Shapiro Wilk，Jarque Bera和Lilliefors

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