07-数据的统计分析与描述.ppt

* * 1、给出数组data的频数表的命令为： [N,X]=hist(data,k) 此命令将区间[min(data),max(data)]分为k个小区间（缺省为10），返回数组data落在每一个小区间的频数N和每一个小区间的中点X. 2、描绘数组data的频数直方图的命令为： hist(data,k) 4.4 频 数 直 方 图 的 描 绘 * * 1、正态总体的参数估计 设总体服从正态分布，则其点估计和区间估计可同时由以下命令获得： [muhat,sigmahat,muci,sigmaci] = normfit(X,alpha) 此命令在显著性水平alpha下估计数据X的参数（alpha缺省时设定为0.05），返回值muhat是X的均值的点估计值，sigmahat是标准差的点估计值, muci是均值的区间估计,sigmaci是标准差的区间估计. 4.5 参数估计 * * 2、其它分布的参数估计 有两种处理办法: 一.取容量充分大的样本（n50），按中心极限定理，它近似地 服从正态分布； 二.使用Matlab工具箱中具有特定分布总体的估计命令. （1）[muhat, muci] = expfit(X,alpha)----- 在显著性水平alpha下，求指数分布的数据X的均值的点估计及其区间估计. （2）[lambdahat, lambdaci] = poissfit(X,alpha)----- 在显著性水平alpha下，求泊松分布的数据X 的参数的点估计及其区间估计. （3）[phat, pci] = weibfit(X,alpha)----- 在显著性水平alpha下，求Weibull分布的数据X 的参数的点估计及其区间估计. * * 参数估计与区间估计 * * 【例】 * * * * 在总体服从正态分布的情况下，可用以下命令进行假设检验. 1、总体方差sigma2已知时，总体均值的检验使用 z-检验 [h,sig,ci] = ztest(x,m,sigma,alpha,tail) 检验数据 x 的关于均值的某一假设是否成立，其中sigma 为已知方差， alpha 为显著性水平，究竟检验什么假设取决于 tail 的取值： tail = 0，检验假设“x 的均值等于 m ” tail = 1，检验假设“x 的均值大于 m ” tail =-1，检验假设“x 的均值小于 m ” tail的缺省值为 0， alpha的缺省值为 0.05. 返回值 h 为一个布尔值，h=1 表示可以拒绝假设，h=0 表示不可以拒绝假设，sig 为假设成立的概率，ci 为均值的 1-alpha 置信区间. 4.6 假设检验 * * 例7 Matlab统计工具箱中的数据文件gas.mat.中提供了美国1993年一月份和二月份的汽油平均价格（price1,price2分别是一，二月份的油价，单位为美分），它是容量为20的双样本.假设一月份油价的标准偏差是一加仑四分币（?=4），试检验一月份油价的均值是否等于115. 解 作假设：m = 115. 首先取出数据，用以下命令： load gas 然后用以下命令检验 [h,sig,ci] = ztest(price1,115,4) 返回：h = 0，sig = 0.8668，ci = [113.3970 116.9030]. 检验结果: 1. 布尔变量h=0, 表示不拒绝零假设. 说明提出的假设均值115 是合理的. 2. sig-值为0.8668, 远超过0.5, 不能拒绝零假设 3. 95%的置信区间为[113.4, 116.9], 它完全包括115, 且精度很 高. . To MATLAB（liti7） * * 例 某车间用一台包装机包装葡萄糖，包得的袋装糖重是一个随机变量，它服从正态分布。当机器正常时，其均值为0.5公斤，标准差为0.015。某日开工后检验包装机是否正常，随机地抽取所包装的糖9袋，称得净重为（公斤） 0.497, 0.506, 0.518, 0.524, 0.498, 0.511, 0.52, 0.515, 0.512 问机器是否正常？ 解：总体μ和σ已知，该问题是当σ2为已知时，在α=0.05水平下，根据样本值判断μ=0.5还是μ≠0.5。为此提出假设： 原假设：H0: μ=μ0=0.5 备择假设：H1: μ≠