数理统计与随机过程ch9-方差课题.ppt

R程序运行结果 R源程序 juice-data.frame(X = c(76,67,81,56,51, 82,69,96,59,70, 68,59,67,54,42,63,56,64,58,37), A = gl(4, 5), B = gl(5, 1) ) juice.aov-aov(X~A+B, data=juice) summary(juice.aov) Df Sum Sq Mean Sq F value Pr(F) A 3 1182.95 394.32 10.722 0.0010329 ** B 4 1947.50 486.87 13.239 0.0002342 *** Residuals 12 441.30 36.77 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 解: 按题意需检验假设(2.25), (2.26)。 ， 的值已算出载于上表。现在 r=4，s=5。由(2.27)得到： 方差分析表如下： 由于F0.05(3,12)=3.4910.72,，F0.05(4,12)=3.2613.24，得： 拒绝H01及H02，即，认为不同时间下颗粒状物含量的均值有显著差异；也认为不同地点下颗粒状物含量的均值有显著差异。 即，认为时间和地点对颗粒物的含量影响均为显著。 三、小结 1. 双因素等重复试验的方差分析步骤 2. 双因素无重复试验的方差分析步骤 (3) 研究统计特性; (2) 分解平方和; (4) 确定拒绝域. (1) 建立数学模型; (3) 研究统计特性; (2) 分解平方和; (4) 确定拒绝域. (1) 建立数学模型; * * * * * * * * * * * * * * * * * * * * * * * * * * * * 记号 总平均 双因素试验方差分析的数学模型 检验假设 2. 研究统计特性; 检验步骤 1. 分解平方和; 3. 确定拒绝域. 1.分解平方和 总偏差平方和(总变差) 误　差 平方和 因素 A 的 效应平方和 因素 B 的 效应平方和 因素A,B的交 互效应平方和 2.研究统计特性 自由度 数学期望 3.确定拒绝域 表9.9　双因素试验的方差分析表 方差来源 平方和 自由度 均　方 F　比 因素　A 因素　B 交互作用 误　　差 总　　和 记 * 对GMAT数据得出方差分析表 例1　一火箭用四种燃料,三种推进器作射程试验. 每种燃料与每种推进器的组合各发射火箭两次,得 射程如下（以海里计）. 在MATLAB中求解 函数:anova2 格式:p=anova2(x,reps) 说明:执行平衡的双因素试验的方差分析来比较x 中两个或多个列或行的均值.不同列的数据代表某 一因素的差异,不同行的数据代表另一因素的差异. 如果每行列对有多于一个的观察值,则变量reps指 出每一单元观察点的数目,每一单元包含reps行. Matlab源程序: a=[58.2,56.2,65.3;52.6,41.2,60.8; 49.1,54.1,51.6;42.8,50.5,48.4; 60.1,70.9,39.2;58.3,73.2,40.7; 75.8,58.2,48.7;71.5,51.0,41.4]; p=anova2(a,2) R源程序： juice-data.frame( X = c(58.2, 52.6,56.2,41.2,65.3,60.8,49.1,42.8,54.1,50.5,51.6,48.4, 60.1,58.3,70.9,73.2,39.2,40.7,75.8,71.5,58.2,51.0,48.7,41.4), A = gl(4, 6), B = gl(3, 2) ) juice.aov-aov(X~A*B, data=juice) summary(juice.aov) R程序输出结果 Df Sum Sq Mean Sq F value Pr(F) A 3 261.68 87.23 4.4174 0.025969 * B 2 370.98 185.49 9.3939 0.003506 ** A:B 6 1768.69 294.78 14.9288 6.151e-05