区间估计不仅仅给出参数的近似取值,还给出了取该值的误差范围。要求参数的区间估计,首先要求出该参数的点估计,然后构造一个含有该参数的随机变量,并根据一定的置信水平求该估计值的误差范围。
在用mle函数进行最大似然估计时,还有几种语法格式,例如:
[phat, pci] = mle('dist', data):返回最大似然估计和95%置信区间。
[phat, pci] = mle('dist', data, alpha):返回指定分布的最大似然估计值和100(1-alpha)置信区间。
[phat, pci] = mle('dist', data, alpha, p1):该形式仅用于二项分布,其中p1为试验次数。
这几种语法格式不仅可以返回参数的估计量,还可以返回参数取值的置信区间。
从一批灯泡中随机地取5只灯泡进行寿命试验,测得寿命如下(以小时计)
1050 1100 1120 1250 1280
这里假设灯泡寿命服从正态分布,求灯泡寿命平均值的95%置信区间。
code.matlab
>> X=[1050 1100 1120 1250 1280];
>> [p,ci]=mle('norm',X,0.05)
p =
1.0e+003 *
1.1600 0.0892
ci =
1.0e+003 *
1.0818 0.0339
1.2382 0.1445
所以,灯泡寿命平均值的95%置信区间为[1081.8 1238.2]。