t distribution

x <- rnorm(1000, mean = 5, sd=5)
hist(x)

图片alt

qqnorm(x,col="red",pch=19)
qqline(x,col="black")

图片alt

sapply(1:100,function(i)sample(x,10,replace = T)) |>
    colMeans() |> mean()
# 4.99705893538348

sapply(1:100,function(i)sample(x,10,replace = T)) |>
    colMeans() |> sd()
# 1.5179249240205