计算
m <- 1000
G <- numeric(m)
for(i in 1:1000){
  n <- 100
  z <- sort(rbinom(n,1,0.1))
  S=0
  for(j in 1:n){
    S=S+(2*j-n-1)*z[j]
  }
  G <- S/(n^2*mean(z))
}

hist(G,probability=TRUE)
根据上述代码生成的G画的直方图为什么画出来的频率直方图为什么纵坐标一直大于1