HT <- 144:191 NUM <- c(2,3,6,17,19,56,125,219,463,915,1609,2649,4550,7214,11005,16081, 22098,29903,39048,48312,57703,66639,73332,78051,79829,77866,73767,66321, 57993,48410,39081,29967,22055,15810,10875,7309,4596,2726,1519,939,462, 224,128,50,31,14,5,4) X <- rep(HT,NUM) RNGkind("Mersenne-Twister",normal.kind="Inversion") set.seed(1) V <- rep(0,1000) UV <- rep(0,1000) for (i in 1:1000) { s10 <- sample(X,10) m10 <- mean(s10) V[i] <- sum((s10-m10)^2)/10 UV[i] <- sum((s10-m10)^2)/9 } layout(c(1,2)) hist(V,main="サイズ10の標本の分散の分布",xlim=c(0,80)) hist(UV,main="サイズ10の標本の不偏分散の分布",xlim=c(0,80)) print(c(mean(V),mean(UV)))