矩陣中的Kolmogorov-smirnov檢驗

Question

r代碼

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)