r - 矩阵中的 Kolmogorov-smirnov 检验

Question

代码

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)