嵌套量具 R&R 研究

jfwts 2022年 9月 2日00:05:37侠骨藏经评论1字数 17442阅读58分8秒阅读模式

嵌套量具 R&R 研究 概述

概述

当每个操作员不能测量所有部件时,可使用 嵌套量具 R&R 研究 评估测量系统中的变异。

例如,您有 3 名操作员和 15 个部件。操作员 A 测量部件 1–5 两次,操作员 B 测量部件 6–10 两次,操作员 C 测量部件 11–15 两次。每个部件对于操作员而言都是唯一的,因此,不能有多名操作员测量同一部件。

在何处可找到此分析

要执行嵌套量具 R&R 研究,请选择统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套)

什么情况下使用备择分析

  • 要在具有交叉数据的情况下分析测量系统,请使用 交叉量具 R&R 研究
  • 要在具有其他因子或非平衡设计的情况下分析测量系统,请使用 扩展量具 R&R 研究

嵌套量具 R&R 研究 的数据注意事项

数据注意事项

要确保结果有效,请在收集数据、执行分析和解释结果时注意以下准则。

应按随机顺序收集数据
因为部件嵌套在操作员中,所以数据无法完全随机化。但是,要尽可能地随机化数据,请随机选择操作员,随机选择部件,然后指示操作员按随机顺序测量他们的部件。
请选择表示过程变异的实际或预期极差的部件。
从整个过程范围选择部件可以增加对部件间变异进行良好估计的可能性。例如,不要测量连续部件、来自单个班次或单个生产线的部件或来自一批拒绝品的部件。
部件嵌套在操作员中
当一个因子的每个水平仅与其他因子的一个水平一起出现时,这两个因子就是嵌套的。例如,如果两位操作员测量两组不同但相似的部件,部件会嵌套在操作员下面,并用“部件 (操作员)”来表示。
有关嵌套因子的详细信息,请转到量具 R&R 研究的因子类型。
操作员和部件必须是随机的
当某因子具有许多个可能的水平时,该因子为随机因子,但是在数据中仅包含这些水平的随机样本。
例如,当为研究选择的部件表示生产过程中的所有可能的部件时,部件是随机因子。
有关随机因子的详细信息,请转到量具 R&R 研究的因子类型。

嵌套量具 R&R 研究示例

示例

某工程师想要监控陶瓷组件的抗冲击强度。工程师随机选择 30 个表示过程变异预期极差的样本,并将 10 个随机样本提供给 3 名随机选择的操作员。这 3 名操作员每人测量 10 个不同样本的温度,每个样本测量 2 次,总共测量 60 次。每个部件(样本)对于操作员而言都是唯一的;不会有两名操作员测量同一部件。

由于此测量嵌套在操作员内,因此,此工程师通过执行嵌套量具 R&R 研究来评估因测量系统导致的测量值的变异性。

  1. 打开样本数据 陶瓷部件.MTW.
  2. 选择统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套)
  3. 部件号或批号中,输入部件
  4. 操作员中,输入操作员
  5. 测量数据中,输入响应
  6. 单击确定

解释结果

在方差分析表中,操作员的 p 值为 0.773。由于 p 值较大,工程师无法否定原假设,因此得出如下结论:强度的平均测量值可能与执行测量的操作员无关。但是,部件(操作员)的 p 值为 0.000,小于 0.05。嵌套在每名操作员内的不同部件的平均测量值显著不同。

使用“研究变异 %”将测量系统变异和总体变异进行比较。总体量具 R&R 等于研究变异的 27.86%,部件间变异等于 97.15%。根据应用场合,总体量具 R&R 变异可能是可接受的。有关更多信息,请转到我的测量系统是否可以接受?。

同样,此测量系统可以区分 5 个不同的类别。此结果意味着测量系统可以区分部件。按照 AIAG 的要求,您需要至少 5 个可区分类别才能得到满足要求的测量系统。有关更多信息,请转到可区分类别数的使用。

这些图形还提供有关测量系统的如下信息:

  • 在变异分量图中,大部分变异是用部件间变异解释的。
  • 在“R 控制图(按操作员)”中,所有的数据都受控制,这表明三名操作员都按照一致的方式进行测量。
  • 在“Xbar 控制图(按操作员)”中,多个点位于控制限外部。因此,许多变异是因部件间的差异所致。
  • “按部件”图显示部件之间的差异很大。
  • 在“按操作员”图中,每名操作员的测量值按照大致相同的量变化。而且,部件平均值的变化率很小。尽管始终存在一定的变异,但是数据表明操作员按照相似的方式测量部件。

量具 R&R 研究 - 嵌套方差分析

响应 的量具 R&R (嵌套) 来源 自由度 SS MS F P 操作员 2 2.618 1.30922 0.2594 0.773 部件 (操作员) 27 136.285 5.04758 34.5709 0.000 重复性 30 4.380 0.14601 合计 59 143.283
方差分量 方差分量 来源 方差分量 贡献率 合计量具 R&R 0.14601 5.62 重复性 0.14601 5.62 再现性 0.00000 0.00 部件间 2.45079 94.38 合计变异 2.59679 100.00
量具评估 研究变异 %研究变 来源 标准差(SD) (6 × SD) 异 (%SV) 合计量具 R&R 0.38211 2.29265 23.71 重复性 0.38211 2.29265 23.71 再现性 0.00000 0.00000 0.00 部件间 1.56550 9.39300 97.15 合计变异 1.61146 9.66874 100.00

可区分的类别数 = 5

响应 的量具 R&R (嵌套) 报告

嵌套量具 R&R 研究

嵌套量具 R&R 研究 输入数据

输入数据

统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套)

完成以下步骤指定量具 R&R 研究数据。

  1. 部件号或批号中,输入包含部件或批次名称或编号的列。
  2. 操作员中,输入包含操作员名称或编号的列。
  3. 测量数据中,输入包含观测到的测量值的列。

 

在此工作表中,部件包含部件号标识符,操作员包含操作员标识符,响应包含每个部件的测量值。每个部件对于操作员而言都是唯一的;不会有两名操作员测量同一部件。
C1 C2 C3
部件 操作员 响应
1 A 15.4257
1 A 16.8677
2 A 15.5018
2 A 15.1628
6 B 13.1025
6 B 15.5494

嵌套量具 R&R 研究 指定研究信息

指定研究信息

统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套) > 量具信息

输入要用于识别此研究的信息。

量具名称
输入量具标识符。
研究日期
输入进行研究的日期。
报表人
输入报告研究信息的人员的姓名。
量具公差
输入量具公差。量具公差表示量具所具有的辨别力(也称为测量增量)。有关详细信息,请转到什么是量具公差?。
其他
输入要记录的有关此研究的其他信息。

嵌套量具 R&R 研究指定选项

指定选项

统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套) > 选项
研究变异
输入一个常量,将其乘以标准差以估计测量值中的变异。默认乘数为 6,这是包含 99.73% 的测量值所需标准差的数量。
注意

如果您只提供一个规格,并且需要包括 99.73% 的测量值,则不需要将 6 更改为 3。对于这些计算,Minitab 将常量除以 2。

过程公差
输入规格限或公差范围,以对测量系统变异和客户规格进行比较。当您输入规格限或公差范围时,Minitab 将在输出中显示“% 公差”列。“% 公差”显示测量系统变异占公差的百分比。
  • 至少输入一个规格限:输入测量值的一个或两个规格限,以对测量系统变异和客户规格进行比较。
    注意

    当您输入至少一个规格限时,Minitab 可以计算将产品误分类的概率。

    • 规格下限:输入规格下限。
    • 规格上限:输入规格上限。
  • 规格上限 - 规格下限:如果在具有规格上下限的情况下,不希望计算误报产品的概率,请选择此选项。输入规格限之间的差值。这也称为公差范围。
历史标准差

输入总变异的已知值,它等于部件间变异与测量系统变异之和。使用下拉菜单可以指定是使用历史标准差估计过程变异,还是使用研究中的部件估计过程变异。

选择使用研究中的部件估计过程变异时,Minitab 将在输出中显示“% 过程”。“% 过程”列显示每个方差分量占过程标准差的百分比。

显示错误分类的概率
如果已在过程公差下指定了至少一个规格或输入了公差范围,则选择此选项显示错误分类的概率。有关详细信息,请转到所有统计量和图形并单击“错误分类概率”。
不显示百分比贡献
指定是否显示贡献百分比。
不显示百分比研究变异
指定是否显示研究变异百分比。
在单独的图页上绘制图表,每页一张图
分 别显示每个图形。默认情况下,Minitab 一起显示 6 个图形。
标题
您可以为量具研究输入自定义标题。

嵌套量具 R&R 研究 指定置信区间

指定置信区间

统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套) > 置信区间

指定是否包括置信区间并指定置信水平和类型。

包括置信区间
在量具 R&R 表中显示置信区间。

置信水平
输入置信区间的置信水平。通常,置信水平为 95% 即可。
置信区间
选择要显示的置信区间或边界类型:

  • 双侧:显示具有置信下限和置信上限的置信区间。
  • 下限:仅显示下限。下限是可能小于参数的值。
  • 上限:仅显示上限。上限是可能大于参数的值。

嵌套量具 R&R 研究的存储统计量

存储统计量

统计 > 质量工具 > 量具研究 > 量具 R&R 研究(嵌套) > 存储

可以将分析统计量保存到工作表中,以便可以将它们用在其他分析、图形和宏中。Minitab 将选定的统计量存储在最后一个数据列后面。存储列的名称以一个数字结尾,如果您将同一个统计量存储多次,结尾的数字会递增。

存储可用的统计量取决于您为分析选择的选项。有关其中任何统计量的详细信息,请转到嵌套量具 R&R 研究 的所有统计量和图形。

方差分量和百分比贡献

方差分量和百分比贡献
Minitab 将这些值存储在 VarComp1VarComp2 列中。
方差分量的置信区间
Minitab 将置信区间值存储在 IVarComp1IVarComp2 中。
百分比贡献的置信区间
Minitab 将置信区间值存储在 IPVarComp1IPVarComp2 中。

标准差

标准差
Minitab 将这些值存储在 StdDev1 列中。
标准差的置信区间
Minitab 将置信区间值存储在 SIStdDev1SIStdDev2 列中。

与公差之比(%)

与公差之比(%)
Minitab 将这些值存储在 Tolerance1 列中。
与公差之比(%)的置信区间
Minitab 将置信区间值存储在 ITolerance1ITolerance2 列中。

测量系统变异及百分比

测量系统变异及百分比
Minitab 将这些值存储在 StudyVar1StudyVar2 列中。
研究变异的置信区间
Minitab 将置信区间值存储在 IStudyVar1IStudyVar2 列中。
百分比研究变异的置信区间
Minitab 将置信区间值存储在 IPStudyVar1IPStudyVar2 列中。

与过程变异之比(%)

与过程变异之比(%)
Minitab 将这些值存储在 Process1 列中。
百分比过程变异的置信区间
Minitab 将置信区间值存储在 IProcess1IProcess2 列中。

可区分类别数

可区分类别数
Minitab 将这些值存储在 Distinct1 列中。
可区分类别数的置信区间
Minitab 将置信区间值存储在 IDistinct1IDistinct2 中。

解释 嵌套量具 R&R 研究 的主要结果

主要结果

完成以下步骤解释嵌套量具 R&R 研究。主要输出包括变异评估、测量和测量变异的图表。

步骤 1:使用方差分析表识别显著因子和交互作用

使用方差分析表识别哪些变异源是显著的。方差分析表的“来源”列中包含以下项:

  • 操作员:来自操作员的变异。
  • 部件(操作员):来自嵌套在每个操作员中的部件的变异。
  • 错误或重复性:非部件或操作员解释的变异。

量具 R&R 研究 - 嵌套方差分析

响应 的量具 R&R (嵌套) 来源 自由度 SS MS F P 操作员 2 2.618 1.30922 0.2594 0.773 部件 (操作员) 27 136.285 5.04758 34.5709 0.000 重复性 30 4.380 0.14601 合计 59 143.283
方差分量 方差分量 来源 方差分量 贡献率 合计量具 R&R 0.14601 5.62 重复性 0.14601 5.62 再现性 0.00000 0.00 部件间 2.45079 94.38 合计变异 2.59679 100.00
量具评估 研究变异 %研究变 来源 标准差(SD) (6 × SD) 异 (%SV) 合计量具 R&R 0.38211 2.29265 23.71 重复性 0.38211 2.29265 23.71 再现性 0.00000 0.00000 0.00 部件间 1.56550 9.39300 97.15 合计变异 1.61146 9.66874 100.00

可区分的类别数 = 5

响应 的量具 R&R (嵌套) 报告

主要结果:P

在此示例中,操作员的 p 值为 0.773。因为 p 值大于 0.05,您将无法否定原假设,并可以断定平均强度测量值可能不取决于进行测量的操作员。但是部件(操作员)的 p 值是 0.000 且小于 0.05。在每个操作员自身嵌套的各个部件的平均测量值显著不同。

步骤 2:评估每个测量误差源的变异性。

使用方差分量 (VarComp) 和贡献率评估每个测量误差源的变异性。来源如下:

  • 合计量具 R&R:是重复性和再现性方差分量之和。
  • 重复性:是指当同一操作员测量同一批次的部件时产生的测量变异性。
  • 再现性:不同操作员测量部件时产生的测量变异性。
  • 部件间:由于部件的不同所引起的测量变异性。
注意

如果操作员只能测量部件一次(如使用破坏性试验),您必须能够假定一个批次中所有部件的相同程度高到足以声称其为相同的部件。如果无法做出此假定,则一个批次中部件之间的变异将掩盖测量住系统变异。

理想情况下,只有很少的变异性应由重复性与再现性导致。而绝大部分变异性都应由部件之间的差异引起。

量具 R&R 研究 - 嵌套方差分析

响应 的量具 R&R (嵌套) 来源 自由度 SS MS F P 操作员 2 2.618 1.30922 0.2594 0.773 部件 (操作员) 27 136.285 5.04758 34.5709 0.000 重复性 30 4.380 0.14601 合计 59 143.283
方差分量 方差分量 来源 方差分量 贡献率 合计量具 R&R 0.14601 5.62 重复性 0.14601 5.62 再现性 0.00000 0.00 部件间 2.45079 94.38 合计变异 2.59679 100.00
量具评估 研究变异 %研究变 来源 标准差(SD) (6 × SD) 异 (%SV) 合计量具 R&R 0.38211 2.29265 23.71 重复性 0.38211 2.29265 23.71 再现性 0.00000 0.00000 0.00 部件间 1.56550 9.39300 97.15 合计变异 1.61146 9.66874 100.00

可区分的类别数 = 5

响应 的量具 R&R (嵌套) 报告

主要结果:方差分量、贡献百分比

合计量具 R&R 的贡献百分比是 5.62%,部件间变异是 94.38%。当部件间变异的贡献百分比很高时,测量系统可以准确区分各个部件。

量具 R&R 研究 - 嵌套方差分析

响应 的量具 R&R (嵌套) 来源 自由度 SS MS F P 操作员 2 2.618 1.30922 0.2594 0.773 部件 (操作员) 27 136.285 5.04758 34.5709 0.000 重复性 30 4.380 0.14601 合计 59 143.283
方差分量 方差分量 来源 方差分量 贡献率 合计量具 R&R 0.14601 5.62 重复性 0.14601 5.62 再现性 0.00000 0.00 部件间 2.45079 94.38 合计变异 2.59679 100.00
量具评估 研究变异 %研究变 来源 标准差(SD) (6 × SD) 异 (%SV) 合计量具 R&R 0.38211 2.29265 23.71 重复性 0.38211 2.29265 23.71 再现性 0.00000 0.00000 0.00 部件间 1.56550 9.39300 97.15 合计变异 1.61146 9.66874 100.00

可区分的类别数 = 5

响应 的量具 R&R (嵌套) 报告

主要结果:研究变异百分比

使用研究变异百分比(% 研究变异)可以将测量系统变异与总体变异进行比较。% 研究变异使用过程变异,过程变异等于 6 乘以过程标准差。如果输入公差值,Minitab 将显示“% 公差”列,如果输入历史标准差,Minitab 将显示“% 过程”列。

根据 AIAG 的原则,如果测量系统的变异小于过程变异的 10%,则测量系统是可以接受的。合计量具 R&R 是研究变异的 23.71%。根据应用程序的情况,合计量具 R&R 变异可能是可以接受的。改进测量系统的纠正措施可能包括对操作员进行培训或者使用更好的量具。有关详细信息,请转到我的测量系统是否可以接受?。

 

嵌套量具 R&R 研究
主要结果:变异分量图

变异图的分量显示与测量误差源的变异。如果您输入公差值,Minitab 将显示“% 公差”条;如果您输入历史标准差,Minitab 将显示“% 过程”条。

此图显示大部分变异来自于部件间的变异,说明测量系统的大部分变异是由部件间的差异导致的。

 

步骤 3:检查图形以获得有关量具研究的更多信息

量具 R&R 图形还提供有关测量系统的如下信息:

变异分量图
显示最大的变异分量是否为部件之间的变异。
在可接受的测量系统中,最大的变异分量为部件之间的变异。
按操作员显示的 R 控制图
显示是否有任何点位于控制上限上方。
如果操作员以一致的方式测量部件,则点将落在控制限内。
按操作员显示的 Xbar 控制图
显示大多数点是否落在控制限外部。
您为量具 R&R 研究选择的部件应当表示典型的部件间变异性。因此,部件平均值之间的变异应当更大,图形中也应当显示大多数点落在控制限外部。
按部件(操作员)的测量值图形
显示每位操作员对于每个部件的多个测量值是否比较接近,这表明了部件间变异与部件内变异。
如果每个部件的多个测量值比较接近,则表明由同一位操作员针对同一个部件得到的测量值之间的变异很小。
操作员测量值图
显示操作员之间的差异是否小于部件之间的差值。
跨不同操作员的水平直线表示每位操作员的测量值均值相似。理想情况下,每位操作员的测量值之间的差值应相同。
嵌套量具 R&R 研究

嵌套量具 R&R 研究 的所有统计量和图形

所有统计量和图形

为嵌套量具 R&R 研究提供的每个统计量和图形查找定义和解释指南。

自由度

每个 SS(平方和)的自由度 (DF)。通常,DF 用于测量计算每个 SS 时可用的信息量。

SS

平方和 (SS) 是各个平方距离之和,用来度量来自不同源的变异性。合计 SS 表明来自总体平均值的数据中的变异数量。SS 操作员指示每位操作员的平均测量值和总体平均值之间的变异量。

SS 合计 = SS 操作员 + SS 部件(操作员)+ SS 重复性

MS

均方 (MS) 是来自各个源的数据中的变异性。MS 考虑了不同源的水平数或可取值的个数不同在计算中的影响。

MS = SS/每一种变异源的 DF 值

F

该统计量用于确定操作员或部件(操作员)是否会对测量结果产生显著影响。

F 统计量越大,该因子在响应或测量变量的变异性中所起的作用就越大。

P

如果原假设成立,P 值就是获得至少与从样本计算的值一样极端的检验统计量(如 F 统计量)的概率。

解释

使用方差分析表中的 p 值可以确定平均测量值之间是否存在显著差异。

如果 p 值较小,则说明所有操作员或部件(操作员)均值均相同的假定可能不为真。

要确定平均测量值之间是否存在显著差异,可将 p 值与显著性水平(表示为 α 或 alpha)进行比较,以评估原假设。原假设声明组均值均相等。通常,显著性水平为 0.05 即可。显著性水平为 0.05 时,表示实际不存在差异时断定存在差异的风险为 5%。

P 值 ≤ α:至少一个均值存在显著差异
如果 p 值小于或等于显著性水平,则可否定原假设并得出至少有一个均值与其他均值显著不同的结论。例如,至少一位操作员的测量结果不同。
P 值 > α:均值之间不存在显著差异
如果 p 值大于显著性水平,则无法否定原假设,因为您没有足够的证据断定总体均值存在差异。例如,您无法得出操作员的测量结果不同的结论。
但是,也不能得出均值相同的结论。可能存在差异,但您的检验功效不够大,检测不到差异。

方差分量

方差分量是针对 ANOVA 表中每个源的预估方差分量。

解释

使用方差分量评估每个测量误差源的变异性。

对于可接受的测量系统来说,最大的变异分量为部件之间的变异。如果重复性和再现性对于变异量贡献很大,您需要研究问题的来源并采取纠正措施。

方差分量贡献率

% 贡献是每个方差分量占总体变异性的百分比。它等于每个源的方差分量除以总变异,再乘以 100。

解释

使用贡献率评估每个测量误差源的变异性。

对于可接受的测量系统来说,最大的变异分量为部件之间的变异。如果重复性和再现性对于变异量贡献很大,您需要研究问题的来源并采取纠正措施。

标准差(SD)

标准差 (SD) 是每个变异来源的标准差。此标准差等于该源方差分量的平方根。

标准差是一种测量变异的便利方式,因为标准差具有与部件测量和公差相同的单位。

研究变异 (6 * SD)

研究变异等于每个变异来源的标准差与 6 或您在研究变异中指定的乘数相乘。

通常,过程变异会定义为“6s”,其中“s”是作为总体标准差(表示为 σ 或 sigma)估计值的标准差。当数据正态分布时,接近 99.73% 的数据处于均值的 6 个标准差内。要定义不同的数据百分比,请使用其他标准差乘数。例如,如果要知道 99% 的数据处于哪个范围,可使用乘数 5.15 而非默认乘数 6。

%研究变异 (%SV)

% 研究变异等于每个变异来源的研究变异除以总研究变异,然后再乘以 100。

% 研究变异是该源的计算方差分量 (VarComp) 的平方根。因此,VarComp 的 % 贡献之和等于 100,但 % 研究方差值之和不等于 100。

解释

利用 %研究变异将测量系统变异与总体变异做比较。如果利用测量系统评估过程改进情况,如降低部件之间的变异性,则 %研究变异可以对测量精度进行较好地估计。如果要评估测量系统的能力以评估部件与规格的比较,则 %公差是合适的度量。

%公差 (SV/Toler)

%公差作为每个来源的研究变异进行计算,通过除以过程公差并乘以 100 计算得出。

如果输入公差,Minitab 则会计算 %公差,它将测量系统变异与规格相比。

解释

利用 %公差可相对于规格评估部件。如果利用测量系统评估过程改进情况,如降低部件之间的变异性,则适合使用 %研究变异度量。

% 过程 (SV/Proc)

如果您输入历史标准差但使用研究中的部件估计过程变异,则 Minitab 会计算 % 过程。% 过程将测量系统变异与历史过程变异进行比较。% 过程等于每个源的研究变异除以历史过程变异,再乘以 100。默认情况下,过程变异等于 6 乘以历史标准差。

如果您使用历史标准差估计过程变异,Minitab 将不显示 % 过程,因为 % 过程与 % 研究变异相同。

95% 置信区间

95% 置信区间 (95% CI) 是可能包含每个测量错误度量真实值的值范围。

Minitab 提供了方差分量的置信区间、方差分量的贡献百分比、标准差、研究变异、研究变异百分比、公差百分比和可区分类别数。

解释

因为数据样本是随机的,所以两个量具研究不可能得出相同的置信区间。但是,如果多次重复进行您的研究,则特定百分比的包含未知真实测量的结果置信区间将发生错误。这些包含参数的置信区间的百分比是区间的置信水平。

例如,具有 95% 的置信水平,置信区间便可以有 95% 的置信度包含真实值。置信区间有助于评估结果的实际意义。使用专业知识确定置信区间是否包含对您有实际意义的值。如果置信区间范围太广而不实用,则请增大样本大小。

假定 VarComp 的重复性是 0.044727,对应的 95% CI 是 (0.035, 0.060)。从该数据计算得出的估计重复性变异是 0.044727。置信区间 0.035 到 0.060 包含真实重复性变异的置信度为 95%。

可区分类别数

可区分类别数是在量具 R&R 研究中使用的度量,用于识别测量系统检测测量特征中的差异的能力。可区分类别数表示跨越产品变异极差的非重叠置信区间数,由您选择的样本数定义。可区分类别数还表示您的测量系统可识别的过程数据中的分组数。

解释

由汽车工业行动组织 (AIAG) 出版的《测量系统分析手册》1声明了 5 个或更多个类别,用于指示可接受的测量系统。如果可区分类别数小于 5,则测量系统可能不具备足够的分辨力。

通常,当可区分类别数小于 2 时,测量系统将没有控制过程的值,因为系统无法对部件进行区分。当可区分类别数是 2 时,仅可将部件划分为两组,如高和低。当可区分类别数是 3 时,可以将部件划分为 3 组,如高、中和低。

有关详细信息,请转到使用可区分类别数。

误分类概率(P)

当您指定至少一个规格限时,Minitab 可以计算将产品误分类的概率。因为量具的变异,部件的测量值不会总是等于部件的真实值。测量值和真实值之间的差异造成了对部件误分类的可能性。

Minitab 可以计算误分类的联合概率和条件概率。

联合概率
如果以前对部件的可接受情况不了解,请使用联合概率。例如,您正在生产线中抽样,但并不了解具体每个部件合格与否。您可能会做出两种误分类:

  • 可能某部件是不合格的,而您接受了它。
  • 可能某部件是合格的,而您拒绝了它。
条件概率
如果以前对部件的可接受情况有所了解,请使用条件概率。例如,您正在从一堆返工品中或从即将作为合格产品发货的一批产品中抽样。您可能会做出两种误分类:

  • 可能会接受某个从一堆需要返工的不合格产品中抽样的部件(也称为误接受)。
  • 可能会拒绝某个从一堆即将发货的合格产品中抽样的部件(也称为误拒绝)。

解释

由三名操作员测量十个部件,每个部件测量三次。下图显示了与规格限制相比较,测量值的分布情况。一般而言,生产过程拥有越多变异且生产越多接近于规格限制的部件,误分类的概率就越高。

嵌套量具 R&R 研究

Gage R&R Study - ANOVA Method

Two-Way ANOVA Table With Interaction

Source DF SS MS F P Part 9 88.3619 9.81799 492.291 0.000 Operator 2 3.1673 1.58363 79.406 0.000 Part * Operator 18 0.3590 0.01994 0.434 0.974 Repeatability 60 2.7589 0.04598 Total 89 94.6471

α to remove interaction term = 0.05

Two-Way ANOVA Table Without Interaction

Source DF SS MS F P Part 9 88.3619 9.81799 245.614 0.000 Operator 2 3.1673 1.58363 39.617 0.000 Repeatability 78 3.1179 0.03997 Total 89 94.6471

Gage R&R

%Contribution Source VarComp (of VarComp) Total Gage R&R 0.09143 7.76 Repeatability 0.03997 3.39 Reproducibility 0.05146 4.37 Operator 0.05146 4.37 Part-To-Part 1.08645 92.24 Total Variation 1.17788 100.00

Process tolerance = 2.5

Study Var %Study Var %Tolerance Source StdDev (SD) (6 × SD) (%SV) (SV/Toler) Total Gage R&R 0.30237 1.81423 27.86 72.57 Repeatability 0.19993 1.19960 18.42 47.98 Reproducibility 0.22684 1.36103 20.90 54.44 Operator 0.22684 1.36103 20.90 54.44 Part-To-Part 1.04233 6.25396 96.04 250.16 Total Variation 1.08530 6.51180 100.00 260.47

Number of Distinct Categories = 4

Probabilities of Misclassification

Joint Probability Part is bad and is accepted 0.037 Part is good and is rejected 0.055
Conditional Probability False Accept 0.151 False Reject 0.073

Gage R&R for Measurement

接受一个不合格部件的联合概率是 0.037。拒绝一个合格部件的联合概率是 0.055。

误接受的条件概率(也即在重新检验过程中接受实际不符合规格的部件)是 0.151。误拒绝的条件概率(也即在重新检验过程中拒绝实际符合规格的部件)是 0.073。

 

变异分量图

变异分量控制图是量具 R&R 研究结果的图形总结。

图中所表示的变异源包括:

  • 合计量具 R&R:重复性方差分量和再现性方差分量之和。
  • 重复性:同一操作员测量同一批次中的部件时的测量变异性。
  • 再现性:不同操作员测量部件时的测量变异性。
  • 部件间:由于部件的不同所引起的测量变异性。

解释

各彩色条形图分别表示:

% 贡献
% 贡献是每个方差分量占总变异的百分比。它等于将每个源的方差分量除以总变异,再乘以 100。
% 研究变异
% 研究变异是每个源的研究变异所占的百分比。它等于每个源的研究变异除以总研究变异,然后再乘以 100。
% 公差
% 公差将测量系统变异与规格进行比较。它等于每个源的研究变异除以过程公差,再乘以 100。
当您指定过程公差范围或规格限值时,Minitab 将计算此值。
% 过程
% 过程将测量系统变异与过程总变异进行比较。它等于每个源的研究变异除以历史过程变异,再乘以 100。
当您指定历史标准差并选择使用研究中的部件估计过程变异时,Minitab 将计算此值。

在可接受的测量系统中,最大的变异分量为部件之间的变异。

嵌套量具 R&R 研究

R 控制图

R 控制图是一种极差控制图,显示操作员的一致性。

R 控制图包含以下元素:

标绘点
对于每位操作员,每个部件的最大值和最小值之间的差值。R 控制图上的点由操作员标绘,因此您可以看出每位操作员测量的一致程度。
中心线 (Rbar)
过程的总平均值(即所有样本极差的平均值)。
控制限(LCL 和 UCL)
预计的样本极差变异量。Minitab 使用样本内的变异计算控制限。
注意

如果每位操作员对每个部件测量 9 次或更多次,Minitab 显示 S 控制图(而非 R 控制图)。

解释

小的平均极差表示测量系统的变异性低。高于控制上限 (UCL) 的点表示操作员对部件的测量不一致。UCL 的计算包括每位操作员对部件的测量次数以及部件之间的变异。如果操作员对部件的测量一致,那么与研究变异相比,最高测量值与最低测量值之间的极差就很小,且各个点都应该处于受控状态。

嵌套量具 R&R 研究

Xbar 控制图

Xbar 控制图将部件间变异与重复性分量进行比较。

Xbar 控制图包含以下元素。

标绘点
每个部件的平均测量值,由每位操作员标绘。
中心线 (Xbar)
所有操作员测量的所有部件测量值的整体平均值。
控制限(LCL 和 UCL)
控制限基于每个平均值中的重复性估计和测量次数。

解释

为量具 R&R 研究选择的部件应当代表所有可能的部件。因此,该图应当表示期望在部件均值间观察到的变异比重复性变异多。

通常,该图的控制限较窄,其中的许多失控点表示具有低变异的测量系统。

嵌套量具 R&R 研究

部件(操作员)对比图

“按部件(操作员)”图形显示在研究过程中所测量并按部件排列的所有测量值。此图显示因子水平之间的差异。通常,量具 R&R 研究按部件和按操作员排列测量值。但是,使用扩展的量具 R&R 研究,可以在图中绘制其他因子。

在此图中,点表示测量值,而带十字标的圆形符号表示均值。连接线连接每个因子水平的平均测量值。

注意

如果每个水平的观测值大于 9,Minitab 将显示箱线图而非单值图。

解释

如果每个部件的多个测量值的变化尽可能小(一个部件的点接近另一个部件的点),则表示测量系统的变异小。而且,部件的平均测量值的变异足够大时才会显示部件有差异和表示整个过程极差。

嵌套量具 R&R 研究

“按操作员”图

“按操作员”图形显示研究过程中所测量并按操作员排列的所有测量值。此图显示因子水平之间的差异。通常,量具 R&R 研究按部件和按操作员排列测量值。但是,使用扩展的量具 R&R 研究,可以在图中绘制其他因子。

注意

如果每个操作员大于 9 个观测值,Minitab 将显示箱线图而非单值图。

解释

跨操作员的直线水平线表示每位操作员的测量值均值是相似的。理想情况下,每位操作员的测量值之间的差值应相同。

嵌套量具 R&R 研究

1 汽车工业行动组织 (AIAG) (2010)。测量系统分析参考手册,第 4 版。克莱斯勒、福特、通用汽车供方质量要求特别工作组。

嵌套量具 R&R 研究 中的方差分析表的方法和公式

方差分析表

请选择您所选的方法或公式。

平方和

平方和 (SS) 是各个平方距离之和,用来度量来自不同源的变异性。

SS操作员
嵌套量具 R&R 研究
说明
a 部件数
n 仿行数
嵌套量具 R&R 研究
每个操作员的均值
嵌套量具 R&R 研究
总均值
SS部件(操作员)
嵌套量具 R&R 研究
说明
n 仿行数
嵌套量具 R&R 研究
操作员 j 内的每个部件 i 的均值
嵌套量具 R&R 研究
操作员 j 的均值
SS重复性
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
每个观测值
嵌套量具 R&R 研究
操作员 j 内的每个部件 i 的均值
SS合计
SS合计 = SS操作员 + SS部件(操作员) + SS重复性

自由度

每个 SS(平方和)的自由度 (DF)。通常,DF 用于测量计算每个 SS 时可用的信息量。

DF操作员
嵌套量具 R&R 研究
说明
b 操作员数
DF部件(操作员)
嵌套量具 R&R 研究
说明
a 每个操作员测量的部件数
b 操作员数
DF重复性
嵌套量具 R&R 研究
说明
a 部件数
b 操作员数
n 仿行数
DF合计
嵌套量具 R&R 研究
说明
a 部件数
b 操作员数
n 仿行数

均方

均方 (MS) 是来自各个源的数据中的变异性。MS 考虑了不同源的水平数或可取值的个数不同在计算中的影响。

MS操作员
嵌套量具 R&R 研究
MS 部件(操作员)
嵌套量具 R&R 研究
MS 重复性
嵌套量具 R&R 研究

F

F 统计量用于确定操作员或部件(操作员)是否会对测量结果产生显著影响。

F操作员
嵌套量具 R&R 研究
F部件(操作员)
嵌套量具 R&R 研究

p 值

如果原假设成立,P 值就是获得至少与从样本计算的值一样极端的检验统计量(如 F 统计量)的概率。

嵌套量具 R&R 研究 中可区分类别数的方法和公式

可区分类别数

请选择您所选的方法或公式。

可区分类别数

可区分类别数表示要跨越产品变异极差的非重叠置信区间数。您也可以将其视为测量系统可识别的过程数据中的分组数。

嵌套量具 R&R 研究

 

然后,Minitab 会截断该值的尾数,除非该值小于 1。如果小于 1,Minitab 会将可区分类别数设为等于 1。

置信区间

假设 L 和 U 是量具方差和总方差比率的下限和上限,则可区分类别数的下限和上限为:

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

注意

L 和 U 必须介于 (0, 1) 之间。如果 L 和 U 超出该范围,则可区分类别数的下限和上限将缺失。

嵌套量具 R&R 研究 中的量具 R&R 表的方法和公式

量具 R&R 表

请选择您所选的方法或公式。

VarComp

每个源贡献的方差组件。

重复性
可重复性 = MS 可重复性
再现性
嵌套量具 R&R 研究
说明
a 部件数
b 操作员数
n 重复数
合计量具 R&R
总量根 R&R = 可重复性 = 可重复性
部件间
嵌套量具 R&R 研究
说明
n 重复数
合计变异
总变化 = 瓦尔康普总量级 R&R = VarCompPart 到零件

%贡献

% 贡献是每个方差分量占总体变异性的百分比。它等于每个源的方差分量除以总变异,再乘以 100。贡献率通过每个方差分量除以总变异然后乘以 100 计算得出。此列中的百分比增加到 100。

嵌套量具 R&R 研究

 

标准差

标准差 (SD) 是每个变异来源的标准差。此标准差等于该源方差分量的平方根。

在指定将使用历史标准差估计过程变异时,Minitab 执行以下操作:

  • 如果 嵌套量具 R&R 研究则总标准差为 σ,嵌套量具 R&R 研究
  • 否则,Minitab 使用数据估计总标准差和 嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究 历史标准差
嵌套量具 R&R 研究 根据数据计算的总量具标准差
嵌套量具 R&R 研究 部件间标准差

研究变异

研究变异等于每个变异来源的标准差与 6 或您在研究变异中指定的乘数相乘。

通常,过程变异会定义为“6s”,其中“s”是作为总体标准差(表示为 σ 或 sigma)估计值的标准差。当数据正态分布时,接近 99.73% 的数据处于均值的 6 个标准差内。要定义不同的数据百分比,请使用其他标准差乘数。例如,如果要知道 99% 的数据处于哪个范围,可使用乘数 5.15 而非默认乘数 6。

% 研究变异和置信区间

% 研究变异等于每个变异来源的研究变异除以总研究变异,然后再乘以 100。

% 研究变异是该来源的计算方差分量 (VarComp) 的平方根。因此,VarComp 的 % 贡献之和等于 100,但 % 研究方差值之和不等于 100。

嵌套量具 R&R 研究

 

置信区间

如果 L 和 U 是贡献百分比的下限和上限,则研究变异相应百分比的置信区间为:

嵌套量具 R&R 研究

 

%公差和 CI

%公差是每个分量的公差百分比。

如果给定公差(规格上限 - 规格下限),则 %公差将通过指定公差除以每个分量的“研究变异”计算得出。

嵌套量具 R&R 研究

 

如果仅给定一个规格限值,则公差百分比是单侧公差除以每个分量“研究变异”的二分之一。单侧公差是将所有测量的平均值减去给定规格限值的绝对值。

Minitab 仅当您在“选项”子对话框中输入过程公差(规格上限 - 规格下限)或某个规格限值时才显示此值。

置信区间

如果 L 和 U 是方差分量的上下限,则相应公差百分比的置信区间是:

嵌套量具 R&R 研究

 

说明
k k 是研究常量,其默认值是 6

%过程和 CI

%过程是过程变异的百分比,通过估计标准差与历史标准差的比率计算得出。

嵌套量具 R&R 研究

 

仅当在“选项”子对话框中输入历史标准差时,Minitab 才会显示此值。

置信区间

如果 L 和 U 是方差分量的上下限,则相应过程变异的置信区间是:

嵌套量具 R&R 研究

 

说明
k k 是研究常量,其默认值是 6

嵌套量具 R&R 研究 中误分类概率的方法和公式。

误分类概率

请选择您所选的方法和公式。

误分类概率(P)

当至少输入一个规格限时,Minitab 会将误分类概率作为联合概率和条件概率进行计算。

联合概率

可能某部件是不合格的,而您接受了它:

嵌套量具 R&R 研究

 

可能某部件是合格的,而您拒绝了它:

嵌套量具 R&R 研究

 

条件概率

可能某部件是不合格的,而您接受了它(误接受):

嵌套量具 R&R 研究

 

可能某部件是合格的,而您拒绝了它(误拒绝):

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

F(X,Y) 是二变量正态随机向量 (X,Y)T 的累积分布函数 (CDF):

均值,μ = (θ,θ)T

嵌套量具 R&R 研究

 

F(X) 和 F(Y) 是相应边际 CDF。

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

嵌套量具 R&R 研究 的置信区间中的方差分量的方法和公式

置信区间中的方差分量

请选择您所选的方法和公式。

量具 R&R 置信区间的常见表示法和规则

对于所有方差分量,方差分量的下限和上限不能为负值。如果使用公式计算的值为负,则它们设置为零。

对于介于 0 和 1 之间的所有比率,下限和上限也应该介于 0 和 1 之间。如果边界超出范围,则相应地设置为 0 或 1。

Notation

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
the α *100 percentile of the chi-square distribution with nq degrees of freedom
F=(nq, n=) the α *100 percentile of the F distribution with nq and nγ degrees of freedom
J the number of operators
I the number of parts
K the number of replicates

对于自由度:

运算符: n2+J +1

零件: n1=J(1+I)

复制: n4=IJ(K=1)

系统操作器 = S12

MSPart*操作员 = S22

MS 复制 = S32

再现性方差分量置信区间

Minitab 计算精确的 (1-α) *100% 置信区间的上下限。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

公式

嵌套量具 R&R 研究
嵌套量具 R&R 研究

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

再现性(或操作员)方差分量置信区间

Minitab 使用修改的大样本 (MLS) 方法、近似 (1–α) *100% 置信区间的上下限。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

公式

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

总量具方差分量置信区间

Minitab 使用修改的大样本 (MLS) 法计算近似 (1–α) *100% 置信区间的上下限。要计算单侧置信边界,请将 H 和 G 中的 α/2 替换为 α。

公式

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

部件间方差分量置信区间

Minitab 使用修改的大样本 (MLS) 方法计算近似 (1–α) *100% 置信区间的上下限。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

公式

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

总方差分量置信区间

Minitab 使用修改的大样本 (MLS) 方法计算近似 (1–α) *100% 置信区间的上下限。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

公式

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

嵌套量具 R&R 研究 的置信区间中的方差比的方法和公式

置信区间中的方差比

请选择您所选的方法或公式。

重复性方差和总方差比值的置信区间

可以采用两种计算方法。第一种方法,Minitab 使用修改的大样本 (MLS) 方法计算上下限。如果某条件在计算过程中未得到满足,则 Minitab 将使用 Satterthwaite 近似。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

MLS 方法

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

使用 MLS 方法的上下限存在的两个条件是:

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

如果这两个条件未得到满足,则 Minitab 无法使用此方法构建上下限。Minitab 将使用 Satterthwaite 近似计算上下限。

Satterthwaite 近似

保持此上下限公式不变,并对 L 和 U 进行如下定义:

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
自由度为 nq 的卡方分布的第 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

再现性方差和总方差比值的置信区间

可以采用两种计算方法。第一种方法,Minitab 使用修改的大样本 (MLS) 方法计算上下限。如果某条件在计算过程中未得到满足,则 Minitab 将使用备择近似。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

MLS 方法

通过解二次方程计算近似 (1–α) *100% 置信区间的上下限。

下限
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
上限
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究

第二种方法

如果 B2– 4AC < 0,则二次方程无解。这种情况下,Minitab 将使用第二种方法评估置信区间。通过以下方式计算近似 (1–α) * 100% 置信区间:

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

部件方差和总方差比值的置信区间

可以采用两种计算方法。第一种方法,Minitab 使用修改的大样本 (MLS) 方法计算上下限。如果某条件在计算过程中未得到满足,则 Minitab 将使用备择近似。要计算单侧置信范围,请将 H 和 G 中的 α/2 替换为 α。

MLS 方法

通过解二次方程计算近似 (1–α) * 100% 置信区间的上下限。

下限
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
上限
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究

第二种方法

如果 B2– 4AC < 0,则二次方程无解。这种情况下,Minitab 将使用第二种方法评估置信区间。通过以下方式计算近似 (1–α) * 100% 置信区间的上下限:

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

嵌套量具 R&R 研究

 

表示法

嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
嵌套量具 R&R 研究
说明
嵌套量具 R&R 研究
带有 nq 自由度的卡方分布 α * 100 个百分位数
J 操作员数
I 部件数
K 仿行数

量具方差和总方差比值的置信区间

下限 = 1 –(部件方差和总方差比值的置信区间下限)

上限 = 1 –(部件方差和总方差比值的置信区间上限)

  • 我的微信
  • 微信扫一扫
  • weinxin
  • 我的微信公众号
  • 微信扫一扫
  • weinxin
jfwts
  • 本文由 发表于 2022年 9月 2日00:05:37
  • 转载请务必保留本文链接:https://www.yingxz.cn/26.html
匿名

发表评论

匿名网友 填写信息

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen: