食品安全的定量評(píng)估

1、食品安全問(wèn)題的定量評(píng)估摘要本題是一個(gè)食品安全定量評(píng)估的問(wèn)題。我們組根據(jù)網(wǎng)上公布的數(shù)據(jù)和統(tǒng)計(jì)年鑒，先選出幾個(gè)比較突出的化學(xué)污染物（如鎘、汞、砷、有機(jī)氯、農(nóng)藥DDT），找出中國(guó)人均食物日攝入量及種類(lèi)，之后選取深圳的網(wǎng)上公布的數(shù)據(jù)分析研究深圳的食品安全問(wèn)題。我們組提出兩個(gè)問(wèn)題：一是如何確定污染物各自之間的比重及人均日攝入總量，二是化學(xué)污染物的攝入量對(duì)人體的危害程度的評(píng)價(jià)。針對(duì)問(wèn)題一建立模型一：運(yùn)用層次分析，結(jié)合CAC的食物分類(lèi)標(biāo)準(zhǔn)將化學(xué)污染物分成3類(lèi)，運(yùn)用熵權(quán)法求權(quán)重求出各個(gè)污染物各自之間的比重及在日食入量中占的比重；模型二是運(yùn)用概率評(píng)估模型進(jìn)行分析的，我們知道對(duì)總?cè)巳簛?lái)說(shuō)，其每日攝入食物量和自身體

2、重都是服從正態(tài)分布的。我們用MATLAB軟件中normrnd函數(shù)分別模擬出一組服從正態(tài)分布的日攝入和體重的隨機(jī)數(shù)，由某種化學(xué)污染物攝入量和安全攝入量得到某種化學(xué)污染物對(duì)人體安全危害指數(shù)；然后調(diào)用深圳市2001年全市食品安全檢查數(shù)據(jù)，并對(duì)其進(jìn)行檢驗(yàn)，得出當(dāng)?shù)厥称肺：χ笖?shù)表如下圖，并對(duì)其做出合理的分析。 化學(xué)污染物名稱濃度危害指數(shù)中毒程度鎘0.010241.2安全無(wú)毒汞0.008128.9正常砷0.03556.2安全無(wú)毒有機(jī)氯0.0196.5安全無(wú)毒農(nóng)藥DDT0.0211.4正常最后我們對(duì)模型優(yōu)缺點(diǎn)進(jìn)行分析，并對(duì)模型中存在的問(wèn)題進(jìn)一步改進(jìn)。關(guān)鍵詞：CAC 層次分析，熵權(quán)法，概率評(píng)估，危害指數(shù) 危害

3、程度一、問(wèn)題重述俗話說(shuō)的好“安居樂(lè)業(yè)”，其實(shí)社會(huì)也是如此，如果必需品都沒(méi)有安全保證，那社會(huì)又怎么會(huì)和諧！近年來(lái)，地溝油、硫磺姜、回爐面包、三鹿牛奶等食品安全事件頻發(fā)，食品安全問(wèn)題依然是當(dāng)下中國(guó)最熱門(mén)的一幕丑劇。當(dāng)瘦肉精事件尚未平息，上海染色饅頭又像一重磅炸彈引爆全國(guó)；隨后接棒的則有溫州毒饅頭、廣東墨汁粉條、遼寧毒豆芽等。常年風(fēng)起云涌的食品安全事件，像一雙雙粗重的鋼琴手，一次次肆虐地撩撥人們敏感的神經(jīng)。為此，人大常委會(huì)建議：把食品安全作為國(guó)家安全組成部分，其重要性不亞于金融安全、糧食安全、能源安全、生態(tài)安全；加大對(duì)違反食品安全行為的懲處力度，以法治的方式維護(hù)公民的“食品安全權(quán)”；推進(jìn)食品，尤其是

4、農(nóng)產(chǎn)品的規(guī)模化、產(chǎn)業(yè)化生產(chǎn)經(jīng)驗(yàn)，從根本上保證食品質(zhì)量；參照科學(xué)依據(jù)，結(jié)合中國(guó)國(guó)情，制定出全國(guó)統(tǒng)一的食品安全檢驗(yàn)和認(rèn)定標(biāo)準(zhǔn)；充分發(fā)揮社會(huì)組織作用，強(qiáng)化食品生產(chǎn)行業(yè)協(xié)會(huì)的作用，加強(qiáng)行業(yè)自律；加強(qiáng)宣傳力度，將涉及到食品安全方面的法律法規(guī)匯編成冊(cè)，針對(duì)從業(yè)者流動(dòng)性強(qiáng)等特點(diǎn)，加大宣傳力度并從嚴(yán)考核；充分發(fā)揮社會(huì)監(jiān)督的作用，鼓勵(lì)全社會(huì)共同參與食品安全監(jiān)督等。食品安全問(wèn)題內(nèi)容廣泛、復(fù)雜，為了突出主要削弱次要問(wèn)題，我們小組選擇食品中污染物的分布作為研究食品安全的一個(gè)突破口，建立數(shù)學(xué)模型，利用互聯(lián)網(wǎng)數(shù)據(jù)，定量評(píng)估食品安全問(wèn)題。二、問(wèn)題分析建立數(shù)學(xué)模型即是抽象實(shí)際問(wèn)題，由于關(guān)于食品安全問(wèn)題的方面又太過(guò)繁雜，我們組

5、僅選擇其中一個(gè)方面即食品中污染物的分布進(jìn)行研究，為了定量評(píng)估題目中涉及到的食品安全問(wèn)題。1）問(wèn)題的背景“民以食為天，食以安為先”，食品衛(wèi)生安全關(guān)系國(guó)計(jì)民生，關(guān)系經(jīng)濟(jì)繁榮和社會(huì)穩(wěn)定。近年來(lái)，隨著人民生活水平的不斷提高，食品安全問(wèn)題更加受到人民的關(guān)注。我國(guó)是人口大國(guó)，食品消費(fèi)量大，而糧食生產(chǎn)加工分散，因而食品衛(wèi)生安全監(jiān)控難度較大。學(xué)習(xí)借鑒歐美國(guó)家食品衛(wèi)生安全監(jiān)控的做法，建立化學(xué)污染在食品中的安全模型，用以處理調(diào)查或檢測(cè)數(shù)據(jù)，對(duì)公共食品衛(wèi)生安全做出評(píng)估應(yīng)該是一條可行的途徑。歐美已經(jīng)建立了有關(guān)模型并且開(kāi)始使用，用時(shí)間證明了其可行性。但是模型的建立和使用還存在著很多困難和問(wèn)題。在建立食品安全模型時(shí)需要設(shè)

6、計(jì)合理的周陽(yáng)調(diào)查方案或者數(shù)據(jù)的搜集，并由抽樣數(shù)據(jù)建立比較準(zhǔn)確的人群食品攝入量模型?；瘜W(xué)污染物分布模型的建立面臨的問(wèn)題包括：怎樣利用抽樣率很低的隨機(jī)抽樣數(shù)據(jù)建立模型；怎么由隨機(jī)變量取值大于某一數(shù)值的部分樣本數(shù)據(jù)再加上其他可以利用的信息估計(jì)出這個(gè)隨機(jī)變量的整體分布。危害評(píng)估模型就是利用模型一的結(jié)果對(duì)全國(guó)、某地區(qū)、某類(lèi)食品的安全狀況作出評(píng)價(jià)。2）問(wèn)題分析通過(guò)查閱相關(guān)資料，我們了解到食品中對(duì)人體健康造成威脅的因素大致有三類(lèi)：物理因素，化學(xué)因素，生物因素。其中，物理因素可以通過(guò)良好的生產(chǎn)操作加以避免，因此基本不做討論；生物因素涉及到外部條件，以及該生物危害人體的繁殖和變異過(guò)程，具有較高的復(fù)雜性，也不易做

7、研究；所以，我們主要分析化學(xué)因素?；瘜W(xué)因素主要包括重金屬、有機(jī)農(nóng)藥、化肥、食品添加劑等，他們會(huì)通過(guò)各種復(fù)雜的途徑先污染食物，進(jìn)而危害人體健康，根據(jù)所污染所產(chǎn)生的危害的嚴(yán)重性程度劃分等級(jí)，確定食品中化學(xué)污染物的危害指數(shù)，進(jìn)而完成對(duì)食品安全的定量評(píng)估。為此，我們需要解決一下兩個(gè)方面的問(wèn)題：?jiǎn)栴}一：查閱相關(guān)資料，根據(jù)所得數(shù)據(jù)求出各種食品中不同污染物的污染水平，及人均日攝入量中各化學(xué)污染物的含量。問(wèn)題二：計(jì)算出各類(lèi)化學(xué)元素的人均攝入量,結(jié)合問(wèn)題一得到的分類(lèi)后各類(lèi)食物中所含的化學(xué)污染物的含量，怎樣計(jì)算出各類(lèi)化學(xué)污染物的人均攝入量和人體每日安全攝入量，如何根據(jù)這些數(shù)據(jù)確定出該化學(xué)污染物人均攝入量與危害指數(shù)

8、的關(guān)系；并最終確定某種化學(xué)污染物實(shí)際測(cè)定值與食品安全的關(guān)系。三、模型假設(shè)3.1假設(shè)各化學(xué)污染物之間不會(huì)相互作用；3.2假設(shè)被評(píng)估人群在調(diào)查期內(nèi)的體重是一個(gè)定值；3.3假設(shè)被評(píng)估人群的體重服從正態(tài)分布；3.4假設(shè)被評(píng)估人群每天攝入量是服從正態(tài)分布的；3.5假設(shè)化學(xué)殘留物在食物中的各項(xiàng)指標(biāo)不會(huì)隨時(shí)間改變；3.6被評(píng)估人群中不同個(gè)體對(duì)化學(xué)污染物的耐受劑量與體重成正比。四、符號(hào)說(shuō)明國(guó)際食品法典委員會(huì)食品與飼料分類(lèi)標(biāo)準(zhǔn)指標(biāo)組成的矩陣矩陣列向量歸一化后的矩陣人均每日消耗食物量第j個(gè)指標(biāo)的熵定義第j個(gè)指標(biāo)的熵權(quán)值第i種食物的攝入量第i種食物中鎘的濃度 第i種中食物汞的濃度第i種食物中砷的濃度第i種食物中有機(jī)

9、氯的濃度第i種食物中農(nóng)藥DDT的濃度個(gè)體的體重食品中某種化學(xué)污染物的人群攝入量食品安全危害指數(shù)安全系數(shù)五、模型的建立與求解5.1模型的建立與求解首先，建立對(duì)食品中某種化學(xué)污染物攝入量的評(píng)估模型，得到食品中該化學(xué)污染物的攝入量評(píng)估值；其次，建立食品安全危害指數(shù)模型，將實(shí)際食品中化學(xué)污染物濃度的測(cè)定量與通過(guò)模型得到的危害指數(shù)等級(jí)表中的各級(jí)濃度范圍進(jìn)行比較，最終可定量評(píng)估出食品的安全程度。5.2問(wèn)題一模型的建立我們參考國(guó)際食品法典委員會(huì)食品與飼料分類(lèi)標(biāo)準(zhǔn)對(duì)食物進(jìn)行分類(lèi)。食品法典分類(lèi)的標(biāo)準(zhǔn)是按4個(gè)層次來(lái)進(jìn)行食品和動(dòng)物飼料分類(lèi)的。首先，按照原料來(lái)源及是否經(jīng)過(guò)加工分為A、B、C、D、E5個(gè)等級(jí)。這5個(gè)等級(jí)

10、分別是A級(jí)植物來(lái)源的初級(jí)食品、B級(jí)動(dòng)物來(lái)源的初級(jí)食品、C級(jí)初級(jí)飼料、D級(jí)植物來(lái)源的加工食品、E級(jí)動(dòng)物來(lái)源的加工食品。各大類(lèi)型中根據(jù)污染物的相似性再分為不同的組，分類(lèi)標(biāo)準(zhǔn)的每個(gè)組中都包含各種推薦指標(biāo)，如潛在污染物殘留量、作物和動(dòng)物在正常情況下的消費(fèi)量以及應(yīng)用農(nóng)藥最大殘留限量的商品比例等。所以我們結(jié)合CAC分類(lèi)標(biāo)準(zhǔn)將部分食物進(jìn)行分類(lèi)，并分別研究各類(lèi)食品中的不同化學(xué)污染物。由于該污染物的種類(lèi)繁多，所以我們挑選對(duì)人體危害性較大的砷、汞、鎘，有機(jī)氯和農(nóng)藥DDT進(jìn)行研究。通過(guò)運(yùn)用層次分析法將各種食品分成三大類(lèi)，即植物性食物、動(dòng)物性食物、其他類(lèi)食物，得到如下層次結(jié)構(gòu)圖：食物中化學(xué)污染物的比重植物性動(dòng)物性其他

11、糧食蔬菜水果肉類(lèi)魚(yú)類(lèi)蛋奶類(lèi)奶粉零食飲料、茶葉奶類(lèi)果蔬罐頭水果魚(yú)肉罐頭魚(yú)類(lèi)皮蛋蛋類(lèi) 依據(jù)層次分析法確定指標(biāo)分層體系，然后采用熵權(quán)法計(jì)算食品中各化學(xué)污染物所占的比重。熵權(quán)法的基本原理及模型：熵是系統(tǒng)無(wú)序程度的度量，它表示從一組不確定事物中提取信息量的多少，可用于度量已知數(shù)據(jù)所包含的有效信息量和確定權(quán)重，在食品安全評(píng)價(jià)中可以采用。通過(guò)熵的計(jì)算確定權(quán)重，即根據(jù)各項(xiàng)監(jiān)測(cè)指標(biāo)值的差異程度u，確定各指標(biāo)的權(quán)重。1） 數(shù)據(jù)矩陣歸一化。設(shè)有m個(gè)評(píng)價(jià)指標(biāo)，n個(gè)評(píng)價(jià)對(duì)象，其原始數(shù)據(jù)矩陣為：式中（i=1，2,，m；j=1,2，n）為第i個(gè)指標(biāo)的第j個(gè)采樣值。矩陣歸一化：2） 在m個(gè)評(píng)價(jià)指標(biāo)n個(gè)評(píng)價(jià)對(duì)象的評(píng)價(jià)估計(jì)問(wèn)題中

12、，第j個(gè)指標(biāo)的熵的定義為： （i=1，2，m）當(dāng)=0時(shí)，令定義第j個(gè)指標(biāo)的熵后，則其熵權(quán)為： 其中 5.12問(wèn)題一模型的求解：層次分析法中對(duì)于第四層的分類(lèi)做一些說(shuō)明，由于對(duì)于不同的化學(xué)污染物，有的調(diào)查數(shù)據(jù)不全包含所有指標(biāo)，因此在計(jì)算比重時(shí)采用化學(xué)污染物含量與分類(lèi)總量的比值作為指標(biāo)的權(quán)重(平均值與樣本數(shù)的乘積再求和與總樣本數(shù)的比)。上述層次結(jié)構(gòu)圖中，我們以汞、鎘、砷、有機(jī)氯和DDT共有四層，其中第四層對(duì)第三層的權(quán)重為其加權(quán)平均值，根據(jù)附表中的一、二、三、四、五表中數(shù)據(jù)，整理化簡(jiǎn)，然后運(yùn)用熵權(quán)法，求出各層次每個(gè)成員的熵權(quán)值。是各種化學(xué)污染物在三類(lèi)食物中所占的比重組成的矩陣如下：第j個(gè)指標(biāo)的熵值矩陣

13、：五種化學(xué)污染物的權(quán)重矩陣：用熵權(quán)法求得的鎘、汞、砷、有機(jī)氯、農(nóng)藥在植物性食物、動(dòng)物性食物、其他食物中的比重如下表：化學(xué)污染物的比重表（1）化學(xué)污染物名稱植物性食物動(dòng)物性食物其他食物鎘的含量比0.4675710.2254550.306973汞的含量比0.6118480.3458610.042291砷的含量比0.1642390.1999660.635796有機(jī)氯的含量比0、164910.1253670.709724農(nóng)藥DDT的含量比0、0271040.1874820.785414再次運(yùn)用熵權(quán)法求鎘、汞、砷、有機(jī)氯、農(nóng)藥DDT的權(quán)重比，結(jié)果如下：食物中各污染物的比重表（2）名稱鎘汞砷有機(jī)氯農(nóng)藥DD

14、T食物中比值0.0336990.2235270.144570.2240260.374177從表（1）中我們可以看到化學(xué)污染物鎘在各類(lèi)食物中比較均衡；汞在植食性食物中含量比較大，而在其他食物中含量比很??；砷在植食性和動(dòng)食性含量比相當(dāng)，都不大，而在其他食物中含量比很大；有機(jī)氯和農(nóng)藥DDT在其他食物中含量比遠(yuǎn)高于植食性和動(dòng)食性食物。從表（2）中我們可以看出，就整體食用食物來(lái)說(shuō)，有機(jī)氯和農(nóng)藥DDT占比重較大，汞和砷所占比重一般，而鎘占最小。根據(jù)我們分析的結(jié)果，希望有關(guān)部門(mén)加強(qiáng)對(duì)這方面的注意力度 5.2問(wèn)題二模型的建立與求解5.21問(wèn)題二模型的建立：由于模型一僅僅是對(duì)食品中化學(xué)污染物分布情況作分析，未能

15、得出該類(lèi)化學(xué)污染物對(duì)人體的危害程度，因此我們用概率評(píng)估模型對(duì)此做一分析。由于食品中化學(xué)污染物的危害作用與其進(jìn)入人體的絕對(duì)量有關(guān)，因此評(píng)價(jià)食品安全以人體對(duì)化學(xué)污染物的實(shí)際攝入量與其安全攝入量比較更為科學(xué)合理，在這樣一種理論背景下，導(dǎo)出可以用來(lái)評(píng)價(jià)食品中某種化學(xué)污染物濃度對(duì)人體影響的食品安全危害指數(shù)：其中，是某種化學(xué)污染物的安全攝入量，我們經(jīng)查有關(guān)資料得出。由于以上模型主要針對(duì)于成人人群，考慮到不同人群的生理特性存在差異，引入一個(gè)安全系數(shù)，且可以由美國(guó)食品質(zhì)量安全法確定，進(jìn)而得到危害指數(shù)的最終模型，如下： 其中我們經(jīng)查資料得出安全系數(shù)FQPA=100，校正系數(shù)R=1。參照美國(guó)環(huán)保署(EPA)模式，

16、即攝入量的理論評(píng)估模型，對(duì)所攝入的食品中某種化學(xué)污染物的攝入量進(jìn)行理論評(píng)估，即：鎘的日攝入量：-(1)汞的日攝入量：-(2) 砷的日攝入量：-(3)有機(jī)氯的日攝入量：-(4)農(nóng)藥DDT：-(5)根據(jù)模型假設(shè)，我們分別將個(gè)體食品日攝入量和體重作為兩個(gè)獨(dú)立分布的總體和，且總體和均服從正態(tài)分布。在實(shí)際應(yīng)用中，通過(guò)大量的統(tǒng)計(jì)數(shù)據(jù)，分別得到服從正態(tài)分布的總體和的特征參數(shù)，即總體的方差和期望,以及總體的方差和期望。數(shù)據(jù)模擬模型:MATLAB的統(tǒng)計(jì)工具箱中，提供了通用的隨機(jī)數(shù)產(chǎn)生函數(shù)和特定分布的隨機(jī)數(shù)產(chǎn)生函數(shù)。我們用normrnd函數(shù)模擬出一組符合正態(tài)分布的隨機(jī)數(shù)。如下：個(gè)體食品日攝入量矩陣： -（6）和總

17、體的矩陣，如下： =-(7)將得到的矩陣(6)和(7)分別代到上式(1)，(2)，(3)，(4)，(5)即分別可得到每種化學(xué)污染物的日攝入量EXP的概率分布，從而可以計(jì)算一系列統(tǒng)計(jì)量，如：平均值，P50，P90，P99等作為被評(píng)估的化學(xué)污染物攝入估計(jì)值。我們用Matlab繪出化學(xué)污染物日攝入量概率密度曲線圖和人的體重概率密度曲線圖，由于有五種化學(xué)污染物要處理，我們?cè)俅蝺H貼出鎘的概率密度分布圖：食品中鎘的概率密度圖鎘的分布圖人體重的概率密度分布圖人體重的分布圖5.23化學(xué)污染物濃度與危害指數(shù)對(duì)應(yīng)關(guān)系的確定:針對(duì)五種具有代表性的化學(xué)污染物，再由模型中兩個(gè)子模型分別求出個(gè)體的化學(xué)污染物日攝入量矩陣，

18、對(duì)所得矩陣中的所有元素進(jìn)行升序排列，因?yàn)榈诎俜治粩?shù)可以保證的人是在相應(yīng)的安全裕度以內(nèi)，所以，我們以作為被評(píng)估人群的化學(xué)污染物日攝入量估計(jì)值。利用得到的五種化學(xué)污染物的日攝入量的估計(jì)值與相應(yīng)的標(biāo)準(zhǔn)值進(jìn)行比較，得出相應(yīng)的食品安全的危害指數(shù)。應(yīng)用結(jié)果如下表： 化學(xué)污染物名稱濃度危險(xiǎn)指數(shù)鎘0-0.00850-1010-5050-100>0.085>100汞0-0.00280-1010-5050-100>0.028>100砷0-0.0570-1010-5050-100>0.57>100有機(jī)氯0-0.0290-1010-5050-100>0.29>100農(nóng)藥

19、DDT0-0.0180-1010-5050-100>0.18>100根據(jù)有關(guān)網(wǎng)上資料得出危害指數(shù)與中毒程度對(duì)應(yīng)關(guān)系表，如下：危害指數(shù)與中毒程度對(duì)應(yīng)表危害指數(shù)0-1010-5050-100>100中毒程度安全無(wú)毒正常輕度中毒中毒根據(jù)網(wǎng)上數(shù)據(jù)調(diào)查有關(guān)資料得出某地的食品化學(xué)污染物含量得出鎘、汞、砷、有機(jī)氯、農(nóng)藥DDT危害指數(shù)，如下：化學(xué)污染物名稱濃度危害指數(shù)中毒程度鎘0.010241.2安全無(wú)毒汞0.008128.9正常砷0.03556.2安全無(wú)毒有機(jī)氯0.0196.5安全無(wú)毒農(nóng)藥DDT0.0211.4正常由該模型數(shù)據(jù)可得出該地區(qū)化學(xué)污染物都處于正常范圍。六、模型的評(píng)價(jià)與推廣模型一

20、中的熵是系統(tǒng)無(wú)序程度的度量，它表示從一組不確定事物中提取信息量的多少，可用于度量已知數(shù)據(jù)所包含的有效信息量和確定權(quán)重，在食品安全評(píng)價(jià)中可以采用。熵作為不確定行的度量，可直接根據(jù)候選方案各指標(biāo)值構(gòu)成的判斷矩陣來(lái)計(jì)算評(píng)價(jià)指標(biāo)的熵權(quán)。熵權(quán)法克服了常規(guī)方法評(píng)價(jià)過(guò)程中的主觀化和缺少科學(xué)合理計(jì)算方法的不足，很好地解決了評(píng)價(jià)過(guò)程中定性指標(biāo)難于客觀量化的問(wèn)題，從而使工程中的評(píng)價(jià)過(guò)程更加公正、科學(xué)。該模型思路清晰明確，可以通過(guò)計(jì)算機(jī)軟件來(lái)簡(jiǎn)化評(píng)標(biāo)決策中的有關(guān)計(jì)算與分析過(guò)程，提高評(píng)標(biāo)效率。由于根據(jù)網(wǎng)上公布的數(shù)據(jù)找出統(tǒng)計(jì)局關(guān)于人均日食物攝入的量及種類(lèi)，具有科學(xué)性；對(duì)指標(biāo)進(jìn)行檢測(cè)的食物種類(lèi)均是日常生活中的食物，具有一

21、般性；模型一僅考慮了化學(xué)污染物在食物中的影響，就有局限性，可是得出的數(shù)據(jù)還是具有一定的使用價(jià)值。模型二的應(yīng)用過(guò)程中，根據(jù)人體每日攝入量和人口體重正態(tài)分布的特征，用數(shù)學(xué)軟件工具模擬出所研究人群的每日攝入量和體重?cái)?shù)據(jù)。通過(guò)對(duì)不同模擬人模型結(jié)果的做出方差統(tǒng)計(jì)作圖，得到既滿足穩(wěn)定性要求，規(guī)模又最小的模擬人數(shù)300。.沒(méi)有大量的試驗(yàn)數(shù)據(jù)，不能得出相應(yīng)的各種化學(xué)污染物之間的相關(guān)系數(shù)。利用單一危害指數(shù)衡量化學(xué)污染物對(duì)食品安全造成的危害，沒(méi)有區(qū)分劇毒、重毒、低毒化學(xué)污染物生物學(xué)效應(yīng)的不同。只考慮了經(jīng)食物進(jìn)入體內(nèi)的化學(xué)污染物而沒(méi)有考慮其他途徑，如接觸等；進(jìn)而也沒(méi)有區(qū)分不同途徑進(jìn)入人體的化學(xué)污染物產(chǎn)生的差異。沒(méi)有

22、考慮食物中的化學(xué)污染物殘留值可能會(huì)隨時(shí)間而變化。不同化學(xué)污染物作用機(jī)理可能相差很遠(yuǎn)，作用器官也不一定相同，這些因素都會(huì)影響相對(duì)系數(shù)。 七、參考文獻(xiàn)1 方積乾，醫(yī)用統(tǒng)計(jì)學(xué)與電腦實(shí)驗(yàn)，上海：上?？茖W(xué)技術(shù)出版社1988。 2 郭祖超，醫(yī)學(xué)統(tǒng)計(jì)學(xué)，北京：人民軍醫(yī)出版社，1999。 3 孫振球，醫(yī)學(xué)統(tǒng)計(jì)學(xué)，北京：人民衛(wèi)生出版社，2002。4 邢航，多因素方差分析中數(shù)學(xué)模型的建立與檢驗(yàn)方法，電大理工，第2期 總第235期：1-6，2008。5 羅祎，論食品安全暴露評(píng)估模擬模型，食品科學(xué)， 1005- 9989(2007)02- 0021- 04：1-4,2007。6 劉淼，關(guān)于二維隨機(jī)概率密度函數(shù)算法的一

23、點(diǎn)補(bǔ)充，中央民族大學(xué)學(xué)報(bào)，2006，15（3）：246-247.7 虞學(xué)軍，復(fù)合因素作用下人體功能狀態(tài)數(shù)學(xué)模型的生理學(xué)分析，航天醫(yī)學(xué)與醫(yī)學(xué)工程，1999年第3期第12卷：1-5,1999。八、附錄表一網(wǎng)上2001年鎘的數(shù)據(jù)食品種類(lèi)國(guó)標(biāo)試驗(yàn)數(shù)平均數(shù)中位數(shù)糧食0.27370.02540.0598蔬菜0.056490.02180.0518水果0.033810.0070.0164肉類(lèi)0.17120.02690.0668蛋類(lèi)0.052780.01060.0269魚(yú)類(lèi)0.13100.01630.0419奶類(lèi)2430.01410.038奶粉350.00740.0179飲料2160.00850.0239零食2

24、190.02380.0553罐頭2320.01090.0278表二網(wǎng)上2001年汞的數(shù)據(jù)食品種類(lèi)國(guó)標(biāo)試驗(yàn)數(shù)平均數(shù)中位數(shù)糧食0.024310.00730.0135蔬菜0.012980.03520.0814水果0.01330.00090.0014肉類(lèi)0.05230.00290.0108皮蛋0.03370.01850.0367蛋類(lèi)0.0350.00820.0166魚(yú)類(lèi)0.370.00110.002奶類(lèi)0.0850.00340.0063奶粉0.01170.0020.0045魚(yú)肉罐頭0.3100.00570.0125果蔬罐頭150.00290.0048零食170.00370.0079果汁飲料130.00

25、240.0034表三網(wǎng)上2001年砷的數(shù)據(jù)食品種類(lèi)國(guó)標(biāo)試驗(yàn)數(shù)平均數(shù)中位數(shù)糧食0.77590.0580.11蔬菜0.56640.0550.101水果0.53910.0360.066肉類(lèi)0.57230.0580.111蛋類(lèi)780.0670.128魚(yú)類(lèi)0.52020.0480.091奶類(lèi)0.53230.0850.199奶粉0.2380.0240.04飲料1.62530.0290.047零食0.5900.170.542罐頭0.57590.0450.1表四網(wǎng)上2001年有機(jī)氯的數(shù)據(jù)食品種類(lèi)國(guó)標(biāo)試驗(yàn)數(shù)平均數(shù)中位數(shù)糧食0.33980.001790.00329蔬菜0.26750.006610.00955水果0

26、.23750.001990.00406肉類(lèi)0.43950.002760.00487魚(yú)類(lèi)0.52810.001340.0023奶粉250.010230.01731奶2110.000690.00144植物油1980.000810.00127茶葉0.5150.000020.0004表五網(wǎng)上2001年農(nóng)藥DDT食品種類(lèi)國(guó)標(biāo)試驗(yàn)數(shù)平均數(shù)中位數(shù)糧食0.24070.002560.00421蔬菜0.16750.004390.00887水果0.13750.004090.0083肉類(lèi)0.23960.002980.00271魚(yú)類(lèi)0.52610.005930.01741奶粉250.004260.00541奶0.121

27、10.001190.00163植物油0.51980.000430.00044茶葉0.2150.000060.00006表六網(wǎng)上2001年指標(biāo)分層測(cè)量數(shù)據(jù)等級(jí)組別商品名污染物分類(lèi)鎘汞砷有機(jī)氯農(nóng)藥DDT植物性糧食蔬菜水果水果果蔬罐頭動(dòng)物性肉類(lèi)魚(yú)類(lèi)魚(yú)類(lèi)魚(yú)肉罐頭蛋奶類(lèi)皮蛋蛋類(lèi)奶類(lèi)其他奶粉零食（植物油）飲料（茶葉）表七熵權(quán)法鎘在三類(lèi)食物中的比重植物性食物動(dòng)物性食物其他食物Xij0.4686350.4852790.2238080.4022140.2940540.2570770.1291510.2206750.519117Kij0.4686350.4852750.2238080.4022140.294052

28、0.2570760.1291510.2206730.519116wj0.4686350.4852750.223808ej0.4022140.2940520.257076表八熵權(quán)法汞在三類(lèi)食物中的比重植物性食物動(dòng)物性食物其他食物Xij0.1658150.1288660.2469140.7995460.169120.2962960.0346390.7020030.45679Kij0.1658150.1288670.2469140.7995440.1691220.2962950.0346390.7020110.456789Wj0.5400460.7400010.968208Ej0.6118480.3

29、458610.042291表九熵權(quán)法砷在三類(lèi)食物中的比重植物性動(dòng)物性其他食物Xij0.3892620.3093350.2157280.3691280.2560010.2606720.2416110.4346620.523599Kij0.3892620.3093360.2157280.3691280.2560020.2606720.2416110.4346630.5236Wj0.9815450.9775310.928558Ej0.1642390.1999660.635796表十熵權(quán)法有機(jī)氯在山中食物中的比重植物性食物動(dòng)物性食物其他食物Xij0.1722810.57620.9249550.6361890.2797490.0732370.191530.144050.001808Kij0.1722810.5762010.9249550.6361890.2797490.0732370.19153