计数型msa：kappa分析中“期望计数”的计算方法

《计数型msa：kappa分析中“期望计数”的计算方法》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、计数型MSA：Kappa分析中“期望计数”的计算方法我们仍用MSA手册中的例子来做说明。两评价人A和B分别对随机抽取的50个零件进行测量，对每个零件的测量随机地重复了3次。设定1表示可接受的决定（即评为1类），0表示不可接受的决定（即评为0类）。测量结果用以下交叉表列出。科恩（Cohen）给出的Kappa计算公式为：MSA手册中定义P0为对角栏框中（观测）计数占总计数的比例，Pe为对角栏框中期望计数占总计数的比例。（第三版中文MSA手册中翻译有误，这里是按英文版翻译过来的）（观测）计数容易理解，如对角栏框中的44表示A和B都评

2、为0类的测量次数，97表示A和B都评为1类的测量次数，而6则表示A评为0类但B评为1类的测量次数，3则表示A评为1类但B评为0类的测量次数。因此，，P0也就是评价人A和B在测量中实际一致性的比率。那期望计数是指什么呢？假如这两位评价人都是任意地（猜测）将50个零件判定为可接受或不可接受，其结果也会达到一定程度的一致，不过这里的一致是偶然达到的，这种偶然达到的一致性比率称为偶然一致性比率，也就是Kappa计算中的Pe。当评价人A与B随机地作评价时，此两人行动一定是独立的，互不影响，从而两人同评为0类的概率为P0+×P+0，两人同

3、评为1类的概率为P1+×P+1。（这里P0+表示A评价人评为0类的测量次数占总测量次数的比率，P1+表示A评价人评为1类的测量次数占总测量次数的比率；同理P+0、P+1分别表示B评价人评为0类的测量次数占总测量次数的比率，和B评价人评为1类的测量次数占总测量次数的比率。这时偶然一致性比率Pe=P0+×P+0+P1+×P+1。因此，这里15.7和68.7分别就是两评价人同时评为0类和1类的期望计数，分别为对应单元格行总计数乘以列总计数除以总计数（即总测量次数）所得。同样31.3和34.3这两个期望计数也是按同样方式得出的。可以看

4、出，期望计数就是指两评价人因偶然性而使其判定可能达到一致的期望值或预期值。在计算kappa时，分子分母同时减去偶然一致比率Pe目的就是为了消除因偶然性造成的评价的一致性，使kappa测量一致性的程度更有效用。计数型测量系统研究的假设试验分析法(2009-07-0801:42:56)标签：测量系统分析msa计数型杂谈分类：测量系统分析MSA假设试验分析法是评估计量型测量系统风险的分析方法之一。该方法对测量系统的评价包括两个方面：一方面运用Kappa对评价人之间的一致性和评价人与标准之间的一致性进行衡量，另一方面就是对测量系统的有

5、效性进行评价。所以说，Kappa分析仅仅是计数型测量系统分析的一部分工作，要判定计数型测量系统是否合适还要对测量系统的有效性进行评价。MSA手册中通过例子的形式已经说明了如何进行计数型测量系统有效性的分析，但对于有效性的置信区间如何计算以及漏发警报率和误发警报率如何计算都没有明示的说明。我们仍然分析MSA手册中的例子（见MSA中文第三版P127表12）。通过表12记录的研究数据，我们知道参考值显示合格的零件有34件，不合格的有16件。其中，39件是A、B、C三个评价者三次测量都判定与参考值一致的；评价者A有8个零件的判定结果与

6、参考值不一致，其中5次将合格零件判定为不合格，3次将不合格零件判定为合格；评价者B有5个零件的判定结果与参考值不一致，其中2次将合格零件判定为不合格，3次将不合格零件判定为合格；评价者C则有10个零件的判定结果与参考值不一致，其中9次将合格零件判定为不合格，6次将不合格零件判定为合格。（这里要注意到判定结果与参考值不一致的零件数不一定等于将合格零件判定为不合格以及将不合格零件判定为合格的总次数，因为每个评价者分别对50个零件各测试了3次，在判定结果不一致的零件中，有1次的结果不一致的，也有2次的结果不一致的）。根据上述原始统计

7、信息，我们可以分别计算出A、B、C三个评价者测量的有效性和整个测量系统的有效性。那么我们可以算出在这个例子中：实际上根据随机抽样检验所计算出的有效性是需要估计置信区间的，按95%的置信度来估计，评价者测量的有效性和测量系统的有效性置信区间的上下限值可依据Beta分布的分位数来确定。置信上限的计算公式为：Betainv(1-α/2,B+1,A-B)置信下限的计算公式为：1-Betainv(1-α/2,A-B+1,B)在置信度为95%时，显著水平α=5%,因此1-α/2=0.975，上述公式中A为总零件数，B为判定结果与参考值一致

8、的零件数。我们通过MicrosoftExcel可以计算出评价者A、B、C测量有效性以及测量系统有效性的置信上限和下限。评价者A的95%UCL=Betainv(0.975,42+1,50-42)=93%评价者A的95%LCL=1-Betainv(0.975,50-42+1,42