我正在进行一项小型模拟研究,以评估两个正态性检验的好坏。我的计划是生成许多不太多的正态分布样本,并确定每个测试有多少次拒绝正态性的零假设。
我目前有的(不完整)代码如下:
现在我想计算每个测试中小于0.05的p值的比例。请问我该如何做?如果这是一个新手问题,我很抱歉,因为我自己对R也很陌生。
谢谢。
我目前有的(不完整)代码如下:
library(nortest)
y<-replicate(10000,{
x<-rnorm(50)
ad.test(x)$p.value
ks.test(x,y=pnorm)$p.value
}
)
现在我想计算每个测试中小于0.05的p值的比例。请问我该如何做?如果这是一个新手问题,我很抱歉,因为我自己对R也很陌生。
谢谢。