基于聯系數的食品安全衛生檢測結果排序和影響因素分析

崔銳利 ，孫齊蕊，孫愛峰

1.白城市食品藥品監督管理局食品化妝品稽查分局，吉林白城 137000；2.吉林大學公共衛生學院，吉林長春 130021；3.白城市衛生急救中心，吉林白城 137000

聯系數是集對分析理論中的一個結構函數[1]，是反 映系統中確定-不確定性分析的一個重要工具。作為食品衛生監督重要組成部分的食品衛生檢測結果，可以用定性二分類（合格-不合格）方法表示，其本身就充滿了確定-不確定因素，因為隨著人們對檢測食品中有害物質的認識不同，合格標準在不同時期可能發生變化，因此可以應用聯系數進行進一步分析。為此，作者將某市2010—2012年及其各季度和各檢測指標的食品安全衛生檢測結果運用聯系數進行分析，在對不同年度不同季度食品安全衛生檢測結果進行排序的同時，找出影響該市食品安全衛生監測結果的主要因素，為順利開展食品衛生監督工作奠定堅實的基礎。現將結果報告如下。

1 資料與方法

1.1 一般資料

以文獻[2]提供的某市2010—2012年食品安全衛生檢測超標情況為基本資料，計算得到不同年度不同季度及其不同指標的合計合格率（%）和不合格率（%），具體數值見表1。由于數據來自文獻，真實可靠。

表1 不同年度不同季度和不同檢測指標檢測結果（%）

1.2 統計方法

對于表1資料，應采用同-反型二元聯系數表示，并計算各聯系數的有效值[3-5]，判定各年（季）度食品安全衛生檢測結果的排序結論，判定各評價指標對食品安全衛生檢測結果的影響程度。

第一步，依據表1給出數據中的合格率作為二元聯系數的同部，不合格率作為反部，寫出各年（季）度各評價指標的歸一化聯系數umn=amn±cmnj，式中m表示年（季）度，n表示評價指標。

第二步，計算非賦權年（季）度歸一化聯系數um。由于同一年（季）度的抽檢樣本數量相同，所以聯系數計算公式為 Um=（am1±cm1j）+（am2±m2j）+（am3±cm3j）， 經過歸一化處理得到非賦權年（季）度聯系數um。

第三步，計算非賦權評價指標歸一化聯系數un。將表1數據的最后一列各指標的合計合格率作為同-反型二元聯系數的同部，不合格率作為反部，得到非賦權評價指標歸一化聯系數un=an±cnj。

第四步，判定結果。令j=-1計算各聯系數的有效值[3-4]，根據數值越大排序越靠前的準則判定各年（季）度食品安全衛生檢測結果的排序，根據數值越大影響程度越低的準則判定各評價指標對食品安全衛生檢測結果的影響程度。

2 結果

2.1 各年（季）度各評價指標歸一化聯系數的表示

依據方法中的第一步，各年（季）度各評價指標歸一化聯系數見表2。

表2 各年（季）度各評價指標歸一化聯系數

2.2 不同年度（季度）食品安全衛生檢測結果的排序

按照方法中的第二步和第四步，計算不同年（季）度非賦權歸一化聯系數及其有效值，按照大小順序分別對年度和季度排序，見表3。

表3 不同年（季）度非賦權歸一化聯系數及其有效值

2.3 不同評價指標影響程度的比較

計算表2中最后一列給出的不同評價指標非賦權歸一化聯系數的有效值，按照從小到大的順序判定各評價指標對食品安全情況的影響程度，見表4。

表4 不同評價指標對食品安全情況影響程度的比較

3 討論

聯系數是集對分析的一個特征函數和重要的數學分析工具，可用來全面和深入地研究兩個集合的確定性和不確定性關系的聯系、可變和轉化，在衛生監督中已經得到一定應用[3-4,6-9]。但上述文獻均從不同角度運用三元聯系數或多元聯系數進行不同的統計分析，獲得較為理想的結果，但尚未見有用同反型聯系數對食品安全衛生檢測結果進行排序和影響因素分析的文獻，因此作者運用集對分析的基本理論，建立歸一化同反型聯系數并計算其有效值，做以上分析，以指導食品衛生監督工作，保證食品安全，維護和保障人群健康水平。

從不同年度（季度）食品安全衛生檢測結果的分析看，2010—2012年的歸一化聯系數分別為u2010=99.31±0.69j、u2011=97.82±2.18j和 u2012=98.90±+1.10j， 其有效聯系值分別為98.62、95.64和97.80，依據判定準則可知不同年度食品安全衛生檢測結果的排序為2010年＞2012年＞2011年；第一季度至第四季度的歸一化聯系數分別為 u1=98.48±1.52j、u2=98.07±1.93j、u3=99.56±0.44j和 u4=98.64±1.36j，其有效聯系值分別為 96.96、96.14、99.12 和97.28，依據判定準則可知不同季度食品安全衛生檢測結果的排序為第三季度＞第四季度＞第一季度＞第二季度，與秩和比法的分析結論完全相同，但對不同年度的排序與文獻[2]的結論不同。

從不同檢測指標影響該市食品安全情況的評價看，微生物的歸一化聯系數為98.45±1.55j，重金屬為99.14±0.86j，添加劑為 98.33±1.67j，其有效聯系值分別為96.90、98.28和96.66，依據判定準則可知各指標的影響程度排序為添加劑＞微生物＞重金屬，與秩和比法[10]的分析結論完全相同，所以應重視添加劑和微生物對食品安全情況的影響，并重點加強食品衛生監測。

在計算各評價指標歸一化聯系數時，其實質是計算合計（平均）檢測合格率和不合格率。這時不能簡單地由不同年度或不同季度計算得到的率相加后求平均。因為不同年度或不同季度的檢測樣本數量是不同的，所以必須按照率的定義，計算不同年度或不同季度總的合格樣本數量和不合格樣本數量，之后去除3個年度采樣的總數量得到合計（平均）合格率和不合格率，據此來建立二元聯系數。而計算非賦權年（季）度歸一化聯系數時，由于同一年（季）度的抽檢樣本數量相同，所以可以按照將幾個率相加后，取其平均值作為幾個率的平均 （合計）率。這是應用此方法進行排序和影響因素分析應該注意的問題，以免出現錯誤結論，誤導衛生監督實踐。

綜上所述，該文所得結論與秩和比法完全相同，唯一缺陷是不能進行多重比較，但原理更為明了，計算更為簡單，操作性較強，值得在食品衛生監測結果排序和影響因素分析中得到推廣應用。

[1]蔣云良，趙克勤，劉以安，等.信息處理集對分析[M].北京：清華大學出版社，2015.

[2]范宇暉，孟彩霞，王艷菊，等.2010年-2012年深圳市食品衛生安全檢測結果分析[J].臨床醫藥實踐，2014，23（4）：409-414.

[3]汪麗偉，孫齊蕊，孫愛峰.聯系數在公共場所衛生監督質量綜合評價中的應用[J].中國衛生產業，2014，11（21）：62-63.

[4]史景明，孫愛峰.聯系數在食品衛生質量合格率排序和潛在發展趨勢分析中的應用 [J].中國醫藥指南，2013，11（33）：590-591.

[5]張紹林，孫愛峰.聯系數有效值在抗菌藥物配伍治療布魯菌病效果排序中的應用[J].中國醫藥指南，2012，10（31）：373-374.

[6]吳瑩，孫愛峰，孫齊蕊.集對分析在食品衛生監督工作質量分類和排序中的應用[J].中國衛生產業,2014,11(21):61,63.

[7]邱昭君，孫愛峰.聯系數在不同地區衛生監督質量聚類分析中的應用[J].中國醫藥指南,2011,9(36):490-492.

[8]夏淼，孫愛峰.集對分析在職業衛生監督質量綜合評價中的應用價值[J].中國衛生產業，2015，12(3):10-11.

[9]張杰，孫愛峰，李雪云.聯系數在不同年度不同種類食品檢測總合格率聚類分析中的應用[J].中國醫藥指南，2012，10(19):392-393.