Z比分?jǐn)?shù)評(píng)定法在大腸菌群能力驗(yàn)證評(píng)定中的應(yīng)用

◎ 安鐵鵬，張 莉，屠錦娣，周爭艷，候國峰，陳 江，商 帥

（陜西科儀陽光檢測技術(shù)服務(wù)有限公司，陜西 西安 710016）

能力驗(yàn)證是認(rèn)證管理機(jī)構(gòu)評(píng)價(jià)和監(jiān)督實(shí)驗(yàn)室技術(shù)能力的重要依據(jù)之一，是判斷和監(jiān)控實(shí)驗(yàn)室技術(shù)能力的持續(xù)狀況、識(shí)別實(shí)驗(yàn)室間差異的一種有效手段。實(shí)驗(yàn)室間比對(duì)的方法對(duì)數(shù)據(jù)分析提出了挑戰(zhàn)，雖大多數(shù)實(shí)驗(yàn)室間比對(duì)可提供單峰和近似對(duì)稱的數(shù)據(jù)，但大多數(shù)能力驗(yàn)證數(shù)據(jù)集包含一定比例的偏離數(shù)據(jù)主體的結(jié)果，這種離群結(jié)果可能具有高變異性，從而使得均值和標(biāo)準(zhǔn)差值的傳統(tǒng)統(tǒng)計(jì)方法變得不可靠[1]。將本次提供的4組樣品分別分發(fā)到123家實(shí)驗(yàn)室進(jìn)行定量檢測，將結(jié)果匯總后進(jìn)行對(duì)數(shù)轉(zhuǎn)換，分析采用GB/T 28043—2019中的穩(wěn)健技術(shù)（參見附錄C）。避免了極端值（即離群結(jié)果）對(duì)統(tǒng)計(jì)結(jié)果的影響，從而可以比較客觀、科學(xué)地評(píng)價(jià)每一個(gè)實(shí)驗(yàn)室的結(jié)果。本次選取一組共39個(gè)樣本數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

1 材料與方法

1.1 儀器與試劑

高壓蒸汽滅菌鍋、生化培養(yǎng)箱、顯微鏡；大腸桿菌標(biāo)準(zhǔn)菌株ATCC25922、乳糖蛋白胨培養(yǎng)液、伊紅美藍(lán)培養(yǎng)基。

1.2 試驗(yàn)方法

按照《生活飲用水標(biāo)準(zhǔn)檢驗(yàn)方法 微生物指標(biāo)》（GB/T 5750.12—2006）中總大腸菌群（多管發(fā)酵法）進(jìn)行樣品前處理。本次能力驗(yàn)證為定量分析（MPN法），根據(jù)總大腸菌群陽性的管數(shù)，查MPN檢索表3，報(bào)告每100 mL水樣中的總大腸菌群最可能數(shù)（Most Probable Number，MPN）值[2]。

1.3 統(tǒng)計(jì)方法及評(píng)價(jià)

對(duì)本次能力驗(yàn)證計(jì)劃檢測項(xiàng)目實(shí)驗(yàn)室的檢測結(jié)果經(jīng)對(duì)數(shù)轉(zhuǎn)換，統(tǒng)計(jì)分析采用GB/T 28043—2019中的穩(wěn)健技術(shù)（參見附錄C），先計(jì)算中位數(shù)med（x），再根據(jù)中位數(shù)計(jì)算中位絕對(duì)離差MADe（x），采用迭代算法A計(jì)算數(shù)據(jù)的穩(wěn)健平均值c和穩(wěn)健標(biāo)準(zhǔn)差b，直到穩(wěn)健均值和標(biāo)準(zhǔn)差的第3個(gè)有效數(shù)字連續(xù)兩次不再發(fā)生變化[1]。根據(jù)公式δ=1.5b求得能力評(píng)定標(biāo)準(zhǔn)差δ值，最后根據(jù)以下公式計(jì)算Z值：

式中：x-實(shí)驗(yàn)室測試結(jié)果；c-穩(wěn)健平均值；δ-能力評(píng)定標(biāo)準(zhǔn)差。

應(yīng)用Z值評(píng)定準(zhǔn)則如下：|Z|≤2，表明“滿意”，無需采取進(jìn)一步措施；2＜|Z|＜3，表明“有問題”（可疑值）；|Z|≥3，表明“不滿意”（離群值）[3]。

2 結(jié)果統(tǒng)計(jì)與分析

2.1 實(shí)驗(yàn)室測定結(jié)果及對(duì)數(shù)值

對(duì)匯總的檢測結(jié)果按照數(shù)據(jù)大小進(jìn)行排序，并求log（10），表示為x（1），x（2）…x（p）；具體見表 1。

表1 大腸菌群統(tǒng)計(jì)結(jié)果對(duì)數(shù)值表

2.2 中位數(shù)med（x）的計(jì)算

將測定數(shù)據(jù)按遞增順序排列[4]，表示為x（1），x（2）…x（p）。

2.2.1 中位絕對(duì)離差MADe（x）的計(jì)算

中位絕對(duì)離差MADe（x）用于估計(jì)正態(tài)分布數(shù)據(jù)的總體標(biāo)準(zhǔn)差，能夠極大抵消離群值得影響，先計(jì)算絕對(duì)偏差di（i=1,…,p）：

再根據(jù)公式MADe（x）=1.483med（d）計(jì)算中位絕對(duì)離差。

2.2.2 穩(wěn)健平均值c和穩(wěn)健標(biāo)準(zhǔn)差b的計(jì)算

穩(wěn)健平均值c和穩(wěn)健標(biāo)準(zhǔn)差b的初始值計(jì)算：

根據(jù)以下步驟更新穩(wěn)健平均值c和穩(wěn)健標(biāo)準(zhǔn)差b：

對(duì)每個(gè)x（i）（i=1,2,…p），計(jì)算：

再由下式計(jì)算穩(wěn)健平均值c和穩(wěn)健標(biāo)準(zhǔn)差b

穩(wěn)健估計(jì)值c和b可由迭代計(jì)算得出，即根據(jù)（6）～（9）公式多次更新c和b的值，直到該過程收斂，即穩(wěn)健均值和標(biāo)準(zhǔn)差的第3個(gè)有效數(shù)字連續(xù)兩次不再發(fā)生變化。

2.2.3Z值計(jì)算

式中：x-實(shí)驗(yàn)室測試結(jié)果對(duì)數(shù)值；c-穩(wěn)健平均值；δ-能力評(píng)定標(biāo)準(zhǔn)差。

凍干粉中大腸菌群穩(wěn)定均值、標(biāo)準(zhǔn)差和Z 值計(jì)算詳見表2。

表2 凍干粉中大腸菌群穩(wěn)定均值、標(biāo)準(zhǔn)差和Z值計(jì)算表

（續(xù)表2）

2.3 直方圖顯示

為了更直觀的展現(xiàn)能力驗(yàn)證統(tǒng)計(jì)結(jié)果，可采用直方圖表示有關(guān)統(tǒng)計(jì)量。Z比分?jǐn)?shù)序列直方圖按照Z比分?jǐn)?shù)大小順序顯示出每個(gè)實(shí)驗(yàn)室的Z比分?jǐn)?shù)值，并標(biāo)有實(shí)驗(yàn)室的編號(hào)。根據(jù)代表實(shí)驗(yàn)室的柱狀圖可以清晰辨別出實(shí)驗(yàn)室結(jié)果及所在位置，見圖1。

圖1 Z比分?jǐn)?shù)直方圖

本組39家實(shí)驗(yàn)室，最終通過數(shù)據(jù)統(tǒng)計(jì)有36家實(shí)驗(yàn)室結(jié)果滿意，2家結(jié)果可疑，1家結(jié)果不滿意。

3 結(jié)論

大腸菌群是一類在一定培養(yǎng)條件下能發(fā)酵乳糖、產(chǎn)酸產(chǎn)氣的需氧和兼性厭氧菌的革蘭氏陰性無芽胞桿菌[5]。本次大腸菌群檢測能力驗(yàn)證共有123家實(shí)驗(yàn)室參加，結(jié)果不滿意的實(shí)驗(yàn)室21個(gè)，占17.1%；結(jié)果可疑的實(shí)驗(yàn)室8個(gè)，占6.5%。

對(duì)于Z比分?jǐn)?shù)統(tǒng)計(jì)數(shù)據(jù)分析，采用GB/T 28043—2019中的穩(wěn)健技術(shù)（參見附錄C）迭代算法，避免了極端值（即離群結(jié)果）對(duì)統(tǒng)計(jì)結(jié)果的影響，從而可以比較客觀、科學(xué)地評(píng)價(jià)每一個(gè)實(shí)驗(yàn)室的結(jié)果。