深圳蔬菜質(zhì)量評(píng)價(jià)與抽樣優(yōu)化模型的算法實(shí)現(xiàn)

深圳蔬菜質(zhì)量評(píng)價(jià)與抽樣優(yōu)化模型的算法實(shí)現(xiàn)摘要蔬菜銷售行業(yè)在國(guó)民經(jīng)濟(jì)中占有不可忽視的地位，而蔬菜安全則是關(guān)乎國(guó)民健康和國(guó)家發(fā)展的重中之重。本文以深圳蔬菜抽樣體系為研究對(duì)象，在MATLAB、SPSS、SQL等軟件的數(shù)據(jù)處理基礎(chǔ)上，運(yùn)用指數(shù)平滑、模糊層次分析、多維灰聚類、風(fēng)險(xiǎn)矩陣等多種算法，從安全的總指標(biāo)出發(fā)，完成了對(duì)深圳蔬菜抽樣方法的綜合評(píng)價(jià)與優(yōu)化。同時(shí)，我們建立線性回歸、需求預(yù)測(cè)、曲線擬合等通用模型，實(shí)現(xiàn)全部問(wèn)題的解答。針對(duì)問(wèn)題一，讀取整合處理數(shù)據(jù)；將蔬菜分為八大類填充入庫(kù),計(jì)算各類蔬菜各月銷售比重；通過(guò)多維尺度分析法將八類蔬菜聚合為五類；用樣本平移法、指數(shù)平滑法，擬合填充缺失月份數(shù)據(jù)；用EXCEL繪制各月各類蔬菜銷量折線圖和全年各類蔬菜銷量均值比重餅圖；計(jì)算農(nóng)貿(mào)市場(chǎng)、批發(fā)市場(chǎng)、商場(chǎng)超市、生產(chǎn)基地各渠道銷售量所占比重；根據(jù)交易記錄，計(jì)算新加渠道（零售、電商）的銷售比重(6.8%，0.3%)；分析圖表得出結(jié)論：綠色蔬菜和瓜豆類蔬菜冬季波動(dòng)最為明顯，全年銷量以綠色蔬菜為甚(56.5%)，蔥蒜類銷量最次(0.4%)，選擇農(nóng)貿(mào)市場(chǎng)(47.0%)和批發(fā)市場(chǎng)(32.6%)購(gòu)買蔬菜的市民最多，因此后續(xù)研究應(yīng)當(dāng)特別注意上述類別、季節(jié)和銷售渠道的蔬菜品質(zhì)。針對(duì)問(wèn)題二，在SQL中計(jì)算農(nóng)貿(mào)市場(chǎng)(5.450%)、批發(fā)市場(chǎng)(2.750%)、商場(chǎng)超市(8.630%)、生產(chǎn)基地(1.480%)四大銷售渠道的抽檢不合格率（零售和電商的不合格率缺失）；對(duì)各渠道進(jìn)行打分構(gòu)造優(yōu)先關(guān)系矩陣，再用MATLAB計(jì)算六大渠道風(fēng)險(xiǎn)權(quán)重；將前四類渠道的權(quán)重與不合格率做線性擬合，算得零售(7.270%)、電商(6.240%)不合格率；利用公式p=i=06piqi=4.804%計(jì)算平均不合格率，以及忽略某渠道后剩余渠道不合格率的改變值?pi；分析知針對(duì)問(wèn)題三，從三方面分析：對(duì)于鮮菜中風(fēng)險(xiǎn)物質(zhì)成分，采用多維灰聚類法分析，構(gòu)造了含有五個(gè)指標(biāo)的權(quán)重矩陣來(lái)判斷風(fēng)險(xiǎn)物質(zhì)的灰度特征，確定了21個(gè)中高風(fēng)險(xiǎn)的鮮菜樣本。對(duì)樣本中出現(xiàn)的頻數(shù)較高以及危險(xiǎn)性較高、超標(biāo)嚴(yán)重的分析物質(zhì)進(jìn)行匯總分析。由此得到深圳市在檢測(cè)物質(zhì)的食品安全方面仍存在一定的風(fēng)險(xiǎn)的結(jié)論；針對(duì)地域因素，通過(guò)計(jì)算各省樣本的抽查合格率對(duì)出口省份進(jìn)行安全性排名，運(yùn)用K-聚類分析法綜合考慮多種因素將可輸出省份分為3類，可見(jiàn)該市進(jìn)口蔬菜中存在風(fēng)險(xiǎn)較高者；結(jié)合深圳市民的飲食習(xí)慣和反季節(jié)蔬菜的特性，得出結(jié)論：深圳市面臨反季節(jié)鮮菜以及本地受污染蔬菜風(fēng)險(xiǎn)。針對(duì)問(wèn)題四，按照簡(jiǎn)單隨機(jī)抽樣原則估計(jì)樣本量，蔬菜抽樣總量為全年11400批次，每月950批次。根據(jù)各轄區(qū)蔬菜樣品不合格率及原蔬菜抽樣比例，規(guī)定各轄區(qū)每月抽樣任務(wù)；針對(duì)不同季節(jié)制定八大類蔬菜抽樣方案；根據(jù)問(wèn)題一二得到的六類渠道銷量及不合格率，規(guī)定抽樣場(chǎng)所的比例，采用不同抽樣方式，農(nóng)貿(mào)市場(chǎng)、生產(chǎn)基地、批發(fā)市場(chǎng)采用分層抽樣及系統(tǒng)抽樣法，零售攤販采用分級(jí)抽樣，電商渠道對(duì)貨源和倉(cāng)庫(kù)中的產(chǎn)品以及物流過(guò)程、配送過(guò)程中抽樣。針對(duì)問(wèn)題五，綜合前幾題的結(jié)論，就深圳蔬菜安全問(wèn)題向市政府提出建議。關(guān)鍵詞：深圳蔬菜模糊層次分析多維尺度分析多維灰聚類指數(shù)平滑K-聚類分析風(fēng)險(xiǎn)矩陣

問(wèn)題重述問(wèn)題背景隨著時(shí)代發(fā)展，人類生活得到了極大的便利。但與此同時(shí)，技術(shù)手段的不完善、食品抽檢方案的不健全、非法分子管制措施的相對(duì)缺失，都使得食品安全問(wèn)題日趨嚴(yán)重。眾所周知，民以食為天。蔬菜作為極富營(yíng)養(yǎng)價(jià)值的天然食品，其監(jiān)督自然不容小覷。而深圳作為近40年前改革開(kāi)放的窗口，以及擁有2000多萬(wàn)人口廣東第二大城市，其新鮮蔬菜的消費(fèi)以輸入型為主，輸入渠道、生產(chǎn)地區(qū)、氣候條件、管理模式、菜品種類等都具有多樣化的特點(diǎn)。因此，它的食品安全監(jiān)管也更應(yīng)當(dāng)提上日程。建立合理有效的深圳蔬菜安全評(píng)估模型和抽檢方案，將對(duì)人類健康和社會(huì)穩(wěn)定提供極大保障。問(wèn)題提出上述諸多因素都對(duì)傳統(tǒng)的蔬菜安全監(jiān)管工作方法構(gòu)成新的挑戰(zhàn)。附件1是深圳市2017年食用農(nóng)產(chǎn)品監(jiān)測(cè)抽檢方案。附件2是深圳市2017年1-3月食用農(nóng)產(chǎn)品監(jiān)測(cè)抽檢詳細(xì)計(jì)劃。附件3是近兩年全國(guó)蔬菜及其制品抽樣檢驗(yàn)不合格列表。我們的任務(wù)是根據(jù)背景材料和收集的數(shù)據(jù)，建立數(shù)學(xué)模型回答以下問(wèn)題：?jiǎn)栴}一：收集數(shù)據(jù)，調(diào)查、分析深圳市各季節(jié)新鮮蔬菜消費(fèi)種類及不同輸入渠道的比例。問(wèn)題二：評(píng)估檢驗(yàn)中忽略某些輸入渠道的抽樣導(dǎo)致的食品安全風(fēng)險(xiǎn)。問(wèn)題三：基于附件評(píng)估深圳新鮮蔬菜食品安全風(fēng)險(xiǎn)。問(wèn)題四：根據(jù)科學(xué)的統(tǒng)計(jì)推斷原理，設(shè)計(jì)一套深圳市新鮮蔬菜抽樣檢驗(yàn)方案，以滿足不同層面監(jiān)測(cè)工作的需要，提高食品安全防控的效率與質(zhì)量。問(wèn)題五：基于你的研究結(jié)果給深圳市政府寫一封建議信。問(wèn)題分析問(wèn)題一要求我們根據(jù)附件和自行搜集的資料，提取各個(gè)季度深圳各類新鮮蔬菜的消費(fèi)情況比例，以及各類蔬菜輸入渠道的比例。針對(duì)該問(wèn)題，我們首先根據(jù)近兩年深圳蔬菜質(zhì)檢結(jié)果信息，使用SQL對(duì)搜得數(shù)據(jù)進(jìn)行預(yù)處理和簡(jiǎn)單篩查，得出各月份各類蔬菜銷售情況表格。再通過(guò)相關(guān)存儲(chǔ)過(guò)程的建立，計(jì)算各類銷售渠道占比。最終，使用樣本平移法對(duì)所得數(shù)據(jù)進(jìn)行檢驗(yàn)，使模型更具合理性。問(wèn)題二要求我們分析各輸入渠道的可忽略性極其影響。針對(duì)該問(wèn)題，我們依然利用深圳蔬菜質(zhì)檢報(bào)告，基于問(wèn)題一的結(jié)論，通過(guò)大量數(shù)據(jù)在SQL數(shù)據(jù)庫(kù)中歸類分析，計(jì)算各渠道蔬菜產(chǎn)品的比重和不合格率。再根據(jù)運(yùn)籌學(xué)中的模糊層次分析法，計(jì)算各渠道的安全風(fēng)險(xiǎn)權(quán)重。接著，通過(guò)已有渠道信息散點(diǎn)圖的EXCEL函數(shù)擬合，預(yù)測(cè)新加渠道的不合格率。最后，比較刪去特定渠道對(duì)總不合格率的影響，得出風(fēng)險(xiǎn)評(píng)估結(jié)果。問(wèn)題三要求我們基于附件評(píng)估深圳市的新鮮蔬菜食品安全風(fēng)險(xiǎn)?；谏钲谑絮r菜的輸入性特征，我們將影響新鮮蔬菜食品安全的因數(shù)主要分為三類：鮮菜中存在的危險(xiǎn)物質(zhì)帶來(lái)的物質(zhì)性風(fēng)險(xiǎn)、生產(chǎn)管理不規(guī)范或污染嚴(yán)重區(qū)域生產(chǎn)的鮮菜產(chǎn)品帶來(lái)的地域性風(fēng)險(xiǎn)、生產(chǎn)方式不同以及本地飲食習(xí)慣不同帶來(lái)的其他風(fēng)險(xiǎn)。運(yùn)用灰度聚類法對(duì)未檢測(cè)危險(xiǎn)物質(zhì)進(jìn)行分類并且確定具有高度及中度風(fēng)險(xiǎn)的物質(zhì)，并確定其比例，以對(duì)危險(xiǎn)物質(zhì)帶來(lái)的風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià)。對(duì)于已檢測(cè)的物質(zhì)，我們采用分類匯總和樣本描述的方法，以得到已檢測(cè)物質(zhì)中出現(xiàn)漏檢錯(cuò)檢可能存在的風(fēng)險(xiǎn)。對(duì)于地域性風(fēng)險(xiǎn)，我們先通過(guò)樣本的檢驗(yàn)不合格率對(duì)于全國(guó)可能出口的省市進(jìn)行劃分，并通過(guò)K-聚類分析結(jié)合鐵路距離、平均運(yùn)價(jià)、季節(jié)性、是否有出口記錄等特性對(duì)出口省市進(jìn)行分類，得出主要出口、次要出口、以及極少出口的各個(gè)省市，結(jié)合檢驗(yàn)不合格率對(duì)于鮮菜產(chǎn)品帶來(lái)的地域性風(fēng)險(xiǎn)進(jìn)行評(píng)估。對(duì)于生產(chǎn)方式不同以及本地飲食習(xí)慣不同帶來(lái)的其他風(fēng)險(xiǎn)，結(jié)合實(shí)際數(shù)據(jù)以及反季節(jié)蔬菜的特性對(duì)該類風(fēng)險(xiǎn)進(jìn)行定性分析，得到其他風(fēng)險(xiǎn)的評(píng)價(jià)。問(wèn)題四要求我們根據(jù)科學(xué)的統(tǒng)計(jì)推斷原理，設(shè)計(jì)一套滿足不同層面監(jiān)測(cè)工作的需要的深圳市新鮮蔬菜抽樣方案，并提高食品安全防控的效率與質(zhì)量。針對(duì)該問(wèn)題，我們首先分析了現(xiàn)有抽樣方案的不足之處，之后采用簡(jiǎn)單隨機(jī)抽樣的方法估計(jì)樣本量，依據(jù)各轄區(qū)新鮮蔬菜樣品不合格率及原新鮮蔬菜抽樣比例，為各轄區(qū)分配具體抽樣任務(wù)。根據(jù)問(wèn)題一、問(wèn)題二所得結(jié)論，確定不同抽樣渠道抽樣比例以及不同季節(jié)蔬菜抽樣比例。不同渠道的抽樣采用分層抽樣、多階段抽樣、系統(tǒng)抽樣、簡(jiǎn)單隨機(jī)抽樣相結(jié)合的方式，以更加全面準(zhǔn)確地覆蓋區(qū)域的情況，提高評(píng)價(jià)精度。問(wèn)題模型假設(shè)假設(shè)所收集數(shù)據(jù)皆為可靠數(shù)據(jù)，不存在錯(cuò)誤信息。假設(shè)質(zhì)檢報(bào)告樣本中各類蔬菜和輸入渠道比例與實(shí)際市場(chǎng)蔬菜銷售和比例完全相同。假設(shè)每一年相同時(shí)間蔬菜消費(fèi)種類、輸入渠道比例相同，市民飲食習(xí)慣不變，忽略不同年份天氣、氣候、政策等影響可能帶來(lái)的差異。假設(shè)模型涉及的所有蔬菜所屬分類都包含于本文所列八大類中，不存在模糊分類的情況，蔬菜與類型一一對(duì)應(yīng)。假設(shè)各城市間鐵路長(zhǎng)度即為蔬菜運(yùn)輸距離，且各類蔬菜生產(chǎn)省份到深圳的距離均以該省省會(huì)城市計(jì)算，運(yùn)價(jià)以普通運(yùn)輸價(jià)格估計(jì)。假設(shè)附件中收集的全國(guó)數(shù)據(jù)中的不合格率與全國(guó)鮮菜出口樣品中的不合格率相同假設(shè)鮮菜樣品中存在的所有危險(xiǎn)物質(zhì)均包括在全國(guó)的檢驗(yàn)報(bào)告中。假設(shè)2015年的全國(guó)農(nóng)村統(tǒng)計(jì)年鑒中有關(guān)于鮮菜產(chǎn)量的數(shù)據(jù)具有平常性和代表性。假定各個(gè)省份的蔬菜抽樣比例相同。文中不同渠道運(yùn)輸?shù)臏囟纫约跋鄬?duì)濕度都大致相同，而在蔬菜的運(yùn)輸流通環(huán)節(jié)中，運(yùn)輸?shù)臅r(shí)間越長(zhǎng)，蔬菜變質(zhì)的安全風(fēng)險(xiǎn)越高。在蔬菜生產(chǎn)環(huán)節(jié)中，各個(gè)銷售渠道各自獲得的蔬菜的品質(zhì)大體一致。符號(hào)說(shuō)明符號(hào)含義符號(hào)含義不安全性評(píng)價(jià)指標(biāo)平均不合格率農(nóng)藥殘留第渠道所占比例新鮮程度n最小樣本量非法添加樣本總體數(shù)農(nóng)貿(mào)市場(chǎng)安全風(fēng)險(xiǎn)u正態(tài)變量值批發(fā)市場(chǎng)安全風(fēng)險(xiǎn)不合格率商場(chǎng)超市安全風(fēng)險(xiǎn)δ容許誤差生產(chǎn)基地安全風(fēng)險(xiǎn)各因素相對(duì)不安全性的優(yōu)先關(guān)系零售安全風(fēng)險(xiǎn)各銷售渠道相對(duì)因素的優(yōu)先關(guān)系電商安全風(fēng)險(xiǎn)對(duì)應(yīng)的模糊一致陣元素層對(duì)層的權(quán)重向量第渠道的蔬菜不合格率層對(duì)因素的權(quán)重向量總權(quán)重向量模型建立與求解5.1問(wèn)題一的模型建立與求解5.1.1問(wèn)題一的分析本題要求我們收集調(diào)查數(shù)據(jù)，并分析深圳各季蔬菜消費(fèi)種類及輸入渠道比例。針對(duì)這一問(wèn)題，我們給出如下求解思路和具體步驟：步驟一：根據(jù)深圳市市場(chǎng)和監(jiān)督管理委員會(huì)給出的官方數(shù)據(jù)（深圳市蔬菜質(zhì)量安全監(jiān)測(cè)結(jié)果），導(dǎo)入SQL數(shù)據(jù)庫(kù)進(jìn)行預(yù)處理，剔除無(wú)用數(shù)據(jù)和不詳數(shù)據(jù)。步驟二：在已有數(shù)據(jù)中添加蔬菜的八大分類信息列，通過(guò)查詢和更新語(yǔ)句進(jìn)行完全填充。步驟三：利用多維尺度分析法將蔬菜八大類細(xì)分為五類。步驟四：使用樣本平移法對(duì)缺失數(shù)據(jù)進(jìn)行估計(jì)，再用指數(shù)平滑法將數(shù)據(jù)擬合。步驟五：建立各季節(jié)、各月份、各抽樣場(chǎng)所（即輸入渠道）抽檢情況視圖，并編寫存儲(chǔ)過(guò)程分別計(jì)算相應(yīng)時(shí)段、相應(yīng)渠道各類別蔬菜占比。步驟六：將所得結(jié)果導(dǎo)入EXCEL，繪制消費(fèi)折線圖和蔬菜類別餅圖。5.1.2菜品樣本分類的確定與數(shù)據(jù)填充查閱原題所給的附件一，易知蔬菜可分為鱗莖類、葉菜類、蕓苔屬類、瓜類、茄果類、豆類、根莖類、水生類和其他小類。為了數(shù)據(jù)處理的便捷和結(jié)果的合理性，我們對(duì)上述種類進(jìn)行了修正，得到根菜類、鮮豆類、茄果瓜菜類、蔥蒜類、嫩莖葉花菜類、水生蔬菜類、薯芋類和野生蔬菜類八大類，具體信息見(jiàn)下表。表1.1深圳市蔬菜銷售種類[1]種類名稱根菜類白蘿卜、青蘿卜、胡蘿卜等鮮豆類豌豆、豆芽、荷蘭豆等茄果瓜菜類茄子、番茄、甜椒、黃瓜等蔥蒜類蒜苗、蔥等嫩莖葉花菜類大白菜、上海青、筍、西蘭花等水生蔬菜類蓮藕、茭白等薯芋類山藥、姜、芋頭等野生蔬菜類枸杞、香椿等至此，深圳蔬菜消費(fèi)種類信息已在數(shù)據(jù)庫(kù)中補(bǔ)充完畢。5.1.3基于多維尺度分析方法的蔬菜分類細(xì)分考慮到八個(gè)種類的規(guī)律性和走勢(shì)的復(fù)雜性，在后續(xù)處理中可能很難對(duì)其進(jìn)行分析，在此我們使用SPSS中的多維尺度分析法對(duì)蔬菜類別做進(jìn)一步細(xì)分。經(jīng)過(guò)多方資料的收集，我們認(rèn)為影響蔬菜種類相似度的因素有如下幾種：營(yíng)養(yǎng)成分：蛋白質(zhì)、脂肪、淀粉、水、無(wú)機(jī)鹽和維生素六大類；食用器官：根、莖、葉、花、果實(shí)和種子六種；生產(chǎn)地域：各高原、平原、江河流域等；可能危害相似度：農(nóng)藥殘留、非法添加、微生物、食品添加劑等。以上述指標(biāo)為細(xì)則，我們對(duì)已經(jīng)分好的八類蔬菜進(jìn)行打分，詳見(jiàn)下表：表1.2蔬菜相似度評(píng)分表根菜類鮮豆類茄果瓜菜類蔥蒜類嫩莖葉花菜類水生蔬菜類薯芋類野生蔬菜類根菜類0鮮豆類2.5170茄果瓜菜類21.50蔥蒜類2.252.5830嫩莖葉花菜類2.252.252.752.920水生蔬菜類2.771.52.711.670.720薯芋類1.3332.633.082.83.170野生蔬菜類3.253.423.32.82.051.673.50基于矩陣元素分值越低越相近的原則，此表格可以在概念圖中得以直觀呈現(xiàn)：圖1.1蔬菜種類的二維尺度空間概念圖本文采取的多維尺度分析法并不同于以圖表形式輸出結(jié)果的其他方法，我們得到的是一系列的數(shù)值說(shuō)明：FormatrixStress=.13195RSQ=.88836其中Stress值（應(yīng)力值）為0.13說(shuō)明本次擬合水平（信度）一般，RSQ值為0.89>0.60說(shuō)明本次擬合通過(guò)檢驗(yàn)，其效度是可以接受的。綜上，我們可將新鮮蔬菜聚合為瓜豆類、蔥蒜類、綠色蔬菜類、野生蔬菜類和地下作物類五種：瓜豆類：茄果瓜菜類+鮮豆類；解釋：二者營(yíng)養(yǎng)成分相似，主要都為糖類和淀粉；且二者的食用部分都是果實(shí)。蔥蒜類：蔥蒜類；綠色蔬菜類：嫩莖葉花菜類+水生蔬菜類；解釋：二者營(yíng)養(yǎng)成分都是維生素和礦物質(zhì)，且食用部分都以根莖葉為主。（4）野生蔬菜類：野生蔬菜類；（5）地下作物類：薯芋類+根菜類。解釋：二者食用部分都以根莖為主。至此，蔬菜細(xì)分完成。5.1.4缺失數(shù)據(jù)處理初始數(shù)據(jù)信息量龐大，且由多個(gè)文檔整合而成。考慮到后續(xù)數(shù)據(jù)處理的便捷性需求，我們將其導(dǎo)入數(shù)據(jù)庫(kù)中，通過(guò)查詢語(yǔ)句剔除空行和無(wú)效列，并增加菜品分類屬性列和不同時(shí)段不同場(chǎng)所抽檢視圖，以備后續(xù)計(jì)算之需。具體過(guò)程參見(jiàn)附錄。通過(guò)軟件對(duì)數(shù)據(jù)屬性的歸類，我們發(fā)現(xiàn)2016年部分月份的數(shù)據(jù)處于缺失狀態(tài)，且特定區(qū)段的數(shù)據(jù)規(guī)律性欠佳。我們收集到的數(shù)據(jù)是深圳市市場(chǎng)和質(zhì)量監(jiān)督管理委員會(huì)發(fā)布的2016年7月至2017年2月的《深圳市農(nóng)產(chǎn)品質(zhì)量安全豬肉和蔬菜例行監(jiān)測(cè)結(jié)果》。對(duì)于缺失數(shù)據(jù)，我們采用了移動(dòng)平均法和指數(shù)平滑法，最終，我們采用指數(shù)平滑法得到結(jié)果。樣本平移法[2]我們首先采用了一種根據(jù)真實(shí)數(shù)據(jù)構(gòu)造采樣點(diǎn)的方法來(lái)解決數(shù)據(jù)不完整的問(wèn)題。其基本思想類似于連續(xù)函數(shù)中的插值，在采樣環(huán)節(jié)中，所采集的樣本數(shù)據(jù)是離散的，因此和連續(xù)函數(shù)不同，構(gòu)造出的為樣本數(shù)據(jù)也應(yīng)該是離散的，并且統(tǒng)計(jì)上服從相似形狀的分布，稱之為“偽樣本空間”。偽樣本空間應(yīng)該滿足：樣本標(biāo)志值（被測(cè)量值）的概率密度函數(shù)和真實(shí)樣本空間樣本標(biāo)志值具有相似的分布，即概率密度函數(shù)的形式不發(fā)生變化。采樣數(shù)據(jù)是抽樣率很低的隨機(jī)抽樣數(shù)據(jù)，而且由于蔬菜的季節(jié)性、區(qū)域性、多樣性、差異性，抽樣無(wú)法得到詳細(xì)、完整的數(shù)據(jù)，缺少了三至六月的數(shù)據(jù)，而蔬菜消費(fèi)種類的比例隨月份有明顯的變化，這樣會(huì)對(duì)各季節(jié)新鮮蔬菜消費(fèi)比例的評(píng)估造成不良影響。所以，針對(duì)三至六月的數(shù)據(jù)分析采用了樣本平移法。我們將廣州江南果菜批發(fā)市場(chǎng)的每月蔬菜銷售情況作為參考，假設(shè)其每月蔬菜種類銷售變化情況反映深圳市的蔬菜種類銷售情況。原因是廣州江南果菜批發(fā)市場(chǎng)是中國(guó)乃至東南亞地區(qū)最大的果菜集散地之一，蔬菜交易量一直穩(wěn)居全國(guó)第一，每天的蔬菜成交量達(dá)1000萬(wàn)公斤，占廣州市蔬菜上市量70%，具有參考性。根據(jù)廣州江南果菜市場(chǎng)2015年5月、2016年3月、4月、6月、7月蔬菜銷售情況分析，可以總結(jié)得到：三月，冬儲(chǔ)菜的交易步入尾聲，新的蔬菜尚未批量上市，市場(chǎng)供應(yīng)的蔬菜主要以冬儲(chǔ)菜為主，其中銷量比較大的為大白菜、蘿卜、土豆。海南的反季節(jié)瓜豆進(jìn)入盛產(chǎn)期，三月還沒(méi)到本地和北方瓜豆的上市期。四月，市場(chǎng)銷售量較大的以供應(yīng)到飯?zhí)玫拇箢愂卟藶橹鳎渲兴脑落N量排名前三的蔬菜有大白菜、土豆、青瓜，其中，青瓜的供應(yīng)量明顯增多，銷量提高了三成左右。在換季階段，市民普遍食用的蔬菜由葉菜類蔬菜逐漸向瓜豆類蔬菜轉(zhuǎn)移。海南瓜豆類上市量減少，到月下旬只有較少的上市量，本地瓜豆在下旬開(kāi)始少量上市。五月，市場(chǎng)又開(kāi)始進(jìn)入瓜豆類蔬菜的銷售旺季。本月銷售量最大的蔬菜為土豆、胡蘿卜、蘿卜、洋蔥頭、大白菜、椰菜、萵筍、大肉姜、青瓜、南瓜。五月中下旬，高溫、多雨、潮濕的天氣不利于葉菜類生長(zhǎng)，菜葉容易腐爛，葉菜類供應(yīng)減少。海南瓜豆類上市量將進(jìn)一步減少，到月底估計(jì)只有微量的海南蔬菜上市，來(lái)自北方的蔬菜和本地瓜豆類蔬菜，供應(yīng)量逐漸增加，成為市場(chǎng)銷售的主流蔬菜。瓜豆類、茄果類蔬菜供應(yīng)量越來(lái)越多，葉菜類蔬菜供應(yīng)量越來(lái)越少。新產(chǎn)洋蔥和新產(chǎn)的土豆也開(kāi)始大量上市。六月，銷量排名前三的蔬菜有紹菜、椰菜、西紅柿，銷量環(huán)比上漲。北方露天蔬菜進(jìn)入采收旺季，以瓜豆類為主的蔬菜供應(yīng)量持續(xù)增加，瓜豆類蔬菜價(jià)廉物美，是市民最熱愛(ài)的選擇。葉類蔬菜水分較大，在適宜的溫度和日照下生長(zhǎng)周期縮短，市場(chǎng)供應(yīng)量增加。七月，銷量排名前三的蔬菜有土豆、紹菜、西紅柿，土豆銷售量大幅增加，但椰菜、白蘿卜、青瓜等排名較前的蔬菜在本月銷售量均有所下降，原因在于現(xiàn)在市場(chǎng)銷售的外地蔬菜主要來(lái)自于北方，受到北方地區(qū)暴雨天氣的影響，供應(yīng)到江南市場(chǎng)的外地蔬菜有所減少。產(chǎn)地來(lái)自于湖北的蔬菜供應(yīng)量下降尤為明顯，如大白菜、椰菜、白蘿卜等，銷量在本月有兩成左右的降幅，北方雨水天氣是本月江南市場(chǎng)外地蔬菜銷量減少的主要原因。據(jù)此，我們得到深圳市一年蔬菜銷量估計(jì)值。指數(shù)平滑法對(duì)于缺失值的預(yù)測(cè)，我們又采用移動(dòng)平均法結(jié)合簡(jiǎn)單季節(jié)性的修正對(duì)缺失的幾個(gè)月的數(shù)據(jù)進(jìn)行擬合，移動(dòng)平均法則不考慮較遠(yuǎn)期的數(shù)據(jù)。并在加權(quán)移動(dòng)平均法中給予近期資料更大的權(quán)重；而指數(shù)平滑法則兼容了全期平均和移動(dòng)平均所長(zhǎng)，不舍棄過(guò)去的數(shù)據(jù)，但是僅給予逐漸減弱的影響程度，即隨著數(shù)據(jù)的遠(yuǎn)離，賦予逐漸收斂為零的權(quán)數(shù)。簡(jiǎn)單季節(jié)性修正模型適用于沒(méi)有趨勢(shì)并且季節(jié)性影響隨時(shí)間變動(dòng)保持恒定的序列。其平滑參數(shù)是水平和季節(jié)。簡(jiǎn)單季節(jié)性指數(shù)平滑法與ARIMA模型極為相似，包含零階自回歸、一階差分、一階季節(jié)性差分和一階、階和階移動(dòng)平均數(shù)。根據(jù)指數(shù)平滑理論，我們得出五類蔬菜模型的平穩(wěn)方，以及總模型的平穩(wěn)方，如下圖所示：表1.3各類蔬菜模型統(tǒng)計(jì)量綜合模型模型擬合統(tǒng)計(jì)量平穩(wěn)的R方瓜豆類-模型_1.896蔥蒜類-模型_2.868地下作物類-模型_3.860綠色蔬菜類-模型_4.689野生蔬菜類-模型_5.795總模型平穩(wěn)的R方.860以之前兩年所有數(shù)據(jù)為基礎(chǔ)，我們擬合出2016年3月至2016年6月各類蔬菜銷量波動(dòng)值，以折線圖給出：圖1.2各類蔬菜模型統(tǒng)計(jì)量綜合因變量為季節(jié)性數(shù)據(jù)，因此，采用平穩(wěn)的R方來(lái)描述擬合程度更加合適。另外，由表1.3可知，在模型擬合和模型統(tǒng)計(jì)量中，平穩(wěn)的R方均大于0.65，所以擬合的程度較好。從圖1.2也可以看出，觀測(cè)值與擬合值具有良好的一致性。至此，我們已完成了所有缺失數(shù)據(jù)的合理填補(bǔ)。5.1.5深圳市各季節(jié)新鮮蔬菜消費(fèi)情況及種類計(jì)算在分類完成的前提下，我們針對(duì)問(wèn)題一的題干要求，首先在數(shù)據(jù)庫(kù)中建立各月份、各季度的抽檢結(jié)果視圖。接下來(lái)，我們通過(guò)一系列存儲(chǔ)過(guò)程的建立，對(duì)原始信息進(jìn)行抽取，分別計(jì)算相應(yīng)時(shí)段的蔬菜種類。最終，將所有數(shù)據(jù)導(dǎo)入EXCEL，整合四個(gè)季度的數(shù)據(jù)，求得全年均值。前文已經(jīng)提及，考慮到2016年3月至6月相關(guān)數(shù)據(jù)的不完善，我們需要對(duì)此進(jìn)行彌補(bǔ)和修正。綜合近兩年的同類信息，我們向原數(shù)據(jù)中輸入大量基數(shù)樣本，并使用指數(shù)平滑法對(duì)新樣本進(jìn)行擬合，獲得相應(yīng)月份的蔬菜抽樣比例。2016年7月至2017年2月各月蔬菜消費(fèi)比例情況如下表所示：表1.4全年每月各類蔬菜消費(fèi)比例表瓜豆類蔥蒜類綠色蔬菜類野生蔬菜類地下作物類1月31.7%0.7%58.5%2.1%6.9%2月34.9%1.0%55.1%1.7%7.2%3月36.5%0.1%56.5%0.6%6.3%4月36.6%0.2%57.6%0.4%5.2%5月38.4%0.3%54.9%0.4%6.0%6月38.0%0.2%54.8%3.1%3.9%7月38.5%0.0%52.4%3.8%5.3%8月38.0%1.8%50.5%4.4%5.3%9月32.2%0.6%59.6%4.1%3.5%10月29.7%0.0%63.3%3.8%3.1%11月32.0%0.0%57.9%5.0%5.1%12月31.2%0.3%57.3%4.7%6.3%若按季度統(tǒng)計(jì)，可簡(jiǎn)化上表如下：表1.5全年各季度各類蔬菜消費(fèi)比例表瓜豆類蔥蒜類綠色蔬菜類野生蔬菜類地下作物類總計(jì)第一季度34.39%0.60%56.74%1.47%6.80%1第二季度37.67%0.23%55.77%1.30%5.03%1第三季度36.23%0.80%54.17%4.10%4.70%1第四季度31.00%0.10%59.56%4.50%4.84%1將表1.5繪制為折線圖，可得如下結(jié)果：圖1.3全年每月各類蔬菜消費(fèi)比例折線圖由上圖可知：綠色蔬菜于秋冬季節(jié)進(jìn)入銷售旺季，而春夏季銷量略減；相反地，瓜豆類蔬菜在秋冬季節(jié)銷售情況并不理想，而春夏季略有上升；其余蔬菜基數(shù)較小，此處不作分析。我們可以綜合深圳的季節(jié)特性對(duì)這一情況作解釋：每年春初本地蔬菜大量上市，而瓜豆類蔬菜未進(jìn)入上市期——因此綠色蔬菜銷量較高，而瓜豆類蔬菜稍低；4-5月份春夏之際，隨著季節(jié)的轉(zhuǎn)變，市民的蔬菜選擇也在改變，即從葉菜類蔬菜向瓜豆類蔬菜的轉(zhuǎn)移——因此綠色蔬菜銷量緩慢下降，而瓜豆類銷量緩慢上升；11月左右，天氣轉(zhuǎn)涼，葉菜類開(kāi)始大量上市——因此綠色蔬菜銷量出現(xiàn)高峰。2016年7月至2017年2月全年蔬菜消費(fèi)平均比例情況如下表所示：圖1.4全年蔬菜消費(fèi)比例均值餅圖由上圖可知，深圳人的蔬菜飲食以綠色蔬菜類為主，瓜豆類次之，蔥蒜類最少。從深圳人的習(xí)慣看，蔬菜谷類在整個(gè)飲食結(jié)構(gòu)中占比較小，而肉類居多。而在蔬菜大類內(nèi)，綠色蔬菜和瓜豆顯然是多數(shù)人的選擇，蔥蒜多做調(diào)味品，因此比重最低。5.1.6深圳市新鮮蔬菜輸入渠道比例計(jì)算根據(jù)題干要求，我們首先按照原數(shù)據(jù)給出的農(nóng)貿(mào)市場(chǎng)、批發(fā)市場(chǎng)、商場(chǎng)超市、生產(chǎn)基地四類輸入渠道（即蔬菜來(lái)源）對(duì)抽檢信息進(jìn)行分類，并在數(shù)據(jù)庫(kù)中分別創(chuàng)建視圖。其后，我們?cè)诮⒚麨椤扒辣壤钡拇鎯?chǔ)過(guò)程，對(duì)上述四類渠道的占用比例進(jìn)行統(tǒng)一計(jì)算，得到以下初始比例圖表：表1.6原有四類渠道所占比重農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地50.6%35.1%9.6%4.7%由于題干中提示“新鮮蔬菜消費(fèi)屬性是輸入型的”，且近年來(lái)蔬菜輸入渠道日趨多樣化，我們?cè)诓殚喍喾劫Y料后決定在表中增設(shè)零售和電商兩個(gè)渠道。零售：根據(jù)某一特定年份內(nèi)批發(fā)市場(chǎng)和零售商販的相對(duì)比例，可計(jì)算出零售渠道的銷售比重。電商：據(jù)統(tǒng)計(jì)，廣東省蔬菜行業(yè)市場(chǎng)規(guī)模為3950億元，而2015年廣東省農(nóng)產(chǎn)品電商交易額約129億元。其中，生鮮農(nóng)產(chǎn)品（果蔬肉類）電商交易額占農(nóng)產(chǎn)品總交易額的18.6%，蔬菜在果蔬肉類生鮮農(nóng)產(chǎn)品電商交易額中占比最少，約為1/2。經(jīng)過(guò)上述修正，我們得出更為合理的六類蔬菜銷售渠道比例圖表：表1.7修正后六類渠道所占比重農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地零售電商比例47.0%32.6%8.9%4.4%6.8%0.3%由上圖可知，選擇農(nóng)貿(mào)市場(chǎng)進(jìn)行蔬菜消費(fèi)的市民最多，批發(fā)市場(chǎng)次之，選擇電商（即網(wǎng)購(gòu)蔬菜）的市民最少。5.2問(wèn)題二的模型建立與求解5.2.1問(wèn)題二的分析本題需要我們基于所獲數(shù)據(jù)評(píng)估抽檢方案忽略特定輸入渠道時(shí)導(dǎo)致的安全風(fēng)險(xiǎn)。為解決該問(wèn)題，我們需要建立模型對(duì)缺失渠道的不合格率進(jìn)行分析，具體步驟如下：步驟一：從上題數(shù)據(jù)出發(fā)，在數(shù)據(jù)庫(kù)中編寫存儲(chǔ)過(guò)程計(jì)算各類已知渠道的不合格率。步驟二：通過(guò)各類因素的補(bǔ)充，利用模糊層次分析法，在MATLAB中實(shí)現(xiàn)各輸入渠道風(fēng)險(xiǎn)權(quán)重的求解。步驟三：將所得結(jié)果導(dǎo)入EXCEL，對(duì)各渠道的權(quán)重和不合格率進(jìn)行描點(diǎn)和函數(shù)擬合，計(jì)算新增渠道的不合格率。步驟四：綜合考慮各渠道所占比例，分析將其忽略所導(dǎo)致的安全風(fēng)險(xiǎn)。步驟五：使用故障樹法進(jìn)行輔助分析。5.2.2已知渠道不合格率的計(jì)算在上題建立四類渠道視圖的基礎(chǔ)上，我們?cè)跀?shù)據(jù)庫(kù)中繼續(xù)編寫名為“合格率”的存儲(chǔ)過(guò)程，對(duì)每一渠道的不合格率進(jìn)行計(jì)算，結(jié)果如下圖:表2.1原有四類渠道不合格率農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地5.45%2.75%8.63%1.48%5.2.3基于模糊層次分析法的安全評(píng)估模型建立問(wèn)題一考慮了零售攤販和電商銷售渠道，而這兩個(gè)渠道的不合格率無(wú)法從已有數(shù)據(jù)中獲取。針對(duì)這一問(wèn)題，本文采用模糊層次分析法[3]，計(jì)算六類渠道的風(fēng)險(xiǎn)權(quán)重。首先我們構(gòu)建層次結(jié)構(gòu)圖，如下所示（相關(guān)符號(hào)定義已在前文給出）:圖2.1食品安全風(fēng)險(xiǎn)層次結(jié)構(gòu)圖不安全性評(píng)價(jià)指標(biāo)（）不安全性評(píng)價(jià)指標(biāo)（）農(nóng)藥殘留（）新鮮程度（）非法添加（）農(nóng)貿(mào)市場(chǎng)安全風(fēng)險(xiǎn)批發(fā)市場(chǎng)安全風(fēng)險(xiǎn)商場(chǎng)超市安全風(fēng)險(xiǎn)生產(chǎn)基地安全風(fēng)險(xiǎn)零售安全風(fēng)險(xiǎn)電商安全風(fēng)險(xiǎn)目標(biāo)層準(zhǔn)則層方案層我們將影響各渠道安全風(fēng)險(xiǎn)的因素分為農(nóng)藥殘留（）、新鮮程度（）和非法添加（）三大類。同時(shí)，我們將六大銷售渠道的安全風(fēng)險(xiǎn)分為農(nóng)貿(mào)市場(chǎng)安全風(fēng)險(xiǎn)（）、批發(fā)市場(chǎng)安全風(fēng)險(xiǎn)（）、商場(chǎng)超市安全風(fēng)險(xiǎn)（）、生產(chǎn)基地安全風(fēng)險(xiǎn)（）、零售安全風(fēng)險(xiǎn)（）和電商安全風(fēng)險(xiǎn)（）。而目標(biāo)則為生鮮蔬菜不安全性的評(píng)價(jià)指標(biāo)（）。設(shè)為各層各元素之間優(yōu)先關(guān)系矩陣。其中的取值規(guī)范見(jiàn)下表：表2.2優(yōu)先關(guān)系矩陣元素標(biāo)度含義標(biāo)度含義0因素不如因素重要0.5因素與因素同樣重要1因素比因素重要顯然，優(yōu)先關(guān)系矩陣為模糊互補(bǔ)陣。我們對(duì)其按行求和，求得以下數(shù)組，即為。設(shè)為模糊一致陣。以為元素構(gòu)造，通過(guò)算式實(shí)現(xiàn)。設(shè)為每一層各元素相對(duì)上一層的重要度，即權(quán)重。我們通過(guò)計(jì)算矩陣各行元素之和與所有元素和相除的商，對(duì)向量進(jìn)行賦值：。分別輸入層、層優(yōu)先關(guān)系矩陣如下：表2.3層優(yōu)先關(guān)系矩陣0.51000.50110.5表2.4層優(yōu)先關(guān)系矩陣0.51011000.50100110.51110000.50001010.50110110.5、層同理，此處不予列出。利用上述公式，我們對(duì)4個(gè)矩陣分別進(jìn)行模糊一致化、行歸一化處理，得到權(quán)重向量、、、。最后，我們通過(guò)算得層元素對(duì)層的權(quán)重向量：。將此結(jié)果在EXCEL直觀呈現(xiàn)如下：表2.5六類渠道安全風(fēng)險(xiǎn)權(quán)重表農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地零售電商權(quán)重0.15590.13120.22990.09720.20220.18365.2.4基于函數(shù)擬合的不合格率計(jì)算結(jié)合表2.1，將上述權(quán)重和前四類渠道的不合格率在折線圖中描點(diǎn)，并擬合為一次函數(shù)，見(jiàn)下圖：圖2.2四類渠道安全風(fēng)險(xiǎn)權(quán)重與不合格率相關(guān)曲線根據(jù)曲線方程，我們將零售渠道的風(fēng)險(xiǎn)權(quán)重（）和電商渠道的風(fēng)險(xiǎn)權(quán)重（）分別代入，求得二者不合格率和。至此，我們已獲得全部所需數(shù)據(jù)。5.2.5忽略特定渠道的安全風(fēng)險(xiǎn)計(jì)算模型基于之前所有準(zhǔn)備工作，我們現(xiàn)在可以直接建立模型求解問(wèn)題二。根據(jù)六大新鮮蔬菜銷售渠道各自的不合格率，和表1.7提供的渠道比例，可求得平均不合格率?？紤]到刪除某一渠道的抽樣，將會(huì)導(dǎo)致平均不合格率或高或低的改變，產(chǎn)生新的。而改變值能夠反映該渠道的影響力大小。（1）當(dāng)時(shí)，該渠道樣本的缺失將導(dǎo)致平均不合格率的升高。越大，說(shuō)明被忽略的不合格樣本越多，忽略該渠道導(dǎo)致的安全風(fēng)險(xiǎn)越大。（2）當(dāng)時(shí)，該渠道樣本的缺失將導(dǎo)致平均不合格率的降低。越小，說(shuō)明該渠道本身的不合格樣本越多，忽略該渠道導(dǎo)致的安全風(fēng)險(xiǎn)將顯著降低。我們將全部信息整合如下：表2.6六類渠道綜合信息和可忽略性特征變量農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地零售電商均值比例0.470.3260.0890.0440.0680.003--不合格率(%)5.4502.7508.6301.4807.2706.2404.804權(quán)重0.15590.13120.22990.09720.20220.1836--剩余渠道平均不合格率(%)4.2325.7984.4314.9574.6244.800--平均不合格率改變值(%)-0.5730.994-0.3740.153-0.180-0.004--由上表知，忽略農(nóng)貿(mào)市場(chǎng)樣本所減少的平均不合格率最多，為，因此，忽略農(nóng)貿(mào)市場(chǎng)渠道導(dǎo)致的安全風(fēng)險(xiǎn)最大。商場(chǎng)超市次之，為。而批發(fā)市場(chǎng)和生產(chǎn)基地的，說(shuō)明這兩個(gè)渠道本身的不合格率很低，將這它們忽略是可行的。因此，為了使蔬菜安全得到更優(yōu)保障，應(yīng)當(dāng)適當(dāng)增加農(nóng)貿(mào)市場(chǎng)和商場(chǎng)超市的樣本容量，并減少批發(fā)市場(chǎng)和生產(chǎn)基地所抽樣本容量。5.2.5.2故障樹分析法的不同渠道相對(duì)風(fēng)險(xiǎn)排序在進(jìn)行新鮮果蔬的渠道風(fēng)險(xiǎn)比較時(shí)，主要考慮的因素為蔬菜在不同渠道的供應(yīng)鏈中受到的二次污染。我們將不同渠道的蔬菜發(fā)生的安全風(fēng)險(xiǎn)的概率用故障樹法進(jìn)行表示，并進(jìn)行比較。故障樹如圖所示。圖2.3故障樹導(dǎo)致蔬菜質(zhì)量風(fēng)險(xiǎn)的因素分為三種：原產(chǎn)品質(zhì)量風(fēng)險(xiǎn)運(yùn)輸階段風(fēng)險(xiǎn)出售階段風(fēng)險(xiǎn)其中運(yùn)輸階段風(fēng)險(xiǎn)受到運(yùn)輸時(shí)長(zhǎng)以及是否采用冷鏈物流的影響，而出售階段風(fēng)險(xiǎn)主要受到蔬菜從存儲(chǔ)到出售的時(shí)間的影響。其中，我們假定在蔬菜生產(chǎn)環(huán)節(jié)中，各個(gè)銷售渠道各自獲得的蔬菜的品質(zhì)大體一致，且由于深圳市蔬菜的進(jìn)口特性，不同渠道的獲取蔬菜的原產(chǎn)品質(zhì)存在一定的差異。在考慮蔬菜運(yùn)輸環(huán)節(jié)時(shí)，我們認(rèn)為新鮮果蔬在運(yùn)輸過(guò)程中會(huì)受溫度、相對(duì)濕度、運(yùn)輸時(shí)間等因素影響，致其新鮮程降低，甚至?xí)霈F(xiàn)果蔬腐爛變質(zhì)的問(wèn)題，這對(duì)于果蔬的安全風(fēng)險(xiǎn)有較大的影響。本文中不同渠道運(yùn)輸?shù)臏囟纫约跋鄬?duì)濕度都大致相同，而在蔬菜的運(yùn)輸流通環(huán)節(jié)中，運(yùn)輸?shù)臅r(shí)間越長(zhǎng)，蔬菜變質(zhì)的安全風(fēng)險(xiǎn)越高。與此同時(shí)，在流通方式的選擇方面，不同的流通方式也會(huì)較大程度的影響蔬菜的品質(zhì)，是否采用冷鏈物流對(duì)于蔬菜的質(zhì)量有著較為直接的影響。而在不同銷售渠道的蔬菜存儲(chǔ)的環(huán)節(jié)當(dāng)中，我們將蔬菜在不同運(yùn)輸渠道在運(yùn)輸環(huán)節(jié)當(dāng)中的質(zhì)量變化以細(xì)菌和微生物的增長(zhǎng)水平來(lái)表示。細(xì)菌和微生物的增長(zhǎng)一般以分為四個(gè)階段：遲緩期、對(duì)數(shù)期、穩(wěn)定期、衰亡期。蔬菜在運(yùn)輸過(guò)程中，我們假定細(xì)菌一直處于對(duì)數(shù)期階段。在對(duì)數(shù)期間，細(xì)菌的增長(zhǎng)隨時(shí)間成指數(shù)函數(shù)變化，時(shí)間運(yùn)輸時(shí)間越長(zhǎng)對(duì)于蔬菜的影響越大。而冷鏈物流（低溫運(yùn)輸）在一定程度上抑制了細(xì)菌的分裂。所以蔬菜在運(yùn)輸過(guò)程中的風(fēng)險(xiǎn)主要由以上兩點(diǎn)構(gòu)成（1.時(shí)間；2.是否冷鏈）。圖2.4蔬菜運(yùn)輸方式我國(guó)蔬菜的運(yùn)輸方式主要分為以下三類，自營(yíng)配送，共同配送和第三方配送。其中自營(yíng)配送是一種傳統(tǒng)類型的蔬菜配送方式，零售點(diǎn)和農(nóng)貿(mào)市場(chǎng)較多采用這種方式進(jìn)行蔬菜配送。這種模式反應(yīng)較為快速靈活，農(nóng)戶可以較好地控制配送活動(dòng)，但是其運(yùn)送量一般較小且較為分散，運(yùn)輸方式也較為不統(tǒng)一，存在管理上的缺陷。與此同時(shí)，這種配送方式較為簡(jiǎn)陋，一般不采用冷鏈物流，其過(guò)程中的貨物質(zhì)量也不能得到較好地保證。共同配送主要是由以一家或多家蔬菜生產(chǎn)基地或蔬菜批發(fā)市場(chǎng)為中心組成的新配送公司，或從不同的地方分散集貨，統(tǒng)一對(duì)多家用戶進(jìn)行配送，按照用戶的需求將貨物配送到各個(gè)指定的接貨點(diǎn)。共同配送可以提高配送效率和運(yùn)送的裝載率，但與此同時(shí)，由于配送由多個(gè)主體共同結(jié)合，加大了管理的難度，難免會(huì)出現(xiàn)時(shí)間延長(zhǎng)等情況。蔬菜的第三方配送主要是由超市以及電商把自己需要完成的配送業(yè)務(wù)通過(guò)專業(yè)的第三方配送機(jī)構(gòu)進(jìn)行配送的運(yùn)作模式。在這種高度專業(yè)的配送服務(wù)能夠較好地降低配送損耗，多采用冷鏈配送的方式來(lái)保證蔬菜的質(zhì)量。在存儲(chǔ)時(shí)長(zhǎng)方面，由于零售和農(nóng)貿(mào)市場(chǎng)、生產(chǎn)基地、批發(fā)市場(chǎng)的蔬菜采購(gòu)量和流動(dòng)量都較大，蔬菜一般較早就被售出，儲(chǔ)存時(shí)間一般都較短。而在超市中，蔬菜的銷售一般較慢，且會(huì)在冷庫(kù)中長(zhǎng)期保存，儲(chǔ)存時(shí)間較長(zhǎng)。而電商銷售的蔬菜的儲(chǔ)存時(shí)間則更加長(zhǎng)。由于底事件中各單元的不合格率根據(jù)現(xiàn)有的數(shù)據(jù)無(wú)法進(jìn)行具體的表示，我們只能根據(jù)風(fēng)險(xiǎn)樹法中底事件中估計(jì)的相對(duì)不合格率對(duì)于不同渠道的安全風(fēng)險(xiǎn)進(jìn)行定性的相對(duì)比較。根據(jù)上文中的分析，我們得出了以下底事件概率。進(jìn)貨質(zhì)量不良（X1）沒(méi)有進(jìn)行冷鏈運(yùn)輸（X2）運(yùn)輸時(shí)長(zhǎng)較長(zhǎng)(X3)銷售滯留（X4）結(jié)果農(nóng)貿(mào)市場(chǎng)0.1331.0000.0630.0670.2625批發(fā)市場(chǎng)0.1330.0000.1880.0670.2商場(chǎng)超市0.2000.0000.1880.2000.4生產(chǎn)基地0.0670.0000.2500.2000.266667零售0.2001.0000.1250.2000.525電商0.2670.0000.1880.2670.533333運(yùn)用故障樹算法對(duì)于頂事件概率進(jìn)行計(jì)算，得到以上結(jié)果，我們簡(jiǎn)單做一總結(jié)：風(fēng)險(xiǎn)程度：電商＞零售＞商場(chǎng)超市＞生產(chǎn)基地＞農(nóng)貿(mào)市場(chǎng)＞批發(fā)市場(chǎng)。而在不合格率分析中，我們得到的結(jié)果是：風(fēng)險(xiǎn)程度：商場(chǎng)超市＞零售＞電商＞農(nóng)貿(mào)市場(chǎng)＞批發(fā)市場(chǎng)＞生產(chǎn)基地。兩種方法所得結(jié)果略有不同，可能的原因如下：首先，合格率與風(fēng)險(xiǎn)的定義略有差異，本方法通過(guò)故障樹計(jì)算的的蔬菜風(fēng)險(xiǎn)概率描述的是蔬菜在購(gòu)買、運(yùn)輸與存儲(chǔ)三個(gè)過(guò)程當(dāng)中蔬菜存在食品安全風(fēng)險(xiǎn)的相對(duì)排序，安全風(fēng)險(xiǎn)并不等同于合格率，部分蔬菜發(fā)生安全風(fēng)險(xiǎn)，但仍然處于合格狀態(tài)。其次，生產(chǎn)基地和電商的風(fēng)險(xiǎn)排名要高于其不合格率排名。由于生產(chǎn)基地運(yùn)輸時(shí)長(zhǎng)較高，所以故障樹分析中存在較大的比例份額。而實(shí)際中，生產(chǎn)基地多采用大批量的運(yùn)送方法，在運(yùn)輸過(guò)程中由于擠壓更容易發(fā)生變質(zhì)。而電商作為一種新型的方式在管理上還存在一定的不足。最后，故障樹分析法由于缺乏相對(duì)準(zhǔn)確的數(shù)據(jù)，我們采用了不同蔬菜流通渠道對(duì)于故障樹的底事件的相對(duì)風(fēng)險(xiǎn)排名來(lái)確定各渠道之間的相對(duì)風(fēng)險(xiǎn)進(jìn)行輔助評(píng)價(jià)排序，在流通環(huán)節(jié)中，還存在售賣方式不同（如按斤出售或按箱出售）等因素沒(méi)有考慮，所以結(jié)論也存在相應(yīng)的局限性。至此完成問(wèn)題二的解答。5.3問(wèn)題三的模型建立與求解5.3.1問(wèn)題三的分析本題要求我們基于附件評(píng)估深圳新鮮蔬菜食品安全風(fēng)險(xiǎn)。附件一為2017年深圳市食用農(nóng)產(chǎn)品質(zhì)量安全監(jiān)測(cè)抽檢方案，附件二為深圳1-3月的分區(qū)分渠道的抽檢計(jì)劃，附件三為全國(guó)不合格蔬菜及其制品的匯總報(bào)表，其中包涵了不同的地區(qū)和鮮菜樣品的詳細(xì)數(shù)據(jù)。針對(duì)這一問(wèn)題，我們給出如下求解思路和具體步驟：步驟一：首先我們對(duì)于深圳市在風(fēng)險(xiǎn)物質(zhì)方面存在的食品安全風(fēng)險(xiǎn)進(jìn)行分析。根據(jù)附件三中全國(guó)不合格蔬菜及其制品的匯總報(bào)表以及附件一中深圳市食用農(nóng)產(chǎn)品質(zhì)量安全監(jiān)測(cè)抽檢方案，對(duì)于鮮菜制品中的抽樣項(xiàng)目進(jìn)行分類，將鮮菜的農(nóng)產(chǎn)品安全檢測(cè)物質(zhì)分為兩種：1.抽樣方案中未提及的風(fēng)險(xiǎn)物質(zhì)，2.抽樣報(bào)告中已提及的風(fēng)險(xiǎn)物質(zhì)。對(duì)于抽樣方案中未提及的風(fēng)險(xiǎn)物質(zhì)，我們首先采用定性分析，對(duì)于分析的風(fēng)險(xiǎn)物質(zhì)進(jìn)行初步的定性，接著使用多維灰聚類法對(duì)于風(fēng)險(xiǎn)物質(zhì)進(jìn)行定量的分析，得到深圳市在抽樣方案中未提及的風(fēng)險(xiǎn)物質(zhì)的缺失檢驗(yàn)風(fēng)險(xiǎn)。對(duì)于抽樣報(bào)告中已提及的風(fēng)險(xiǎn)物質(zhì)，我們也要考慮到漏檢和缺檢對(duì)于深圳市食品安全的影響程度，分類匯總之后得到深圳市檢驗(yàn)得樣品鮮菜當(dāng)中主要存在的不合格項(xiàng)目以及超標(biāo)情況較為嚴(yán)重的不合格項(xiàng)目，提出在可能存在的漏檢或缺檢的情況下對(duì)于深圳市食品安全的影響。步驟二：接著我們對(duì)于深圳市在選擇鮮菜進(jìn)口產(chǎn)地方面的食品安全風(fēng)險(xiǎn)進(jìn)行分析。首先根據(jù)附件三中全國(guó)不合格蔬菜及其制品的匯總報(bào)表結(jié)合中國(guó)農(nóng)村統(tǒng)計(jì)年鑒2015對(duì)于鮮菜安全程度省排名進(jìn)行確定，并使用K-聚類分析法對(duì)深圳市從不同省份采購(gòu)的鮮菜量的估計(jì)比較將全國(guó)31個(gè)可能的蔬菜出口省份進(jìn)行聚類，根據(jù)以上情況結(jié)合風(fēng)險(xiǎn)的比率提出深圳市現(xiàn)行進(jìn)口產(chǎn)地對(duì)于深圳市食品安全可能帶來(lái)的風(fēng)險(xiǎn)。步驟三：最后我們對(duì)于深圳市在選擇鮮菜種類方面的食品安全風(fēng)險(xiǎn)進(jìn)行分析。首先經(jīng)過(guò)調(diào)查對(duì)于反季節(jié)蔬菜進(jìn)行分類，分析了反季節(jié)蔬菜對(duì)于深圳市食品安全的影響。接著結(jié)合附件數(shù)據(jù)對(duì)于深圳市的市民主要喜愛(ài)的蔬菜種類進(jìn)行分析，結(jié)合地域特性進(jìn)行風(fēng)險(xiǎn)分析。5.3.2鮮菜中存在的危險(xiǎn)物質(zhì)帶來(lái)的物質(zhì)性風(fēng)險(xiǎn)5.3.2.1對(duì)于抽樣標(biāo)準(zhǔn)中沒(méi)有提及的有害物質(zhì)在附件一的抽樣方法的說(shuō)明當(dāng)中，提到了蔬菜、水果和食用菌抽樣按《農(nóng)藥殘留分析樣本的采樣方法》（NY/T789-2004）規(guī)定執(zhí)行。該標(biāo)準(zhǔn)中，檢測(cè)了包括滅幼脲、多菌靈、吡蟲啉、甲氨基阿維菌素苯甲酸鹽、烯酰嗎啉、蟲螨腈、咪鮮胺、嘧菌酯、二甲戊樂(lè)靈、噻蟲嗪、氟啶脲等58種農(nóng)藥殘留在內(nèi)的多少種農(nóng)藥殘留，然而在全國(guó)的蔬菜抽樣調(diào)查當(dāng)中，檢查出了除了抽樣調(diào)查之外的非法添加、污染物等危險(xiǎn)物質(zhì)?？紤]到未檢查物質(zhì)的風(fēng)險(xiǎn)性對(duì)于深圳市鮮菜質(zhì)量的潛在風(fēng)險(xiǎn)，我們使用定性分析以及多維灰聚類法對(duì)于未檢測(cè)組織的風(fēng)險(xiǎn)程度進(jìn)行評(píng)估。一、風(fēng)險(xiǎn)評(píng)估中的定性分析方法下面對(duì)鮮菜制品中風(fēng)險(xiǎn)較大的風(fēng)險(xiǎn)項(xiàng)目進(jìn)行進(jìn)行定性分析：非法添加新鮮蔬菜中的非法添加物質(zhì)主要是熒光增白劑以及6-芐基腺嘌呤。熒光增白劑是一類化學(xué)物質(zhì)，它的主要作用是對(duì)蔬菜進(jìn)行增白，使用熒光增白劑加入食品如蘑菇、面粉、腐竹中可以掩蓋其發(fā)黃及不新鮮的品相。由于熒光增白物質(zhì)能夠促進(jìn)細(xì)胞對(duì)于放射線的吸收、促進(jìn)細(xì)胞產(chǎn)生變異。所以熒光增白劑被列為潛在的致癌因素之一。6-芐基腺嘌呤即植物細(xì)胞分裂素，主要用于促進(jìn)細(xì)胞分裂、促進(jìn)果實(shí)生長(zhǎng)、促進(jìn)細(xì)胞增大、促進(jìn)種子發(fā)芽等。在鮮菜種植中主要被用于豆芽種植。有調(diào)查顯示長(zhǎng)期食用添加了6-芐基腺嘌呤的食品，會(huì)損害人類的生理功能以及內(nèi)分泌系統(tǒng)，因此國(guó)家禁止在生產(chǎn)中使用6-芐基腺嘌呤。微生物蔬菜中的微生物超標(biāo)主要以大腸桿菌為主。大腸桿菌一般作為腸道菌的正常組成部分，但存在一部分的致病大腸桿菌，對(duì)人類的消化系統(tǒng)帶來(lái)強(qiáng)烈的傷害。在距離養(yǎng)牛場(chǎng)較近的蔬菜種植場(chǎng)中的蔬菜容易受到大腸桿菌的污染。污染物蔬菜中的污染物主要為重金屬，主要以汞、鎘、鉛為主。蔬菜中的重金屬主要來(lái)源于農(nóng)藥的殘留以及種植土地的污染。重金屬在人體內(nèi)具有累積效應(yīng)，當(dāng)重金屬累積到一定程度之后，會(huì)引發(fā)人體的一系列疾病，所以國(guó)家對(duì)于蔬菜中的重金屬含量有著嚴(yán)格的規(guī)定。其他物質(zhì)蔬菜中還存在一些其他物質(zhì)例如生物毒素赭曲霉毒素A以及亞硝酸鹽等。這些物質(zhì)通常是由于蔬菜制品存儲(chǔ)不當(dāng)引起。超標(biāo)少量即對(duì)人體會(huì)造成極大的影響，但一般樣品出現(xiàn)較少，這里不做過(guò)多討論。二、風(fēng)險(xiǎn)評(píng)估中的多維灰聚類法[4]影響鮮菜食品安全的未檢測(cè)因素有很多，其中含有大量的未知信息，具有一定的灰性。在缺少樣本數(shù)據(jù)的情況下，為了最大程度實(shí)現(xiàn)對(duì)于因素的量化與白化，我們選用多維灰度聚類法對(duì)于數(shù)據(jù)進(jìn)行統(tǒng)計(jì)與分析。設(shè)定一定的評(píng)判標(biāo)準(zhǔn)，將樣本分為高中低三類進(jìn)行聚類評(píng)估。最終得到各個(gè)樣本的權(quán)系數(shù)，對(duì)樣本得出相應(yīng)的灰度判定。以下為其基本原理：設(shè)有個(gè)評(píng)估樣本，項(xiàng)指標(biāo)，個(gè)不同的判別灰類，樣本關(guān)于指標(biāo)的量化評(píng)價(jià)為，指標(biāo)在灰類的白化權(quán)函數(shù)為，為指標(biāo)對(duì)評(píng)估目標(biāo)所占重，且，因此，樣本在灰類的綜合權(quán)系數(shù)為（１）稱為樣本綜合權(quán)系數(shù)矩陣的行向量，由的最大值確定樣本的所屬灰類。構(gòu)建歸一化矩陣對(duì)于不同的物質(zhì)以及不同的國(guó)家要求標(biāo)準(zhǔn)，我們應(yīng)采用不同的方法對(duì)于數(shù)據(jù)進(jìn)行歸一化。對(duì)于非法添加物質(zhì)，我們采用檢測(cè)值減去限量值的方法，對(duì)于限量為不得出現(xiàn)取零。對(duì)于微生物、污染物、農(nóng)殘以及其他物質(zhì)采用檢測(cè)物質(zhì)超標(biāo)倍數(shù)的方法進(jìn)行歸一化，得到的歸一化矩陣如下表（部分）。表3.1超標(biāo)物質(zhì)歸一化矩陣食品添加劑污染物非法添加農(nóng)殘微生物1.001.000.081.001.001.001.000.071.001.001.001.000.041.001.001.001.000.751.001.001.001.000.571.001.001.001.000.021.001.001.001.000.061.001.001.001.000.071.001.001.001.000.201.001.001.001.000.031.001.001.001.000.021.001.001.001.000.021.001.00確定參數(shù)的正負(fù)極性這5個(gè)參數(shù)數(shù)值隨現(xiàn)場(chǎng)的安全風(fēng)險(xiǎn)程度成正比，是正極指標(biāo)。確定風(fēng)險(xiǎn)程度區(qū)間表3.2超標(biāo)物質(zhì)風(fēng)險(xiǎn)程度區(qū)間考慮因數(shù)食品添加劑污染物非法添加農(nóng)殘微生物150.004.000.2050.00160.007.002.000.065.003.531.001.00011主要根據(jù)國(guó)家現(xiàn)行的食品安全質(zhì)量評(píng)估標(biāo)注確定。構(gòu)造白化權(quán)函數(shù)及權(quán)系數(shù)（以污染物為例）高風(fēng)險(xiǎn)：f中風(fēng)險(xiǎn)：f低風(fēng)險(xiǎn)：f其他物質(zhì)與此同理。計(jì)算權(quán)系數(shù)矩陣將5個(gè)參數(shù)處理后的數(shù)值帶入公式，可得各自權(quán)系數(shù)矩陣。。確定權(quán)重值根據(jù)以上定性分析中所提及的鮮菜樣品中的風(fēng)險(xiǎn)內(nèi)容，我們將鮮菜作物中的相關(guān)參數(shù)劃分不同的等級(jí)，并認(rèn)為此五項(xiàng)參數(shù)是造成鮮菜質(zhì)量不合格的全部原因。將3個(gè)權(quán)系數(shù)矩陣的對(duì)應(yīng)值進(jìn)行加權(quán)平均，分別側(cè)重考慮了非法添加、污染物（重金屬元素）以及微生物和其他物質(zhì)得到了不合格鮮菜樣品的綜合權(quán)系數(shù)矩陣。由綜合權(quán)系數(shù)矩陣中各不合格樣品行向量的最大值，判斷各樣品所屬灰類?；翌惙植紙D如下圖所示。并列出具有較高風(fēng)險(xiǎn)的鮮菜的詳細(xì)信息表。表3.3分別側(cè)重考慮非法添加、污染物（重金屬元素）以及微生物和其他物質(zhì)的灰類分布圖表3.4中風(fēng)險(xiǎn)以及高風(fēng)險(xiǎn)蔬菜的危險(xiǎn)物質(zhì)表（部分）蔬菜名稱不合格項(xiàng)目分類不合格項(xiàng)目檢測(cè)數(shù)值法規(guī)限制風(fēng)險(xiǎn)等級(jí)小豆芽非法添加6-芐基腺嘌呤（6-BA）0.75mg/kg0.0（mg/kg）高竹蓀污染物鎘2.1mg/kg≤0.5mg/kg高莧菜污染物鎘0.38mg/kg≤0.2mg/kg高蔬菜（生姜）污染物鉛1.4mg/kg≤0.1mg/kg高木耳菜污染物鎘0.33mg/kg≤0.2mg/kg高韭菜污染物汞0.034mg/kg≤0.01mg/kg高馬鈴薯污染物汞0.039mg/kg≤0.01mg/kg高飯豆(雀蛋豆)生物毒素赭曲霉毒素A14.3μg/kg≤5.0μg/kg高黃豆芽非法添加6-芐基腺嘌呤（6-BA）0.2mg/kg不得檢出中黃豆芽非法添加6-芐基腺嘌呤0.071mg/kg不得檢出中茄子污染物鎘（以Cd計(jì)）0.2mg/kg≤0.05mg/kg中雙孢菇污染物總汞0.2mg/kg≤0.1mg/kg中嫩生姜污染物亞硝酸鹽(以NaNO2計(jì))39.4mg/kg≤20mg/kg中結(jié)論：深圳市中主要存在的危險(xiǎn)物質(zhì)為非法添加中的6-芐基腺嘌呤，以及污染物中的鎳、鎘、汞等重金屬物質(zhì)，微生物中的大腸桿菌，生物毒素中的赭曲霉毒素A，以及污染物中的亞硝酸鹽，這些物質(zhì)成分由于不屬于深圳市食用農(nóng)產(chǎn)品質(zhì)量安全監(jiān)測(cè)抽檢方案中，且其對(duì)于該市的鮮菜安全存在較大的影響。5.3.2.2對(duì)于抽樣標(biāo)準(zhǔn)中已提及的有害物質(zhì)一、被檢出頻率較高的物質(zhì)表3.5不合格項(xiàng)目判定結(jié)果不合格項(xiàng)目頻率不合格農(nóng)殘硫丹11水胺硫磷4辛硫磷11對(duì)于抽樣標(biāo)準(zhǔn)中已提及的有害物質(zhì)，我們首先對(duì)于其被檢測(cè)的最終頻數(shù)按大小進(jìn)行分類，選取了其中最高的三類物質(zhì)，分別為硫丹、水胺硫磷與辛硫磷，在全國(guó)此類農(nóng)藥殘留超標(biāo)較為嚴(yán)重的情況下，深圳市同樣面臨著這一類物質(zhì)的安全影響。在某類物質(zhì)大量的超標(biāo)情況下，容易出現(xiàn)對(duì)含該類風(fēng)險(xiǎn)物質(zhì)漏檢的情況。二、毒性較高的幾類農(nóng)殘表3.6（LD50）毒性評(píng)級(jí)[5]危險(xiǎn)物質(zhì)毒性評(píng)級(jí)（LD50）危險(xiǎn)物質(zhì)毒性評(píng)級(jí)（LD50）苯醚甲環(huán)唑低毒氧樂(lè)果中毒多菌靈微毒乙酰甲胺磷低毒二氧化硫低毒毒死蜱低毒氟蟲腈高毒代森錳鋅微毒克百威劇毒滴滴涕中毒樂(lè)果中毒敵敵畏高毒硫丹劇毒啶蟲脒中毒氯氰菊酯高毒腐霉利微毒水胺硫磷劇毒甲氨基阿維菌素苯甲酸鹽低毒涕滅威中毒甲胺磷劇毒辛硫磷低毒甲基異柳磷劇毒滅蠅胺低毒附件三的全國(guó)不合格蔬菜及其制品的匯總報(bào)表中，鮮菜被檢驗(yàn)出的危險(xiǎn)物質(zhì)中存在克百威、硫丹、水胺硫磷、甲胺磷、甲基異柳磷等幾類劇毒物質(zhì)以及氟蟲腈等高毒性農(nóng)藥，國(guó)家的標(biāo)準(zhǔn)中對(duì)這幾類農(nóng)藥有著嚴(yán)格的控制，劇毒和高毒農(nóng)藥殘留的含量被嚴(yán)格地控制。由于深圳市鮮菜的輸入性特性，無(wú)法對(duì)原產(chǎn)地的種植情況進(jìn)行檢查，而一旦出現(xiàn)漏檢劇毒農(nóng)藥輸入鮮菜的情況，后果將是比較嚴(yán)重的。三、超標(biāo)程度較高的幾類農(nóng)殘?jiān)诓缓细駱颖局?，存在超?biāo)上千倍的重金屬以及農(nóng)藥殘留硫丹的情況，對(duì)于鮮菜產(chǎn)品主要采用抽查的方法進(jìn)行檢查，而蔬菜等農(nóng)產(chǎn)品極易出現(xiàn)檢查成分波動(dòng)的情況，在沒(méi)有制定合理的抽樣方法的前提下，深圳市仍然具有被高超標(biāo)鮮菜制品影響的潛在威脅。5.3.3生產(chǎn)管理不規(guī)范或污染嚴(yán)重區(qū)域生產(chǎn)的鮮菜產(chǎn)品帶來(lái)的地域性風(fēng)險(xiǎn)5.3.3.1鮮菜安全程度省排名的確定根據(jù)題目所提供的附件3之中的全國(guó)的抽樣檢測(cè)的不合格報(bào)告數(shù)據(jù)，我們選擇了其中的鮮菜有關(guān)的數(shù)據(jù)進(jìn)行了重點(diǎn)分析，經(jīng)過(guò)分類統(tǒng)計(jì)得到了其中的各個(gè)省份的被抽樣的不合格個(gè)數(shù)以及各個(gè)省份的蔬菜年產(chǎn)量。[6]用不合格頻數(shù)除以蔬菜量估計(jì)出各省市的不合格鮮菜產(chǎn)品占全部蔬菜產(chǎn)量的比值，并根據(jù)該數(shù)據(jù)從小到大得到大致的蔬菜合格率的省排名。表3.7根據(jù)鮮菜不合格率得出的安全程度的省排名（由低風(fēng)險(xiǎn)到高風(fēng)險(xiǎn)）省份不合格頻數(shù)百分比蔬菜年產(chǎn)量排名安徽41.2908369北京10.3236216福建165.16170618甘肅123.871481.920廣東92.9031306.911廣西51.6131104.68貴州309.677847.723海南61.935239.421河北72.2581220.410河南41.2901745.83黑龍江30.968265.713湖北41.29011454湖南165.1611283.714吉林134.194214.627江蘇123.8711354.912江西20.645563.35遼寧72.258492.115內(nèi)蒙古165.161265.726寧夏6721.613117.331青海72.25850.530山東20.6451832.91山西268.387252.829陜西103.22649017上海30.968132.219四川51.61312766天津41.29089.925西藏10.32323.924新疆185.806296.728云南41.290900.87浙江10.323619.12重慶185.806681.822合計(jì)3101005.3.3.2深圳市從不同省份采購(gòu)的鮮菜量的估計(jì)比較對(duì)于不同的省份，我們綜合考慮了一個(gè)省份的出口蔬菜量、鐵路距離、該地運(yùn)送往深圳的工作日、出口的蔬菜種類的種類系數(shù)、該地運(yùn)輸需要的估計(jì)工作日、該地主要的生產(chǎn)時(shí)間、是否有大量其他糧食輸出深圳記錄等因素對(duì)于全國(guó)可能的蔬菜蔬菜區(qū)域進(jìn)行分類。（忽略臺(tái)灣、香港、澳門）表3.7數(shù)據(jù)的詳細(xì)說(shuō)明表列1說(shuō)明單位省份的出口蔬菜量以中國(guó)農(nóng)村統(tǒng)計(jì)年鑒2015的數(shù)據(jù)估計(jì)噸/每年鐵路距離以最近鐵路距離估計(jì)（到省會(huì)）千米工作日以正常季節(jié)正常速度速度估計(jì)天種類得分將該地區(qū)的是否主要生產(chǎn)某類產(chǎn)品的01矩陣與深圳市對(duì)各類產(chǎn)品的需求比例的列向量相乘得到的結(jié)果數(shù)值生產(chǎn)季節(jié)主要分為冬春季節(jié)以及夏秋季節(jié)分類單位是否對(duì)深圳大量出口其他產(chǎn)品以其他糧食制品對(duì)深圳的出口程度進(jìn)行模擬擬合0="無(wú)出口"，1="少量出口"，2="正常出口"，3="大量出口"在對(duì)于影響因素進(jìn)行處理之后，我們運(yùn)用spss軟件對(duì)31個(gè)考慮地區(qū)進(jìn)行分類K均值聚類分析，k均值聚類分析的過(guò)程如下：1、隨機(jī)選取k個(gè)聚類質(zhì)心點(diǎn)（clustercentroids）為。2、重復(fù)下面過(guò)程直到收斂對(duì)于每一個(gè)樣例i，計(jì)算其應(yīng)該屬于的類

對(duì)于每一個(gè)類j，重新計(jì)算該類的質(zhì)心

根據(jù)以上原理以及影響因子的相關(guān)性分析，我們將影響因子的程度的重要程度按照運(yùn)價(jià)、蔬菜產(chǎn)量、輸出需要工作日、種類得分、是否有大量蔬菜輸出深圳記錄、季節(jié)性的順序排序。將影響因素輸入spss程序后，對(duì)聚類的分類個(gè)數(shù)進(jìn)行多次嘗試，經(jīng)過(guò)多次嘗試后將我們的聚類個(gè)數(shù)分為三類。此時(shí)得到模擬的結(jié)果較為滿意，與文獻(xiàn)中提到的主要場(chǎng)地基本相同。其中聚類類別1中的各省份與深圳市的聯(lián)系密切程度一般，預(yù)計(jì)出口程度一般，省份主要為北部以及西北部省份，估計(jì)影響因素主要為地域原因帶來(lái)的運(yùn)輸成本以及運(yùn)輸時(shí)間的帶來(lái)的可能產(chǎn)品變質(zhì)。聚類類別2中的各省份預(yù)測(cè)出口程度較高，主要為南方的近距離的相似氣候以及擁有相近的蔬菜作物的省份以及高產(chǎn)量、低運(yùn)價(jià)輸出的部分北方省份。聚類類別3中的各省份預(yù)計(jì)出口程度較低，主要由于地處遙遠(yuǎn)造成的高額的運(yùn)價(jià)或過(guò)低的產(chǎn)量帶來(lái)的相關(guān)影響。以下為迭代記錄以及分類：表3.8K-均值聚類算法的迭代記錄表迭代歷史記錄迭代聚類中心內(nèi)的更改12311129.654782.913983.196210.16673.035150.2053175.62375.995276.446442.19375.995.0005.000.000.000表3.9基于多種因素的深圳市與各個(gè)省份之間的蔬菜預(yù)計(jì)出口程度聚類結(jié)果類別號(hào)列2列3列4類別號(hào)頻率頻率百分比頻率百分比1北京18.3332安徽16.667甘肅18.333福建16.667遼寧18.333廣東16.667內(nèi)蒙古18.333廣西16.667寧夏18.333貴州16.667青海18.333海南16.667山西18.333河北16.667陜西18.333河南16.700上海18.333湖北16.667四川18.333湖南16.667天津18.333江蘇16.667重慶18.333江西16.667合計(jì)12100.000山東16.6673黑龍江125.000云南16.667吉林125.000浙江16.667西藏125.000合計(jì)15100.000新疆125.000合計(jì)41005.3.3.3省份產(chǎn)品的威脅在鮮菜食品質(zhì)量安全排名省份的后13位省份（比例系數(shù)大于0.02）福建，上海，甘肅，海南，重慶，貴州，西藏，天津，內(nèi)蒙古，吉林，新疆，山西，青海，寧夏中，屬于深圳市鮮菜來(lái)源主要省份的有福建、貴州以及海南，對(duì)于這些省份，深圳市受到的生產(chǎn)地區(qū)帶來(lái)的風(fēng)險(xiǎn)比較大。5.3.4生產(chǎn)方式不同以及本地飲食習(xí)慣不同帶來(lái)的其他風(fēng)險(xiǎn)進(jìn)行分析5.3.4.1反季節(jié)一、我國(guó)的反季節(jié)蔬菜的基本類型第一類是于我國(guó)的低緯度地區(qū)生長(zhǎng)的，在冬日保持一定的溫度與光照，使得蔬菜仍然能夠正常生長(zhǎng)。這一類蔬菜的主要來(lái)源為廣東、云南、海南等較熱帶的區(qū)域，于冬季或春季銷往我國(guó)北部地區(qū)。第二類屬于正常季節(jié)冷藏，并于反季節(jié)銷售的類型，這類作物在收獲季節(jié)就已貯藏好，比如常見(jiàn)的土豆等適合冷藏保存的蔬菜作物，在反季節(jié)進(jìn)行銷售。第三類則是大棚蔬菜。前兩種都不算是反季節(jié)，對(duì)于生產(chǎn)地而言都是順應(yīng)季節(jié)的。二、對(duì)于深圳市反季節(jié)蔬菜的分析深圳市位于我國(guó)的低緯度地區(qū)，其本身就有種植第一類反季節(jié)蔬菜，因此深圳市并不受到外地第一類反季節(jié)蔬菜的影響。第二類反季節(jié)作物重要考慮的是冷藏對(duì)于蔬菜帶來(lái)的影響。蔬菜在冷藏時(shí)溫度較低,蔬菜的呼吸作用減弱,有機(jī)物分解的慢,故能較長(zhǎng)時(shí)間地保持新鮮。但與此同時(shí)，低溫會(huì)抑制果蔬的酵素活動(dòng)，從而使殘毒無(wú)法分解。在冷藏的過(guò)程中，蔬菜容易出現(xiàn)亞硝酸鹽的累積，并且在冷藏過(guò)程中，也會(huì)出現(xiàn)一定程度的腐敗，容易導(dǎo)致蔬菜出現(xiàn)生物毒素的積累。深圳市屬于亞熱帶地區(qū)，氣候干燥濕潤(rùn)，不利于新鮮蔬菜的保存，因此蔬菜冷藏是一種較為常見(jiàn)的現(xiàn)象。所以該類反季節(jié)蔬菜會(huì)對(duì)深圳市的蔬菜質(zhì)量造成不利的影響。第三類考慮大棚蔬菜的影響。利用大棚進(jìn)行蔬菜種植過(guò)程中，會(huì)有氮?dú)狻喯跛釟怏w和二氧化碳等有毒有害氣體產(chǎn)生，而在大棚蔬菜的種植當(dāng)中，也容易出現(xiàn)過(guò)度施肥導(dǎo)致的各類影響。氮?dú)鈺?huì)對(duì)蔬菜帶來(lái)一定的腐蝕性，導(dǎo)致蔬菜根部或是葉片出現(xiàn)斑點(diǎn)和潰爛，引起鮮菜質(zhì)量的下降。當(dāng)氨態(tài)氮肥使用過(guò)度時(shí)，會(huì)導(dǎo)致亞硝酸氣體產(chǎn)生，從而導(dǎo)致土壤呈酸性，在冬季進(jìn)行種植時(shí)，由于通風(fēng)的不足，也容易出現(xiàn)有毒物質(zhì)的累積。我國(guó)大棚蔬菜的使用很廣泛，考慮到深圳市的鮮菜產(chǎn)品的輸入型特性，大棚蔬菜也會(huì)帶來(lái)很大的問(wèn)題。[7]5.3.4.2深圳市民對(duì)蔬菜偏愛(ài)程度考慮到深圳市的飲食習(xí)慣，并通過(guò)深圳市農(nóng)貿(mào)市場(chǎng)提供的近一年的各蔬菜種類的銷量排名以及深圳市有關(guān)蔬菜價(jià)格的變化情況。深圳市對(duì)于本地綠色蔬菜的需求還是較為強(qiáng)烈的。但深圳本地土質(zhì)污染嚴(yán)重。土質(zhì)將直接影響到鮮菜生產(chǎn)的質(zhì)量，所以深圳市還面臨著本地污染土質(zhì)出產(chǎn)的蔬菜的風(fēng)險(xiǎn)。5.4問(wèn)題四的模型建立與求解5.4.1問(wèn)題四的分析本題要求根據(jù)科學(xué)的統(tǒng)計(jì)推斷原理，設(shè)計(jì)一套深圳市新鮮蔬菜抽樣方案，滿足不同層面監(jiān)測(cè)工作的需要，提高食品安全防控的效率與質(zhì)量。步驟一：估計(jì)抽樣方案的樣本量。步驟二：規(guī)定各轄區(qū)具體抽樣任務(wù)。步驟三：針對(duì)不同輸入渠道，確定抽樣比例。步驟四：針對(duì)不同季節(jié)，針對(duì)不同蔬菜類型，確定抽樣比例，提高安全防控的效率。步驟五：針對(duì)不同輸入渠道，設(shè)計(jì)抽樣方法。步驟六：規(guī)定各轄區(qū)設(shè)計(jì)抽樣方案需遵循的規(guī)則。5.4.2原方案分析深圳市2017年食用農(nóng)產(chǎn)品監(jiān)測(cè)抽檢方案規(guī)定了蔬菜的抽樣環(huán)節(jié)、抽樣品種和數(shù)量。其中，現(xiàn)有抽檢環(huán)節(jié)包括農(nóng)批市場(chǎng)、農(nóng)貿(mào)市場(chǎng)、超市商店、配送企業(yè)。農(nóng)產(chǎn)品抽樣總量為全年31200批次，分12個(gè)月完成，每月2600批次。其中，蔬菜的抽樣品種以鱗莖類、葉菜類、蕓苔屬類、瓜類、茄果類、豆類、根莖類、水生類八大類蔬菜為主，其它類蔬菜為輔?，F(xiàn)有方案沒(méi)有基于科學(xué)的統(tǒng)計(jì)推斷原理規(guī)定樣本量；忽略了電商、零售攤販等輸入渠道；沒(méi)有規(guī)定各轄區(qū)每月蔬菜具體抽樣任務(wù)；沒(méi)有對(duì)季節(jié)、蔬菜種類做出具體抽樣比例及方法規(guī)定；無(wú)法在一定置信區(qū)間內(nèi)推斷總體的合格率。對(duì)此，我們?cè)O(shè)計(jì)一套深圳市新鮮蔬菜抽樣方案，綜合考慮深圳市消費(fèi)的蔬菜輸入渠道、種類、銷量，根據(jù)統(tǒng)計(jì)推斷原理，確定合理的樣本量、科學(xué)的樣本量分配，明確抽樣品種、數(shù)量、抽樣方法以及判定原則?，F(xiàn)代抽樣檢驗(yàn)方法即統(tǒng)計(jì)抽樣檢驗(yàn)方法，與傳統(tǒng)的、不科學(xué)的抽樣檢驗(yàn)方法區(qū)別，是建立在概率統(tǒng)計(jì)理論基礎(chǔ)上的，抽樣檢驗(yàn)是利用批或過(guò)程中隨機(jī)抽取的樣本，對(duì)批或過(guò)程的質(zhì)量進(jìn)行檢驗(yàn)[8]，抽樣過(guò)程的關(guān)鍵是采用什么樣的抽樣方式來(lái)保證抽查樣品的代表性。對(duì)于新鮮蔬菜，既要有代表性，能反映當(dāng)?shù)匦迈r蔬菜質(zhì)量整體水平以及不同來(lái)源、不同輸入渠道蔬菜的安全狀況，也要考慮典型性[9]。5.4.3深圳市新鮮蔬菜抽樣方案5.4.3.1樣本量估計(jì)[9]首先，按照隨機(jī)抽樣原則估計(jì)樣本量，用公式①計(jì)算。①公式①中n0為最小樣本量，為總體數(shù)，uα22是對(duì)應(yīng)標(biāo)準(zhǔn)正態(tài)分布尾部面積之和為α的正態(tài)變量值，是之前調(diào)查的不合格率，容許誤差為δ②δ取值不同，計(jì)算的樣本量不同且差異很大。為避免人力物力的浪費(fèi)，需要慎重對(duì)待樣本含量估計(jì)中δ取值的問(wèn)題。不合格率取平均不合格率4.8%，δ取2％，α取0.05，即置信概率為95％，則uα22=1.96，利用公式②計(jì)算得最小樣本量n=865。由于不能完全做到隨機(jī)抽樣，為保證樣本具有代表性，可將樣本量增加10％，n=865新鮮蔬菜抽樣總量為全年11400批次，分12個(gè)月完成，每月950批次。5.4.3.2各轄區(qū)具體抽樣任務(wù)根據(jù)各轄區(qū)新鮮蔬菜樣品不合格率及各轄區(qū)原新鮮蔬菜抽樣比例，提高原不合格率高的轄區(qū)抽樣比例，提高食品安全防控的質(zhì)量，規(guī)定各轄區(qū)每月具體任務(wù)數(shù)量如下：坪山區(qū)50批次、鹽田區(qū)60批次、坪山新區(qū)50批次、南山區(qū)100批次、福田區(qū)100批次、羅湖區(qū)90批次、寶安區(qū)130批次、光明新區(qū)55批次、龍華區(qū)30批次、大鵬新區(qū)45批次、龍華新區(qū)70批次、龍崗區(qū)170批次。表4.1深圳市各轄區(qū)每月新鮮蔬菜抽樣規(guī)定抽樣區(qū)域不合格率原抽樣比例現(xiàn)抽樣比例現(xiàn)抽樣批次坪山區(qū)12.52.18%5%50鹽田區(qū)7.35.23%7%60坪山新區(qū)5.34.63%5%50南山區(qū)5.310.18%10%100福田區(qū)5.39.77%10%100羅湖區(qū)5.29.99%10%90寶安區(qū)5.315.32%14%130光明新區(qū)4.16.67%6%55龍華區(qū)3.43.16%3%30大鵬新區(qū)3.25.09%5%45龍華新區(qū)2.77.02%7%70龍崗區(qū)2.620.77%18%1705.4.3.3抽樣場(chǎng)所安排抽樣需覆蓋深圳市各區(qū)所有蔬菜輸入渠道，除原先已抽樣場(chǎng)所農(nóng)貿(mào)市場(chǎng)、批發(fā)市場(chǎng)、商場(chǎng)超市、生產(chǎn)基地，新增加零售攤販以及電商渠道。根據(jù)問(wèn)題二中得到的六類渠道銷量比例以及不合格率，規(guī)定抽樣場(chǎng)所的比例如表4.2。表4.2深圳市各渠道每月新鮮蔬菜樣本量分配抽樣環(huán)節(jié)農(nóng)貿(mào)市場(chǎng)批發(fā)市場(chǎng)商場(chǎng)超市生產(chǎn)基地零售電商抽樣數(shù)量/批次5001801501010010抽樣比例（%）52.62016.71.111.11.1（各渠道抽樣相對(duì)比例=渠道銷量比例×不合格率）5.4.3.4抽樣蔬菜品種規(guī)定抽樣主要針對(duì)根菜類、鮮豆類、茄果瓜菜類、蔥蒜類、嫩莖葉花菜類、水生蔬菜類、薯芋類、野生蔬菜類八大類蔬菜。由于蔬菜具有季節(jié)性、區(qū)域性、多樣性，蔬菜的抽樣比例，應(yīng)該針對(duì)不同的季節(jié)月份制定相應(yīng)的抽樣方案。例如，春季應(yīng)該相對(duì)增加水生蔬菜類的抽樣，相對(duì)減少薯芋類的抽樣，提高食品安全防控的效率。夏季應(yīng)該加強(qiáng)茄果瓜菜類的抽查力度，嫩莖葉花菜類的抽樣數(shù)量相對(duì)減少。表4.3深圳市各渠道每月新鮮蔬菜樣本量分配春（3-5月）夏（6-8月）秋（9-11月）冬（12-2月）根菜類5.07%4.20%2.83%5.34%鮮豆類8.17%7.83%5.90%8.31%茄果瓜菜類28.53%29.87%25.41%24.33%蔥蒜類0.20%0.67%0.20%0.67%嫩莖葉花菜類54.60%51.17%59.65%55.31%水生蔬菜類2.27%1.60%0.63%1.74%薯芋類0.77%0.97%1.07%1.47%野生蔬菜類0.40%3.70%4.30%2.84%5.4.3.5抽樣方法對(duì)不同渠道進(jìn)行抽樣，采用分層抽樣、多階段抽樣、系統(tǒng)抽樣、簡(jiǎn)單隨機(jī)抽樣相結(jié)合的方式，可以更加科學(xué)全面地得到整個(gè)地區(qū)的情況，可以提高評(píng)價(jià)精度；對(duì)變化性和不確定性進(jìn)行全面分析，其可靠性可以從分析過(guò)程和結(jié)果中得到[10]。對(duì)于不同的輸入渠道，采用多種抽樣方式。農(nóng)貿(mào)市場(chǎng)和生產(chǎn)基地根據(jù)銷售量分配樣本量，將同類蔬菜分別編號(hào)，用系統(tǒng)抽樣法抽取樣本。系統(tǒng)抽樣，是將蔬菜編號(hào)按一定次序排列，然后按相等的間隔抽取樣本單位。例如，某一類蔬菜有125份，抽取5份樣本，則間隔為125/5=25，即從1號(hào)樣品開(kāi)始，每隔25號(hào)抽取一份樣品。批發(fā)市場(chǎng)批發(fā)市場(chǎng)采用分層抽樣的方式。將批發(fā)市場(chǎng)按銷售量分層。例如，將銷售量較大的布吉農(nóng)產(chǎn)品批發(fā)市場(chǎng)、海吉星農(nóng)產(chǎn)品批發(fā)市場(chǎng)作為第一層，寶安農(nóng)產(chǎn)品批發(fā)市場(chǎng)、福田農(nóng)產(chǎn)品批發(fā)市場(chǎng)作為第二層，銷量較小的作為第三層，按照銷量比例，分配樣本量，采用系統(tǒng)抽樣法在每一層中抽樣。商場(chǎng)超市因商場(chǎng)超市總體數(shù)量有限，并且可以獲得所有商場(chǎng)超市的名錄，所以可以采用分層抽樣的方法。例如，將商場(chǎng)超市分為大型商場(chǎng)、小型超市兩層，在各層內(nèi)隨機(jī)抽樣，將所有的商場(chǎng)和超市編號(hào)，用系統(tǒng)抽樣法，隨機(jī)抽取目標(biāo)樣本量。樣本量根據(jù)每一層所占比例按比例分配，若大型商場(chǎng)蔬菜量占60%，則抽樣數(shù)量為150*60%=90批次，小型超市抽樣數(shù)量為150-90=60批次。零售攤販因零售攤販數(shù)量多、無(wú)法獲得全部零售攤販名錄，所以采用分級(jí)抽樣的方法，即先隨機(jī)抽取一級(jí)單元，再?gòu)某榈玫囊患?jí)單元中隨機(jī)抽取范圍更小的二級(jí)單元。例如，要抽取30個(gè)零售攤販，先將抽樣轄區(qū)分成30個(gè)街區(qū)，從中隨機(jī)抽取6個(gè)街區(qū)，每個(gè)街區(qū)再隨機(jī)抽取零售攤販。電商渠道電商渠道是之前抽查過(guò)程忽略的抽樣渠道，所以應(yīng)該加強(qiáng)抽樣力度，因?yàn)閺南M(fèi)者在電商上購(gòu)買新鮮蔬菜與蔬菜交付到消費(fèi)者這段時(shí)間的溫度與保質(zhì)期控制會(huì)受到物流的影響，進(jìn)而影響蔬菜的新鮮程度、潔凈程度，所以，除了對(duì)電商的貨源和倉(cāng)庫(kù)中的產(chǎn)品抽樣，應(yīng)該增加物流過(guò)程、配送過(guò)程中的抽樣，確保最終交付到消費(fèi)者的蔬菜食品安全得到保障。其他規(guī)定各轄區(qū)制定抽樣計(jì)劃須遵循以下原則：（1）各轄區(qū)抽樣渠道需覆蓋全部六類抽樣渠道，按照農(nóng)貿(mào)市場(chǎng)站總抽樣量50%，批發(fā)市場(chǎng)占總抽樣量的20%，其余在批發(fā)市場(chǎng)、商場(chǎng)超市、零售、電商完成，其中，應(yīng)對(duì)之前沒(méi)有抽取的零售攤販與電商重點(diǎn)抽查。（2）抽樣位點(diǎn)需覆蓋轄區(qū)內(nèi)所有蔬菜經(jīng)營(yíng)場(chǎng)所，每個(gè)經(jīng)營(yíng)為點(diǎn)本年度內(nèi)需要覆蓋到位，對(duì)于不合格率較大的商場(chǎng)超市、零售攤販、農(nóng)貿(mào)市場(chǎng)、電商渠道，適當(dāng)增加抽樣品種數(shù)量和頻次。（3）按照每月銷量（見(jiàn)表1.4）進(jìn)行抽樣比例分配。對(duì)于每月主銷蔬菜與反季節(jié)蔬菜應(yīng)該重點(diǎn)抽查與二次抽查。（4）針對(duì)市民喜愛(ài)的蔬菜，應(yīng)該二次抽查，將漏查帶來(lái)的食品安全風(fēng)險(xiǎn)降到最低。5.4.3.6結(jié)果分析通過(guò)以上抽樣方案，每月抽取950批次，通過(guò)不合格率評(píng)價(jià)深圳市蔬菜質(zhì)量。置信概率95%，允許誤差不超過(guò)2%，若檢測(cè)結(jié)果的不合格率小于平均不合格率4.8%，則本次方案的樣本量設(shè)計(jì)滿足要求，否則，說(shuō)明深圳蔬菜質(zhì)量安全狀況比預(yù)期的安全狀況差，需要加大監(jiān)測(cè)力度,增加樣本量。5.5問(wèn)題五的解答尊敬的深圳市政府領(lǐng)導(dǎo)：您好！我們是來(lái)自外省某高校的學(xué)生。關(guān)于貴市蔬菜安全問(wèn)題，我們希望能夠發(fā)表一些看法和建議。若您能撥冗詳讀這封信，我們將不勝感激。作為廣東第二大開(kāi)放城市和珠江三角洲的重要城市，貴市無(wú)時(shí)無(wú)刻不發(fā)揮著經(jīng)濟(jì)中心和綜合樞紐的作用。經(jīng)濟(jì)發(fā)展與社會(huì)繁榮相互促進(jìn)。全市2000多萬(wàn)人口當(dāng)中，就有1000萬(wàn)左右為流動(dòng)人口。另外，據(jù)可靠消息顯示，貴市的蔬菜有90%以上來(lái)自外地。在如此背景之下，食品安全尤其是蔬菜安全，就更應(yīng)當(dāng)引起重視。近年來(lái)，越來(lái)越多的人開(kāi)始擔(dān)心蔬菜存在安全隱患。去年10月，貴市食品藥品監(jiān)管局還曾發(fā)布關(guān)于試行市民免費(fèi)快速檢測(cè)蔬菜受理點(diǎn)的公示。出于個(gè)人興趣，以及作為當(dāng)代大學(xué)生的使命感和責(zé)任感，我們花了將近一周時(shí)間對(duì)貴市蔬菜安全情況做了一系列調(diào)查。根據(jù)質(zhì)檢局和其他渠道提供的抽樣監(jiān)測(cè)數(shù)據(jù)，我們發(fā)現(xiàn)了整個(gè)蔬菜抽檢體系的許多問(wèn)題，現(xiàn)將我們的發(fā)現(xiàn)列舉如下：蔬菜消費(fèi)的季節(jié)性與渠道差異性。1.根據(jù)我們的統(tǒng)計(jì)結(jié)果，貴市市民喜食綠葉菜、茄果瓜菜和豆類蔬菜，少食蔥蒜菜類，這種蔬菜飲食結(jié)構(gòu)與全國(guó)其余地區(qū)基本無(wú)異。2.貴市蔬菜消費(fèi)具有明顯的季節(jié)性。早春季節(jié)，葉菜銷量較大；春夏之際，天氣轉(zhuǎn)熱，葉菜消費(fèi)逐漸被瓜豆蔬菜所取代；秋冬季節(jié)，葉菜銷量反彈明顯。相比而言，其他菜類則一直處于銷量低谷。3.市民的蔬菜來(lái)源主要是當(dāng)?shù)剞r(nóng)貿(mào)市場(chǎng)和批發(fā)市場(chǎng)，選擇超市的極少，選擇電商的最少。究其原因，市民更加信賴前兩個(gè)蔬菜輸入渠道的安全性，而大多數(shù)人也對(duì)網(wǎng)購(gòu)蔬菜的新鮮程度和有害與否存疑。因此，貴政府應(yīng)當(dāng)充分考慮當(dāng)?shù)厥忻竦娘嬍沉?xí)慣，加大對(duì)農(nóng)貿(mào)批發(fā)市場(chǎng)蔬菜銷售攤點(diǎn)的監(jiān)管和懲治力度，并制定一系列與季節(jié)波動(dòng)相關(guān)的抽檢方案。農(nóng)貿(mào)市場(chǎng)和超市蔬菜更具威脅。1.據(jù)統(tǒng)計(jì)，在超市蔬菜的不合格率竟達(dá)8.63%之高，從側(cè)面看，這也是少有市民選擇這一渠道的原因。2.相反，生產(chǎn)基地和批發(fā)市場(chǎng)的蔬菜較為安全，不合格率基本處于1-2%區(qū)間內(nèi)。據(jù)此我們可以不精確地得出“越靠近生產(chǎn)源頭的蔬菜越安全”的結(jié)論。3.若在抽檢過(guò)程中去除農(nóng)貿(mào)市場(chǎng)這一場(chǎng)所，則抽檢結(jié)果的可靠性將驟降。換句話說(shuō)，忽略農(nóng)貿(mào)市場(chǎng)的蔬菜檢驗(yàn)，對(duì)全市蔬菜安全的影響是致命的。綜上，我們建議貴市完善法律法規(guī)，從源頭遏制蔬菜安全威脅。另外，在每月例行的“豬肉產(chǎn)品和蔬菜質(zhì)量安全監(jiān)測(cè)”中，建議增加農(nóng)貿(mào)市場(chǎng)和超市的樣本數(shù)，這將有利于更多蔬菜安全問(wèn)題的排查。三、進(jìn)行鮮菜檢驗(yàn)時(shí)應(yīng)擴(kuò)大檢驗(yàn)種類并考慮地域性和本地特性。1．除常見(jiàn)的農(nóng)藥殘留可能對(duì)于鮮菜的危險(xiǎn)較大之外，鮮菜當(dāng)中也會(huì)隱藏一些在例如重金屬元素以及6-芐基腺嘌呤等對(duì)食品安全有重大影響的“影子殺手”，在鮮菜的檢查當(dāng)中，這也是一項(xiàng)不容忽視的因素。2.中國(guó)幅員遼闊，鮮菜的場(chǎng)地眾多，然而不同的區(qū)域以及不同的自然環(huán)境生產(chǎn)出的鮮菜質(zhì)量不同，對(duì)于進(jìn)口鮮菜的選擇，我們應(yīng)該“擇優(yōu)錄取”，并加強(qiáng)對(duì)不合格率較高的省份的抽查力度。3.反季節(jié)種植是一種新型種植模式，這種新型的種植方式方便了我們的日常生活，但與此同時(shí)，反季節(jié)鮮菜由于其特殊的種植模式也會(huì)帶來(lái)不同的安全問(wèn)題，因此我們應(yīng)該“因物制宜”，選擇合適的檢驗(yàn)方式。四、應(yīng)以零售攤販和電商為重點(diǎn)抽查對(duì)象，關(guān)注當(dāng)季、反季蔬菜和市民偏愛(ài)的蔬菜。針對(duì)貴市的蔬菜問(wèn)題，建立一套合理可靠的抽檢體系極為必要。但從我們調(diào)查統(tǒng)計(jì)的數(shù)據(jù)來(lái)看，貴市現(xiàn)有的抽樣方案科學(xué)性欠佳，忽略了電商、零售攤販等輸入渠道，對(duì)各轄區(qū)每月蔬菜具體抽樣任務(wù)、不同季節(jié)蔬菜種類抽樣比例及方法沒(méi)有具體規(guī)定。因此，我們提議設(shè)計(jì)一套深圳市新鮮蔬菜抽樣方案，綜合考慮深圳市消費(fèi)的蔬菜輸入渠道、種類、銷量，明確抽樣品種、數(shù)量、抽樣方法以及判定原則。值得注意的是，在新的蔬菜抽檢方案中應(yīng)當(dāng)對(duì)之前忽視的零售攤販與電商重點(diǎn)抽查；對(duì)于每月主銷蔬菜與反季節(jié)蔬菜應(yīng)該重點(diǎn)抽查與二次抽查；針對(duì)市民喜愛(ài)的蔬菜，應(yīng)該二次抽查，確保蔬菜安全和市民健康。以上即為我們以自己的淺陋學(xué)識(shí)提出的幾點(diǎn)建議。最后我們想說(shuō)，一個(gè)城市的發(fā)展好壞，關(guān)鍵在于民心所向。倘若要使人民能夠更好地在深圳生存發(fā)展，就絕不能忽視群眾健康和食品安全的重要性。衷心希望您能考慮我們的建議，并祝深圳發(fā)展越來(lái)越好。此致敬禮！XXX2017.4.26模型檢驗(yàn)、評(píng)價(jià)與推廣6.1模型的檢驗(yàn)6.1.1對(duì)于K-聚類法中的各個(gè)影響因素的相關(guān)性與重要程度排序的檢驗(yàn)表6.1關(guān)于鐵路距離以及運(yùn)價(jià)之間的正態(tài)性以及相關(guān)性檢驗(yàn)相關(guān)性檢驗(yàn)鐵路距離運(yùn)價(jià)鐵路距離Pearson相關(guān)性1.871**顯著性（雙側(cè)）.000N3131運(yùn)價(jià)Pearson相關(guān)性.871**1顯著性（雙側(cè)）.000N3131正態(tài)性檢驗(yàn)Kolmogorov-SmirnovaShapiro-Wilk統(tǒng)計(jì)量dfSig.統(tǒng)計(jì)量dfSig.鐵路距離.11131.200*.96531.388運(yùn)價(jià).17131.022.93931.079考慮到鐵路距離與運(yùn)價(jià)之間可能存在的相關(guān)關(guān)系，我們首先對(duì)鐵路距離與運(yùn)價(jià)樣本進(jìn)行正態(tài)性檢驗(yàn)，由于樣本數(shù)小于200，所以我們查看Shapiro-Wilk數(shù)值，P值均大于0.05，所以接受原假設(shè)，數(shù)據(jù)符合正態(tài)分布。接著對(duì)數(shù)據(jù)進(jìn)行相關(guān)性檢驗(yàn)，鐵路距離與運(yùn)價(jià)間的pearson相關(guān)性系數(shù)為0.871，P值小于0.05，所以鐵路距離與運(yùn)價(jià)之間強(qiáng)相關(guān)。因此影響因素中將運(yùn)價(jià)作為較重要的因子，同時(shí)將鐵路距離的重要性后移是正確的。同理得出，其他因數(shù)之間是否存在相關(guān)關(guān)系以及其相關(guān)程度。KMO和Bartlett的檢驗(yàn)取樣足夠度的Kaiser-Meyer-Olkin度量。.786Bartlett的球形度檢驗(yàn)近似卡方110.481df15Sig..000經(jīng)過(guò)Bartlett的球形度檢驗(yàn)，P值小于0.05,說(shuō)明各變量之間具有一定的相關(guān)性，因子分析有效。6.1.2平滑指數(shù)法所得估計(jì)結(jié)果的檢驗(yàn)?zāi)Ｐ徒y(tǒng)計(jì)量模型Ljung-BoxQ(18)統(tǒng)計(jì)量DFSig.模型14.09916.591在平滑指數(shù)法預(yù)測(cè)時(shí)，我們選用了簡(jiǎn)單季節(jié)性模型，而不是非季節(jié)性模型。檢驗(yàn)中LBQ值=0.591>0.05，則一個(gè)或多個(gè)滯后的自相關(guān)可能顯著不同于零，說(shuō)明在這段時(shí)間內(nèi)各個(gè)值并不是獨(dú)立和隨機(jī)的,該結(jié)果符合簡(jiǎn)單季節(jié)性的假設(shè)6.2模型的評(píng)價(jià) 6.2.1模型的優(yōu)點(diǎn)1．對(duì)于風(fēng)險(xiǎn)物質(zhì)的評(píng)價(jià)，采用多維灰聚類評(píng)估方法，其為數(shù)學(xué)層次分析為基礎(chǔ)，結(jié)合相關(guān)專業(yè)知識(shí)和實(shí)踐經(jīng)驗(yàn)的一種結(jié)合主觀與客觀的綜合評(píng)價(jià)方法，常用于食品方面的風(fēng)險(xiǎn)分析，對(duì)于存在一定灰度的問(wèn)題有著良好的解決功能。2．模型中采用了多種聚類方法，根據(jù)不同的需要，選擇了最適合的聚類模型，如多維尺度分析和K-聚類分析，并對(duì)聚類模型進(jìn)行了檢驗(yàn)，制作了直觀的分類圖形。3．采用