假设检验方差分析回归分析.doc

《数理统计》课内实验报告 学生姓名: 张学阳 1009300132 及 学 号： 学 院: 理学院 班 级: 数学101 课程名称： 数理统计 实验题目: 假设检验、方差分析、回归分析 指导教师 姓名及职称： 王亮红 讲师 郭新辰 教授 韩 玉 副教授 2012年11月13日 目 录 一、实验目的 1 二、实验内容 1 三、实验要点及说明 2 四、实现方法 2 五、实验结果 3 六、源程序清单 4 七、思考及总结 5 一、实验目的 1. 掌握假设检验、方差分析、回归分析的概念； 2. 学会利用Matlab软件实现对实验数据的假设检验、方差分析、回归分析等统计分析方法。 二、实验内容 1. 某种零件的尺寸方差为，对一批这类零件检查6件得尺寸数据（单位：毫米）为： 32.56 29.66 31.64 30.00 21.87 31.03 设零件尺寸服从正态分布，问这批零件的平均尺寸能否认为是32.50毫米（显著性水平）。 2. 按照规定，每100克罐头番茄汁中，维生素C的含量不得少于21毫克，现从某厂生产的一批罐头中抽取17个，测得维生素C的含量（单位：毫克）如下： 22 21 20 23 21 19 15 13 16 23 17 20 29 18 22 16 25 已知维生素C的含量服从正态分布，试检验这批罐头的维生素含量是否合格（显著性水平）。 4．为寻求适应某地区的高产油菜品种，今选了五种不同品种进行 试验，每一品种在四块试验田上得到在每一块田上的亩产量如下： 品种 田块 A1 A2 A3 A4 A5 1 256 244 250 288 206 2 222 300 277 280 212 3 280 290 230 315 220 4 298 275 322 259 212 问题：试分析不同品种的平均亩产量是否有显著差异？若存在显著性差异，哪个品种的亩产量高？并给出参数的估计值。 (均值、方差、置信区间) 三、实验要点及说明 对于问题1和问题2首先判断属于标准差已知还是未知的检验，然后根据检验统计量判断是否落入拒绝域，再判断是否合理。 四、实现方法 1.总体，待检验的原假设与备择假设分别为 这是一个双边检验问题，检验的拒绝域为 取显著水平查表知，具体数据可由Matlab求解。 2.这是一个关于均值的单边的假设检验问题。原假设是 由于标准差未知，故采用检验，其拒绝域为 显著性水平，可以查表和通过Matlab计算得到相应的数值。 4.若取，由Matlab求解得到值为0.0162，由于值小于，故拒绝原假设。下面是参数估计 各水平均值的最大似然估计 通过Matlab求解得到264,、277.25、269.75、285.5、212.5 因此A4产量最高 方差的无偏估计为 置信区间为 求解得到 [234.503, 293.497] [247.753, 306.747] [240.253, 299.247] [256.003, 314.997] [183.003, 241.997] 五、实验结果 问题一结果： H =1 sig = 1.2923e-011 结果H =1，说明在0.05的水平下，可拒绝原假设，即这批零件的尺寸不能认为是32.50毫米。 问题二结果： h =0 sig =0.1581 结果表明，h=0，即在显著水平为0.025的情况下，不能拒绝原假设，认为这批罐头的维生素含量合格。 问题四结果： P = 0.0162 图一 anova1函数的数据结构结果图形 图二 anova1函数的箱形图 由于P的值小于0.05，故拒绝原假设，认为各台机器生产的薄板厚度有显著的差异。 六、源程序清单 问题一源程序： X=[32.56 29.66 31.64 30.00 21.87 31.03]; [H,sig]=ztest(X,32.5,1.1,0.05,0) 问题二源程序： Y=[22 21 20 23 21 19 15 13 16 23 17 20 29 18 22 16 25]; [h,sig]=ttest(Y,21,0.025,-1) 问题四源程序： Z=[256 244 250 288 206; 222 300 277 280 212; 280 290 230 315 220; 298 275 322 259 212]; P=anova1(Z) mean(Z) mean(Z)-sqrt(766.1)*2.1314/2 mean(Z)+sqrt(