【数学建模算法】(30)数据的统计描述和分析(下)

2.参数估计

利用样本对总体进行统计推断的一类问题是参数估计,即假定已知总体的分布,通常是,估计有关的参数,如。参数估计分点估计和区间估计两种。

2.1.点估计

点估计是用样本统计量确定总体参数的一个数值。评价估计优劣的标准有无偏性、最小方差性、有效性等,估计的方法有矩法、极大似然法等。

最常用的是对总体均值和方差(或标准差)做点估计。如果抛开一切数学观点,从纯直接的角度看,一组自然、合理的点估计显然是(在字母上加表示它的估计值)

2.2.区间估计

点估计虽然给出了待估参数的一个数值,却没有告诉我们这个估计值的精度和可信程度。一般地,总体的待估参数记作(如),由样本算出的 的估计量记作,人们常希望给出一个区间,使 以一定的概率落在此区间内。若有:

则称为的置信区间,分别称为置信下限和置信上限,称为置信概率或置信水平, 称为显著性水平。

给出的置信水平为的置信区间,称为 的区间估计。置信区间越小,估计的精度越高;置信水平越大,估计的可信程度越高。但是这两个指标显然是矛盾的,通常是在一定的置信水平下使置信区间尽量小。通俗地说,区间估计给出了点估计的误差范围。

2.3.参数估计的Matlab实现

Matlab 统计工具箱中,有专门计算总体均值、标准差的点估计和区间估计的函数。对于正态总体,命令是:

[mu,sigma,muci,sigmaci]=normfit(x,alpha)

其中 x 为样本(数组或矩阵),alpha 为显著性水平 (alpha 缺省时设定为 0.05),返回总体均值和标准差 的点估计 mu 和 sigma,及总体均值和标准差 的区间估计muci 和 sigmaci。当 x 为矩阵时,x 的每一列作为一个样本。

Matlab 统计工具箱中还提供了一些具有特定分布总体的区间估计的命令,如expfit,poissfit,gamfit,你可以从这些字头猜出它们用于哪个分布,具体用法参见帮助系统。

3.假设检验

统计推断的另一类重要问题是假设检验问题。在总体的分布函数完全未知或只知其形式但不知其参数的情况,为了推断总体的某些性质,提出某些关于总体的假设。例如,提出总体服从泊松分布的假设,又如对于正态总体提出数学期望等于的假设等。假设检验就是根据样本对所提出的假设做出判断:是接受还是拒绝。这就是所谓的假设检验问题。

3.1.单个总体均值的检验

假设检验的种类:
双边检验:
右边检验:
左边检验:

3.1.1.已知,关于 的检验( 检验)

在 Matlab 中 检验法由函数 ztest 来实现,命令为:

[h,p,ci]=ztest(x,mu,sigma,alpha,tail)

其中输入参数 x 是样本,mu 是中的,sigma是总体标准差, alpha 是显著性水平 (alpha 缺省时设定为 0.05),tail 是对备选假设的选择:为时用tail=0(可缺省);为时用tail=1;为时用tail=-1。输出参数h=0表示接受,h=1表示拒绝,p表示在假设下样本均值出现的概率,p越小越值得怀疑,ci是的置信区间。

例3 某车间用一台包装机包装糖果。包得的袋装糖重是一个随机变量,它服从正态分布。当机器正常时,其均值为 0.5 公斤,标准差为 0.015 公斤。某日开工后为检验包装机是否正常,随机地抽取它所包装的糖 9 袋,称得净重为(公斤):
0.497 0.506 0.518 0.524 0.498 0.511 0.520 0.515 0.512
问机器是否正常?

解:总体已知,,未知。于是提出假设和
Matlab实现如下:

x=[0.497 0.506 0.518 0.524 0.498...
0.511 0.520 0.515 0.512];
[h,p,ci]=ztest(x,0.5,0.015)

求得,说明在 0.05 的水平下,可拒绝原假设,即认为这天包装机工作不正常。

3.1.2.未知,关于 的检验(检验)

在 Matlab 中 t 检验法由函数 ttest 来实现,命令为:

[h,p,ci]=ttest(x,mu,alpha,tail)

例4 某种电子元件的寿命(以小时计)服从正态分布,均未知.现得 16 只元件的寿命如下:
159 280 101 212 224 379 179 264
222 362 168 250 149 260 485 170
问是否有理由认为元件的平均寿命大于225(小时)?

解:按题意需检验

取。Matlab 实现如下:

x=[159 280 101 212 224 379 179 264 ...
222 362 168 250 149 260 485 170];
[h,p,ci]=ttest(x,225,0.05,1)

求得,说明在显著水平为 0.05 的情况下,不能拒绝原假设,认为元件的平均寿命不大于 225 小时。

3.2.两个正态总体均值差的检验(检验)

还可以用 t 检验法检验具有相同方差的 2 个正态总体均值差的假设。在 Matlab 中由函数 ttest2 实现,命令为:

[h,p,ci]=ttest2(x,y,alpha,tail)

与上面的 ttest 相比,不同处只在于输入的是两个样本 x,y(长度不一定相同),而不是一个样本和它的总体均值;tail 的用法与 ttest 相似,可参看帮助系统。

例5 在平炉上进行一项试验以确定改变操作方法的建议是否会增加钢的得率,试验是在同一平炉上进行的。每炼一炉钢时除操作方法外,其它条件都可能做到相同。先用标准方法炼一炉,然后用建议的新方法炼一炉,以后交换进行,各炼了 10 炉,其得率分别为:

标准方法 78.1 72.4 76.2 74.3 77.4 78.4 76.0 75.6 76.7 77.3
新方法 79.1 81.0 77.3 79.1 80.0 79.1 79.1 77.3 80.2 82.1

设这两个样本相互独立且分别来自正态总体和,均未知,问建议的新方法能否提高得率?(取。)
(1)需要假设检验

(2)Matlab实现

x=[78.1 72.4 76.2 74.3 77.4 78.4 76.0 75.6 76.7 77.3];
y=[79.1 81.0 77.3 79.1 80.0 79.1 79.1 77.3 80.2 82.1];
[h,p,ci]=ttest2(x,y,0.05,-1)

求得:
。表明在的显著水平下,可以拒绝原假设,即认为建议的新操作方法较原方法优。
注:ttest2 既可以做方差相等的,又可以做方差不相等的假设检验,使用格式为:

h = ttest2(x,y,alpha,tail, 'unequal' )

3.3.分布拟合检验——检验法

假设为:

:总体 的分布函数为,
:总体的分布函数不是。

在用下述检验法检验假设时,若在假设下的形式已知,但其参数值未知,这时需要先用极大似然估计法估计参数,然后作检验。
检验法的基本思想如下:将随机试验可能结果的全体分为个互不相容的事件:
。于是在假设下,我们可以计算(或),。在次实验中,事件出现的频率与往往有差异,但一般来说,若为真,且试验的次数又甚多时,则这种差异不应该很大。基于这种想法,皮尔逊使用:

作为假设检验的统计量。并证明了以下定理:

定理:若充分大,则当为真时(不论中的分布属什么分布),上述统计量总是服从于自由度为的分布,其中是被估计的参数的个数。

于是,若在假设下算得统计量有:

则在显著性水平下拒绝,否则就接受。

注意:使用检验法时,要求样本容量不小于50,以及每个都不小于 5,而且最好是在 5 以上。否则应适当地合并,以满足这个要求。

例6 下面列出了84 个伊特拉斯坎人男子的头颅的最大宽度(mm),试检验这些数据是否来自正态总体(取)。
141 148 132 138 154 142 150 146 155 158
150 140 147 148 144 150 149 145 149 158
143 141 144 144 126 140 144 142 141 140
145 135 147 146 141 136 140 146 142 137
148 154 137 139 143 140 131 143 141 149
148 135 148 152 143 144 141 143 147 146
150 132 142 142 143 153 149 146 149 138
142 149 142 137 134 144 146 147 140 142
140 137 152 145

编写Matlab程序:

clc
x=[141 148 132 138 154 142 150 146 155 158 ...
150 140 147 148 144 150 149 145 149 158 ...
143 141 144 144 126 140 144 142 141 140 ...
145 135 147 146 141 136 140 146 142 137 ...
148 154 137 139 143 140 131 143 141 149 ...
148 135 148 152 143 144 141 143 147 146 ...
150 132 142 142 143 153 149 146 149 138 ...
142 149 142 137 134 144 146 147 140 142 ...
140 137 152 145];
mm=minmax(x) %求数据中的最小数和最大数
hist(x,8) %画直方图
fi=[length(find(x<135)),...
    length(find(x>=135&x<138)),...
    length(find(x>=138&x<142)),...
    length(find(x>=142&x<146)),...
    length(find(x>=146&x<150)),...
    length(find(x>=150&x<154)),...
    length(find(x>=154))] %各区间上出现的频数
mu=mean(x),sigma=std(x) %均值和标准差
fendian=[135,138,142,146,150,154] %区间的分点
p0=normcdf(fendian,mu,sigma) %分点处分布函数的值
p1=diff(p0) %中间各区间的概率
p=[p0(1),p1,1-p0(6)] %所有区间的概率
chi=(fi-84*p).^2./(84*p)
chisum=sum(chi) %皮尔逊统计量的值
x_a=chi2inv(0.9,4) %chi2分布的0.9分位数

得皮尔逊统计量chisum=2.2654,故在水平0.1下接受,,即认为数据来自正态分布总体。

你可能感兴趣的:(【数学建模算法】(30)数据的统计描述和分析(下))