JMP常用工具整理PPT学习教案

1、会计学1JMP常用工具整理常用工具整理认识数据的分布第1页/共66页认识数据的分布第2页/共66页n柏拉图展示的是根据出现频率进行排序后的数据类别，例如：返工的原因。第3页/共66页JMP里的柏拉图n请打开QUALITY CONTROL.jmp，里面的Flows是瑕疵的类型，右侧记录了瑕疵“发生的时刻”。n这组数据没有直接提供瑕疵计数。所以要用“X分组”来区隔。第4页/共66页基准（可调）JMP里的柏拉图nJMP会以第一个“发生时刻”值的瑕疵多少排列顺序作为不同时刻分隔的基准，并可自行调整。n于是，可以观察各种瑕疵在不同“作业时刻”的具体数量分布。第5页/共66页打开文件 ，看I-MR单值 移

2、动全距图控制图看的就是“异常点”第6页/共66页n因I-MR数据基于正态分布，所以“移动极差”只有超限点控制图看的就是“异常点”第7页/共66页n打开文件 SPC-X.jmp，看X-Bar 分组均值范围图n若无GROUP列，则在此输入组内个数X-Bar 分组均值范围图第8页/共66页X-Bar 分组均值范围图第9页/共66页打开P-CHART.jmp，反映了一段时间内被检验产品的不良数，其中分母不尽一致。当然，若分母一致，可自行输入样本常数。第10页/共66页数据包含以下特征:数据正态5个一组的SPC基本正常Max=85.752，Min=56.749若将LSL=57、USL=80以及LSL=5

3、7、USL=86分别输入看流程能力分组若非连续采集（比如每周分组若非连续采集（比如每周4 4算一组等），可在上方算一组等），可在上方“分组依据分组依据”里选类似里选类似SubgroupSubgroup（此时（此时SubgroupSubgroup算属性）。算属性）。JMP缺省显示缺省显示长期能力长期能力，尽管显示出的是，尽管显示出的是Cpk第11页/共66页若需单独显示传统意义上的若需单独显示传统意义上的短期能力短期能力，则需点选最下方的分组大小（缺省，则需点选最下方的分组大小（缺省5 5个一组）。个一组）。此时，前面的此时，前面的Cpkpk变成了真正的变成了真正的PpkPpk第12页/共66页

4、数据Camshaft.mtw (Minitab 工作表)n凸轮轴由机床削切凸轮轴由机床削切n数据在第数据在第3列列(Supp2)n规格是规格是 600 5n凸轮轴生产的子群大小为凸轮轴生产的子群大小为1n流程能力如何？流程能力如何？第13页/共66页流程能力分析步骤nSPC稳定性检查，稳定性检查，对异常点作出处理对异常点作出处理。n正态检定。如果数正态检定。如果数据非正态，需做拟据非正态，需做拟合转换。合转换。n执行流程能力分析执行流程能力分析。所以，建议所以，建议JMP先从先从“分析分析分布分布”开始逐步推进开始逐步推进分析。分析。第14页/共66页长期长期Sigma水平水平短期短期Sigm

5、a水平水平JMP里的“Z基准”反映出流程的西格玛水平n注意，这里的长短期有注意，这里的长短期有1.5 的水平位移。看来，是按的水平位移。看来，是按照照MOTOROLA的观点去表现的。的观点去表现的。第15页/共66页n 特性不清楚的分布特性不清楚的分布利用利用BOX-COX转换，转换成正态分布转换，转换成正态分布案例分析案例分析: 打开文件打开文件BOX-COX.JMP，其中的，其中的D3数据，规格数据，规格为为LSL=0.06，USL=0.15，请计算流程能力。，请计算流程能力。D3Percent0.140.130.120.110.100.090.080.070.0699.999959080

6、70605040302010510.1Mean = 0.05所以我们不能否定零假设。数据不能提供足够的证据否定平均强度等于2.85磅。单样本 Z: force mu = 2.85 与 2.85 的检验假定标准差 = 0.100492 平均值变量 N 平均值 标准差 标准误 95% 置信区间 Z Pforce 49 2.8463 0.1005 0.0144 (2.8182, 2.8745) -0.26 0.798第24页/共66页我们首先使用Shapiro-Wilk W检验评估正态分布。统计 基本统计量 正态性检验正态检验值0.68060.05，服从正态。第25页/共66页单样本 T: Thic

7、kness mu = 3 与 3 的检验 平均值变量 N 平均值 标准差 标准误 95% 置信区间 T PThickness 18 3.00294 0.00310 0.00073 (3.00140, 3.00448) 4.04 0.001P值d) 计算P值：分析 分布 均值检验e) 因为p=0.001 以X拟合Y 勾选 “假定等方差”第28页/共66页n首先按“不等方差”做等方差检验。正态看F检验，非正态看Levene检验，从它们各自的p值是否大于0.05，确定是否“等方差”。n如果“等方差”，则按“均值/方差分析/合并的t”；若“不等方差”，则“t检验”。P=0.1830.05，均值相等JM

8、P输出第29页/共66页利用软件检验两个相关的样本是否相等n两个相关的样本必须是随机抽取的n每个抽样总体都应该大致呈正态分布统计 基本统计量 配对 t第30页/共66页一轮胎公司认为他们新生产的轮胎的里程数较竞争者的有提高。选择了12部车，用新轮胎跑1000哩，再用竞争者的轮胎跑1000哩 。假定里程的差异服从正态分布。第31页/共66页第32页/共66页或直接看“Welch检验”第33页/共66页打开n因为数据没有进行堆栈，先运行“表堆叠”。n接着再针对堆叠表，进入“分析以X拟合Y” ；或“分析拟合模型”（可看残差） 。第34页/共66页n首先通过JMP的“分布”栏做各组数据的SPC稳定性研

9、究与正态检定。n进入“分析以X拟合Y”。n接着还是做等方差检验（“不等方差”），步骤同前面的“双样本检验”。n如果“等方差”，则看“均值/方差分析”；pF”即p值。JMP的3水平以上一元ANOVA第35页/共66页R2 =0.50712一元ANOVA的数据分析（“分析以X拟合Y” ）npF”即p值。1第36页/共66页一元ANOVA的数据分析（“分析拟合模型” ）n首先通过“分布”做各组数据的SPC稳定性研究与正态检定。n如果各组数据都是等方差，也可以直接运行JMP的“分析拟合模型”，同样可得判定各组数据是否“相等”的p值和因子显著性的R方值。12第37页/共66页n在有些实验中，我们发现当改

10、变其他因子的不同水平时，一个因子的水平的主效果有所改变。在这种情况下因子间具有交互作用。n两条直线不平行表示存在交互作用。浓度1浓度2催化剂 132385424催化剂 2交互作用1-15550454035302520c co on ns si is st te en nc ce e平平均均值值-11catalyzerY Yi ie el ld d 交交互互作作用用图图数据平均值第38页/共66页n论社会贡献，“品质”和“智商”有无交互作用？n在质量检查工作，性格沉稳与否和注意力之间是否存在交互作用？n收入水平和奖励水平之间呢？第39页/共66页多元ANOVA分析目的：提供分析具有显著交互作用的

11、2因子实验的机会 切换到工作表 Montint.jmp输出变量: Yield（望大）输入变量: Temperature温度 (Low, Med, High) Catalyst Solution催化剂浓度 (Low, Med, High)nJMP的多元方差分析依然是“分析拟合模型”，且“特质”依旧“标准最小二乘法”。nJMP的多元方差分析不区分“平衡ANOVA”和“一般线性模型GLM”。n将所有因子作“完全析因”，以观察所有交互影响。第40页/共66页多元ANOVA分析n看“效应检验”，2因子交互影响的p=0.029 Catalyst Temp*Catalyst。n各因子对Y较有利的适宜水平分别

12、是？第41页/共66页多元ANOVA分析路径第42页/共66页在在Minitab中中, 可通过以下两种方法得到一般线性回归模型可通过以下两种方法得到一般线性回归模型(最佳拟合线最佳拟合线)：统计 回归 回归统计 回归 拟合线图第43页/共66页R2 = 87.7%P-value = 0.000P-value 0.05表示在统计上回归关系显著回归关系所表达的回归关系所表达的Y的变异的的变异的87.7% 在统计上是显著的。在统计上是显著的。R2 =0.877截距和系数的截距和系数的p线性拟合线性拟合”的的“标绘残差标绘残差”第46页/共66页随时间随机随时间随机正态性正态性正态且均值为正态且均值为

13、0随机随机JMP的残差图表 看拟合直线两看拟合直线两侧侧 的点是否对称，的点是否对称， 判断正态与否。判断正态与否。 非正态则一边非正态则一边数数 据多，说明流程据多，说明流程 噪音大。噪音大。观察时序。若图形不观察时序。若图形不随机，说明或许存在随机，说明或许存在与与“时间时间”相关的因相关的因子在影响着子在影响着Y。第47页/共66页n发现发现“预测值预测值-残差残差”图有图有特殊形态，可点特殊形态，可点“二元拟二元拟合合”里的里的“特殊拟合特殊拟合”。JMP的残差图表随机随机 希望散点对称、无特希望散点对称、无特殊殊 形态地散布。从而反形态地散布。从而反映映 出残差的方差为一恒出残差的方

14、差为一恒定定 的常数。的常数。 否则，分析数据背后否则，分析数据背后有有 否异常？流程的噪音否异常？流程的噪音干干 扰大吗？是否有将不扰大吗？是否有将不同同 组的数据混放？组的数据混放？ 为求恒方差，则需对为求恒方差，则需对数数 据进行据进行转换处理：转换处理： 曲线型升阶；曲线型升阶； 喇叭型用对数转喇叭型用对数转换；换； 椭 圆 型 用椭 圆 型 用Freeman-Freeman- Tukey Tukey转换。转换。第48页/共66页n拟合置信曲线拟合置信曲线=CI（深绿色带）n单值置信曲线单值置信曲线=PI（浅绿色带）第49页/共66页n运行运行“分析分析拟合模型拟合模型”。n选择选择“

15、刻画器刻画器”，最下方会，最下方会出现出现“预测刻画器预测刻画器”。n将红线将红线十字坐标十字坐标自由移动，自由移动，可见按回归公式下可见按回归公式下Y随着随着X变变化而出现化而出现拟合值（红字）拟合值（红字），下面括弧内是该下面括弧内是该拟合值的置拟合值的置信区间信区间，及，及CI。n双击红字双击红字X，出现输入框，输，出现输入框，输入给定入给定X=1.15，Y=91.473，CI（90.947，91.999）。）。给定X的置信区间CI第50页/共66页给定X的预测区间PIn选择选择“保存列保存列预测公式预测公式”和和“保存保存列列单值置信限公式单值置信限公式” ，在数据表里，在数据表里会增

16、加三列数据。会增加三列数据。n在数据表里增加在数据表里增加1列在最后，双击列在最后，双击X空格，填入给定空格，填入给定X=1.15，边，边上会出现上会出现Y拟合值拟合值=91.473，以及，以及PI预测区间（预测区间（89.130，93.815），）， 比比CI（90.947，91.999） 要宽。要宽。第51页/共66页打开Reactor 32 Runs.jmp，该例可同时作为DOE分析和历史数据的多元回归分析。n对于百分比数值的对于百分比数值的Y，如果其，如果其分母相同，可只针对分子做连分母相同，可只针对分子做连续性回归分析。若分母不同，续性回归分析。若分母不同，可在可在“特质特质”里选里

17、选“广义线性广义线性模型模型”直接分析即可。直接分析即可。n如果数据量不大，且凭专业经如果数据量不大，且凭专业经验断定不存在验断定不存在3次以上交叉，次以上交叉，可直接选可直接选“析因次数析因次数”（阶数（阶数在下方可输入，缺省在下方可输入，缺省2阶）。阶）。但但“完全析因完全析因”下的下的R-Sq肯肯定定=1。n需要手写输入的多阶交互因子需要手写输入的多阶交互因子可选中相应几个因子后按可选中相应几个因子后按“交交叉叉”添加。添加。多元回归分析第52页/共66页n出现出现Lenth PSE，无论是多元，无论是多元回归分析还是在回归分析还是在DOE分析里都分析里都需要模式缩减。直至出现需要模式缩

18、减。直至出现“拟拟合汇总合汇总”里的里的“均方根误差均方根误差”为止。为止。多元回归分析第53页/共66页n在在“预测刻画器预测刻画器”的的“设置意愿设置意愿”设置完成，必须再执行设置完成，必须再执行“最大化意最大化意愿愿”才能出现目标才能出现目标Y下的因子预测值。下的因子预测值。n从从“模拟器模拟器”可仿真观察最终实现的可能波动范围，自然越窄越好。可仿真观察最终实现的可能波动范围，自然越窄越好。多元回归分析下的因子预测第54页/共66页JMPJMP用用“以以X X拟合拟合Y Y”做逻辑回做逻辑回归。请注意检查每个归。请注意检查每个X X和和Y Y的的属性。属性。JMPJMP数据表里的数据表里

19、的Y Y“审核出错审核出错”的基准时间缺省是以第一的基准时间缺省是以第一个个Y Y值为准，如果要像值为准，如果要像JMPJMP那那样缺省以样缺省以1 1作为基准事件，作为基准事件，则需对则需对Y Y列右键点选列右键点选“列信列信息息”里的里的“列属性列属性”里的里的“值排序值排序”，将右框里的，将右框里的“1 1”上移即可。上移即可。JMP里的逻辑回归第55页/共66页PP 值，说明模值，说明模型整体显著成立型整体显著成立。PP值排序”处理，此处X2“材料”将内部资深-1设定为缺省基准水平。练习第62页/共66页根据p值是否写出逻辑回归函式写出公式：lnp/(1-p) = -21.2 + 0.0155lnp/(1-p) = -21.2 + 0.0155报价报价 - 1.127 - 1.127MaterialMaterial（1-1-内部专家内部专家）如有2种以上的“材料”，会有相对“基础事件”p（“失败”）的另一函式。参考因子（Material）水平：1，即“内部资深”，这个是可调的。 = ep1-p-21.2+0.0155报价-1.127Material(1) =p1-p1e21.2-0.0155报价+1.127M (1) PP 值，说明模值，说明模型整体显著成立型整体显著成立。PP 值，说明模型的截距值，说明模型的截距和系数