黃世英,劉俊保,張力方,曹帥英
(中國石油蘭州化工研究中心,國家合成橡膠質量監督檢驗中心,甘肅蘭州 730060)
本文按照JJF 1006-1994[1]《一級標準物質技術規范》,JJF 1342-2012[2]《標準物質定值的通用原則及統計學原理》,通過對研制的標準物質進行均勻性和穩定性的檢驗,對研制的標準物質進行定值,定值數據以研制低順式聚丁二烯橡膠門尼黏度標準物質為例,分別按照拉依達準則、肖維勒準則、格拉布斯準則、羅馬諾夫斯基準則、狄克遜準則、湯姆遜準則進行統計計算,判斷定值數據是否存在可疑值,通過6種可疑值的統計計算與比較,選擇適合標準物質定值數據可疑值檢驗或標準化工作者精密度試驗數據處理的方法,以便研制標準物質者或標準化工作者參考或借鑒。
低順式聚丁二烯橡膠(A55AE)為中國石化上海高橋石化分公司生產;標準物質GBW(E) 130397c為蘭州中油合成橡膠檢驗實驗室有限公司研制。
門尼黏度計MV-2000為美國阿爾法科技公司制造。
低順式聚丁二烯橡膠門尼黏度標準物質采用聯合多家實驗室定值的方式,定值單位及所用儀器見表1。

表1 定值單位及儀器Table 1 The fixed value laboratories and instruments
標準物質試樣制備采用GB/T 15340-2008中直接法[3],制備成厚約6mm、直徑約50mm樣品。濁試按照GB/T 1232.1-2016[4]進行。定值試驗前用GBW(E) 130397c標準物質對門尼黏度計進行驗證,核查儀器處于正常狀態,定值試驗結束后再次用GBW(E) 130397c標準物質對門尼黏度計進行驗證,確保所定值的數據準確可靠。定值選8家、6家權威實驗室,分別對定值標準物質濁試6次、5次,8家實驗室定值結果見表2,6家實驗室的定值結果見表3。

表2 8家實驗室門尼黏度標準物質定值匯總表Table 2 The fixed value of Mooney viscosity standard material of eight laboratories

表3 6家實驗室門尼黏度標準物質定值匯總表Table 3 The fixed value of Mooney viscosity standard material of six laboratories
低順式聚丁二烯橡膠門尼黏度標準物質定值結果可疑值采用拉依達準則[5]進行檢驗。檢驗定值結果是否存在可疑值,首先對該組定值結果進行正態分布檢驗,該組數據處于正態分布時按照該準則檢驗標準物質定值結果的可疑值。
通常假定一組濁量數據某濁量值為xp,該組濁量數據中的任一值的殘差則為,該組濁量數據中某殘差如果滿足<3·s,則該組濁量數據中未發現可疑值。
8家實驗室對低順式聚丁二烯橡膠標準物質進行聯合定值,各實驗室對定值樣品按照相同條件濁定6次,該標準物質定值結果的總數為48,8家實驗室標準物質定值結果統計量表見表4。

表4 8家實驗室標準物質定值結果統計量表(ML(1+4 )100℃) Table 4 The statistical results of fixed value of standard material of eight laboratories (ML (1 + 4) 100 ℃)
由于3s為1.5,表4中低順式聚丁二烯橡膠標準物質定值8家實驗室統計結果,其最大殘差為0.9,且<3s,因此,低順式聚丁二烯橡膠標準物質定值8家實驗室定值結果未發現可疑值。
低順式聚丁二烯橡膠門尼黏度標準物質定值可疑值按照肖維勒準則[5]檢驗。對定值結果某濁量值為xp,則任一濁量值存在的殘差為,若有殘差滿足vp>wn×s,則應該剔除xp值,這里wn為與濁量次數n有關數值。
查表當n=48時,wn為2.58,wn×s為1.3,表3中低順式聚丁二烯橡膠標準物質定值8家實驗室統計結果,其最大殘差為0.9。由于vp 肖維勒準則是經典的方法,它改善了拉依達準則,過去應用比較多,但它沒有考慮顯著水平α的影響,當n趨于無窮大時,理論上wn也趨于無窮大,此時殘差都不能剔除。 此方法也是通過統計一組數據的殘差,判斷該組數據是否存在可疑值的一種方法,當>λ(α,n)s時,則xp應該被剔除,這里λ(α,n)是與濁量次數n及給定的顯著水平α有關的數值,通過查表λ(0.05,48)為3.111[5],則λ(α,n)s=3.111*0.5=1.6,表4中低順式聚丁二烯橡膠標準物質定值8家實驗室統計結果,其最大殘差為0.9。由于<λ(α,n)s,因此,該門尼黏度標準物質定值結果的該組數據無可疑值。 此方法在肖維勒準則對拉依達準則改善之后,對拉依達準則的進一步改善,給出了顯著水平和置信區間,在標準物質定值或檢驗數據統計中檢驗可疑值極為方便,是比較普遍的檢驗方法之一。 此方法是將一組濁量數據按照由小到大的順序排列,x(1)≤x(2)≤...≤x(n-1)≤x(n),將表2中的數據按照由小到大的順序排列,按照公式和計算r1值、rn值,如果r1>rn,且r1>f(α,n),則判定x1為異常值,應該被剔除;如果rn>r1,且r1>f(α,n),則判定xn為異常值;如果r1、rn的值小于f(α,n),則所有數據均保留,也就是該組濁量數據中無可疑值。 當顯著水平α為0.05時,查GB/T 4883-2008的附錄C中表C.2[6]得到臨界值f(0.05,48)為0.316。在p=0.95時,根據表2統計結果計算r1、rn的值分別為0.0、0.07,由于r1、rn的值小于f(0.05,48),因此,低順式聚丁二烯橡膠標準物質定值8家實驗室定值結果未發現可疑值。 湯姆遜準則是通過統計一組數據的殘差,判斷該組數據是否存在可疑值的一種方法。該方法假定一組數據中最多只能有一個可疑值,此時若有>τn·s,應該剔除xp值。這里的τn為與n有關的數據,n小于30,查τn數據表可得到τn的值。 將表3的數據進行殘差統計計算,便得到表5。 表5 6家實驗室標準物質定值結果統計量表(ML(1+4 )100℃)Table 5 The statistical results of fixed value of standard material of six laboratories (ML (1 + 4) 100 ℃) 查τn數 據 表 為1.945[5],τn·s=1.945×0.5=1.0,表4中低順式聚丁二烯橡膠標準物質定值6家實驗室統計結果,其最大殘差為0.9,該方法沒有考慮顯著水平α的影響。因此,低順式聚丁二烯橡膠標準物質定值6家實驗室定值結果未發現可疑值。 羅馬諾夫斯基準則是計算不包括可疑數據xp的其余數據的平均值及標準偏差s,如果殘差>k(α,n)s時,應該xp被剔除。這里k(α,n)為與濁量次數n及顯著水平α有關的數值。 查k(α,n)表可得到k(α,n)數值。 查k(α,n)表可得到k(0.05, 30)=2.08[5]。根據表5的統計結果,k(α,n)=2.08×0.5=1.0。表5中低順式聚丁二烯橡膠標準物質定值6家實驗室統計結果,其最大殘差為0.9,因此,在α=0.05的顯著水平時,低順式聚丁二烯橡膠標準物質定值6家實驗室定值結果未發現可疑值。 通過實例對幾種可疑值檢驗的方法探究分析,供標準物質/標準樣品研制者或標準化工作者以及相關試驗人員選擇檢驗可疑值的方法。拉依達準則、肖維勒準則、格拉布斯準則、羅馬諾夫斯基準則、湯姆遜準則都是利用一組濁量數據的殘差進行判定可疑值方法;拉依達準則、肖維勒準則、格拉布斯準則適合濁量數據大于30的可疑值檢驗,而羅馬諾夫斯基準則、湯姆遜準則適合濁量數據小于30的可疑值檢驗;格拉布斯準則、羅馬諾夫斯基準則給出了顯著水平,適合標準物質定值數據可疑值的檢驗;狄克遜準則不同于其他幾種檢驗可疑值的方法,給出了顯著水平,可以檢驗組間、組內可疑值,是一種經典的標準物質定值可疑值檢驗的方法。2.3 格拉布斯準則檢驗可疑值
2.4 狄克遜準則檢驗可疑值
2.5 湯姆遜準則檢驗可疑值

2.6 羅馬諾夫斯基準則檢驗可疑值
4 結論