Minitab两因素方差分析续[专业技术]

1、使用双因子方差分析 (ANOVA) 过程可在存在两 个固定因子时固定因子时检验总体平均值的相等性。此过程 要求因子水平每一组合的观测值数必须相同（平 衡）。 仅当需要拟合可加性模型（Fit additive model） （无交互作用项的模型）时，其中一个或这两个 因子才可以为随机值。 双因子方差分析过程不支持多重比较。 注：如果数据平衡，且您需要检查涉及随机因子 的交互作用，那么可以使用统计 方差分析 平 衡方差分析。如果需要使用多重比较对平均值进 行比较，或者如果数据不平衡，那么可以使用统 计 方差分析 一般线性模型。 注解注解1：关于平衡两因素和平衡设计方差分析的区别：关于平衡两因素和平

2、衡设计方差分析的区别 1高级教学 注解注解2：关于平均值分析：关于平均值分析 平均值分析的英文缩写 ANOM 是看上去像方差分 析的英文缩写 ANOVA，平均值分析可检验总体 平均值的相等性。 Minitab 显示的图形类似于控制图，该图显示因子因子 的每个水平的平均值的每个水平的平均值如何与总体平均值（也称为 总均值）进行比较。Minitab 对与总体平均值显著 不同的平均值进行标记。因此，平均值分析可以 说明水平平均值何时不同以及差异是什么水平平均值何时不同以及差异是什么。 通过方差分析，如果可以假定响应大致按正态分 布，那么可以使用平均值分析。另外，当响应由 比率（二项数据）和计数（Po

3、isson 数据）组成 时，可以使用特殊的平均值分析版本。使用二项 数据时，样本数量 (n) 必须为常数。 2高级教学 均值分析图示例 Minutes Strength 181510 321321321 2 0 -2 效效应应 -1.578 1.578 0 181510 7 6 5 M Mi in nu ut te es s 平平均均值值 5.300 7.145 6.222 321 8 6 4 2 S St tr re en ng gt th h 平平均均值值 5.300 7.145 6.222 D De en ns si it ty y 的的双双因因子子正正态态平平均均值值分分析析 Alph

4、a = 0.05 交交互互效效应应 M Mi in nu ut te es s 的的主主效效应应S St tr re en ng gt th h 的的主主效效应应 3高级教学 图例分析 使用平均值分析的主效应图可检验“每个因子的 水平平均值等于指定 a 水平时的总体平均值”这 一假设。Minitab 为双因子设计中的每个因子显示 一个主效应图。主效应图显示： 标绘点 每个因子水平中的样本平均值。 中心线（绿色） 总体平均值。 决策的上限和下限（红色） 用来检验此假设。 Minitab 查找位于决策限之外的样本平均值，并用 红色符号对其进行标记。 4高级教学 如果样本平均值超出决策限，那么可以否

5、 定“平均值等于总体平均值”这一假设。 如果样本平均值未超出决策限，那么不能 否定“平均值等于总体平均值”这一假设。 5高级教学 注解注解3：等方差检验：等方差检验 Bonferroni 置信区间 Bonferroni 置信区间使用全族误差率。假设 该过程的全族置信水平为 95%。全族误差 率等于 1 - 置信水平 = 1 - 0.95 = 0.05。 Bonferroni 法通过将全族误差率分割在各个 区间之中。假设有六个区间。将每个区间 的单个误差给定为 0.05 / 6 = 0.00833，计 算单个置信水平 1 - 0.0083 = 0.9917。由 于置信水平较大 (0.9917)，

6、因此单个区间 通常相当宽。这种方法使得一个或多个置 信区间不能覆盖其相关总体标准差的概率 最多为 0.05。 6高级教学 与单元（配对因素）对应的总体标准差的 点估计值是指该单元中观测值的样本标准 差。一个单元至少要有两个观测值来计算 样本标准差。如果没有，那么该单元的点 估计值在输出中为空白。 标准差的置信区间以卡方分布为基础。此 分布为非对称，因此，置信区间也是非对 称的。 7高级教学 95% 标准差 Bonferroni 置信区间 方法 类型 经验 N 下限 标准差 上限 1 0 4 2.80384 5.88784 40.4990 1 1 4 1.84435 3.87298 26.640

7、0 2 0 4 2.26721 4.76095 32.7478 2 1 4 1.98261 4.16333 28.6371 3 0 4 2.88359 6.05530 41.6509 3 1 4 2.42820 5.09902 35.0732 示例示例 8高级教学 示例注解：示例注解： 标准差的 Bonferroni 置信区间显示以下内容： 公路类型：第一个因子。 经验：第二个因子。 N：单元中的观测值数。例如，在六个因子水平 组合的每一单元中有四个观测值。 下限和上限：为每个 sigma给定的 95.0% 置信区 间时的下端点值和上端点值。每个区间提供对应 单元的总体标准差的一个估计值。例如

8、，区间 (2.80384, 40.4990) 为公路类型 = 1 和经验 = 0 估计总体标准差。根据此区间， sigma介于 2.80384 与 40.4990 之间。 9高级教学 注解注解4：minitab方差齐性检验方差齐性检验 Minitab 显示了用于判断方差是否相等的两种检验 的结果：Bartlett 检验和 Levene 检验。在两种检 验中，原假设 (Ho) 是考虑的总体方差（或等效的 总体标准差）相等，备择假设 (H1) 指并非所有的 方差都相等。 检验的选项取决于分布属性： 当数据来自正态分布时使用 Bartlett 检验。对于 偏离正态性的情况，Bartlett 检验的功

9、能并不强大。 当数据来自连续但不一定正态的分布时，请使用 Levene 检验。 10高级教学 注解注解5：主效应图：主效应图 将主效应图与方差分析一起关联使用。当 平均响应值跨因子水平而更改时，主效应 随即出现。使用此图 检查每个因子的水平平均值 比较多个因子的水平平均值 11高级教学 具有多个因子时，主效应图将是最佳选择。可以 将水平平均值中的更改进行比较，以查看哪些因 子对响应（反应变量）的效应最大。某一因子的 不同水平对响应具有不同效应时，便会出现主效 应。对于有两个水平的因子，可能会发现一个水 平会提高平均值，而另一个水平则不然。这种差 异就是主效应。 Minitab 通过绘制每个因子

10、水平的平均响应值创建 主效应图。以线连接每个因子水平的各个点。 Minitab 还在总体平均值处绘制了一条参考线。查 看此线可以确定对某个因子是否存在主效应。 12高级教学 当线为水平时（与 x 轴平行），则不存在 主效应。因子的每个水平以相同的方式影 响响应，响应平均值在所有因子水平中相 同。 当线不水平时（与 x 轴不平行），则存在 主效应。不同因子水平对响应的影响不同。 标绘点之间垂直位置的差异越大（线与 X 轴不平行的程度越大），主效应的量值就 越大。 13高级教学 Minitab 方差分析基础方差分析基础 寻找因素与反应变量关系式的方法论 一元配置分散分析(DATA形态为 Stack

11、 的时候) 一元配置分散分析(DATA形态为 Unstack 的时候) 二元配置分散分析 平均分析 平衡方差分析(在各水准反复相同的时候) 一般线型模型 支份分散分析 检定方差的同一性 区间 Plot 主效果 Plot 交互效果 Plot 14高级教学 Minitab One Way ANOVA(One Way ANOVA(单因素方差分析单因素方差分析) ) 因子为一个, 反复数为对所有水准不相同也可, Radom实验。 在数据为一个 Col中以 Stack 形态保存时使用。 Response:指定反应变量 Factor:指定说明变量(要因) Comparisons:检定多重比较 Store

12、residuals:保存残差 Store fits:保存水准平均值 DF:自由图(Degree of Freedom) SS:乘方的和(Sum of Square) MS:不偏分散(Mean of Square) F:F-概率值 P:P-value(留意概率) 留意水准比 p-value 大则有影响。 即水准间有差。 (级区间有变动) - 上面的 p值大于 0.05，故没有影响。 EXH_AOV.MTW (先需要检定 RESPONSE值的 正态性) 15高级教学 Graphs. Dotplots / Boxplots 图象输出 option Residual Plots:对残差提供多样的 pl

13、ot - 残差只有随正态性时，它的结果值才能 判断为正确。 存在各范围间的重叠区间各点呈现直线状态时，意味着正态性 Minitab One Way ANOVA(One Way ANOVA(单因素方差分析单因素方差分析) ) 16高级教学 当数据按水准类别指定在 Col 时使用(Unstack 形态) 剩余事项与 Stack 情况相同 Responses:指定按各水准别 有反应值的Col Minitab One Way ANOVA(Unstacked) One Way ANOVA(Unstacked) 17高级教学 因子为 2个，把因子各水准的组合全部Radom实施的实验。 数据应为 Stack

14、 形态。 Response:实验结果数据 Row factor:B因子 Column factor:A因子 Store residuals:保存残差 Fit additive model:选择交互作用的有无 Lake与 Interaction 的 p值 大于 0.05，故不会 引起效果。 Suppleme的 p值 小于 0.05，故 Suppleme 的 水准间有差。 看左图可知道 Suppleme 的平均间有差。 看左图可知道 Lake 的平均间没有差。 EXH_AOV.MTW Minitab Two-way ANOVATwo-way ANOVA（两因素方差分析两因素方差分析） 18高级教学

15、 用 Graph 来显示因子的平均值，检讨因子的哪个水准有影响 -方差分析是对水平间有无差别的分析 -平均分析是对全体平均与各水平平均间有无差 别的分析 Response:反应(结果)值 Distribution of Data:资料的分布形态 -Normal:正态分布, Factor 1:因子水准 Col (单因素) Factor 2:因子水准第二 Col (两因素) -Binomial:二项分布 -Poisson:Poisson分布 Alpha level:留意水准 脱离管理线则有影响 用两个因子的交互作用效果 Main Effect:主要因 Minutes 的 3水平(值=18)时有影响

16、 Strength 的 3水平(值=3)时有影响 EXH_AOV.MTW Minitab Analysis of MeansAnalysis of Means（均数分析）（均数分析） 19高级教学 Minitab Balanced ANOVABalanced ANOVA（平衡设计方差分析）（平衡设计方差分析） 所有单元的观察个数相同时使用 Response:反应变量数据 Model:指定需分析的因子 Random factors:指定变量因子 Probtype|Calculat的标记为考虑交互作用 效果的计算实施. EXH_AOV.MTW Probtype, Calculat, Probtyp

17、e*Calculat 等比留意水准(0.05) 小，故判断为 各因子的水准间存在散布的差。 Engineer 为变量因子故无统计意义。 20高级教学 Minitab Test for Equal VariancesTest for Equal Variances（等方差检验）（等方差检验） 检定2总体以上的方差是否一致 - 原假设 : 所有水平的方差一致 - 对立假设 : 至少一个以上的方差不一样 正态分布数据时:Bartletts Test 包括正态分布的连续性数据时:Levenes Test 因 p-value 比留意水准(0.05)大，故选择归属假设，即所有水平的方差一致。 EXH_AO

18、V.MTW 21高级教学 Minitab Main Effects PlotMain Effects Plot（主效应图）（主效应图） 对主效应的水平间差异比较 Responses:指定反应值 Factors:指定因子 Base plots on:指定plot基准 Supplement 在2水平时值特大。 Lake在各水准间无太大的变动。 EXH_AOV.MTW 22高级教学 Minitab Interactions PlotInteractions Plot（交互效应图）（交互效应图） 交互作用的水平间差异比较 Display full interaction plot matrix: 作成

19、为 matrix 可知道按 Field 水准变更的 Variety 各水准的 变动及平均值。 -平均是 Variety 4,6水准比别的水准小。 -变动是 Variety 2 水准比别的水准大。 -水准间 Cross 角度越大，交互作用效果就 越大。 ALFALFA.MTW 23高级教学 DOE(实验设计实验设计) Minitab 24高级教学 正交设计 多因子试验与正交表多因子试验与正交表 无交互作用情况下的正交设计无交互作用情况下的正交设计 有交互作用情况下的正交设计有交互作用情况下的正交设计 裂区法裂区法 多指标的数据分析多指标的数据分析 Minitab 25高级教学 多因子试验与正交表

20、 Minitab 26高级教学 Minitab 27高级教学 Minitab 28高级教学 Minitab 29高级教学 交互作用 Minitab 30高级教学 Minitab 31高级教学 Minitab 32高级教学 Minitab 33高级教学 正交表 Minitab 34高级教学 Minitab 35高级教学 Minitab 36高级教学 Minitab 37高级教学 Minitab 38高级教学 Minitab 39高级教学 无交互作用情况下的正交设计 Minitab 40高级教学 Minitab 41高级教学 Minitab 42高级教学 Minitab 43高级教学 Minita

21、b 44高级教学 Minitab 45高级教学 Minitab 46高级教学 数据的直观分析 Minitab 47高级教学 Minitab 48高级教学 Minitab 49高级教学 Minitab 50高级教学 Minitab 51高级教学 Minitab 52高级教学 数据的方差分析 Minitab 53高级教学 Minitab 54高级教学 Minitab 55高级教学 Minitab 56高级教学 Minitab 57高级教学 Minitab 58高级教学 Minitab 59高级教学 Minitab 60高级教学 Minitab 61高级教学 Minitab 62高级教学 Minit

22、ab 63高级教学 Minitab 64高级教学 Minitab 65高级教学 Minitab 实验设计基础实验设计基础 如何实施实验如何选取数据, 如何解释才能以最少的实验次数 迅速获得最大的信息量的计划方法. 实验的成败，只有把以往的经验或者理论性、 技术性知识等的原有技术与 依照实验计划法的知识结合起来才有可能. Create Factorial Design:析因实验设计 Define Custom Factorial Design:在变更当前的 实验设计而再指定时使用。 Analyze Factorial Design:得出实验分析结果 Factorial Plot:主效果, 交互效

23、果 plot 作成 Contour/Surface(Wireframe)Plots:展现实验的 反应表面 Overlaid Contour Plot:以视觉性展示多个反应 变量的妥协领域 Response Optimizer:寻找满足目标值因子的 最佳组合 Factorial:析因设计 RS Design:反应曲面设计 Mixture Design:混合试验设计 Modify Design:修正试验设计 Display Design:实验计划后生成的内容通过 Worksheet 可见 66高级教学 Minitab 创建田口设计创建田口设计 使用“创建田口设计”可以在工作表中设置田口正交表设计。 设计的每行都指定 一个用于试验游程的因子级别的组合。 田口设计用于稳健参数设计（试验设计的一种），其中的主要目标是在调整（或 保持）目标的同时，找出使响应变异最小化的因子设置。田口设计为设计在各种 条件下始终运行优良的产品提供了一种强大而高效的方法。