ペアの数

  • 分散はばらつき:\sum_{i,j} (x_i-x_j)^2をペア数で割ったもの(の半分の値)が知りたい
  • 有限個の要素のときペア数はn(n-1)なので…
X <- c(rnorm(10),runif(10))
var(X)
sum((outer(X,X,"-"))^2)
sum((outer(X,X,"-"))^2)/(n*(n-1))/2