这个是文献截图beta(2,2),n=20时,用KS法检验样本,模拟1000次,当α《0.05时记为1,求和后得出功效。
我用的程序是
> library(fBasics)
> library(akima)
> n=10000
> c4=numeric(n)
> for (i in 1:n){
+ a1<-rbeta(20,2,2)
+ b4<-ks.test(a1,'pnorm')
+ c4
=ifelse(b4$p.value<0.05,1,0)
+ sum(c4==1)
+ }
> d4<-c4
> sum(d4)*100/n
[1] 100
>
得出的功效确实100%?
这是怎么一回事?我用文献中其他几种方法(SW法,JB法)得出的结果基本和文献时一致的,可是KS法却完全不一样,为什么?求助大家!!!!!!!!!!!!!!!!