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

崔銳利 ，孫齊蕊，孫愛峰

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

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

1 資料與方法

1.1 一般資料

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

表1 不同年度不同季度和不同檢測指標(biāo)檢測結(jié)果（%）

1.2 統(tǒng)計(jì)方法

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

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

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

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

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

2 結(jié)果

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

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

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

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

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

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

2.3 不同評價(jià)指標(biāo)影響程度的比較

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

表4 不同評價(jià)指標(biāo)對食品安全情況影響程度的比較

3 討論

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

從不同年度（季度）食品安全衛(wèi)生檢測結(jié)果的分析看，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ù)判定準(zhǔn)則可知不同年度食品安全衛(wèi)生檢測結(jié)果的排序?yàn)?010年＞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ù)判定準(zhǔn)則可知不同季度食品安全衛(wèi)生檢測結(jié)果的排序?yàn)榈谌径龋镜谒募径龋镜谝患径龋镜诙径?，與秩和比法的分析結(jié)論完全相同，但對不同年度的排序與文獻(xiàn)[2]的結(jié)論不同。

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

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

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

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

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

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

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

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

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

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

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

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