Kolmogorov-smirnov测试矩阵

Question

r code

B = 10
CC <- numeric(B)
for(j in 1:10)
  N = rnorm(1000)
  print(ks.test(N,pnorm)$statistic)
  CC[j] <- ks.test(N,pnorm)$statistic
  print(CC[j])
print(CC)

如果你在r中运行它，你将获得一个包含9个零的矩阵，并且只有最后一个术语是实际的数值。有人可以解释为什么会发生这种情况，我怎样才能得到其他指数的数值。

Answer 1

可能是因为缺少{}。试试这个：

B = 10
CC <- numeric(B)
for(j in 1:10) {
  N = rnorm(1000)
  print(ks.test(N,pnorm)$statistic)
  CC[j] <- ks.test(N,pnorm)$statistic
  print(CC[j])
}
print(CC)