我是R的新人。如果你能帮助我那会很棒。我的问题如下。我有这段代码:
N <- 1000 # number of simulations
n <- 30
prob <- 0.001
Y2 <- rep(0, times=N)
for(i in 1:N) {
X <- rbinom(n, size=1, prob)
#binomial experiment
Y2[i] <- sum(X)
}
table(Y2)
library(dgof)
ks.test(Y2, "pnorm", mean=prob*n, sd=sqrt(n*prob*(1-prob)))
我想将这个代码生成10次,所以我将获得10个不同的p值和D.此外,参数n将不是固定的,但它将是一个值向量(50,100,500,1000,5000) )。所以我想要一张像照片一样的ks.test矩阵。 THX。enter image description here