食品質(zhì)量安全——抽檢數(shù)據(jù)分析

1、食品質(zhì)量安全抽檢數(shù)據(jù)分析目錄目錄問題的重述問題的重述一一問題的分析問題的分析二二模型的假設(shè)模型的假設(shè)三三符號約定符號約定四四模型的建立與求解模型的建立與求解五五問題的分析問題的分析六六模型的評價與推廣模型的評價與推廣七七一、問題的重述一、問題的重述1、評價評價食品安全食品安全情況變化情況變化趨勢趨勢2、找出與找出與食品質(zhì)量食品質(zhì)量相關(guān)因素相關(guān)因素3、改進改進食品抽檢食品抽檢方法方法1、評價食品安全情況變化趨勢、評價食品安全情況變化趨勢逐年統(tǒng)計各項指標合格率逐年統(tǒng)計各項指標合格率逐年統(tǒng)計各項指標超標百分比逐年統(tǒng)計各項指標超標百分比橫向、縱向比較分析橫向、縱向比較分析2、找出與食品質(zhì)量相關(guān)因素、找

2、出與食品質(zhì)量相關(guān)因素橫向?qū)Ρ葯M向?qū)Ρ纫蛩胤治鲆蛩胤治鰈 生產(chǎn)環(huán)節(jié)生產(chǎn)環(huán)節(jié)l 流通環(huán)節(jié)流通環(huán)節(jié)l 餐飲消費餐飲消費環(huán)節(jié)環(huán)節(jié)環(huán)節(jié)環(huán)節(jié)l Logestic回回歸歸分析模型分析模型l 虛擬變量虛擬變量l 比較比較Wald統(tǒng)統(tǒng)計值計值l 食品產(chǎn)地食品產(chǎn)地l 省內(nèi)省內(nèi)l 省外省外l 檢測季節(jié)檢測季節(jié)l 食品檢測地點食品檢測地點l 包裝包裝3、改進食品抽檢方法、改進食品抽檢方法1、統(tǒng)計處理上一次抽檢結(jié)、統(tǒng)計處理上一次抽檢結(jié)果果.3、調(diào)整下一、調(diào)整下一次次每種每種食品的檢食品的檢測測次數(shù)次數(shù)2、根據(jù)評價根據(jù)評價結(jié)果，結(jié)果，考慮產(chǎn)地考慮產(chǎn)地季節(jié)等季節(jié)等因素因素.以抽檢結(jié)果的方差和平均值以抽檢結(jié)果的方差和平均值為

3、參數(shù)的評價函數(shù)模型為參數(shù)的評價函數(shù)模型三、模型假設(shè)三、模型假設(shè)1實際深圳市市場監(jiān)督管理局網(wǎng)站上的所有檢測數(shù)據(jù)都準確可靠實際深圳市市場監(jiān)督管理局網(wǎng)站上的所有檢測數(shù)據(jù)都準確可靠2食品都可以歸類到食品分類系統(tǒng)當中食品都可以歸類到食品分類系統(tǒng)當中4劃分之后同一地區(qū)的生產(chǎn)廠家到深圳具有相同的運輸距離劃分之后同一地區(qū)的生產(chǎn)廠家到深圳具有相同的運輸距離5對任何產(chǎn)品抽檢一批次的成本相同對任何產(chǎn)品抽檢一批次的成本相同3抽檢指標都可以歸類到抽檢項目分類系統(tǒng)當中抽檢指標都可以歸類到抽檢項目分類系統(tǒng)當中四、符號約定四、符號約定五、模型的建立與求解五、模型的建立與求解5.1 深圳市主要食品安全情況的變化趨勢深圳市主要食

4、品安全情況的變化趨勢5.2 食品安全相關(guān)因素分析食品安全相關(guān)因素分析5.3 抽樣方法的改進抽樣方法的改進5.1.1 檢測指標的分類檢測指標的分類微生物微生物 工業(yè)殘留及工業(yè)殘留及其他其他 食品添食品添加劑加劑 重金屬重金屬 檢測指標檢測指標超標類超標類未達標類未達標類其他類其他類食品的分類食品的分類1.乳及乳制品（簡稱乳制品）；2.脂肪，油和乳化脂肪制品（簡稱食用油）；3.冷飲，飲料及酒（簡稱飲品）；4.水果、蔬菜（包括塊根類）、豆類、食用菌、堅果以及籽類等（簡稱果蔬）；5.可可制品、巧克力和巧克力制品、糖果及甜味料（簡稱甜味料）；6.糧食及糧食制品，包括大米、面粉、雜糧、淀粉等（簡稱糧食）；

5、7.焙烤食品，包括面包、糕點、餅干等（簡稱焙烤食品）；8.肉蛋類及其制品（簡稱肉蛋制品）；9.水產(chǎn)及其制品，包括魚類、甲殼類、貝類、軟體類、棘皮類等水產(chǎn)及其加工制品等（簡稱水產(chǎn)品）；10.調(diào)味品11.其他類，即不包含在前十類中的食品。數(shù)據(jù)的處理數(shù)據(jù)的處理v1）合格率）合格率20102010年年20112011年年20122012年年乳制品乳制品100.00%100.00%99.80%99.80%100.00%100.00%食用油食用油98.09%98.09%98.11%98.11%99.09%99.09%飲品飲品98.54%98.54%97.72%97.72%98.07%98.07%果蔬果蔬9

6、4.59%94.59%92.04%92.04%97.99%97.99%甜味料甜味料97.89%97.89%91.35%91.35%95.61%95.61%糧食糧食99.59%99.59%94.48%94.48%97.04%97.04%焙烤食品焙烤食品99.73%99.73%81.67%81.67%96.91%96.91%肉蛋制品肉蛋制品95.45%95.45%86.16%86.16%93.63%93.63%水產(chǎn)品水產(chǎn)品96.84%96.84%92.86%92.86%94.88%94.88%調(diào)味品調(diào)味品98.60%98.60%91.59%91.59%94.99%94.99%其他類其他類98.01

7、%98.01%96.77%96.77%95.62%95.62%80.0%82.0%84.0%86.0%88.0%90.0%92.0%94.0%96.0%98.0%100.0%乳制品食用油飲品果蔬甜味料糧食焙烤食品肉蛋制品水產(chǎn)品調(diào)味品其他類2010年2011年2012年圖圖 1 2010年年2012年各類主要年各類主要食品食品 抽檢抽檢合格率比較圖合格率比較圖2）超標百分比）超標百分比數(shù)據(jù)的處理數(shù)據(jù)的處理非非微生物：微生物：微生物：微生物：表表 5 三年微生物超標情況（未做對數(shù)處理前）三年微生物超標情況（未做對數(shù)處理前）超標超標數(shù)量數(shù)量超標百分比超標百分比均值均值中位數(shù)中位數(shù)均值的均值的標準誤標

8、準誤標準差標準差2010年年30105961% 2933%46225312011年年294123269% 2567%648111202012年年1254296%2567%559表表 6 部分微生物超標情況特別嚴重的產(chǎn)品一覽表部分微生物超標情況特別嚴重的產(chǎn)品一覽表抽檢日期抽檢日期樣品名稱樣品名稱受檢單位受檢單位產(chǎn)地產(chǎn)地標準要求標準要求檢測結(jié)果檢測結(jié)果2010/9/1波紋干津市市南北特食品實業(yè)有限公司深圳菌落總數(shù)750(cfu/g)8.1106(cfu/g)2011/2/18鹵豆干深圳市南山區(qū)金華城百貨商場日用品經(jīng)營部長沙菌落總數(shù)750(cfu/g)1.9107(cfu/g)2011/9/28香辣

9、卷深圳市佳滿福百貨有限公司筍崗分店長沙菌落總數(shù)750(cfu/g)1.4108(cfu/g)2012/12/3即食煎牛扒深圳市三合春餐飲管理有限公司深圳大腸菌落90(MPN/100g)2.4104(MPN/100g)超標數(shù)量超標數(shù)量均值均值均值的標準誤均值的標準誤標準差標準差微生物微生物（處理前）（處理前）2010年30105961%462.32531.92011年294123269%648.611120.82012年1254296%5.359.2微生物微生物（處理后）（處理后）2010年3071.61%0.07850.42982011年29464.78%0.02580.44252012年12

10、557.71%0.03940.4399重金屬重金屬2010年8270.88%0.91992.60182011年23347.70%0.91214.37442012年3197.33%0.83151.4402食品添加劑食品添加劑2010年41205.11%0.42532.72342011年264200.66%0.11151.81182012年133141.44%0.12871.4844表表 7 三年內(nèi)各類食品微生物三年內(nèi)各類食品微生物、重金屬、食品添加劑重金屬、食品添加劑超標情況超標情況圖圖 2 各各類食品微生物類食品微生物、重金屬、重金屬、食品添加劑食品添加劑三年內(nèi)超標情況三年內(nèi)超標情況比較（微生

11、物為處理后）比較（微生物為處理后）72%271%205%65%348%201%58%197%141%0%50%100%150%200%250%300%350%400%微生物重金屬食品添加劑2010年2011年2012年表表 8 2010年六類主要食品微生物年六類主要食品微生物、重金屬重金屬、食品添加劑超標情況匯總、食品添加劑超標情況匯總2010年年微生物微生物重金屬重金屬食品添加劑食品添加劑食品分類食品分類超標數(shù)量超標百分比超標數(shù)量超標百分比超標數(shù)量超標百分比處理前處理后飲品飲品12100%79.0%0257.3%果蔬果蔬11 284078% 103.4%3 308.3%24296.2%糧食糧

12、食21643%33.9%00焙烤食品焙烤食品001128.0%肉蛋制品肉蛋制品62387%45.8%01900.0%水產(chǎn)品水產(chǎn)品73943%62.6%00表表 9 2011年六類主要食品微生物年六類主要食品微生物、重金屬重金屬、食品添加劑超標情況匯總、食品添加劑超標情況匯總2011年年微生物微生物重金屬重金屬食品添加劑食品添加劑食品分類食品分類超標數(shù)量超標百分比超標數(shù)量超標百分比超標數(shù)量超標百分比處理前處理后飲品飲品6713%36.8%0374.5%果蔬果蔬76 406698%86.1%6 309.2%88264.1%糧食糧食21446%72.6%1 550.0%17247.0%焙烤食品焙烤食

13、品321175%40.1%0 137227.0%肉蛋制品肉蛋制品1495978%55.7%8 553.4%683.3%水產(chǎn)品水產(chǎn)品195220%44.3%2 115.0%0表表 10 2012年六類主要食品微生物年六類主要食品微生物、重金屬重金屬、食品添加劑超標情況匯總、食品添加劑超標情況匯總2012年年微生物微生物重金屬重金屬食品添加劑食品添加劑食品分類食品分類超標數(shù)量超標百分比超標數(shù)量超標百分比超標數(shù)量超標百分比處理前處理后飲品飲品34217% 101.1%03131.1%果蔬果蔬38329%52.4%049149.8%糧食糧食0027146.3%焙烤食品焙烤食品201332%41.0%0

14、31113.2%肉蛋制品肉蛋制品815104%55.6%010205.7%水產(chǎn)品水產(chǎn)品113785%44.0%1 340.0%2134.5%3）獨立樣本）獨立樣本T檢驗檢驗數(shù)據(jù)的處理數(shù)據(jù)的處理獨立樣本獨立樣本T檢驗檢驗表122010年和2011年食品添加劑超標情況比較的顯著性檢驗食品添食品添加劑加劑方差方程的方差方程的Levene檢驗檢驗均值方程的均值方程的t檢驗檢驗FSig.tdfSig.均值差值標準誤差值95%置信區(qū)間下限上限假設(shè)方差相等10.120.0020.1363030.8920.0450.328-0.6020.691假設(shè)方差不相等0.101460.9200.0450.440-0.8

15、410.930年份年份N均值均值標準差標準差均值的標準誤均值的標準誤2010412.05102.72340.425320112642.00651.81180.1115表112010年和2011年食品添加劑超標百分比的均值、標準差比較獨立樣本獨立樣本T檢驗檢驗表142011年和2012年食品添加劑超標情況比較的顯著性檢驗食品添食品添加劑加劑方差方程的方差方程的Levene檢驗檢驗均值方程的均值方程的t檢驗檢驗FSig.tdfSig.均值差值標準誤差值95%置信區(qū)間下限上限假設(shè)方差相等7.9100.0053.2583950.0010.5920.182-0.2350.950假設(shè)方差不相等3.4773

16、150.0010.5920.170-0.2570.927年份年份N均值均值標準差標準差均值的標準誤均值的標準誤20112642.00651.81180.111520121331.41431.48430.1287表132011年和2012年食品添加劑超標百分比的均值、標準差比較表表 15 三年內(nèi)微生物、三年內(nèi)微生物、重金屬重金屬 超標超標情況比較的顯著性檢驗情況比較的顯著性檢驗微生物微生物重金屬重金屬處理前處理后檢驗?zāi)攴輽z驗?zāi)攴?0年、11年11年、12年10年、11年11年、12年10年、11年11年、12年方差的方差的Levene檢驗檢驗F0.0644.4080.1550.1611.1191

17、.244Sig0.8010.0450.6940.6890.2990.276均值的均值的t檢驗檢驗t-0.0851.8340.8081.499-0.4660.582Sig0.9320.0680.4200.1350.6450.645v分析結(jié)果分析結(jié)果 微生物微生物：超標百分比均值穩(wěn)中有降，說明食品中微生物的安全情況逐年略有改善 重金屬重金屬：超標數(shù)量較少，2012年的超標情況稍要好于2011年 食品添加劑食品添加劑：2011年較2010年的超標情況基本沒有太大改變，但2012年食品添加劑的安全情況較2011年有了顯著性的改善。 總的來說，2012年的各項食品安全情況較前一年都有了一定的改善，說明深

18、圳市市場監(jiān)督管理局對食品的定期抽檢還是有其積極意義的。數(shù)據(jù)的處理數(shù)據(jù)的處理5.2 食品安全相關(guān)因素分析食品安全相關(guān)因素分析5.2.2 食品產(chǎn)地食品產(chǎn)地5.2.1 食品各個環(huán)節(jié)安全情況食品各個環(huán)節(jié)安全情況5.2.4 食品銷售地點（即檢測地點）食品銷售地點（即檢測地點）5.2.3 檢測季節(jié)檢測季節(jié)5.2.5 食品包裝食品包裝5.2.6 綜合因素分析綜合因素分析5.2.1 食品各個環(huán)節(jié)安全情況食品各個環(huán)節(jié)安全情況抽檢總數(shù)抽檢總數(shù)合格數(shù)合格數(shù)合格率合格率生產(chǎn)加工環(huán)節(jié)生產(chǎn)加工環(huán)節(jié)77974795.9%流通銷售環(huán)節(jié)流通銷售環(huán)節(jié)1454142497.9%餐飲服務(wù)環(huán)節(jié)餐飲服務(wù)環(huán)節(jié)2537247397.5%表表

19、 16 2010年食品各個環(huán)節(jié)抽檢情況對比表年食品各個環(huán)節(jié)抽檢情況對比表v生產(chǎn)環(huán)節(jié)的抽檢合格率較其他兩個環(huán)節(jié)要稍低一些 可能是由于生產(chǎn)廠家中衛(wèi)生情況較差造成的。5.2.2 食品產(chǎn)地食品產(chǎn)地食品食品產(chǎn)地產(chǎn)地與與食品質(zhì)量關(guān)系食品質(zhì)量關(guān)系Spearman秩秩相關(guān)性分析相關(guān)性分析列聯(lián)相關(guān)進行分析列聯(lián)相關(guān)進行分析廣東省外廣東省外分析方法：分析方法：v Spearman秩相關(guān)性分析秩相關(guān)性分析運輸距離與食品質(zhì)量v 列聯(lián)相關(guān)進行列聯(lián)相關(guān)進行分析分析城市本身與食品質(zhì)量v 產(chǎn)地劃分產(chǎn)地劃分廣東省內(nèi)城市廣東省外省份、自治區(qū)和直轄市廣東省廣東省內(nèi)內(nèi)廣東省內(nèi)廣東省內(nèi)城市城市距離距離抽檢總數(shù)抽檢總數(shù)合格數(shù)合格數(shù)合格率合

20、格率不合格率不合格率深圳深圳0118771169698.48%1.52%東莞東莞691159111896.46%3.54%惠州惠州92.314213796.48%3.52%中山中山12136735897.55%2.45%廣州廣州1431182114897.12%2.88%佛山佛16%3.84%江門江門151949398.94%1.06%珠海珠83%2.17%開平開平233918896.70%3.30%揭西揭西2671239577.24%22.76%普寧普寧28414713793.20%6.80%揭陽揭陽327524994.23%5.77%揭東

21、揭東335918795.60%4.40%汕頭汕頭34130128795.35%4.65%潮州潮州3615252100.00%0.00%廣東省外廣東省外省份省份距離距離抽檢總數(shù)抽檢總數(shù)合格數(shù)合格數(shù)合格率合格率不合格率不合格率海南海南6796565100.00%0.00%廣西廣西73715515197.42%2.58%湖南湖南78330325082.51%17.49%江西江西81116816698.81%1.19%福建福建84166165498.94%1.06%湖北湖北111419919597.99%2.01%安徽安徽1243787798.72%1.28%浙江浙江129725224697.62%2

22、.38%江蘇江00%0.00%上海上海148844744098.43%1.57%云南云00%0.00%河南河南159519619599.49%0.51%重慶重慶168921420495.33%4.67%山東山東186436636599.73%0.27%河北河北1959919098.90%1.10%四川四川199010710396.26%3.74%山西山西199916015999.38%0.62%北京北京216115615196.79%3.21%天津天津21636868100.00%0.00%內(nèi)蒙古內(nèi)蒙古2640192192100.00%0.0

23、0%黑龍江黑龍江3350135135100.00%0.00%v 值df漸進 Sig. Pearson 卡方670.447350.000有效案例中的 N203915.2.2 食品產(chǎn)地食品產(chǎn)地v數(shù)據(jù)處理數(shù)據(jù)處理 廣東省外以省份劃分，產(chǎn)地距離為省會城市到深圳市的高速公路距離 廣東省內(nèi)以城市為依據(jù)進行劃分，產(chǎn)地距離為該城市到深圳市的高速公路距離 以上所有距離通過百度地圖得到 統(tǒng)一剔除檢測數(shù)目小于50的產(chǎn)地v方法方法 Spearman相關(guān)性分析Spearman秩相關(guān)性秩相關(guān)性分析分析結(jié)果結(jié)果省內(nèi)省內(nèi)省外省外全國全國運輸距離合格率相關(guān)系數(shù)-0.2700.5700.468Sig.0.3500.0140.00

24、4N141832是否相關(guān)否是是Spearman秩相關(guān)性分析秩相關(guān)性分析 結(jié)果說明：結(jié)果說明： 運輸距離與食品質(zhì)量的相關(guān)性： 廣東省內(nèi)城市無 廣東省外的各個省份中度相關(guān)性，正相關(guān)性 即運輸距離越遠食品合格率越高，食品質(zhì)量越好。 結(jié)論分析結(jié)論分析 首先，廣東省內(nèi)檢測食品數(shù)量較多 其次，由于運輸距離短的食品流通更為便捷，由此食品種類更為繁雜，監(jiān)管難度更為大。 再次，當運輸距離達到一定長度后，食品運輸方式或會采用空運等更為快速且安全的運輸方式5.2.3 檢測季節(jié)檢測季節(jié)乳制乳制品品食用食用油油飲品飲品果蔬果蔬甜味甜味料料糧食糧食焙烤焙烤食品食品肉蛋肉蛋制品制品水產(chǎn)水產(chǎn)品品調(diào)味調(diào)味品品10-110-29

25、7.6% 98.5% 98.9%99.4%97.5% 96.8%100%10-3100%100%100% 92.3%99.7% 99.7% 85.7%100% 98.5%10-4100%93.9% 87.3% 96.0% 99.4% 97.5% 96.8% 97.2%100%11-198.7% 95.0% 96.7% 89.1%92.2% 66.9% 73.3% 98.0%11-2100% 98.5%100% 94.6% 92.3%100%98.0%11-3100%97.5% 87.6%100%100% 96.3% 95.9% 88.2% 99.0%11-4100% 98.1% 98.6% 9

26、4.3% 88.1% 95.0% 88.8% 86.5%84.5%12-1100%98.8% 99.5% 90.0% 99.2%100% 93.3%12-2100% 99.4% 95.7% 86.5% 97.8% 98.3% 91.3% 88.6% 95.7% 95.8%12-3100%100% 97.1%100%100% 95.8% 98.4% 96.0%12-4100% 97.8% 99.4% 99.6%95.9% 92.4% 95.2% 94.2% 93.5%13-1100%99.5% 95.7%100%100%100% 92.1%13-2100%100%100% 91.7% 93.0%

27、 85.6%86.7%表 32 各類食品分季度合格率匯總表5.2.3 檢測季節(jié)檢測季節(jié)v方法：列聯(lián)分析方法：列聯(lián)分析 某些食品的樣本量較少，或不合格數(shù)量為0，無法進行卡方檢驗 僅對飲品、果蔬、糧食、焙烤食品、肉蛋制品以及水產(chǎn)品這六類主要食品進行檢驗 以一、二、三、四分別表示四個季度，0表示不合格，1表示合格 以果蔬為例，得到列聯(lián)表如下：5.2.3 檢測季節(jié)檢測季節(jié)v 卡方檢驗卡方檢驗值值dfSig.(雙側(cè)雙側(cè))Pearson卡方17.45930.001有效案例中的N2840表表 30 六類主要六類主要食品食品分分季節(jié)季節(jié)列聯(lián)匯總表列聯(lián)匯總表第一季度第一季度第二季度第二季度第三季度第三季度第四季

28、度第四季度合格不合格合格不合格合格不合格合格不合格飲品飲品計數(shù)777 15 254 7 1225 29 988 25 期望774 18 255 6 1225 29 990 23 果蔬果蔬計數(shù)861 33 275 29 518 25 1049 50 期望851 43 289 15 517 26 1046 53 糧食糧食計數(shù)809 47 721 32 1171 16 1156 36 期望828 28 728 25 1148 39 1153 39 焙烤食焙烤食品品計數(shù)476 130 205 17 989 18 491 31 期望556 50 204 19 923 84 479 43 肉蛋制肉蛋制品品

29、計數(shù)426 62 391 36 280 16 766 76 期望443 45 388 40 269 27 746 78 水產(chǎn)品水產(chǎn)品計數(shù)98 2 222 10 373 22 330 18 期望95 5 221 11 376 19 331 17 表表 31 六類主要食品與季節(jié)因素影響的卡方檢驗表六類主要食品與季節(jié)因素影響的卡方檢驗表食品種類食品種類監(jiān)測總數(shù)監(jiān)測總數(shù)卡方值卡方值自由度自由度P值值飲品33200.88130.830果蔬284017.45930.001糧食398829.59630.000焙烤食品2357197.44730.000肉蛋制品205312.53730.006水產(chǎn)品10752.

30、43430.4875.2.3 檢測季節(jié)檢測季節(jié)v逐類分析逐類分析 飲品的季節(jié)因素相差很小，即受季節(jié)影響較小，原因可能是飲品的包裝較好，保質(zhì)期較長，因此不同的季節(jié)對飲品質(zhì)量影響不大 果蔬受季節(jié)影響較大，第一季度質(zhì)量較好，第二季度質(zhì)量略差。原因可能在于第二季度氣溫較高，使得某些預(yù)包裝果蔬類食品更容易變質(zhì)。 糧食受季節(jié)影響也較大，但與蔬菜不同的是，其第一、二季度質(zhì)量較差，第三季度質(zhì)量較好，原因在于第三季度一般為糧食收割季，所以糧食制品質(zhì)量較好。 焙烤食品在列聯(lián)分析中表現(xiàn)出受季節(jié)因素影響顯著，其中，第一季度質(zhì)量較差，第三、四季度質(zhì)量較好。 肉蛋制品和焙烤食品情況類似，第一季度較差，第二、三季度質(zhì)量好。

31、表表 28 飲品等五飲品等五類食品分季度合格率匯總表類食品分季度合格率匯總表飲品飲品果蔬果蔬糧食糧食焙烤食品焙烤食品肉蛋制品肉蛋制品10-298.5%98.9%99.4%97.5%10-3100%92.3%99.7%99.7%85.7%10-493.9%87.3%99.4%97.5%96.8%11-196.7%89.1%92.2%66.9%73.3%11-2100%94.6%100%11-397.5%87.6%100%96.3%95.9%11-498.6%94.3%95.0%88.8%86.5%12-198.8%99.5%99.2%100%93.3%12-295.7%86.5%98.3%91.

32、3%88.6%12-397.1%100%95.8%98.4%96.0%12-499.4%99.6%95.9%92.4%95.2%13-199.5%95.7%100%100%92.1%13-291.7%93.0%85.6%5.2.3 檢測季節(jié)檢測季節(jié)第一季度第一季度第二季度第二季度第三季度第三季度第四季度第四季度飲品飲品略好略好一般一般一般一般一般一般果蔬果蔬較好較好較差較差一般一般一般一般糧食糧食較差較差略差略差很好很好一般一般焙烤食品焙烤食品很差很差略差略差很好很好較好較好肉蛋制品肉蛋制品較差較差較好較好較好較好一般一般5.2.3 檢測季節(jié)檢測季節(jié)v各種食品的質(zhì)量多多少少都與季節(jié)有著各種食品

33、的質(zhì)量多多少少都與季節(jié)有著關(guān)聯(lián)關(guān)聯(lián)v不同食品種類與季節(jié)關(guān)系不盡不同食品種類與季節(jié)關(guān)系不盡相同相同v在加強相關(guān)的監(jiān)管的過程中，需要在加強相關(guān)的監(jiān)管的過程中，需要對癥下藥對癥下藥v根據(jù)每類食品受季節(jié)影響的特點，制定相應(yīng)根據(jù)每類食品受季節(jié)影響的特點，制定相應(yīng)的抽檢方案。的抽檢方案。5.2.4 食品銷售地點（即檢測地點食品銷售地點（即檢測地點）v方法：列聯(lián)分析方法：列聯(lián)分析 深圳市按城區(qū)劃分，共9個城區(qū) 各個城區(qū)食品合格率表35寶 安寶 安區(qū)區(qū)福 田福 田區(qū)區(qū)光明新光明新區(qū)區(qū)龍 崗龍 崗區(qū)區(qū)羅 湖羅 湖區(qū)區(qū)南 山南 山區(qū)區(qū)坪山新坪山新區(qū)區(qū)鹽 田鹽 田區(qū)區(qū)龍華新龍華新區(qū)區(qū)合格數(shù)合格數(shù)1596950175

34、 15158778302122428不合格不合格數(shù)數(shù)8625106425141381合格率合格率94.9%97.4%94.6%95.9%97.2%98.3%94.2%96.8%88.9%5.2.4 食品銷售地點（即檢測地點食品銷售地點（即檢測地點）v 值值df漸進漸進Sig.Pearson卡方31.36980.000有效案例中的N6561v檢驗檢驗結(jié)果結(jié)果5.2.5 食品包裝食品包裝v方法：列聯(lián)分析方法：列聯(lián)分析 食品按包裝方式，分為包裝和散裝 對每一類食品分別進行分析，減少因素間干擾 某些種類食品由于不合格的散裝檢測數(shù)為0，無法進行列聯(lián)分析。 僅對食用油、飲品、果蔬、甜味料、糧食、肉蛋制品、

35、水產(chǎn)品和調(diào)味品進行分析。表表 37 各種類食品的合格各種類食品的合格與否與否 與與不同包裝的計數(shù)與分析表不同包裝的計數(shù)與分析表合格、包合格、包裝裝合格、散合格、散裝裝不合格、不合格、包裝包裝不合格、不合格、散裝散裝Sig.乳制品乳制品132630/食用油食用油232143635.070.00飲品飲品3816501637.630.00果蔬果蔬24963830481299.820.00甜味料甜味料23655320.470.49糧食糧食471181123570.510.48焙烤食品焙烤食品117410/肉蛋制品肉蛋制品1923273632286.990.00水產(chǎn)品水產(chǎn)品1652423320.500.

36、00調(diào)味品調(diào)味品1264862018.330.00其他其他2081932325/總和總和152361120582284243.230.005.2.5 食品包裝食品包裝v分析結(jié)果分析結(jié)果 除了糧食與甜味料外，其他種類食品的合格與否均顯示出與食品包裝相關(guān)，即包裝的食品合格率更高，并且相關(guān)性非常顯著。 總和數(shù)據(jù)也顯示出，食品包裝對食品合格的相關(guān)性頗為顯著。 散裝的食品衛(wèi)生安全情況較為嚴重，需要加強監(jiān)管。5.2.6 綜合因素分析綜合因素分析v 5.2.6 綜合因素分析綜合因素分析表表32食品檢測地點虛擬變量編號表食品檢測地點虛擬變量編號表v j值對值對應(yīng)的應(yīng)的城區(qū)城區(qū)龍崗龍崗區(qū)區(qū)羅湖羅湖區(qū)區(qū)寶寶安安區(qū)

37、區(qū)鹽田鹽田區(qū)區(qū)南山南山區(qū)區(qū)福田福田區(qū)區(qū)光明光明新區(qū)新區(qū)坪山坪山新區(qū)新區(qū)龍華龍華新區(qū)新區(qū)j1234567895.2.6 綜合因素分析綜合因素分析v 表表40 產(chǎn)地、包裝、季節(jié)和檢測地點與食品質(zhì)量的產(chǎn)地、包裝、季節(jié)和檢測地點與食品質(zhì)量的Logistic多多重共線性診斷表重共線性診斷表多重共線性數(shù)據(jù)多重共線性數(shù)據(jù)容忍度VIF季節(jié)0.9111.098產(chǎn)地0.7861.272檢測地0.9961.004包裝0.7451.3425.2.6 綜合因素分析綜合因素分析v 5.2.6 綜合因素分析綜合因素分析 可以采用最大似然估計法（MIE）。最大似然法的思想是先對初始值進行估計，根據(jù)算法變動回歸參數(shù)，進行迭代運

38、算，使得對數(shù)似然值增大，由此估計函數(shù)。v檢驗檢驗Logistic回歸的顯著回歸的顯著性性 用Hosmer 和 Lemeshow 檢驗 零假設(shè)是模型能夠很好的擬合。 當sig.值0.05，因此不能否認零假設(shè)，即模型能夠很好的擬合整體，不存在顯著的差異。5.2.6 綜合因素分析綜合因素分析v 表表 產(chǎn)地、包裝、季節(jié)和檢測地點與食品質(zhì)量的產(chǎn)地、包裝、季節(jié)和檢測地點與食品質(zhì)量的Logistic回歸分析回歸分析結(jié)果結(jié)果Variables in the Equation BS.E.WalddfSig.Exp(B)Step 1a季節(jié) 25.9983.000 季節(jié)(1)-.541.4121.7191.190.

39、582季節(jié)(2).438.457.9201.3381.550季節(jié)(3)1.301.5475.6701.0173.674產(chǎn)地.000.000.4511.5021.000檢測地 24.0688.002 檢測地(1)1.099.5733.6721.0553.001檢測地(2)-.761.3215.6261.018.467檢測地(3)-.639.5881.1801.277.528檢測地(4).837.5752.1201.1452.309檢測地(5).113.424.0711.7911.119檢測地(6).3711.048.1251.7231.449檢測地(7)-.590.656.8101.368.55

40、4檢測地(8)34.4571.357E8.00011.0009.213E14包裝.770.2966.7421.0092.159Constant3.180.495 41.2391.000 24.045 模型的Hosmer 和 Lemeshow 檢驗的結(jié)果如下表： 由上表知，sig.值0.05，因此不能否認原假設(shè)零假設(shè)，即模型能夠很好的擬合整體，不存在顯著的差異。由于數(shù)據(jù)量較大，對3000多個數(shù)據(jù)回歸得到如此結(jié)果，已可說明回歸顯著。Hosmer and Lemeshow TestHosmer and Lemeshow TestStepChi-squaredfSig.19.6468.291Varia

41、bles in the EquationVariables in the Equation BS.E.WalddfSig.Exp(B)Step 3c季節(jié) 25.8103.000 季節(jié)(1)-.563.4111.8781.171.570季節(jié)(2).395.453.7601.3831.485季節(jié)(3)1.281.5455.5261.0193.602檢測地 24.1408.002 檢測地(1)1.094.5733.6441.0562.988檢測地(2)-.763.3215.6471.017.466檢測地(3)-.656.5881.2471.264.519檢測地(4).842.5752.1491.14

42、32.322檢測地(5).112.424.0701.7921.119檢測地(6).3791.048.1311.7181.461檢測地(7)-.567.655.7511.386.567檢測地(8)34.1041.112E8.00011.000 6.471E14包裝.865.26610.5981.0012.374Constant3.208.49442.2181.00024.720 模型的Hosmer 和 Lemeshow 檢驗的結(jié)果如下表： 由上表知，sig.值0.05，因此不能否認原假設(shè)零假設(shè)，即模型能夠很好的擬合整體，不存在顯著的差異。相比之前全部變量都進入建模結(jié)果更好。Hosmer and

43、Lemeshow TestHosmer and Lemeshow TestStepChi-squaredfSig.34.8508.7735.2.6 綜合因素分析綜合因素分析 總體相關(guān)性總體相關(guān)性 根據(jù)wald值判定，季節(jié)，抽檢地，包裝和產(chǎn)品是否合格表現(xiàn)出相關(guān)性。且季節(jié)抽檢地包裝 季節(jié)因素季節(jié)因素: 第四季度與食品質(zhì)量的相關(guān)性顯著，且相較第一季度質(zhì)量要好 不同種類食品的影響情況也需要分開考慮。 檢測地點因素檢測地點因素: 是否是寶安區(qū)食品與食品質(zhì)量的相關(guān)性顯著，且相對龍崗區(qū)要差。 其他城區(qū)的食品質(zhì)量良莠不齊，且分布更為均勻5.2.6 綜合因素分析綜合因素分析 食品包裝食品包裝: 對食品質(zhì)量起顯著

44、影響,有包裝的食品的合格率更高。 食品產(chǎn)地食品產(chǎn)地 在回歸分析中相關(guān)性不強，分析如下: 如前文所述，由于數(shù)據(jù)信息不齊全，在進行回歸分析時需要篩選出信息齊全的數(shù)據(jù)，那么在此篩選中，因為篩選出的數(shù)據(jù)與原始數(shù)據(jù)的信息分布無可避免地會有所出入，食品產(chǎn)地對食品質(zhì)量的影響或許會有所減弱。 相對其他因素而言，有分析中可以看出，從全國角度來看，食品產(chǎn)地的因素影響不大。5.3 抽樣方法的改進抽樣方法的改進改進模型改進模型5.3.1 市場模型的簡化市場模型的簡化5.3.3 產(chǎn)品抽檢方法改進對比測試產(chǎn)品抽檢方法改進對比測試5.3.4 權(quán)重系數(shù)的權(quán)重系數(shù)的，的確定的確定5.3.5 實證實證數(shù)據(jù)對比模型數(shù)據(jù)對比模型5.

45、3.6 實證實證模型結(jié)果模型結(jié)果5.3.1 市場模型的簡化市場模型的簡化市場中只有已知的若干類食品，且每類食品下有若干種食品。市場中只有已知的若干類食品，且每類食品下有若干種食品。同種食品的檢測項目相同且最終合成為一個指標，我們對食品的檢測值同種食品的檢測項目相同且最終合成為一個指標，我們對食品的檢測值進行如下處理，將檢測值統(tǒng)一除以標準值，特別地，如果是微生物檢測進行如下處理，將檢測值統(tǒng)一除以標準值，特別地，如果是微生物檢測項目，則將檢測值與標準值分別取對數(shù)后再做除法，我們稱這個過程為項目，則將檢測值與標準值分別取對數(shù)后再做除法，我們稱這個過程為“指標歸一化指標歸一化”過程。過程。同一批次的食

46、品質(zhì)量相同同一批次的食品質(zhì)量相同改進模型改進模型v目標目標 在總的抽檢次數(shù)不變的情況下提高抽檢的效率，找到更多的不合格產(chǎn)品。v方法方法 合理分配每種產(chǎn)品的抽檢率，提高抽檢效率 抽檢率與兩個因素有關(guān): 之前抽檢批次結(jié)果的方差 抽檢結(jié)果方差越大，說明該種類產(chǎn)品質(zhì)量波動情況較大，則我們將適當增加我們的抽檢率 之前抽檢批次結(jié)果的超標情況 抽檢結(jié)果超標情況嚴重，我們會增加該類食品的抽檢次數(shù)。改進模型改進模型v 5.3.3 產(chǎn)品抽檢方法改進對比測試產(chǎn)品抽檢方法改進對比測試v通過隨機數(shù)生成相關(guān)的數(shù)據(jù)，定義檢出率為抽檢出的不合格產(chǎn)品占市場中所有不合格產(chǎn)品的比例。對于同一份隨機生成的數(shù)據(jù)，在相同的抽檢次數(shù)下比較

47、兩種方法檢出率的大小，從而比較兩種方法的優(yōu)劣。隨機產(chǎn)品質(zhì)量生成參數(shù)表食品種類食品種類市場比重市場比重質(zhì)量隨機函數(shù)質(zhì)量隨機函數(shù)11000normrnd(0.5,0.3,1,1000);22000normrnd(0.6,0.2,1,2000);33000normrnd(0.45,0.4,1,3000);44000normrnd(0.8,0.34,1,4000);55000normrnd(0.7,0.78,1,5000);5.3.3 產(chǎn)品抽檢方法改進對比測試產(chǎn)品抽檢方法改進對比測試v 5.3.3 產(chǎn)品抽檢方法改進對比測試產(chǎn)品抽檢方法改進對比測試v v v 5.3.5 實證數(shù)據(jù)對比模型實證數(shù)據(jù)對比模型

48、v 5.3.5 實證數(shù)據(jù)對比模型實證數(shù)據(jù)對比模型v生成隨機檢測值的方法生成隨機檢測值的方法 1、對于每個批次產(chǎn)品質(zhì)量值，生成一個介于(0,1)的隨機數(shù)。 2、如果該隨機數(shù)大于該類產(chǎn)品在當前季節(jié)的合格率時，則認為該批次產(chǎn)品是不合格產(chǎn)品。并在之前該類食品在該季節(jié)的統(tǒng)計數(shù)據(jù)中隨機找到的一個不合格產(chǎn)品的檢測值，代表該批次產(chǎn)品的質(zhì)量檢測值。 3、如果該隨機數(shù)小于等于該類產(chǎn)品在當前季節(jié)的合格率時，則認為該批次產(chǎn)品是合格產(chǎn)品。同樣由中的說明，產(chǎn)品為合格時，檢測值介于(0,1)之間，因此再次在(0,1)間生成一個隨機數(shù)，作為該批次產(chǎn)品的質(zhì)量檢測值。5.3.6 實證模型結(jié)果實證模型結(jié)果v 5.3.6 實證模型結(jié)果實證模型結(jié)果v v為了測試我們方法的穩(wěn)定性，我們每次模擬為了測試我們方法的穩(wěn)定性，我們每次模擬20年，即年，即240個月的抽檢周期。共運行個月的抽檢周期。共運行1000次。次。v得到平均檢出率提高的均值水平為得到平均檢出率提高的均值水平為10.75%六、六、模型的結(jié)果分析模型的結(jié)果分析6.1 各主要食品總體安全的情況分析各主要食品總體安全的情況分析 6.3 抽檢方法改進對抽檢效率的影響分析抽檢方法改進對抽