基于聯(lián)系數(shù)的食品安全衛(wèi)生檢測結果排序和影響因素分析

崔銳利 ，孫齊蕊，孫愛峰

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

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

1 資料與方法

1.1 一般資料

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

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

1.2 統(tǒng)計方法

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

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

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

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

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

2 結果

2.1 各年（季）度各評價指標歸一化聯(lián)系數(shù)的表示

依據(jù)方法中的第一步，各年（季）度各評價指標歸一化聯(lián)系數(shù)見表2。

表2 各年（季）度各評價指標歸一化聯(lián)系數(shù)

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

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

表3 不同年（季）度非賦權歸一化聯(lián)系數(shù)及其有效值

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

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

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

3 討論

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

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

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

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

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

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

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

[3]汪麗偉，孫齊蕊，孫愛峰.聯(lián)系數(shù)在公共場所衛(wèi)生監(jiān)督質量綜合評價中的應用[J].中國衛(wèi)生產業(yè)，2014，11（21）：62-63.

[4]史景明，孫愛峰.聯(lián)系數(shù)在食品衛(wèi)生質量合格率排序和潛在發(fā)展趨勢分析中的應用 [J].中國醫(yī)藥指南，2013，11（33）：590-591.

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

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

[7]邱昭君，孫愛峰.聯(lián)系數(shù)在不同地區(qū)衛(wèi)生監(jiān)督質量聚類分析中的應用[J].中國醫(yī)藥指南,2011,9(36):490-492.

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

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