r语言 - 矩阵中的柯尔莫哥罗夫-斯米尔诺夫测试

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)