回归分析作业详解.docx

1、回归分析作业详解一、假设检验:1. 某药厂最近研制出一种新的降压药，为了验证其疗效，选择15个高血压病人进行实验。数据表是服药前后的血压值。选用适当的统计方法验证该药是否有效。patient12345678before11513512713010390101104after1091201251301059490100patient9101112131415before10989120113118130120after9090110103100121108解：设血压值总体分布服从XN(,S2)假设检验问题:1 为服药前血压值样本均值，2为服药后血压值样本均值 原假设H0: 1=2 即服药前后血压

2、值相等。 备择假设H1：12 即服药前后血压值不相等。然后采用双边检验的方法，使用Minitab进行检验，得到以下数据： Paired T-Test and CI: X(服药之前）, Y(服药之后） Paired T for X(服药之前） - Y(服药之后） N Mean StDev SE MeanX(服药之前） 15 113.60 14.21 3.67Y(服药之后） 15 106.33 13.04 3.37Difference 15 7.27 7.27 1.8895% lower bound for mean difference: 3.96T-Test of mean differenc

3、e = 0 (vs 0): T-Value = 3.87 P-Value = 0.001由以上数据可以得出：p值=0.001=0.05，拒绝原假设，服用新的降压药后血压值有显著的差异，即该降压药有效。二、方差分析：1、对于硅酸盐水泥的抗折强度，用四种不同的配方方法收集了以下数据：配方法抗折强度13129300028652890232003300297531503280029002985305042600270026002765（1） 检验配方法影响水泥砂浆强度的假设。（2） 选择一种比较方法对均值进行比较。解答：1. （1）原假设：以上四种配方方法下硅酸盐水泥的抗折强度相等。备择假设：以上四种

4、配方方法下硅酸盐水泥的抗折强度存在差异。使用Minitab进行检验，得到以下数据： One-way ANOVA: 抗折强度 versus 配方方法 Source DF SS MS F P配方方法 3 489740 163247 12.73 0.000Error 12 153908 12826Total 15 643648S = 113.3 R-Sq = 76.09% R-Sq(adj) = 70.11% Individual 95% CIs For Mean Based on Pooled StDevLevel N Mean StDev 1 4 2971.0 120.6 2 4 3156.3

5、136.0 3 4 2933.8 108.3 4 4 2666.3 81.0 2600 2800 3000 3200Pooled StDev = 113.3 抗折强度与配方方法的箱线图 P值=0.0000.05，拒绝原假设。可知不同配方方法下硅酸盐水泥的抗折强度存在显著差异。(2) 选择一种比较方法对均值进行比较：方差分析表中ANOVA的计算及箱线图的绘制可以得到各组间差异显著的结论，再对均值进行多重比较的检验。从“统计方差分析单因子（StatANOVAOne way）”进入相关界面后，点击“比较（Comparisons）”。选中Tukey算法，设定整体误差率为5%。点击“确定”后，分析结果如

6、下：Grouping Information Using Tukey Method配方法 N Mean Grouping2 4 3156.3 A1 4 2971.0 A3 4 2933.8 A4 4 2666.3 BMeans that do not share a letter are significantly different.第一组数据Tukey 95% Simultaneous Confidence IntervalsAll Pairwise Comparisons among Levels of 配方方法Individual confidence level = 98.83%配方

7、方法 = 1 subtracted from:配方方法 Lower Center Upper -+-+-+-+-2 -52.6 185.3 423.1 (-*-)3 -275.1 -37.2 200.6 (-*-)4 -542.6 -304.8 -66.9 (-*-) -+-+-+-+- -400 0 400 800第一组数据结论：“配方法1配方法4”。原因分析：配方法2减去1的样本均值为185.3，但总体的均值差不确定，以98.83%的置信度可以判断，总体的均值差在（-52.6,423.1）内，置信区间包含0,即配方法1与配方法2总体的均值差不排除为0，即两者抗折强度平均值无显著差异。配方法

8、3减去1的样本均值为-37.2，但总体的均值差不确定，以98.83%的置信度可以判断，总体的均值差在（-275.1,200.6）内，置信区间包含0,即配方法1与配方法3总体的均值差不排除为0，即两者抗折强度平均值无显著差异。配方法4减去配方法1时的样本均值为-304.8，总体的均值差异98.83%的概率在（-542.6，-66.9）内上下为负，得到配方法1的均值肯定超过配方法4的均值的结论。第二组数据配方法 = 2 subtracted from:配方法 Lower Center Upper -+-+-+-+-3 -460.3 -222.5 15.3 (-*-)4 -727.8 -490.0

9、-252.2 (-*-) -+-+-+-+- -400 0 400 800第二组数据结论：均值“配方法2配方法4” 原因分析：配方法2与配方法3二者的抗折强度的平均值无显著差异，配方法4减去配方法2时的样本均值为-490.0，但总体的均值差异98.83%的概率在（-727.8，-252.2）内上下为负，配方法2的均值肯定超过配方法4的均值， 第三组数据配方法 = 3 subtracted from:配方法 Lower Center Upper -+-+-+-+-4 -505.3 -267.5 -29.7 (-*-) -+-+-+-+- -400 0 400 800第三组结论：均值“配方法3配方

10、法4”分析：配方法4减去配方法3时的样本均值为-267.5，但总体的均值差异98.83%的概率在（-505.3，-29.7）内上下为负，配方法3的均值肯定超过配方法4的均值。 总结论：“配方法1配方法4，配方法2配方法4，配方法3配方法4”，三组均值差异显著，其余各组间无显著差异。2、纺织厂有很多织布机，设每台机器每分钟织出同样的布，为了研究这一假设，随机选取5台织布机并测定它们在不同时间的产量，得出数据：织布机产量114.014.114.214.014.1213.913.813.914.014.0314.114.214.114.013.9413.613.814.013.913.7513.81

11、3.613.913.814.0（1） 说明为什么这是一种随机效应实验。织布机的产量相等吗？（2） 估计织布机间的变异。（3） 估计实验的误差方差。解答： （1）随机选取5台织布机并测定它们在不同时间的产量，则该实验是一次随机效应实验。然后进行产量分析：原假设：产量相等。备择假设：产量不相等。使用Minitab对数据进行分析：One-way ANOVA: 产量 versus 织布机 Source DF SS MS F P织布机 4 0.3416 0.0854 5.77 0.003Error 20 0.2960 0.0148Total 24 0.6376S = 0.1217 R-Sq = 53.5

12、8% R-Sq(adj) = 44.29%结论：p=0.030.05,所以拒绝原假设，每个织布机的产量不相等。（1） 估计织布机间的变异即估计数据方差分量，使用Minitab进行分析：Variance Components % ofSource Var Comp. Total StDev织布机 0.014 48.82 0.119Error 0.015 51.18 0.122Total 0.029 0.170 Expected Mean Squares1 织布机 1.00(2) + 5.00(1)2 Error 1.00(2)结论：织布机之间变异的方差为0.014（3）估计实验的误差方差，数据分

13、析：Variance Components % ofSource Var Comp. Total StDev织布机 0.014 48.82 0.119Error 0.015 51.18 0.122Total 0.029 0.170结论：实验的误差方差为0.0153、电视机厂感兴趣于对彩色显像管四种不同的涂层对显像管的电导率是否有影响。测得电导率的数据如下：涂层电导率1143141150146215214913714331341361321274129127132129（1） 涂层使电导率有差异吗？（2） 估计总均值与处理效应。（3） 计算涂层4的均值的95%区间估计。计算涂层1与涂层4之间的均

14、值差的99%区间估计。（4） 假定现在采用的是涂层4，你会向工厂推荐那种涂层（使电导率最小）？。解答：（1）原假设：彩色显像管四种不同的涂层对显像管的电导率没影响备择假设：彩色显像管四种不同的涂层对显像管的电导率有影响使用minitab对数据进行分析：One-way ANOVA: 电导率 versus 涂层 Source DF SS MS F P涂层 3 844.7 281.6 14.30 0.000Error 12 236.3 19.7Total 15 1080.9S = 4.437 R-Sq = 78.14% R-Sq(adj) = 72.68%分析结果：p值=0.0000.05，拒绝原假

15、设，彩色显像管四种不同的涂层对显像管的电导率有影响 (2)每个涂层的均值1,2,3,4分别如下所示：Individual 95% CIs For Mean Based on Pooled StDevLevel N Mean StDev -+-+-+-+-1 4 145.00 3.92 (-*-)2 4 145.25 6.65 (-*-)3 4 132.25 3.86 (-*-)4 4 129.25 2.06 (-*-) -+-+-+-+- 126.0 133.0 140.0 147.0分析结果: 样本总均值为137.94,根据均值结果和处理效应的计算公式每个涂层的均值减去总均值,得到处理效应1

16、、2、3、4分别为7.06,7.31，-5.69，-8.69。 （3）对涂层4数据进行正态性检验，结果如下图所示： 分析结果：p值=0.3350.05，该样本服从正态分布，且总体方差未知，所以使用单样本T检验求置信区间。结果如下：Variable N Mean StDev SE Mean 95% CI涂层4 4 129.25 2.06 1.03 (125.97, 132.53)得到涂层4的95%CI为（125.97,132.53）涂层1与涂层4之间的均值差的99%区间估计采用单样本t检验，对二者差值d进行正态性检验，结果如下图所示分析结果：p值=0.1620.05,差值服从正态分布，采用单样本

17、t检验，结果如下：Variable N Mean StDev SE Mean 99% CI涂层1-涂层4 4 15.75 2.06 1.03 (9.73, 21.77)结论：涂层1与涂层4之间的均值差的99%的区间估计为（9.73,21.77）（4）根据箱线图得到涂层4均值小且稳健性强。 4、研究化学过程的产率。两个最重要的变量有可能是压强与温度，每一因素选取三个水平，进行有两次重复的实验，产率数据如下： 温度压强200215230低90.490.790.290.290.690.4中90.190.589.990.390.690.1高90.590.890.490.790.990.1（1） 分析这

18、些数据并做出结论。（产率当然越高越好）（2） 论述模型的适合性。（3） 画出主效应及交互作用图。（4） 在什么条件下运行这一过程？解： （1）采用双因素固定效应模型方差分析的方法，在Minitab中对数据进行分析：Two-way ANOVA: 产率 versus 温度, 压强 Source DF SS MS F P温度 2 0.30111 0.150556 8.47 0.009压强 2 0.76778 0.383889 21.59 0.000Interaction 4 0.06889 0.017222 0.97 0.470Error 9 0.16000 0.017778Total 17 1.2

19、9778S = 0.1333 R-Sq = 87.67% R-Sq(adj) = 76.71%分析：温度和压强的p值分别为0.009和0.000，均小于0.05，因此温度和压强对产率的影响显著温度和压强的交互作用p值=0.4700.05，得出温度和压强的交互作用不显著。 （2）双因素固定效应模型方差分析，其一般模型如下： 是总体平均效应；是行因素A第i水平的效应；j是列因素B第j水平效应；是与j之间的交互作用的效应；是随机误差分量。两个因素最初假设是固定的，处理效应规定为跟总平均的偏差，所以同样，交互作用效应是固定的并限定其满足。试验有n次重复，所以共有abn个观察值。将总平方和进行分解可以得

20、到SST= SSA +SSB +SSAB +SSE。对应项除以其自由度可以得到均方的每一比值MSA/MSE ，MSB/MSE ，MSAB/MSE服从F分布，其分子的自由度分别a-1,b-1,(a-1)(b-1)，分母的自由度是ab(n一1)，而临界区域就是F分布的上尾部，检验方法通常概括在一张方差分析表中，可以得到一个方差分析表。可以根据得到的p值判断F值的极端程度，进而得到结论，因此模型合适。（3）主效应如下： 主效应图 交互作用图如下： 交互作用图（4）对上题主效应图分析得：在温度高，压强为215的情况下有最大产率， 在温度高，压强为215的情况下运行这一过程。5、工程师推测金属部件的表面

21、光洁度（越高越好）受进料速度和切割深度的影响。选取三种进料速度和随机选取四种切割深度，每种组合三次重复。数据如下：进料速度切割深度0.150.180.200.250.2074798299646888104607392960.2592989910486104108110888895990.309910410811498991101111029599107（1） 分析这些数据并写出结论。（2） 作出适当的残差图并论述模型的适合性。（3） 求每种进料速度的平均表面光洁度的点估计。（4） 估计切割深度的方差分量。解： （1）采用一般线性模型，在Minitab软件中对题中数据进行分析：分析结果如下： G

22、eneral Linear Model: 表面光洁度 versus 进料速度, 切割深度 Factor Type Levels Values进料速度 fixed 3 0.20, 0.25, 0.30切割深度 fixed 4 0.15, 0.18, 0.20, 0.25Analysis of Variance for 表面光洁度, using Adjusted SS for TestsSource DF Seq SS Adj SS Adj MS F P进料速度 2 3160.50 3160.50 1580.25 55.02 0.000切割深度 3 2125.11 2125.11 708.37 2

23、4.66 0.000进料速度*切割深度 6 557.06 557.06 92.84 3.23 0.018Error 24 689.33 689.33 28.72Total 35 6532.00S = 5.35931 R-Sq = 89.45% R-Sq(adj) = 84.61%结论：进料速度和切割深度以及他们的交互作用p值分别为0.000和0.000和0.018，均小于0.05，所以可知这些因素对于表面光洁度都有显著的影响。（2）残差图如下：模型合适性分析：“残差与观测值顺序”图中可残差值随机地在水平轴上下无规则的波动，即残差值间是相互独立的；“残差与拟合值”图形既不是“漏斗形”也不是“喇叭

24、形”，说明拟合效果很好，残差是等方差的；“正态概率图”中残差值都分布在直线附近，残差是服从正态分布。残差图分析可知：模型拟合性很高（3）每种进料速度的平均表面光洁度的点估计：1 当进料速度为0.20时的分析结果：One-Sample T: 进料速度0.20Variable N Mean StDev SE Mean 95% CI进料速度0.2 12 81.58 14.34 4.14 (72.47, 90.69) 结论：进料速度为0.20的时候平均表面光洁度的点估计为81.58。2 当进料速度为0.25时的分析结果：One-Sample T: 进料速度0.25Variable N Mean StD

25、ev SE Mean 95% CI进料速度0.25 12 97.58 8.01 2.31 (92.50, 102.67)结论：进料速度为0.25的时候平均表面光洁度的点估计为97.58。3 当进料速度为0.30时的分析结果：One-Sample T: 进料速度0.30Variable N Mean StDev SE Mean 95% CI进料速度0.3 12 103.83 6.07 1.75 (99.97, 107.69) 结论：进料速度为0.30的时候平均表面光洁度的点估计为103.83。（4）估计切割深度的方差分量：使用指令显示方差分量进行分析： Variance Components, using Adjusted SSEstimatedSource Value进料速度 123.95切割深度 68.39进料速度*切割深度 21.37Error 28.72结论：切割深度方差分量为68.39。