mean()函数: 用来求均值 样本均值可以用mean(s1)求
对于正态分布N(μ, σ^2)，求样本方差的估计量：
var()函数：用于求方差
先求样本标准差 sigma.hat=sqrt(19*var(s1)/20) 样本方差为 sigma.hat^2
最大似然估计涉及到一个分布，是gumbel distribution(极值分布) 用到的包： library(MASS) library(evd) gumbel.fit=fitdistr(x=s1, densfun = dgumbel, start=list(loc=50, scale=15),lower = c(-Inf, 1e-4)) gumbel.fit