MiniTab16 操作教程
&&&&目录第一章 第二章 minitab概要 管理数据6s第三章第四章 第五章 第六章 第七章 第八章 第九章操作和计算数据使用数据分析和质量工具 基本操作示例 做一个简单分析 高级minitab 质量管理和改善 实验设计0前言<&&&&br />6sminitab 是为质量改善、教育和研究应用领域提供统计软件和服务的先导。是一个很 好的质量管理和质量设计的工具软件，更是持续质量改进的良好工具软件。 minitab 统计软件为质量改善和概率应用提供准确和易用的工具。minitab 被许多世 界一流的公司所采用，包括通用电器、福特汽车、通用汽车、3m、霍尼韦尔、lg、东芝、 诺基亚、以及 six sigma 顾问公司。 作为统计学入门教育方面技术领先的软件包，minitab 也被 4,000 多所高等院校所 采用。 minitab 总部位于state college，pa，usa（美国），在英国和法国设有办事处， 在世界各地拥有分销商。 minitab 包括:
基础和高级统计
回归和方差分析
演示质量的图表
模拟和分布
灵活的数据导入、导出和操纵
spc (statistical process control - 统计过程控制)
doe (design of experiments - 试验设计)
可靠性分析
多变量分析
样本量和幂计算
强大的宏语言 minitab — 系统要求
运行 windows(r) 95/98 或者windows nt(tm) 4，需要16mb 内存。
处理器为 486 或以上的个人计算机；如使用 486 处理器，特别建议采用数学协 处理器。 完全安装需要 35 mb 的硬盘空间
vga 或svga显示器，建议最小分辨率为800 x 600。
需要cd-rom 驱动器
某些功能需要使用鼠标 本 教 材 是为 了 配合 我们 lg 曙光 电 子 six sigma 活 动而 编 写的 ， 主 要着 重 于 对 minitab 软件包的基本操作和运用，没有深入该软件的介绍和运用，更没有跟随新版软件 的升级更新。这些需要使用者自己深入钻研。1第一章minitab 概 要6s1-1 1-2minitab系统操作环境简介 minitab的工作步骤1-3关于minitab的 project文件类型2minitab系统操作环境简介6s工具条在minitab里进行数据分析时，需要使用各种窗口和工具。下面是 minitab环境下各部分的简介。会话窗口 历史窗口数据窗口信息窗口状态条快捷菜单 图形窗口各minitab 窗口作用： 会话窗口（session window）显示诸如统计报表之类的输出文本 。
数据窗口（data windows）在此可以输入、修改数据和查看每个工作表的数据 列。
信息窗口（info window） 概括了每个打开的工作表。可以从下拉列表框中选 择要查看的工作表。
历史窗口（history window ）记录了所用过的命令。
图形窗口（graph windows） 显示各种图形。一次最多只能打开15个图形窗口 。3minitab的工作步骤6s进行一个典型的分析往往要经过很多步骤。以下六章讲的就是与此相关的基 本分析步骤。每章节主要讲了那一步的一些基本、常见的分析内容。以下是各章节的列表及内容简述。章节标题 容主要 内2数据管理 件中导入数据或 成规则数据。在数据窗口中输入和修改数据，从文 将数据以文件形式保存，生3操作和计算数据 在数据窗口中对数据列进行操作，排序、分 组和生成 方程式。 运用数据分析 的分析 从基本的统计到质量管理，运用 一系列4和质量工具5 图形数据方式。生成、修改图形，同时打印并且用各 种格式保存图形。 查看输出文本，并且改变输出文本格式， 打印 和用各种格式保存。6管理会话窗， 生成报告4打开、保存、关闭project文件6s new, 按minitab project, 再按 ok.
open project.
save project■打开一个新的project文件，选择file■打开一个已有的project文件，选择file ■打开一个已有的projectworksheet，选择file
open worksheet.■保存一个project文件，选择file■关闭一个project文件，你必须打开一个新的project文件，或一个已有的project文件，或者退出。 一个minitab project文件包含了所有的工作：收集的数据、输出的文本、图形等。 当保存一个project文件时，所有的东西都保存下来了。 但是一次只能打开一个project文件。project文件里的各种工作表（data、 work- sheet、graph等）可以单独进行处理---可以生成数据、图形，也可以从别的文 件中通过拷贝的形式向一个project文件添加数据和图形。project文件里的大部分窗 口的内容都能够单独保存和打印。也可以忽略某些工作表或图形，这样的话，这 些工作表和图形便不会保存到project文件中去。 每当保存一个project 文件时，保存了以下一些内容： 一. 所有窗口的内容，其中包括： 1. 每个数据窗口（data window）的数据列。 2． 信息窗口（info window）中的常量、矩阵。3．4.会话窗口（session window）和历史窗口（history window）中的所有的文本。 每个图形窗口（graph window）。二. 通过命令 fileproject description 生成的 project 文件注解。 三. 通过命令 editor worksheet descriptions 生成的对每个工作表的注 解。 四. 每个窗口的大小、位置和状态。5minitab的工作步骤在对话框中输入变量6sminitab的对话框和其他软件的对话框没什么两样。其最大的特点就是 变量列表框，出现在很多对话框中。变量列表框一般显示的是一些变量列名、常 量或矩阵。凡是在能够输入变量的文本框中，均可以输入这些列表框中所显示的。 变量列表框的优点是可以准确无误的选择变量名（一般指数据列名）。变量列表框中的内容即 当前工作表的列、矩阵变量文本框只容纳列。当光 标停于其中时，左边的变量 列表框仅显示数据表中的各 列6minitab的工作步骤6s在对话框中选择变量的几种方法用鼠标选择一个变量1. 点击需要填写变量的文本框。2．在变量列表框中用鼠标选择一个变量，然后双击此变量。用鼠标选择多个变量1．点击需要填写变量的文本框。 2．点击变量列表框。 3．点击某个变量，用鼠标拖画数个变量，或者按住ctrl键点击不连续变量。 4．点击 select 按钮。用键盘选择一个变量1．按 tab 键进入需要填写变量的文本框。 2．按 f2 键，使得变量列表框被激活。3．使用上下箭头键，选择变量。 4. 再按 f2 键，将所选择的变量便出现在刚才被激活的文本框中了。7第二章 数据管理(managing data)6s2-1 2-2 2-3数据管理概要 在数据窗口中输入数据 生成规则数据8数据管理概要6s1. 数据保存在工作表中在minitab里，与特定的数据集有关的所有数据都包含在工作表中 。一个 project 文件允许有许多工作表（工作表的个数取决于计算机的内存大小）。 一个工作表可以包含三种数据类型-----数值型（numeric）、文本型（text）和 日期/时间（date/time）型，表现形式为：数据列（columns）、常量（constant）、 矩阵(matrices)。可以在多个窗口中察看数据，但大多时候都是在数据窗口中处理数据列。2. 数 据 的 三 种 类 型 ： 数 值 型 (numeric ) 、 文 本 型 (text) 、 日 期 / 时 间 型 (date/time)数值型数据由数字构成，如：1，2，3，5。 文本型数据由以下几种构成:字母、数字、空格和特殊字母如“test number 4” 日期/时间型数据可以是日期（如jan-1-2000 或者3/13/1999），时间（例如： 09:30:22am),也可以是二者均有（如：5/13/:25 am).minitab 将日期 和时间以数字的形式存储，但是可以选择任意一种格式来显示。3. 与数据有关的三个窗口a.数据窗口（data window） 数据窗口中的内容即为工作表中的数据列。当打开了多个工作表时，每个工 作表都有其自己的数据窗口。9数据管理概要在每个数据窗口里，可以：
察看处于工作表中的各列数据。
输入数值并以多种方式修改。
对数据列进行多种操作，包括：改变格式、字体、名字、 宽度、注解和列的位置。 例: 打开file open worksheet data marketd.mtw6s列号列名字母t表明 是文本列字母d表明是 日期/时间列行号注意：数据窗口不是电子表格 尽管数据窗口有行和列，但它不像excel那样的电子表格，它的单元格数量是有限的。 而且单元格中的数值只能手工输入或者有命令生成，它不能包含有方程式， 它的值不会随其他单元格的值变化而变化。例如：想要c3=c1+c2,首先必须用计 算器计算c1+c2的值。如果改变c1的值，若不用计算器再计算的话，c3的值是不 会改变的。换而言之，c3不随c1变而变。10数据管理概要6sb.信息窗口（info window） 此窗口对当前活动工作表的数据（数据列、常量、矩阵）进行了概括。可以对这 些数据进行修改。如果要察看另外一个工作 表的情况，可以从下拉列 表框中选择一个工作表。如果是文本列，标志为t, 日期列标志为d,其余列为 数值型的列c.会话窗口（session window）1. 2. 3. 4. 5. 6. 打开 file open worksheet data acid.mtw。 选择菜单：data display data， 在 columns, constants, and matrices to display 框中选入 acid1 acid2。 点击 ok 。 选择 window session 便能看见会话窗口中的如下内容：data display row 1 2 3 4 5 acid1 0.123 0.109 0.110 0.109 0.112 acid2 0.109 0.111 0.110 0.110 0.10511在数据窗口中输入数据6s在数据窗口中的一个单元格中输入数据，只要点击一下此单元格，然后输入 数据，再回车。每一列单元格一般表示一 个变量输入数据时，该列的数据会自动格式化。当往一空列输入数据时，minitab会 给该列分配一种数据类型：数值型、文本型或者日期/时间型。如果此数据类型不是数值型，minitab也会在列号后面添加一个标识符加以区分：d表示日期/时间型 数据，t表示文本型数据。12在数据窗口中输入数据6s一、 在一行、一列或者一块单元格中输入数据1) 输入一列数据 a. 点击数据方向箭头使之朝下 b. 输入数据，然后按tab或者enter键移动当前活动单元格。按ctrl+enter 组合键，当前活动单元格便跳到了下一列的顶部。数据方向箭头注意：输入完一个值回车之 后，当前活动单元格往下移 动一格。2) 输入一行数据 a. 点击数据方向箭头使之朝右b.输入数据，然后按tab或者enter键移动当前活动单元格。按ctrl+enter 组合键，当前活动单元格便 跳到了下一行的顶部。数据方向箭头注意：输入完一个值回车之 后，当前活动单元格往右移 动一格。13在数据窗口中输入数据6s3).输入一块数据 a. 选择一块区域，使之高亮度显示。 b. c. 输入数据，当前活动单元格仅仅在所选区域内移动。 如果要取消所选区域，可以按箭头键或者用鼠标点击数据窗口的任何地 方。二、生成规则数据可以按照一定的规则生成数据，如：按顺序生成一简单数据集或按字母顺序 生成一列数据。也可以生成规则的日期/时间型数据。 最常见的就是在一列中按照一定的顺序生成数据。如：从1到100的所有整数， 或者介于1到50之间的所有偶数。尤其是这一列中的 每个数可以重复，或者整个 队列的数据也可以重复。在变量设计分析中，当输入因子水平时，将数据或序列 反复多次是非常有用和简便的。14在数据窗口中输入数据6s1).生成简单的数字集 1. 选择菜单：calc
make patterned data
simple set of numbers.2. 在store patterned data in 的文本框中，输入数据列名（此列既可以是新 的列，也可以是已经存在的列），此列用来保存所生成的数据。 3. 在from first value(起始值）,to last value（最终值）和in steps of （步 长，即增量）文本框中输入数字。 4. 若有必要，也可以在list each value (每个值显示的次数)和list the whole sequence (显示整个编列的次数)文本框中输入一个数字。 2). 生成规则数据的几个例子结 果 from first value 1 10 -0.5 1 1 to last value 100 1 -0.1 4 5 in steps of 1 2 0.1 1 1 list each value 1 1 1 5 2 list the whole sequence 1 1 1 1 2从1到100的所有整数 从10到1的所有偶数 -0.5,-0.4,-0.3,-0.2,-0.1 3445515第三章 操作和计算数据6s3-1 3-2 3-3 3-4 3-3改变列的数据类型和格式 归类数据并生成相应的子工作表 叠加列 对数据进行替换 使用计算器16改变列的数据类型和格式6s进行分析之前，有时需要重新组织和安排输入到minitab里的数据，如： 移动或者删除行和列，将一种数据类型转换成另一种数据类型，或者控制数据 窗口中数据的显示方式。还可以基于原始变量之上再创建新的变量。可以合并两列，创建列的子列。一、改变列的数据类型和格式数据有三种类型：数值型、文本型、日期/时间型。一列只能包含一种类型 的数据。可以给某一空列分配一种数据类型并且改变某一存在列的数据类型。 一旦某一列已经有了一种数据类型的话，还可以为这列指定一种显示格式。 当对此列的格式进行修改时，仅仅是改变了它的显示形式而已，并没有改变其真实值。例如：某一单元格有一个数：1.234567。如果你将其现实形式改成为只显示 两位小数，数据窗口中的单元格会显示1.23，在图形窗口上那一点的显示也会为 1.23，但是所有的计算都仍会使用1.234567，会话窗口的显示形式则由其所使用 的分析命令语句决定，而不是数据窗口中的那种显示格式。1). 改变列的数据类型
给一空列分配一种数据类型 a.选择editor
format column b.选择numeric,text,或data/time改变一非空列的数据类型 a. 选择manip
change data type b. 选择要转换的数据类型.. c. 填完对话框，然后点击ok.17改变列的数据类型和格式6s2). 改变数值型和日期/时间型的格式
改变数值列中数据的小数位数 a. 选择一列或多列。这些列必须是空列（无格式）或者以数值型格式存在。 b. 选择editor
format column
numeric c. 选择fixed decimal width. 在decimal places文本框中， 输入小数的位 数， 然后点击ok. 改变日期/时间列的数据显示形式 单个的日期/时间数值可以是日期， 时间或者二者兼有。入下列数 据 是有效的 格式： 1/1/96 ， 3：04pm, 5/23/pm a. 选择一列或者多列。这些列必须 是空列（无格式）或者以日期/时 间型格式存在。 b. 选择editor
format column
date/time. c. 从date/time column format 文本框中，选择一种格式， 然后点击ok.18归类并且生成相应的子工作表6s二、将数据按不同类别归类并且生成相应的子工作表在一个比较大的数据集里面，经常需要对一组数据进行分析或生成图形。如：在 某项调查中对女性聚焦，或聚焦于某一季度的销售额。对这些你感兴趣的数据，minitab能够另外生成一个子工作表来包含这些数据。 可以基于某种条件来分类，或者对于指定的变量的所有组，分割一个工作表。 1). 按条件分类 例:打开 file open worksheet data market.mtw。 可以基于一个或多个条件对数据分类。如：按两种条件创建一个子集— 季度销售额大于100并且季度广告费用小于15。符合这两个条件的数据行将复到 一个新的工作表，原来的工作表仍然没有变动。 按条件分类 a. 选择manip
subset worksheetminitab自动为新工作表给一 个缺省名，但是可以随意改 动。按condition19归类并且生成相应的子工作表6sb.在specify which rows to include/exclude（指定哪行包括或者不包括）下 面，选择rows that match （匹配的行），然后点击condition(条件），弹出一 对话框，如图：c. 在 condition 一栏，输入条件方程式。例如：想要创建只包含销售额大于100 的子集，只要输入 sales & 100。在此，可以根据条件，选择 一系列数学函数 注意： 如果在subset主对话框选择了exclude rows的话，上面的方程式所生成的子集包 含的数据将会是销售额不大于100的数据。20归类并且生成相应的子工作表6s2). 将一个变量按组分割成工作表 在一个变量里，可以根据所有关键字将一个工作表分割成几组。例如：假设 有一个关于销售额的工作表，有两列，sales这一列为一个财政季度所有的收入；quarter这一列即为季度（1，2，3，4）。如果你是基于变量quarter来分割工作表，minitab会生成四个工作表：第 一季度的所有数据，第二季度的所有数据，等等。原来工作表内容不变 。 将一个变量按组分割成工作表 a. 选择manip
split worksheet如 果 选 中 了 include missing as a by level,假如在by variables 文本框中的任 意列含有空值，minitab会生成另外 一个工作表。b. 在by variables一栏，输入一列或多列，这些列即为分组条件。然后点击 ok.如果选择了多列，对于这些列的每一种组合所得的数据，minitab都会为其 生成一个新的工作表。 例如：假设year这一列有两个数据：；quarter这一列有四个：1， 2，3，4。如果在by variables中选择了year 和quarter，minitab 会生成8个工作表： 1991年每个季度一个工作表,1992年每个季度一个工作表。21叠加列6s有时候需要将两列合并成一列以便用一条命令对其分析。minitab能够快速 容易的将一列数据叠加于另一列之上。你可以将叠加后的数据用另一列存储而 使原来的列保持不便。当叠加列的时候，可以生成一个脚标列标明叠加后的数据是来自于哪一列。脚标还可以用于以后重新分类数据或者生成图形（在图形里数据点因来源于不同的组，显示有所不同）或者将叠加后的数据还原。 叠加列 a. 选择 manip
stack/unstack
stack输入一新列或者已经存在 的列来存储脚标。b. 在 stack the following columns 中，输入想要叠加的列。输入的第一列 将会处于第二列之上，第二列叠加在第三列之上。 c. 在 stock the stacked data in 一栏，输入一新列或者已经存在的列。22对数据进行替换6s可以将一个值转换成另外一个值，或者是将一定范围的值转换成另一个值， 也可以将一种数值型的数据转换成另一数值型的数据（例如：把1.0到1.9之间的所有数变成1），文本转换成文本，文本转换成数字或者数字转换成文本。 把数字型数据编码为文本型数据 a.选择 manip
numeric to textb. 在 code data from columns 一栏，输入一列或者多列。 c. 在 into columns，输入一个或者多个新的或已经存在的列。例如：假设没有叫做 gender的列，当输入“gender”时，便会创建一个名称为 gender 的列。 d. 在位于 original values 下方的第一个文本框中，输入一个数值，或者一定范围的 值（如1:12表示从1到12). e. 在 new 下方的第一文本框中，输入想要转换成的值。 f. 另外，由此还可看出，一次可以转换8个数值。23使用计算器6s使用计算器可以快速的 进行基本的算术运算和复杂的数学函数计算。结果 可以保存在某一列中。使用计算器a. 选择 calc
calculator 。 b. 在 store result in variable 中，输入保存结果的列。 c. 在 expression 中选择变量和各自的函数。也可以通过手动输入表达式。例子：使用计算器对列进行计算 1.例:打开 file open worksheet data
pulse.mtw 2. 选择 calc
calculator 3. 在 store result in variable 中，输入 pulsedif。 4. 在 expression 中，输入算式 pulse2-pulse124使用计算器6sc1 ↓ 1 2 3 4 5 6 7 pulse 1 64 58 62 66 64 74 84…5. 点击 ok 后结果如所示：c2 pulse 2 88 70 76 78 80 84 84…c9 pulse dif 24 12 14 12 16 10 0…25第四章使用数据分析和质量工具6s4-1 4-2 4-3 4-4 4-5 4-6 4-7 4-8概要 描述性统计 置信区间和平均值测试 相关 回归 分散分析 表 质量控制26概要6sminitab提供了许多统计和图表技法来进行数据分析。常用的方法有：
基本统计 回归方差分析实验计划法 管理图 质量工具 可靠度/残差分析 多变量分析 时间序列 表 非参数分析探测数据分析本章将会根据一个例子进行工序能力分析，其中包括：
基本统计——描述性统计，一样本 t 检验，相关性 回 表 归——线性回归 ——交叉表
方差分析——one-way anova
管 理 图——x-r图27描述性统计6sminitab提供了两条命令——display descriptive statistics 和store descriptive statistics, 这两条命令对每一列进行了计算并且将统计结果存储下来了。这些结果可以 从会话窗口中显示出来。 显示描述性统计结果 a. 选择 stat
basic statistics
display descriptive statistics b. 在 variables文本框中，输入想要进行分析统计的数据列. c. 另外，还可以选择by variables(具体作用参见下例）和graph按钮。 一个描述性统计的例子 假如想要对参加脉搏测试的男(sex=1) 女(sex=2)学生进行统计。1.打开文件：pulse.mtw 。 2.选择 stat
basic statistics
display descriptive statistics 。 3.在 variables 一栏，输入 weight 。 4.选中 by variables 并且在文本框中 输入 sex ，然后点击 ok 结果如下：variable weightsex 1 2 sex 1 2n 57 35mean 158.26 123.80median 155.00 122.00trmean 157.61 123.74stdev 18.64 13.37 q1 145.00 115.00 q3 170.00 131.00variable weightse mean 2.47 2.26minimum 123.00 95.00maximum 215.00 150.00 结果解释 在会话窗口中的这段文字对体重进行了描述性的统计。因为选中了by variables 为 sex，所以分别对男生（sex=1)女生（sex=2）进行了描述性统计。从结果可以看出：男 生体重（mean意思为平均值）大于女生体重（mean=123.80)。
注意：1.by variables的意思是对每个关键值（如男和女）分别进行统计计算。 2.graph 按钮的作用是以图形 （直方图 、正态曲线的直方图 、 点 图 、盒 式 图 (boxplot)的方式表示统计结果。28置信区间和平均值测试6sminitab 提供了一些命令计算一样本或两样本的置信区间以及进行平均值检 验。工序能力包括z-test 值、一样本t-test 值、两样本t-test和配对t-test.minitab也 提供了概率估计和概率差别的方法。另外，当你不知道所收集的数据是否服从正态分布时，minitab又提供了计算置信区间和对中央值进行假设检测的方法。下面是一样本t置信区间和平均值假设检验（假设服从正态分布）的一个例子。 计算t置信区间和平均值检验 a. 选择 stat
basic 行单个分析。 c. 以下二者中任选一项：
若计算平均值的置信区间，选中 confidence interval。 statistics
1-sample t b. 在 variables 文本框中，输入包含样本的列，minitab将对样本的每一列进 若进行假设检验，选中 test mean。 d. 若想更进一步，还可以使用下面的选项： 为置信区间指定一个置信度。缺省值为95%。
指定一个原假设检验值。缺省值为：μ=0
通过选择小于(&)、大于(&)、不等于（≠），定义选择性假设。29置信区间和平均值测试6st置信区间的一个例子 假设根据测得的样本数据（脉搏测验），想要得到平均值95%的t置信区间。 1. 打开文件 pulse.mtw。 2. 选择 stat
basic statistics
1-sample t 3. 在 variables 文本框中，输入 pulse1，点击 ok。 session window 输出：t confidence intervals variable pulse1 n 92 mean 72.87 stdev 11.01 se mean 1.15 95.0 % ci (70.59, 75.15) 结果解释 根据输出结果，估计其平均值为72.87，有95%的把握认为真值落在区间（70.95，75.15）。30相关(correlation)6s使用回归对成对变量计算相关系数和相关性。 相关系数是衡量两个变量之 间线性关系的程度。其值介于-1与+1之间。如果一个变量随另外一个变量减少而 递增的话，那么相关系数为负值；相反的，如果这两个变量同时增加或同时减小，那么相关系数为正。计算两列的相关性 a. 选择 stat basic statistics
correlation b. 在 variables 文本框中，输入包含测量数据的列，minitab 便会计算两两 列之间的相关性。31相关(correlation)6s相关性的一个例子 假设要分析学生身高与体重的关系。 1.打开文件 pulse.mtw。2.选择 stat
basic statistics
correlation3.在 variables 文本框中，输入 height 和 weight .点击 ok。 session window 输出：correlations: height, weight pearson correlation of height and weight = 0.785p-value = 0.000 结果解释 根据输出结果，相关值（r=0.785，p-value=0.000）表明身高与体重正相关。32回归(regression)6sstat
regression
regression。response 文本框中,输入包含反应值(y)的列。 predictors 的文本框中,输入预测值(x)变量。回归分析用于找出反应值(y)与一个或者多个预测值（x)之间的关系。 怎样做一个线性回归a. 选择b. 在 c. 在33回归(regression)6s运行一个简单线性回归的例子 假设想要找出体重与身高的关系 1. 打开文件 pulse.mtw2. 选择 stat
regression
regression3. 在 response 文本框中，输入 weight。在 predictors 中，输入 height 。然后点击 ok 。regression analysis: weight versus height the regression equation is weight = - 205 + 5.09 height predictor coef se coef constant -204.74 29.16 height 5.7 t -7.02 12.02 p 0.000 0.000s = 14.79 r-sq = 61.6% r-sq(adj) = 61.2% analysis of variance source df ss ms f p regression 1
144.38 0.000 residual error 90
total 91 51284 unusual observations obs height weight fit se fit residual stresid 9 72.0 195.00 161.87 2.08 33.13 2.26r 25 61.0 140.00 105.86 3.62 34.14 2.38r 40 72.0 215.00 161.87 2.08 53.13 3.63r 84 68.0 110.00 141.50 1.57 -31.50 -2.14r r denotes an observation with a large standardized residual 结果解释 p值等于0，表明体重是身高的一个显著性因子，r2 等于61.6%表明此模型在 反应值中所占的比率。34分散分析(anova)6s分散分析（anova)把比较两样本总体均值的 t-test 扩充到比较多于两个均 值的test。 minitab的anova包括一次、二次或者更为复杂的anova模型、等变量检 验等。当按照单个变量进行分类时，一个一次anova检验总体均值的相等性。 下面我们举一个例子，说明当反应值数据在一列而区分总体的水平值在另一列 时，怎样进行数据分析。如果你在工作表中将每个总体的数据置于不同的列， 可以运用 one-way 命令。 怎样做one-way anova a.选择 stat anova
one-way b.在 response 中，输入包含反应值的列。 c.在 factor 中，输入包含因子水平的列。一个 one-way anova 分析的例子 1.打开文件 pulse.mtw 2.选择 stat
one-way 3.在反应值response中，输入 weight， 在 factor 中，输入 sex。 4.点击 graphs。 5.选中 boxplots of data 和 normal plot of residuals。 点击ok两次后，会话窗内容如下：35分散分析(anova)6sms
f p 90.80 0.000one-way analysis of variance analysis of variance for weight source sex error total df 1 90 91 ss
51284 individual 95% cis for mean based on pooled stdev level 1 2 n mean stdev --+---------+---------+---------+---(--*-) (---*--) --+---------+---------+---------+---pooled stdev = 16.84 120 135 150 16557 158.26 18.64 35 123.80 13.37boxplots of data 和normal plot of residuals 如下：normal probability plot of the residuals(response is weight) 3boxplots of weight by sex(means are indicated by solid circles)2002normal score1 0 -1 -2weight150100-3sex12-40-30-20-100102030405060residual结果解释 分散分析输出结果和盒式图(boxplot)清楚地表明，女性和男性的体重是不一样的。f统计值比较大而p值比较小，表明男女存在统计上显著性差异。女性 体重平均值的95%置信区间介限于 118和130之间但是对于男性而言,95%置信区 间介于155和162之间.36表(table)6s使用 minitab制表功能，可以将数据以表的形式归纳出来,或者对这些数据进 行进一步的分析。也可以进行各种各样的统计，或者进行 χ2 检验，也可以进行 简单或者复杂的相应的分析。交叉表格可以显示一次，二次，多次表格数据。在这些表格里，包含常量、百分数、统计数据，例如平均，标准偏差，最大值以及相关变量。 怎样创建一个统计表格 a. 选择 stat
cross tabulation。 b. 以下二者任选其一： ※ 对于原始数据，在 classification variables 中，输入两到十个包含原始数 据的列。 ※ 对于频率数据，a. 在 classification variables 中，输入两到十个包含分类数据的列。 b.选中 frequencies are in 并且输入包含频率数据的列。37表(table)6s一个显示各列百分比的 二次表的例子假设你要对不同程度吸烟者的数量及各所占比率进行统计。 1.打开文件 pulse.mtw。 2.选择 stat
cross tabulation。 3.在 classification variables 中，输入 smokers 和 activity。 4.选中 column percents 然后点击 ok 。tabulated statistics: smokes, activity rows: smokes 0 1 2 all 100.00 -100.00 columns: activity 1 33.33 66.67 100.00 2 31.15 68.85 100.00 3 23.81 76.19 100.00 all 30.43 69.57 100.00cell contents -% of col 结果解释 行表示smokers 变量：1表示经常吸烟者而2表示不常抽者。列表示水平程度：1=轻 微，2=中等，3=许多（这里的0水平是因为测量者误输入），轻微抽烟者占1/3，抽烟 量最大者占1/4。38质量控制6sminitab提供了多种质量管理方法：管理图、质量计划工具、测量系统分析、 工序能力和可信度。 管理图或者统计工序管理图(spc)可以对工序整个时间的变动进行分析。根据样本数据，这些图表生成了一些总的统计（例如：样本平均值或样本概率）。下面是经常使用的x-r图。从字面意思可知，此图表在一个图形窗口中通常有两 个图：一个是子集平均值的管理图（x )、一个是子集范围的管理图（r），同 时查看这两个图，可以跟踪工序水平和工序变动，还可以看出特殊原因的表现。 子集数据可以由单列或者穿过几列的数行构成。当子集数据尺寸不唯一时， 子集数据最好由单列构成，然后建立另一列作为此子集的标志列。 怎样创建一个x-r图 a. 选择 stat
control charts
xbar-r b.以下二者任选其一 1. 如 果 子 集 处 于 一 列 ， 在 single column 中 输 入 数 据 列 。 在 subgroup size 中，输入子集尺寸或者子集标志列。 2. 如果子集处于行中，在 subgroups across rows of 中输入一系列列。 如果有必要，可以点击 option，可以作如下选择： a. 如果对于μ 和σ 有一个目标值或者根据以前的数据知道这几个参 数，则输入μ 的历史值（总体分布的平均值）和σ （总体分布的 标准偏差） b. 控制 minitab 估计μ 和σ 值的方式。 c. 对于特殊原因做8个检验。 d. 当出现了偏移数据时，使用 box-cox 转换。39质量控制6sx-r图的一个例子 一制造车间制造金属固定件，需要估计其工序能力（目标值：2.4克），力求固件的 重量变动最小。这些固件5个一包。随机选择20包画出x-r图，估计他们的制造工序管理 情况。 1.打开文件 2.选择 fastener.mtw。 stat
controlｃharts
xbar-r。3.在 single column 中，输入weights。在 subgroup size 中，输入5。4.在 historical mean 中输入 2.4，点击 ok 。test results for xbar chart test 1. one point more than 3.00 sigmas from center line. test failed at points: 1040质量控制6sxbar/r chart for weights2.8 2.7 2.6 2.5 2.4 2.3 2.2 2.1 2.0 0 ucl=2.719sample meanmean=2.4lcl=2.081 1 10 20subgroupucl=1.170sample range1.00.5r=0.55350.0lcl=0 结果解释 特殊原因检验表明只有一包的平均重量处于规格下限以下（低于目标值 3个 标准偏差-3σ ）。但是又注意到许多点落在中心线以下，表明此工序生产的固件 经常低于目标值。质量工程师还得进行另外的特殊原因检验或者重估此工序。41第五章 基本操作示例6s5-1 5-2 5-3 5-3 5-4 5-5 5-6 5-7事例背景 打开一个工作表 从键盘输入数据 输入规则数据 描述性统计计算 算术运算 生成散点图 计算相关系数42事例背景6s实验背景 克隆是单个基因细胞从同一个体繁衍而来。研究人员发现一种单性克隆物 （生长快而且质硬的树），这些树某天也许会当作一种能源而取代传统燃料。宾夕法尼亚州立大学的研究人员在两个不同的地方，种植了一些白杨树，一个地方在小溪旁，土质肥沃，灌溉系统好；另一个地方是山地，土壤干燥，呈现 沙性。他们以厘米为单位测量直径，以米为单位测量高度；然后把它烘干，测 量它的重量。研究人员想知道，他们能否根据直径和高度，预测它的重量。在这一章里，你可以知道怎样：一、打开一个工作表 二、输入和编辑数据 三、保存数据 四、基本统计计算 五、数学运算 六、画数据图七、计算相关系数43打开一个工作表6s步骤一、打开一个工作表 当你打开 minitab 时，便开始了一个新而且空的 project。可以通过多种方式向此 project 里添加数据，但是最为普通的方式是打开一个工作表。只不过有一点要注意，打开工作表只是将数据拷贝到此 project 里，对于project 里任何数 据的变化都不会影响到原始数据。 在这一节，使用的数据在文件poplar1.mtw里。 1. 选择文件 file
open worksheet。 2. 找到子目录 data。 3. 点击 poplar1.mtw，打开。 4. 如果未见到数据窗口，按快捷键 ctrl+d。这个工作表包含三个变量，分别为：diameter，height，weight。每个变量 包含15个值——到目前为止只有这么多。44从键盘输入数据6s步骤二、从键盘输入数据 到目前为止只有这么多数据，但是又从现场收集到了5行数据，需要通过键盘来输入。1. 将光标放到第16行。数据窗口应当为下面式样：数据方向箭头2. 确保数据方向箭头向右. 3. 从左到右输入以下数据：1.524.51 1.18 3.17 3.33enter enter enter enter enter2.95.27 2.2 4.93 4.89enter enter enter enter enter0.070.79 0.03 0.44 0.52ctrl ctrl ctrl ctrl ctrl+ + + + +enter enter enter enter enter45输入规则数据6s步骤三、输入规则数据 一般是在数据窗口中通过键盘输入数据，但是当数据按照一定的规则时， 还有一个简便的办法输入。现在要创建一个新的变量来标明所测数据值是来自于土壤肥沃、排水系统好的地方（记为1）还是来自于干燥多沙的地方（记为2）。 这个新的变量叫做site，将包含10个“1”，然后是10个“2”。 1. 选择 calc
make patterned data
simple set of 自动分配一个新的变量给第一个空列。在这里是 c4。 3. 为了标明开始和结束顺序，在 from first value 中输入 1，在 to last value 中输入 ｎumbers。 2. 保存数据于：在 store patterned data in 中输入 site 。minitab2 。list each value 中，输入 10 ，然后点击 ok 。4. 既然有 10个1和 10个2，在 新生成的数据如下：46描述性统计计算6s步骤四、描述性统计计算 minitab提供了许多基本统计方法进行数据分析，例如描述性统计、t检验、 z检验以及相关性。 对于两个不同地方的树，应当分别生成统计表对三个变量diameter、 height 和 weight 进行说明。 1. 选择stat
basic statistics
display descriptive statistics。 2. 在 变 量 列 表 框 中 ， 点 击 diameter 然 后 拖 曳 鼠 标 使 diameter 、 height、weight 高亮度显示，然后点击 select 。 3. 选中by variable，输入 site 。通过 minitab 对不同地方的树在直径、 高度、重量上分别生成统计结果。 4. 点击 graphs。5. 选中 boxplot of data，然后点击ok两次。47描述性统计计算descriptive statistics: diameter, height, weight by site variable diameter height weight site 1 2 1 2 1 2 site 1 2 1 2 1 2 n 10 10 10 10 10 10 mean median 2.598 2.320 3.028 3.250 4.098 4.120 4.255 4.865 0.0 0.399 0.380 minimum 1.060 1.180 1.850 2.200 0. trmean 2.604 3.041 4.175 4.351 0. stdev 0.916 1.284 1.103 1.254 0. q1 2.120 1.488 3.518 2.775 0. q3 3.245 4.053 4.853 5.143 0.6svariable diameter height weightse mean 0.290 0.406 0.349 0.396 0.maximum 4.090 4.770 5.730 5.540 0.boxplots of diameter by site 5 1.0 d ia m e te r 4 3 2 1 1 site 2 0.0 w eightboxplots of weight by site0.51siteboxplots of height by site6 5h eig h t4 3 2site1从盒式图(boxplot)可以判断出，地方2的树比地方1的树大。会话窗口的文本 输出包含以下信息： 在地方2三个变量都显示出比较大的平均值和中央值。而且，变量weight相 对其尺寸来说有一个比较大的标准偏差。在地方2，最小重量只有0.03kg而最大 值是1.11kg，这表明一部分白杨树生长得很好，可是其他的仅仅只是活着而已。4822算术运算6s步骤五、进行算术运算 现在的任务是估计树的重量。根据以前的工作，研究人员发现树的重量和直径平方与高度之积有关。手头有直径和高度的数据，利用计数器计算新的变量。计数器将输入的方程式结果计算出来并保存于指定的变量。 1. 选择 calc
calculator。 2. 将新的变量命名为d2h，表示直径的平方与高之积。在 store resultd2h . 3. 在 expression中，键入 c1**2*c2 ，点击 ok 。in variable 中，输入计算结果显示如下：49生成散点图6sgraph
plot。ok 。步骤六、生成散点图 研究人员得出在重量与d2h之间存在一种关系。通过散点图，可以验证这一点。1. 选择2. 在 y 中，输入 weight . 3. 在 x 中，输入 d2h 然后点击图形窗口输出结果：1.0w eight0.50.0 0 50 d2h 100不规则点由此可以看出，weight和d2h之间存在正的线性关系，也就是说，d2h增加 的话，weight也同样增加。图上出现了一个不正常点：有一棵树只有很小的重 量，但d2h值却相当高。现在将其忽略，但以后还得进一步检查。下一步是计算 这二者的相关性并且将其量化。50计算相关系数6s步骤七、计算相关系数 从散点图知随着d2h增长，weight也同样增长，现在通过计算相关系数找出这两个变量的具体关联性。1. 选择 stat
basic statistics
correlation. 2. 在 variables 中，输入 weight 和 d2h，点击 ok 。会话窗口输出：pearson correlation of weight and d2h = 0.913 p-value = 0.000相关系数是两个变量之间线性关系的度量，其值介于-1和+1之间，此处，相 关系数为正的0.913，接近于1，从而把散点图中所表示的关系量化了。说明变 量 weight和d2h 之间存在很强的线性关系。51第六章 做一个简单分析6s526s背景 宾夕法尼亚州立大学的研究者们种植了成百上千的白杨树并且在一系列受 控条件下让其生长.三年后,抽取一批为样本，以厘米为单位测量其直径、以米 为单位测量其高度、千克为单位称其干重。 前面知道，在干重与d2h之间存在相近的关系。但是到底是一种怎样的关系 呢这些快速生长、质地坚硬的树，将来某一天也许会作为燃料或者化学物质 的一种选择。作为数据分析员，知道了直径和高度的测量值之后，能否可靠的 预测出树的产量。步骤一. 打开文件 poplar2.mtw。 步骤二. 做一个简单的回归分析 在上一个示例结尾，我们知道随着d2h增加，weight也增加。由d2h的值预 测weight的一种方法是使用简单的回归命令。 1. 选择 stat
regression
regression 2. 在 response 中，输入 weight。 3. 在 predictors 中，输入 d2h。不妨还做一系列残差分析图，以便找出 潜在的问题。 4. 点击 graphs。536s5. 在 residuals for plots下，点击 standardized。 6. 在 residuals plots ，点击 histogram of residuals 和 normal plot of residuals。 7. 在 residuals versus the variables下，输入 d2h。 8. 点击 ok 两次。minitab便显示会话窗口中的文本、三个图形。546sregression analysis the regression equation is weight = 0.0196 + 0.00758 d2h predictor coef stdev constant 0.66 d2h 0..0007994 s = 0.1298 r-sq = 83.3% t 0.43 9.49 p 0.673 0.000r-sq(adj) = 82.4%analysis of variance source regression residual error total df 1 18 19 ss 1.1 1.8187 ms 1.8 f 89.99 p 0.000unusual observations obs d2h weight 12 126 1. 0.0700fit 0.9stdev fit 0.4residual st resid 0. x -0.rr denotes an observation with a large standardized residual x denotes an observation whose x value gives it large influence.histogram of the residuals(response is weight)15frequency1050 -4 -3 -2 -1 0 1standardized residual556snormal probability plot of the residuals(response is weight) 2 1normal score0-1-2 -4 -3 -2 -1 0 1standardized residualresiduals versus d2h(response is weight)1standardized residual0-1-2-3-4 0 50不规则点100d2hminitab 显示了两个不规则的数据（12行和15行），进一步分析之前，首 先要检查12行和15行是否包含了有效的数据。566s步骤三 、编辑数据 1. 按键 ctrl+d 2. 找到第一个不规则点。 a. 选择 editorgo to。 b. 在 enter column number or name 中，输入 weight c. 在 enter row number 中，输入 12，点击 ok。在数据窗口中，光 标置于12行 weight 单元格。可以看出，weight 和d2h看上去显得相当大，然后再检查原始记录，发现 12号树生长很好，说明此数据正确。按照同样的方法，找到第15行。检查原始 记录，发现这个数据错误，不是0.07，而是0.7。576s简单的重复一下前面的菜单选择。 1. 首先关闭所有图形窗口。 2. 选择 stat regression regression 。regression analysis the regression equation is weight = 0.0200 + 0.00829 d2h predictor coef stdev t p constant 0.65 1.46 0.160 d2h 0...68 0.000 s = 0.03880 r-sq = 98.5% r-sq(adj) = 98.4% analysis of variance source df regression 1 residual error 18 total 19ss 1.1 1.8379ms 1.5f 1202.89p 0.000unusual observations obs d2h weight fit stdev fit residual st resid 12 126 1.92 0.08 1.39 x 17 107 0.58 0.058 -3.42r r denotes an observation with a large standardized residual x denotes an observation whose x value gives it large influence.如果有一个比较好的模型并且对所有的统计假设都感到满意的话,那么可以 直接抽取一部分作白杨树为样本, 测量直径和高度,然后预测其重量,而不用把 树砍下来再去称量。 由此回归输出结果可以看出,在系数表中,d2h有一个比较高的t值和很小的p 值,有足够的理由说明d2h和weight之间存在很强的关系。 在变量分析表中比较大的f统计值和值比较小的p值将这种关系量化.r的平 方(r-sq)以及r的调和平方(r-sq(adj))的值都大于98%,进一步证明d2h与weight 之间存在很强的线性关系。 在下最后结论之前，再看一下图表。586shistogram of the residuals(response is weight) 10frequency50 -3 -2 -1 0 1 2standardized residualnormal probability plot of the residuals(response is weight) 21normal score0-1-2 -4 -3 -2 -1 0 1 2standardized residualresiduals versus d2h(response is weight) 2 1 0 -1 -2 -3 -4 0 50 100standardized residuald2h从图表residuals versus d2h知道变动并未呈现出常数—此即符合一 个很重要的回归假设模型。在比较大的d2h值处，标准残差也比较大.基于时间 因素,以后章节还要进一步讨论。596s步骤五、画一条合适的回归线 下一步，画一个带回归线的散点图来看看那些测量的数据与回归线的接近程 度。 1. 选择 stat regression fitted line plot。 2. 在 response（y）中，输入 weight。 3. 在 predictor（x）中，输入 d2h，点击 ok。 输出结果如下：606s步骤六、通过brush（刷）图形找特殊点 图形上有一点，其d2h值很大，想知道这到底是一个什么样的点。brush可 以快速方便的找到这一点。 1. 选择 editor brush。 2. 在图形上将手状的图标圈住特殊点。特殊点便会显示出来，同时在数据 窗口里相应的行上做一个标记。刷板 要找的点61第七章 高级minitab6s62
6s背 景 从这些种植的树木中真正能预计到多少可以当作能源的材料呢又怎样将 产量最大化呢 为了将产量最大化，研究人员设计了一个实验来决定地点的不 同（site）和条件的不同（treatment）这两个因素是怎样影响四年生白杨树的 重量的。 他们将树种植在两个地方：site1和site2。site1土壤肥沃、湿润，而 site2干燥多沙。他们将这些树分成四种条件培育。条件1（treatment1）只是 控制（无条件）、条件2（treatment2）：施肥；条件3（treatment3）：灌溉； 条件4（treatment4）：既施肥又灌溉。考虑到各种气候的影响，研究人员在第 一年种植一半的树，第二年种植另外一半，从而获得两分数据。 为了进行分析，需要对样本数据进行统计分析。 步骤一、生成一个新的project。 filenew，点击 minitab project，然后点击 ok。 步骤二、打开一个工作表 1. file open worksheet。 2. 找到子目录data，选择文件 poplar3.mtw。点击 ok。 3. 有两个窗口显示了关于此工作表的信息。此数据窗口包含7个变量：site、year、treatment(实验条件)、diameter（cm）、height、weight、age。646s信息窗口步骤三、生成描述性统计 为了产量最大化，需要聚焦哪个因素对重量的影响最大。先看一下变量 weight的描述性统计。 选取stat basic
statistics display descriptive statistics. 在窗口vaibles 中选取weight.descriptive statistics variable n weight 298 variable weight minimum -99.000mean 1.099 maximum 6.930median 1.640 q1 0.597trmean 1.994 q3 3.455stdev se mean 10.255 0.594看看weight的最小值。很明显，重量不可能出现-99千克。事实上，在此，收集 数据的时候，用-99表示已经死亡的树。 如果把-99留在数据表中，那么所做的分析也是错的。事实上，刚才所做的描述性 统计也受到了一定影响。平均值和中央值人为的偏低而标准偏差人为的偏高。需要将 -99转换成空值。 空值并不影响统计分析的结果。在minitab 中，空值用星号“*”代替。656s步骤四、对数据进行替换 minitab提供了许多数据操作工具。最常用的一个就是code命令，它的作用 是将一个数值用另一种形式的值表示。在这里，要把-99变成*——空值的符号。 1. 选择 manipcode numberic to numeric 2. 在 code data from columns 中，输入 weight。 3. 在 into columns 中，输入 weight。这样将会用新值代替原有的值。 4. 在 original values 中，输入 -99。此即要改变的值。 5. 在 new 一栏，输入星号 * 。然后点击 ok 。在数据窗口变量 weight 栏中，可以看到所有出现 –99 的地方都被星号 * 代替。步骤五、tally数据 在此实验中，每个年龄的树数量各是多少呢使用tally命令可以帮你知道。 1. 选择 stat tables tally。2. 在 variables 中，输入 age。点击 ok。输出结果如下：summary statistics for discrete variables age count 3 147 4 151 n= 298输出结果表明三年树龄的有147棵，四年树龄的有151棵。666s步骤六、根据树龄分割数据表 假设只要对四年生的树进行分析。下面是只对四年生的数据创建新的数据 集的一种技巧。对四年生的树建立一个单独的数据表 1 选择 manip split worksheet。 2. 在 variables 中，输入 age 。点击 ok。 minitab将回根据 age 值分割工作表poplar3.mtw。因为在age列，只有两个可 区分值：3和4。minitab便生成2个新的工作表。一个工作表只包含三年生的，命 名为poplar3.mtw（age=3），另外一个是四年生的，命名为poplar3.mtw （age=4）。将其重新命名： poplar4.mtw（age=4）命名为 4yrolds.mtw。步骤七、用直方图检验正规性 现在对变量weight生成一个直方图 1. 首先将数据窗口4yrolds.mtw激活，使之成为当前窗口。 2. 选择 graph histogram。 3. 在 x 中，输入weight。点击 ok 。676s白杨树的重量近似与正态分布（即有点像钟型）。步骤八、用盒式图(boxplot)根据条件（treatment）变量的不同对weight（重 量）变量进行比较。 现在我们来看一下不同条件下的白杨树重量。盒式图是比较一个变量各个 水平差异的一种很好的图形法。 1.选择 graph boxplot。2.在 y 中，输入 weight。 3.在 x 中，输入 treatment。这样的话，minitab会对各个水平条件 treatment 下的weight分别生成盒式图(boxplot)。4.点击 ok 。686s四分之三四分之一 中央值横贯每个盒子(box)的线条表示数据的中央值或者中间值。上面的顶线和下 面的底线分别表示75%和25%。此盒式图表明在treatment2和treatment4（施肥和 施肥/灌溉）下，树木生长得很好，重量比较大，但是在treatment1和treatment3 （管理和灌溉条件）下，树木长得比较轻。同样，我们再来看看地点的不同对 重量有什么影响。地点1的树是种植在肥沃、排水系统好的土壤里，应当说来比 起种植在地点2（干燥多沙的土壤）中的树要长得好。下面通过盒式图来比较一 下地点的不同对重量有何影响。 5. 选择 edit edit last dialog。 6. 在 x 中，输入 site。 7. 点击 ok。奇怪的是，地点1和地点2树的重量看上去没有明显的不同，虽然每个地方 的散布不同，但是中央值几乎相同。696s步骤九、分散分析 从以上的盒式图已经看出，白杨树的重量在四种条件下明显不同，但是对 于两个地点来说，重量的差别不是很明显。 现在决定使用变动分析看看对于因子site和treatment的不同水准，重量是否 存在统计上的显著性差异。 当有两个或者两个以上的因子时，minitab会给你两种选择：对称anova和 一般线性模式（glm）。对称anova必须是对称设计， 也就是说，对于每个 site 和treatment的组合必须有同样数量的观测值（即： site有两个值，treatment也 必须有两个）。在这里，因为我们的设计不是对称的，所以使用一般线性模式。 1. 选择 stat anova general linear model。2.在 responses 中，输入 weight。下一步，输入适合的glm模式。 3.在 mode中输入 site|treatmnt。这根竖线的意思是告诉 minitab，在此 模式中要包括所有可能的交互作用。此竖线的输入，可以通过按键 shift + ， 或者直接使用感叹号“！”代替。 4.点击 ok 。706sgeneral linear model factor type levels values site fixed 2 1 2 treatmen fixed 4 1 2 3 4 analysis of variance for weight, using adjusted ss for tests source df seq ss adj ss adj ms f p site 1 3.112 2.424 2.424 1.52 0.219 treatmen 3 78.005 78.275 26.092 16.39 0.000 site*treatmen 3 10.509 10.509 3.503 2.20 0.091 error 140 222.873 222.873 1.592 total 147 314.498 unusual observations for weight obs weight fit stdev fit residual st resid 22 0.00 0.200 -2.08r 42 0.67 0.267 -2.20r 43 0.67 0.267 -2.59r 52 0.50 0.250 -2.33r 64 2.89 0.289 -2.08r 69 2.89 0.289 -2.27r 72 5.67 0.33 2.02r r denotes an observation with a large standardized residual.在这个模式中，glm列出了每个因子以及每因子的水平数，然后列出了变 动分析表，最后是不规则数据表。 假设要对此模式中的每个影响进行 f-test。例如：要验证原假设（条件的影 响对于两个地点是相同的—site*treatment的交互作用），用minitab的p值与最常 用的0.05水平值进行比较。因为p-value为0.091，它大于0.05，所以你不能拒绝原 假设。也就是说，不能得出“条件对地点的影响不同”的结论。 现在看看site和treatment的主效果。site的p值是0.219，同样大于0.05，所 以也不能得出“树的重量在两个地点之间存在显著性差异”的结论。而 treatment的p值很小，只有0.000，因此可以说对于不同的条件，树的平均重量存 在显著性差异。 这与以前画的盒式图是相吻合——树的重量在不同条件下不同，而对于两 个地点来说只存在轻微的差别。在作出“条件是影响重量的唯一重要因子”结 论之前，再看看年份对于重量的影响——别忘了，研究人员分两年种植了这批 树。716s步骤十、用盒式图(boxplot)根据年份对重量进行比较 1.选择 graph boxplot。2.在 y 中，输入 weight。 3.在 x 中，输入 year。意思是根据年份分别画出重量的盒式图(boxplot) 。 请看一下data display表中的头两行。iqrange box告诉minitab显示盒式图 (boxplot)的时候，显示25%~75%这一段的百分比范围。outlier symbol的意思 是将所有超出范围的值用星号*表示。现在还想在 iq range box里显示置信度 的盒式图。 4.在 data display 表中的display列，点击第三项。 5.点击 display 右边的下拉箭头，选择 ci box。 6.点击右边的单元格。 7.点击 for each 右边的下拉箭头，选中 graph。意思是叫 minitab在每个 盒式图中包含置信度。缺省的情况下，minitab用垂直方式画盒式图，但是可以 画成水平状。 8.点击 options。9.选中 transpose x and y，点击ok两次。726s2year里面的盒子显示 了95%的置信度101234567weight里面的盒子显示了中央值95%的置信度。盒式图表明种植在第二年的树比 第一年的要重。分两年种植只是为了重复一下数据，为什么年份很重要呢经 过采访研究人员知道，在第一年里他们没有施加除草剂来控制杂草，结果有些 树不是死了就是缺少养分。为了提高其成活率，研究人员决定对第二年的树施 加除草剂。通过分析，得出三个主要结论： 1. 施肥是提高重量的一个行之有效的办法。 2. 在树的成长初期，控制杂草很重要。 3. 虽然白杨树对地点的要求不是很高，为了提高实际产量，给予适当的生 长和营养条件。第二年的树不但重一些，而且分布也比较均匀些。在建议使用 除草剂和施肥之前，还得仔细看看第二年的树。736s步骤十一、快速重复一遍所有的分析。 现对第二年的树重复一遍分析。首先要把四年生的数据分开出来，然后使 用历史窗口和命令行编辑器（command line editor），而不是重复使用一个一 个的对话框。 a.根据年份值将poplar3工作表分割开。 1.确定 poplar3.mt w当前活动工作表。 2.选择 manip split worksheet。 3.在 by variables 中，输入 year.点击 ok。 b.重新命名 4.点击新的工作表poplar3.mtw（year=2），重新命名为year2.mtw。 点击ok。 c.对树龄为2的重复一遍分析 5.打开历史窗口：选择window history。 6.拖窗口的滚动条到命令“histogram”出现。 7.从histogram到最后的内容全部选中。 8.选择edit command line editor。出现一个对话框，包含了刚才 所选中的内容。 9.点击 submit commands。整个分析：重量的直方图、不同条件下重 量的盒式图、不同地点重量的盒式图、变动分析、根据年份重量的盒式图到此 全部做完。74
第八章 质量管理和改善6s766s背景 假设你在汽车制造工厂的一个引擎组装部门工作。其中有一个组件是凸轮 轴，长度规格是600±2mm，长期以来一直存在一个问题，凸轮轴的长度总是超 出规格，由此导致组装线适合度下降、磨损和返工率上升。 作为管理者，需要通过x和r图来监控这些特性。一个月里，对凸轮轴的长 度数据以每5个为一组进行了收集。 步骤一、生成一个新的project。 步骤二、打开一个工作表 1. 选择 fileopen worksheet 2. 在data子目录下选择 camshaft.mtw。点击 ok。这就是一个月里所收集到的数据。注意到第一列，它包含了100个观测值 （20个样本，每样本5个数据，单位是毫米。 步骤三、用r图检查一下范围。 首先生成一个管理图看看样本数据里凸轮轴的长度范围情况。我们希望样 本数据的范围离中心值（估计的平均值）不要太远以免出现大的变动。 1. 选择 stat control charts r。2. 在 single column 中输入 length 。 3. 在 subgroup size 中，输入 5。点击 ok 。776s上限 平均值下限每一点表示每个 班次的范围（最 高值减去最低值）此r图所显示出并没有任何点超出规格。注意到中心线在2.72mm处，比所给 的最大值±2要大得多。说明工序存在很大的变动。 步骤四、用xbar图对特殊原因测试。 生成一个x-bar图看看凸轮轴长度是否存在超出规格的问题。另外，使用 minitab的“eight common runs”命令来找出变动的特殊原因。 1.选择 stat control charts xbar。2 在 single column 中输入 length。 3.在 subgroup size 中，输入 5 。 4.点击 tests 。786s5.选中 perform all eight tests。 6.点击 ok 两次。控制上限平均值每一点表示每个 班次数据的平均值 控制下限此 x-bar 图表明工序超出规格。尤其是有一点在test1时出错，有两点在 test6 时出错。在会话窗口中，说明了这些检验表示的意思。796s7. 选取 window session。test 1. one point more than 3.00 sigmas from center line. test failed at points: 8 test 6. 4 out of 5 points more than 1 sigma from center line (on one side of cl). test failed at points: 12 13此工序有一点离中心线超出 3 sigma，4/5的点离中心线 1 sigma 之外（在 中心线的一边）。 既然问题已经确实存在，下面应当找出原因和解决的办法。不幸的是，抽 样计划并没有对问题精确的在哪里、什么时候发生作出具体的说明，因为每个 班次只抽取了一个样本。好的计划应当是在问题的解决阶段每个班次多抽取几 个样本。在特殊原因找到并解决之后再转换到监控阶段。因此，从这些数据可 以知道能得到什么了。 步骤五、生成正态曲线的直方图。 直方图的正态曲线是一种检查变量分布的有用方法。 1.选择 stat basic statistics
display descriptive statistics。2.在 variables 中，输入 length。 3.点击 graph，出现一个对话框。806s4.选中 histogram of data，with normal curve。 5.点击 ok 两次。一般说来，我们希望像length这样的变量服从正态分布。在这里，直方图 近似于钟形。刚才所生成的直方图当然不是很好的钟形。实际上，从598，599， 601高峰开始呈现出两个独立的、明显不同的分布。 检查记录清单发现凸轮轴供应商有两家，现在开始明白为什么会出现这样 奇怪的直方图了。所以决定从两家供应商获得测定值并对每批数据单独运行x和 r图。每家供应商的数据都保存在supp1 列和supp2列中。816s步骤六、同时显示xbar和r图 1.选择 stat control charts xbar-r2 选择 single column，输入 supp1 。 3 在 subgroup size 中，输入 5 。点击ok 。尽管平均值是599.5，而不是600，supplier1 的平均值和 range 看上去都在管理范 围之内。range的平均值是 1.36mm 。826s另外：对supplier2进行估算 对于supplier2和supplier1一样，使用同样的管理图。操作步骤同上，只不过 在single column中，输入supp2，其余的不变。supplier2 的 x-bar 图反映出有问题。从这个图上，可以看出有两点超出 管理上限。r 图没有显示出工序超出管理，然而，中心线是3.720，几乎是supplier1的3倍。从长远来看，应当等到supplier2生产的凸轮轴在管理范围之 内时才采用supplier1。对于supplier2，应当与其合作，共同减小变动到可接受的水平。836s步骤七、准备工序能力分析 只采用supplier1，使得变动大为减小。从组装线上下来的组件不良数明显 降低，但是问题没有完全消除。决定进行能力分析，看看supplier1是否符合工程规格600±2mm。进行工序能力分析之前，工序必须可控。管理图表明，由于措施得力，工 序现在在管理范围之内。同时还希望凸轮轴的长度服从正态分布。现在由直方 图检验其正态性。 现在看看supplier1的凸轮轴长度的分布。 1.选择 stat basic statistics display descriptive statistics。 2.在 variables中，输入 supp1，点击 ok。此直方图比较令人满意，不象以前那样，再没有多个尖蜂。846s步骤八、 进行工序能力分析。 1.选择 stat qulity tools
capability analysis（normal）2.选中 single column 输入 supp1 。 3.在 subgroup size 中，输入 5 。 4.在 lower spec 中，输入 598 。 5.在 upper spec，输入 602 。 6.点击 options 。7.在 target（adds cpm to table），输入 600 。 8.点击 ok 。856s从上面这个图形可以看出，工序的平均值与目标值不足，中心值落在目标 值的左边。曲线的左边尾部落在了下限以外。因此，有些凸轮轴不符合最低规 格598mm。 我们再看看文字性的统计说明。 cpk 指数表明工序能否在公差范围内进行生产。如果cpk指数为 1，意味着 工序恰好符合规格；小于 1 说明超出规格。理想情况下，希望 cpk 远大于 1， 因为这个指数越大，表明工序越好。supp1 的 cpk指数只有0.9，表明需要减小 变动，靠近目标值。86第九章 实验设计6s876s背 景 假设你在一家化学工厂工作，任务是对生产一种化学产品的化学反应进行 研究。你想通过这个反应提高产量。根据过去的经验已经知道改变温度、气压 和催化剂的种类可以改变反应产量。现在的问题是你的那些同事对于以上那些 因素怎样影响反应的各有一套自己的见解。为了真正的改善，你决定做一个实 验来决定这三个因子中哪一个是真正的影响者。 步骤一、生成一个新的project。 步骤二、生成一个实验设计 选择一个设计 假设你要设计一个实验测试三个因子：时间、温度和催化剂类型。 1. 选择 stat
create factorial design只有两个按钮有效：display available designs和designs。 其他按钮是设置因子水平、特殊 选项、控制会话窗口输出、点击 design按钮之后别的按钮便会生 效。886s2.点击 display available designs。 使用此表可以帮助你选择一个合适的设计，主要基于以下几点： ● 感兴趣的因子数目。● 实验中能进行的实验次数。● 预期的实验解决方案。这里我们对三个因子比较感兴趣，从上面这个表中知道有两种选择： ● resolution iii运行四次的 部分因子设计 ● 运行八次的完全因子设计 三因子两水平设计有 23 个可能的因子组合数。通过选择所有可能的组合设 计（完全因子设计），可以没有混淆地得到结果，也就是说，所有效果或影响 都彼此分离开。也可以通过少做几次实验或组合来得到有意义的结果。 最后决定采用三因子的完全因子设计（完全配置法）比部分因子设计（部 分配置法）要合适。在你的化学工厂里，感兴趣的三个因子—时间、温度和催 化剂种类—不需要很多的开销（经济的和时间耗费方面的）。而且此实验的进 行也不是在生产高峰时候，不会影响生产。如果实验的进行开销很大的话，则 需另择良方。896s3.点击 ok。回到主窗口。 4.选择 2-level factorial（default generators ） 5.在 number of factors 中，选择 3 。 6.点击 designs。 7.在 designs 框中，选择 full factorial。 8 在 number of replicates 中，选择 2 。 9.点击 ok 。 给因子命名并设置因子水平 在两水平因子设计中，将因子设置为两水平—一个低水平、一个高水平。 每次实验都包含因子高低水平的组合。像催化剂这种类别型的数据，实际上没 有高低水平，可以输入文本型的水平值。现将各因子水平设定如下： 因子 温度 气压 低水平 20℃ 1大气压 高水平 40℃ 4大气压催化剂ab1.点击 factors。出现如下对话框。2.将上表对应添入对话框。 3.点击 ok 。906s随机化并保存设计 1.点击 options。出现如下对话框。2.在 base for random data generator 中，输入 9 。 3.确保 store design in worksheet 已经选中。点击 ok 。 4.所有的选择都会生效。这样便生成了一个设计。步骤三、查看生成的数据runorder是实 验时的实验顺 序916s步骤四、收集和输入数据 为了收集实验时的数据，需要创建数据集。 1.在数据窗口中，点击c8输入“yield”。要进行16次实验，需要记录每次 的产量。假设得到了以下16个数据（单位：克）。 2.在“yield”列键入以下数据： 66 ，66，102，98，65，54，107，68，53，66，55，85，108，89，52，63 步骤五、筛选设计 筛选设计的目的是找出影响大的因子。既然已经生成了部分因子设计（部 分配置法）并且收集了一些数据，那么可以对这些结果模式化、生成图形来估 计效果。根据适合的算术模式和两种图形方法可以看出哪个因子对提高反应的 产量影响最大。1.选择 stat doe
analyze factorial design.2.在 responses中，输入 yield 。 3.点击 graphs 。926s4.为了生成两个效果图以便知道哪个效果显著些，选中 normal 和pareto， 使用缺省值α 。 5.两次点击 ok。fractional factorial fit estimated effects and coefficients for yield (coded units) term effect coef stdev coef t p constant 74.81 2.561 29.21 0.000 temp 1.38 0.69 2.561 0.27 0.795 pressure 14.12 7.06 2.561 2.76 0.025 catalyst -30.37 -15.19 2.561 -5.93 0.000 temp*pressure -0.13 -0.06 2.561 -0.02 0.981 temp*catalyst -1.13 -0.56 2.561 -0.22 0.832 pressure*catalyst -13.37 -6.69 2.561 -2.61 0.031 temp*pressure*catalyst -0.13 -0.06 2.561 -0.02 0.981 analysis of variance for yield (coded units) source df seq ss adj ss adj ms f p main effects 3 6.19 .28 0.001 2-way interactions 3 720.69 720.69 240.23 2.29 0.155 3-way interactions 1 0.06 0.06 0.06 0.00 0.981 residual error 8 839.50 839.50 104.94 pure error 8 839.50 839.50 104.94 total 15 6056.44完成了一个完全模式，它包含三个主效果、三个二次交互作用、一个三次交 互作用。在 estimated effects and coefficients 表中，根据p列的值决定哪 个效果显著。因为使用的α 值是0.05，气压、催化剂、气压和催化剂的交互作 用显著，也就是说，他们的p值小于0.05。936s效果图（effect plots) 现在可以使用一般概率图和效果的pareto来看哪个效果影响产量。在一般 概率图中，那些与直线结合得不是很好的点的效果往往很显著。0是近似的平均值， 而1和-1是每边的 标准偏差 。一般概率图将那些低于α 水平的效果标志出来了。这里，气压、催化剂、 气压和催化剂的交互作用显著，此时α =0.10。pareto是另一种有用的工具，它 可以帮你找出哪个效果重要。任何超过这条线 的效果都是显著 的。在 pareto上，显示了效果的绝对值。pareto和一般概率图一样，使用同样 的α 值来区分效果的显著性。从这里再一次看到，压力、催化剂、压力和催化 剂的交互作用很显著。等下要找一个模式，没有温度条件、温度和压力、温度 和催化剂，它们不显著。946s步骤六、找一个精简模式 现在找一个新的模式，筛选出不重要的效果，找出此模式后，画出几个图 形将效果可视化，然后评估此精简模式，作一个残差分析。 1.选择 stat doe
analyze factorial design。 2.点击 terms。3.设置模式 ● 从 include terms in the model up through order中，选择 2。注意 到此举将abc移动到 availble terms 中。 ● 在 selected terms 中，点击 a：temp，然后点击按钮 & 将a：temp移动 到availble terms中。 ● 重复以上操作，将ab、ac移动到 availble terms 中。 4. 点击 ok 。956s5. 点击graphs，将normal、pareto不选中。 6. 选中 histogram,normal plot，residuals versus fits和residuals versus order.点击 ok 。 7. 在主对话框中点击 ok 。966s步骤七、评估此精简模式 会话窗口的输出说明了此模式的一些优点。检查一下p列。评估此模式的一 个好的标准是看p值。如果所有的p值都小于设定的α 值，那么可以说此模式很好。在这里，α 值为0.05。fractional factorial fit estimated effects and coefficients for yield (coded units) term effect coef stdev coef t constant 74.81 2.107 35.51 pressure 14.13 7.06 2.107 3.35 catalyst -30.37 -15.19 2.107 -7.21 pressure*catalyst -13.38 -6.69 2.107 -3.17 analysis of variance for yield (coded units) source df seq ss adj ss main effects 2 8.62 2-way interactions 1 715.56 715.56 residual error 12 852.25 852.25 pure error 12 852.25 852.25 total 15 6056.44p 0.000 0.006 0.000 0.008adj ms .56 71.02 71.02f 31.60 10.08p 0.000 0.008在这个模式中，所有的p值都小于0.05，表明此模式适合于进一步的实验检 验。此模式相当简单而且几乎适合于此模式中所有的条件。残差错误只增加了 一点点。可以通过残差图进一步检查此模式。残差是实际的产量减去预计的产 量。此直方图没有呈现出正态分 布，信息不够，难于下结论，仅仅 凭16个点，难于解释直方图。在零值的两侧，分布很 好，没有明显的增长或减少 的趋势。976s此图相当好此图呈现出一定规则的散布， 但同样仅仅凭16个点很难做出判断。步骤八、结论 显示因子图 现在生成两个图形使效果可视化—主效果图和交互作用图。1.选择 stat doe
factorial plots。2.选中 main effects 然后点击 setup。986s3.在 responses 中，输入 yield 。 4.下一步选择条件：
在 availble 中点击 b：pressure，然后点击 & 按钮。将 b： pressure 移动到 selected 框中。
重复此操作，将 c：catalyst 移动到 selected 框中。点击 ok 。 5.选中 interaction，点击 setup 。 6.重复步骤3、4 7.在主因子图对话框中点击 ok 。 评估此图 首先看一下改变压力（pressure）条件或者使用催化剂a或b的基本效果图。 只包含一个因子的效果图叫做主效果图。所有效果的数值意义都显示在会话窗 口里。此点表示气压 设为高值时的 平均值。 此点表示气压 设为低值时的 平均值。此虚线表示所 有产量的平均 值。此图显示了两个主效果——一个是气压、一个是催化剂。由此可以看出：
压力的主效果在高低设置时有所不同。
催化剂的主效果随类型的不同而有所不同。996s还可以看出：催化剂的类型的主效果比气压的主效果要大：连接催化剂a和 催化剂b的直线的斜率比连接高气压和低气压的直线斜率要小。尽管催化剂类型 对产量的影响比气压的影响要大，但是看其交互作用同样重要。交互作用可以 放大或消除主效果。 为了计算主效果，将高设置或第二水平时的产量值减去低设置或第一水平 时的产量值，得到如下结果： 因子 效果大小 解释气压催化剂14.13-30.37四大气压时产量比一大气压的产量高使用催化剂a的产量比催化剂b的产量高如果因子间没有什么交互作用的话， 此图对于解释改变因子怎样获得最大 收益已经足够了。下一步我们看看交互作用。虽然会话窗口中已经说明了有一 个显著的交互作用，我们再从图形上看看此交互作用到底是怎样的。此处表示产量 的单位。此点表示气压 设为低值、采 用催化剂a时 的产量平均值1006s交互作用图表明了改变一个因子会对另一个因子有什么样的影响。因为交 互作用能够放大或者消除主效果，所以看交互作用是非常重要的。上图所显示 的气压和催化剂的交互作用是两条斜率不同的直线。 在1大气压和4大气压下，催化剂a的产量比催化剂b的大。然而，可以看出， 在4大气压下使用催化剂a和使用催化剂b的产量差别比在1大气压下使用催化剂a 和使用催化剂b的产量差别要大得多。 为了得到最大产量，结果已经表明应当在4大气压下采用催化剂a。101&&&&
09:03:53 10:23:55 10:03:40 09:38:15 09:34:36}