iris

data(morley)
dat <- morley[, 3] + 299000
hist(dat, breaks=seq(299600, 300100, by=20), col=5, ylim=c(0, 17))

# 一様分布する整数の乱数
x <- runif(10)
max(x)

x <- runif(100)
max(x)

n <- 10
x <- runif(n)
((n+1)/n)* max(x) # 不偏推定量