1、食品質量安全抽檢數據分析摘要本文根據文章提出的不同問題，建立相應的數學模型，利用 matlab軟件進行求解，對食品質量進行評價和找規(guī)律以及合理抽檢方 法。針對問題一，對深圳市這三年各主要食品領域微生物、重金屬、 添加劑含量等安全情況的變化趨勢做出定量的綜合評價，首先對數據 處理按季節(jié)分為17個子樣本點的抽樣值進行分析，得到各子樣本食 品安全情況。運用層次分析法決定食品安全單位指標在綜合評價中權 重。可以得到深圳市近三年來的食品質量情況應是明顯提高。年份201120122013食品安全系數0.270.170.095等級III級II級I級針對問題二，我們先通過MATLAB對原始數據進行檢驗，對殘差

2、 向量進行分析，得到了殘差向量分析圖，剔除其中的異常點。運用 MATLAB進行編程，得到各因素的偏回方和：xix1x2x3x4x5x6x7x8Q ( x 106)0.33840.00300.36850.10090.14291.24180.14720.1963根據Q的大小可判斷各因素對食品安全系數的影響程度 ix x x x x x x x63187542結論：在食品質量影響因素中食品產地影響最大，食品加工次之， 季節(jié)影響最小，抽查地點幾乎無影響。針對問題三，根據問題一、二數據結果以定量比較評估的方法 分析了各類影響食品安全的因素及其可能造成的危害性的問題。改進 后的食品抽檢的辦法以主要食品為準

3、則層建立了層次分析（AHP）模 型，對影響食品安檢的危害性因素做出定量分析如問題二個影響因素 大小。由問題一、二結果可以通過建立抽檢模型，即改進后的規(guī)準型 抽樣檢驗模型，并以蒙特卡羅法對抽檢的全過程進行模擬，得到相對 誤差逐漸趨向于0。關鍵詞：方差與回歸分析 殘差向量分析 評價指數 層次分析問題重述“民以食為天”，食品安全關系到千家萬戶的生活與健康。隨著人們對生活 質量的追求和安全意思的提高，食品安全已成為社會關注的熱點，也是政府民 生工程的一個主題。城市食品的來源越來越廣泛，人們消費加工好的食品的比 例也越來越高，因此除食材的生產收獲外，食品的運輸、加工、包裝、貯存、 銷售以及餐飲等每一個環(huán)

4、節(jié)皆可能影響食品的質量與安全。另一方面，食品質 量與安全又是一個專業(yè)性很強的問題，其標準的制定和抽樣檢測及評價都需要 科學有效的方法。深圳是食品抽檢、監(jiān)督最統(tǒng)一、最規(guī)范、最公開的城市之一。 請下載2010年、2011年和2012年深圳市的食品抽檢數據（注意蔬菜、魚類、 雞鴨等抽檢數據的獲?。⒏鶕@些資料來討論：如何評價深圳市這三年各主要食品領域微生物、重金屬、添加劑含量 等安全情況的變化趨勢；從這些數據中能否找出某些規(guī)律性的東西：如食品產地與食品質量的 關系；食品銷售地點（即抽檢地點）與食品質量的關系；季節(jié)因素等能否改進食品抽檢的辦法，使之更科學更有效地反映食品質量狀況且 不過分增加監(jiān)管成

5、本（食品抽檢是需要費用的），例如對于抽檢結果穩(wěn) 定且抽檢頻次過高的食品領域該作怎樣的調整？二、問題分析本問題是在監(jiān)督最統(tǒng)一，最規(guī)范，最公開的城市之一深圳，通過 對食品安全抽查數據，抽樣，合格，不合格以及與近幾年數據比較的 數據比例進行研究，探討使人們追求的食品質量安全與可持續(xù)科學發(fā) 展觀加快吻合的方案。對于市場食品安全的分析，我們知道，和y有關的變量有8個， 研究y與變量x，x，x，*之間的定量關系的問題為多元回歸問1 2 3 8題。又因為許多多元非線性回歸問題都可以化為多元線性回歸問題， 所以對于本問題我們建立了多元線性回歸的數學模型。第一問：由于食品安全生成這一過程涉及到諸多環(huán)節(jié)，所以我們

6、必須考慮 食材的生產收獲外，食品的運輸、加工、包裝、貯存、銷售以及餐飲 等每一個環(huán)節(jié)皆可能影響食品的質量與安全。另一方面，食品質量與 安全又是一個專業(yè)性很強的問題，其標準的制定和抽樣檢測及評價都 需要科學有效的方法。運用統(tǒng)計分析中的方差，分析各個影響因子對 食品安全的影響程度。同時我們確定出哪些因素是主要的，哪些是次 要的。通過表格可以看出，由于抽樣的各種不同方法以及子樣空間的 不同導致的子樣均差以及子樣方差的不同，導致了食品安全系數的標 準不同。首先，考慮抽樣方法的影響。通過對抽樣數據和合格數據的觀察， 我們發(fā)現在不同抽樣方法以及不同子樣本空間下，對食品安全系數有 較大的影響。通過對整個表格

7、的分析我們把抽樣數據分了四個數據梯 度,，分別為2013年的，2012年的，以及2011年和2010年的，然后 列出四組數據，再應用方差分析模型進行分析，判斷出評價深圳市這 三年各主要食品領域微生物、重金屬、添加劑含量等安全情況的變化 趨勢。其次，我們在來考慮諸多環(huán)節(jié)對食品安全系數的影響。我們把數據分 為四個組進行分析，分別為取材環(huán)節(jié)，生產環(huán)節(jié)，流通環(huán)節(jié)，餐飲環(huán) 節(jié)。然后我們同樣應用方差分析模型來對數據進行處理，來判斷各個 環(huán)節(jié)對評價深圳市這三年各主要食品領域微生物、重金屬、添加劑含量等安全情況的變化趨勢。第二問：首先對附表1的數據進行處理，對殘差向量進行，剔除其中的異 常點。然后我們建立了多

8、元線性回歸的數學模型，并采用了最小二乘 法來估計參數。把模型寫成矩陣的形式，化簡整理得其正規(guī)方程組， 通過對正規(guī)方程組的求解，最后得到回歸方程。對于各因素對食品安全的影響程度，由于利用偏回歸平方和Q可i 以衡量每個變量在回歸中所起的作用大小（即影響程度），我們對每 個變量氣的偏回歸平方和Q進行了計算，最后把影響程度由深到淺的 各因素進行了排序。回歸方程的顯著性檢驗：事先我們并不能斷定隨機變量y與一般 變量x，x，x，x之間是否確有線性關系。在求線性回歸方程前1238線性回歸模型只是一種假設，所以在求出線性回歸方程之后，我們需 要對其進行統(tǒng)計檢驗。將總的平方和分解為回歸平方和s回和剩余 平方和S

9、剩，運用F檢驗法來判定線性回歸方程的顯著性?；鼗貧w系數的顯著性檢驗：由于回歸方程顯著并不意味著每個自變 量x，x，x，x對因變量y的影響都是重要的。而我們要找出影1238響食品安全的主要因素，即從回歸方程中剔除那些次要的、可有可無 的變量，這就需要我們對每個變量進行考察。顯然，如果某個變量對 y的作用不顯著，那么在多元線性回歸模型中，它前面的系數P .就可 以取值為零。因此，檢驗因子x.是否顯著等價于檢驗7假設 H： =0。最后再運用T檢驗法來辨別模型中哪些因子是顯著的 及從這些數據中能否找出某些規(guī)律性的東西：如食品產地與食品質量 的關系；食品銷售地點（即抽檢地點）與食品質量的關系；季節(jié)因素

10、等等。第三問：由于回歸系數之間存在相關性，當從原回歸方程中剔除一個變量 時，其他變量，特別是與它密切相關的一些變量的回歸系數就會受到 影響，剔除一個變量后，這個變量對y的影響很大部分轉加到另一個 變量對y的影響上。所以，我們對回歸系數進行一次檢驗后，只能剔 除所有不顯著因子中，值最小的，然后重新建立新的回歸方程，再對 新的回歸系數逐個進行檢驗，直到余下的回歸系數都顯著為止。三、問題假設1、土壤對取材生長狀況基本相同，從而對生成食品的安全系數 的影響相同。2、假設周圍自然環(huán)境因素除去土壤外，其他的對食品安全系數 的生成沒有影響。3、在食品生產過程中，材料都已經轉化成了食品。也就是說產生是完全的，

11、沒有發(fā)生材料的剩余。4、同類食品是相同的，也就是說抽查的食品對生成食品的安全系數影響相同。5、餐飲服務等條件始終是不變的。6、抽查是在隨機下進行，不存在任何的人為干擾。四、參數及符號說明表符號表示的意義單位備注Li第i個子樣點與樣本空間的差值個i = 1.7vi第i個空間的大小個N ij第i個子樣本第j種食品的合格%N .：第i個樣本空間第j子樣本經管理后下一 年的合格率%w.第i個樣本空間中第j子樣本總量個V ij第i個子樣本第j個食品合格率%T管理系數1/月0.1 t 0.5ti (i+1)樣本空間與子樣本之間差距個PI食品安全質量綜合評價指數1Ai準則層項目因子i = 1.4Bi項目層項

12、目因子i = 1.4C.I.一致性指標R.I.平均隨機一次性指標C.R.一致性比例指標Ri為第i年需要處理的不合格食品個i=1.10Hi為預測的i年的食品不合格數個i=1.10表二符號 說明X影響食品質量安全的因素。(i 1,2,.,8)y O各食品安全系數(a -1,2,.,38)_i相互獨立且服從同一正態(tài)分布的隨機變量(，-1,2,.,38)Q i變量X,的偏回歸平方和S 總總平方和S-回回歸平方和S 剩剩余平方和Pi待估計系數參數(i =1,2,.,8)yy的回歸值b_i參數。的最小二乘估計(i-1,2,.,8)食材生長土壤，食材的生產收獲，食品的運輸、加工、包裝、貯存、 銷售以及餐飲八

13、個影響因素五、模型的分析、建立與求解設K取尸個水平，瓦,在水平H下總體x,服從正態(tài)分布/ = r ,這里未知，,4可以互不相同，但假定x,有相同的方差。又設在每個水平H下作了 0次獨立試驗，即從中抽取容量為勺的樣本，記作X ，/=, X服從q2),，= 1,尸,/ = 1,0且相互獨立。將這些數據列成表1 (單因素試 驗數據表)的形式。表1單因素試驗數據表4邑112砒a2%X22f4命&將第，行稱為第，組數據。判斷刀的尸個水平對指標有無顯著影響，相當于要作以 下的假設檢驗H。:巧=/弓=“,.；: pH，,.不全相等由于的取值既受不同水平其的影響，乂受其固定下隨機因素的影響，所以將它 分解為

14、rr夕= 一a，h =，a,=出一i = l,/即=1=1產是總均值，名是水平H對指標的效應由（1）、（2）模型可表為 孔=片+名+與 7=1%N（O,cH）, z = 1,y = 1,.,勺 原假設為（以后略去各選假設）H0 :% = OL之=叫_ =。 1.2統(tǒng)計分析 記 、1 ! E，=二%，元=二/與 勺 J=1n j =l J=1a是第i組數據的組平均值，元是總平均值?？疾烊w數據對元的偏差平方和=1 7=1 經分解可徂 記，=支勺（一元）27=1 3=1 J=1 則 S S（9） S是各組均值對總方差的偏差平方和，稱為組間平方和；3反是各組內的數據對均值偏 差平方和的總和。3反映4

15、不同水平間的差異，3丘則表示在同一水平下隨機誤差的 大小。注意到文（孔-xj2是總體NgS）的樣本方差的呵-1倍，于是有 7=1工（-打偵zW-D 7=1 由;分布的可加性知 TOC o 1-5 h z f r、I ：=i7即s/b2 /(I)且有毛=(77-*(10)對s作進一步分析可得ESa = (r-lK + 叫;(11)J=1當H成立時=(1頑(12)可知若反。成立，5只反映隨機波動，而若H。不成立，那它就還反映了 乂的不同水平的效應。單從數值上看，當H。成立時，由(10)、(12)對于一次試驗應有SE/(n-r)而當H。不成立時這個比值將遠大于1。當丑。成立時，該比值服從白由度昭=-

16、1,n2 = (w- r)的F分布,即F =尸(井_1 /_尸)(13)為檢驗H,給定顯著性水平】，記尸分布的分位數為Fi偵一為一心),檢驗 規(guī)則為F 5M-L(i)時接受丑。，否則拒絕。以上對SA,SE,ST的分析相當于對組間、組內等方差的分析，所以這種假設檢驗方法稱 方必分析。1.3方差分析去將試驗數據按上述分析、計算的結果排成表2的形式，稱為單因素方差分析表 (Matlab中給出的方差分析表)表2單因素方力分析表方差來源平方和自由度均方1 pv分位數概率因素4r-1r-1Fr-l.n-r)Pv誤差島n rnr總和ST即一1最后 列給出大/值的概率Ff相方差分析一般用的顯著性水平是：取1

17、= 0.01,拒絕丑0,稱因素，的影響（或， 各水平的差異）非常顯著；取1 = 0.01,不拒絕H。，但取a = 0.05，拒絕稱因 素刀的影響顯著；取& = 0.05,不拒絕H。，稱因素4無顯著影響。1.4 Matlab 實現Matlab統(tǒng)計工具箱中單因素方差分析的命令是anovaL若各組數據個數相等，稱為均衡數據。若各組數據個數不等，稱非均衡數據。（1）均衡數據處理均衡數據的用法為：p=anoval（x）返回值p是一個概率，當p時接受H。，X為陽X尸的數據矩陣，X的每一列是一個 水平的數據（這里各個水平上的樣本容量nt = m ）o另外，還輸出一個方差表和一個 Box 圖。問題一的模型問題

18、一為：如何評價深圳市這三年各主要食品領域微生物、重金屬、添加劑含量等安全情況的變化趨勢問題的分析與模型的建立對深圳市這三年各主要食品領域微生物、重金屬、添加劑含量等 安全情況的變化趨勢做出定量的綜合評價，首先對數據處理按季 節(jié)分為17個子樣本點的抽檢值進行分析，得到各子樣本食品質量 安全情況。運用層次分析法決定食品安全單位指標在綜合評價中 權重。1）構造各季節(jié)食品安全質量綜合評價的階梯層次結構在構造深圳市食品領域綜合評價模型時，考慮地理位置，人員結 構主要考慮老中青人員比例，經濟結構等4項食品安全指標的綜合效 果。構造其綜合評價的遞階層次結構如下圖：圖：食品綜合評價遞階層次結構2）建立V-A層

19、判斷矩陣建立V-A層判斷矩陣的依據是：考慮A層的四個因素對食品質量 安全綜合評價的重要性時，認為年均指標A1最重要，且A1比微生 物、重金屬、添加劑含量等值重要。而在微生物、重金屬、添加劑含 量等值中，而各類食品都用到添加劑，所以，添加劑比微生物、重金 屬等重要。根據判斷矩陣標度內容及其含義，可得V-A層判斷矩陣A：A1 A2A3A4-1233 一A11/213/23/2A2A =1/32/311A31/32/311 JA43）A-B層判斷矩陣該層的矩陣元素Z確定如下：將每條準則AA2、A3、a4下各指標抽查值用分指數公式（1.1）標準化，其中C是在A.準則下第j個指標抽查值的均值（或最大 值

20、），C是第jJ個抽查指標的最大允許不合格值。對于每條準則下，令 3kj = max3 I j = 1，2,3,4i = 1,2,3,4 計算3 = min3 I j = 1，2,3,4. = 1,2,3,4 氣二（3廣3）/9ljij（1.2）rr =st3 -3 .-sA （ s t =A.（1.3）對于任意的3當 I r 1 1當, -1與3相比， 取3二1st 取廣戲 取 3：二：1/ rt其判斷矩陣元素3由以下確定st外 b = 1/ b=按式（1.1）可得相應于表1.1中各指標抽查值的標準化值由式子（1.2）（1.3）（1.4）得A-B層的判斷矩陣4）用方根法計算出判斷矩陣A , B

21、（）的歸一化特征相量：VA、VB（n）5）計算判斷矩陣B（i）,B（2）,B（3）的最大特征值人 并進行一致性檢驗max一致性指標人=w啤Lmax nwi=1iC.I. = y 其中，n為判斷矩陣的階數n -1CI計算一致性比例C.R.=工其中R.I.為平均一致性指標，當C.R. 0.1 R.I.時判斷矩陣具有好的一直性。6）食品質量的綜合評價指數的計算PI =岌 4 V . V G）j=1 i=1模型的求解及結果以2013年第一至四期食品安全抽樣檢驗產品抽查樣本的32種食 品的抽查值為例，求該抽查樣本的食品質量安全綜合評價指數。根據 題目給出的32種食品抽樣數據值，統(tǒng)計得到下面的平均值表：表

22、1.12013年第一至四期食品安全抽樣檢驗產品指標統(tǒng)計表表項目2010201120122013平均值8.25610.109242.43210.18286微生物平均值8.20750.111883.0250.14688重金屬平均值8.20.10492.10.0675添加劑平均值8.38120.106411.41250.3125按式（1.1）可得相應于表1.1匚尸各指標抽查值的標準化值如表1.2所示。表1.2各指標抽查值的標準化值因子2010201120122013A10.688010.163850.121610.068571A20.683960.167810.151250.055077A30.68

23、3330.157350.1050.025312A40.698440.159610.0706250.11719式子（1.2）（1.3）（1.4）得A-B層的判斷矩陣分別為:B（1）：B1B2B3B417.61578.22949B10.13131111.3843B20.12152111B30.111110.7223711B4B(2)：B1B2B3B417.38667.62379B10.13538111.6134B20.13117111.3763B30.111110.619820.726571B4B(3):B1B2B3B417.19417.91019B10.139111.8059B20.126421

24、11.0899B30.111110.553740.91751B4B(4):B1B2B3B417.724498.3325B10.1294611.27561B20.111110.7839311B30.12001111B44)用方根法計算出判斷矩陣A ,B(i),B(2),B(3)的歸一化特征相量:V =(V , V , V , V ) = (0.461539 0.230769 0.153846 a A 今 少 少0.153846);V二(V，V, 0.0802843) 12V (2) = ( V (2)，V (2)V二(V，V, 0.0802843) 12V (2) = ( V (2)，V (2)

25、，B、BiB20.0724441)V(2)，V(2) ) = (0.72337 0.104756 0.0994301B3B4V=(V，B、 BV=(V，B、 B10.0749515)V=(B0.0814922V，VBB0.08826822)VV=(B0.0814922V，VBB0.08826822)V，V)=0.734622 B3B40.0956177其中，的4個分量分別表示A層的4個因子(A1平均值、A2微生 物平均值、A3、重金屬平均值、A4、添加劑平均值)在綜合評價中的各自的權重。VB(1)的4個分量分別表示只考慮年均值A1的條件下4 個指標各自應占的權重；VB(2)、VB(3)、VB(

26、4)類推。5)計算判斷矩陣B(i),B(2),B(3)的最大特征值人并進行一致性檢max驗人 =況 _L人i= 6.0i2,人2 = 6.08,人3 = 6.09 一致性指標maxnWmaxmaxmaxi=iiC.I. = 其中，n為判斷矩陣的階數n -1CI計算一致性比例C.R.=史土其中R.I.為平均一致性指標，查表R.I.R.I. =0.9，故C.R.=0 0.1C.R.（2）= 0 0.1C.R.（3）= 0 0.1當C.R.0.1時判斷矩陣具有好的一致性。C.R.（4）= 0 0.1的綜合評價指數的計算廣.2610552.266）市場食品質量安匕 PI =岌 4 / . VB ）的綜

27、合評價指數的計算廣.2610552.26j=1 i=1根據深圳市食品安檢局標準處推薦的評價標準（GB3838-2002 食品安檢質量標準以及食品質量分級）為參考依據，見下表 1.3表1.3深圳市食品質量系數（PI）分級質量安全值級別食品質量評價0.80劣V嚴重劣食品這樣就得到了 2013年第一至四期食品安全抽樣檢驗產品抽查樣本的 食品安全的優(yōu)劣狀況，并做出定量的綜合評價。其它年份的食品安全 情況同理可求。求解程序見附錄1。表1.4為求得的各個抽查樣本的食品質量安 全狀況和定量綜合評價。其中C.R.為一致性指標（當C.R.0.1時有好 的一致性）詳細的計算結果見附錄2。2013年第一至四期食品安

28、全抽樣檢驗產品抽查樣本的食品安全情況 的綜合評價可由加權平均法來作出判斷。表1.4 2013年各抽查產品合格評價表抽查產品食品綜合評 價指數C.R.（食品質量安全系數）超市速食面制 品0.2651150.00000000 000.000000000 00.00000008830.0000000000預包裝調味0.2659770.00000000 000.00000008830.00000008830.0000000000預包裝即食小 食品0.2846380.00000008830.000000000 00.00000008830.0000000000一次性紙杯塑 料杯0.2486960.000

29、00000 000.00000008830.00000000 000.0000000883餅干及膨化食 品0.258630.00000000 000.000000000 00.00000000 000.0000000000速凍面制品0.2478460.00000008830.00000008830.00000000 000.0000000000巧克力及其制 品0.2610550.00000000 000.000000000 00.00000000 000.0000000000乳制品及含乳 制品0.3574770.00000000 000.000000000 00.00000000 000.00

30、00000000酒店奶制品0.2901320.00000008830.00000008830.00000000 000.0000000000大米0.3788740.00000000 000.00000008830.00000008830.0000000000食用油0.2380920.00000000 000.00000008830.00000017660.0000000883廣式燒鹵熟肉 制品0.3214150.00000000 000.000000000 00.00000000 000.0000000000煙熏肉制品0.2946080.00000008830.000000000 00.000

31、00000 000.0000000883預包裝茶葉0.2751090.00000000 000.000000000 00.00000008830.0000000883大型企事業(yè)單 位食堂0.8441380.00000000 000.00000008830.00000000 000.0000000000速凍湯圓食品0.2857430.00000008830.000000000 00.00000000 000.0000000000幼兒園食堂食 品0.2705070000.0000000883其中權值七可以反映不同食品質量對整個深圳市食品質量的 影響的不同作用：我們采用目標應用最多的因子實測值與標準

32、值為雙重判定依據的賦 權方法一超標指數法或指數賦權法，即根據各食品的質量分級指數 C來確定權重。其公式為：L(i=1,2, n)乎Iii=1L(i=1,2, n)乎Iii=1I = ioi式中PI為各食品質量安全分級的指數，C為PI對應分級指數 的最大允許值。將我們計算出的PI值代入式1.o5計算得i47PIPI = t a j = 0.317 x x x x x x x 63187542結論：在食品質量影響因素中食品產地影響最大，食品加工次之， 季節(jié)影響最小，抽查地點幾乎無影響。問題三：規(guī)準型抽樣檢驗“民以食為天，食以安為先”，食品安全和衛(wèi)生問題一直是關系國計民生 的重大問題。根據問題一、二

33、數據結果以定量比較評估的方法分析了各類影響 食品安全的因素及其可能造成的危害性的問題。改進后的食品抽檢的辦法以主 要食品為準則層建立了層次分析（AHP）模型，對影響食品安檢的危害性因素 做出定量分析如問題二個影響因素大小。由問題一、二結果可以通過建立抽檢 模型，即改進后的規(guī)準型抽樣檢驗模型，并以蒙特卡羅法對抽檢的全過程進行 模擬，得到相對誤差逐漸趨向于0。就目前“營養(yǎng)強化.（查看更多）面粉”如 何抽檢問題，本題在將各檢測項目與各類食品安全影響因素一一對應的基礎上， 建立逐步優(yōu)化模型，最終得出一次全面檢驗的最小總成本。依據問題一二結果 通過參照計量統(tǒng)計檢驗程序給出了一套完備的連續(xù)多次跟蹤抽檢策略

34、，建立更 加切合實際的多目標規(guī)劃的抽檢策略。最終得到當檢測數量達到100以上時， 該部門進行抽檢工作時可達到檢測可靠性高、成本低、工時少的目的。故針對 至少720個總批次的食品本題最終建議抽檢125批次。綜上，本文較好地解決 了當前食品安全問題中存在的矛盾關系。不僅提高了消費者的滿意度，也兼顧 了生產者的利益，具有較高的社會經濟價值。6模型的評價6.1模型的優(yōu)點本文對于各種因素對于食品質量安全的影響建立了多元線性回 歸模型，全面綜合考慮了各個方面的因素，避免了單一因素分析的不 準確性，得出了合理的數學模型。并且通過各因素的顯著性分析，找 到了影響食品質量安全的主要因素，較符合實際情況，模型可靠

35、，并 且模型相對簡單，利于操作；該方法不僅適用于本題，也適用于其他 方面的數據預測，有實際背景，可運用于實踐，具有廣泛適用性。6.2模型的缺點本文忽略了除了所給因素之外的因素對食品安全的影響，與實際 問題存在偏差。同時是在假設各因素相互獨立的情況下對食品質量的 影響進行分析，可能會導致誤差7模型的改進與推廣應用統(tǒng)計推斷先判斷食品的安全性，又有中心極限定律假設食品的合格服從正 態(tài)分布使其食品安全量化，并且量化成數據。比如超過量化后的數據為不合格， 沒有超過為合格。各個國家的食品安全等級不一樣，反映深圳市也是如此。運 用2010年、2011年和2012年深圳市的食品抽檢數據統(tǒng)計量化后的數據分別作為

36、 各類產品安全等級。生活食品質量標準是指按特定生產方式生產，并經國家有 關的專門機構認定，或許使用綠色食品標志的無污染、無公害、安全、優(yōu)質、 營養(yǎng)型的食品。在許多國家，綠色食品又有著許多相似的名稱和叫法，諸如“生 態(tài)食品”、“自然食品”、“藍色天使食品”、“健康食品”、“有機農業(yè)食品”等。影 響因素有原料生長土壤環(huán)境及食品的生產收獲，此外食品的運輸，加工，包裝， 儲存，銷售，以及餐飲等每個環(huán)節(jié)皆可能影響食品的質量與安全。運用統(tǒng)計分 析中的方差，分析各個影響因子對食品安全運用統(tǒng)計分析中的方差，分析各個 影響因子對食品安全的影響程度的影響程度。最終使人們心中無污染、無公害、 安全、優(yōu)質、營養(yǎng)型的食

37、品與市場食品安全等級數據差極限為零。相反如果人 們認識的食品安全數據比市場標準化低，或沒有很好認識安全食品概念及本質。 我們可以通過媒體，廣告，食品安全講座等。讓人們認識和提高對安全食品實 質。如果人們心目中的標準已經是市場安全等級，若是的話，就可以得出合理 的結論：人們追求更加持續(xù)性，更加科學化的生活質量標準。正好和可持續(xù)科 學發(fā)展觀相吻合。模型中得到最優(yōu)回歸方程的方法是從包含全部變量的回歸方程 中逐次剔除不顯著因子，這種方法是在不顯著因子不多時采用，當不 顯著因子較多時，則工作量將會相當大，因為每剔除一個變量就得重 新計算回歸系數。鑒于以上問題，我們引入了逐步回歸分析的方法，它的基本思想

38、是將因子一個個引入，引入因子的條件是，該因子的偏回歸平方和經 檢驗時顯著的。同時，每引入一個新因子后，要對老因子逐個檢驗， 將偏回歸平方和變?yōu)椴伙@著的因子剔除。這種方法不需要計算偏相關 系數，計算較簡便，并且由于每步都作檢驗，因而保證了最后所得的 方程中所有因子都是顯著的。若回歸方程是擬合好的，就可以進一步利用它來進行預報和控 制。預報問題，用統(tǒng)計數學的語言來說就是一個區(qū)間估計問題。在建 立氣象預報、地震預報、自動控制等數學模型時，都可以用到本文的 模型。8參考文獻馬新民，王逸迅.概率與數理統(tǒng)計M.北京：機械工業(yè)出版社， 2010.劉衛(wèi)國.MATLAB程序設計與應用M.北京：高等教育出版社，2

39、009.茹詩松.回歸分析及其試驗設計M.上海：華東師范大學出版社，1986.陳巧玲，沈曉斌，長江水質污染預測的灰色模型，泉州師范學院學報(自然科學)，2005.11孫靖南，鄒志紅，任廣平.模糊綜合評價在天然水體水質評價中的應用研究,環(huán)境污染治理技術與設備,2005 ,6 (2) :45248.李小燕灰色神經網絡預測模型的優(yōu)化研究, HYPERLINK /p-67589 /p-67589378.html鄧聚龍，灰理論基礎M，武漢:華中科技大學出版社,2002.傅立，灰色系統(tǒng)理論及其運用M，北京:科學技術文獻出版社,1992.劉榮珍，趙軍，模糊評價模型在長江水質評價中的應用，蘭州交通大學學報,20

40、07.1212C.Henry Edwards &David E.penney微分方程及邊值問題，北京：清華大 學出版社，2007附錄：一數據統(tǒng)計名稱抽取 合格 不合格2012.26即食蔬 菜學校學 生餐 娛樂場 所食品 水發(fā)水 產潮式打 冷食品 自制米 面油炸食 品即食水 果快餐食 品城中村 快餐食 品中央廚 房半成 品餐飲火 鍋湯底 學校食 堂面制 品西式菜 肴餐前小 吃食品 餐飲領 域鹽焗 類食品 工業(yè)區(qū) 餐飲 餐飲領 域乳及 其乳制2002000300300030255343406054634340606001001000150150015015006160183830100871320

41、019822001928880100100027270品69654555322001991797366564165623一次性 餐飲具 腌臘肉 制品 餐飲領 域茶葉 2012.25 餅干干果及 堅果食 品預包裝 肉制品動物性水產干制 品、藻類 干制品 交通樞紐 商鋪預包 裝食品8019870189109糕點面包1009010預包裝龜27270苓膏茶葉1201200焙烤咖啡27270及其制品飲料酒47470乳制品及1461460含乳制品火鍋原料20016238食品罐頭50482調味品100964塑料奶瓶橡膠16160奶嘴大米1501500食用植44440物油火鍋底料產品27252腌臘肉25250

42、制品餐飲自2012.242002000. 制飲料餐飲大1001000米餐飲食 用油100982餐飲西式菜肴100973餐飲生食水產1361333品企事業(yè)單位食堂分裝1001000菜肴.飯盒羅湖區(qū).龍崗區(qū)18162餐飲肉丸食用植2012.2380773物油腐竹及淀粉制55532品批發(fā)商批發(fā)環(huán)50491節(jié)糧油士多店預包裝40400食品乳制品及含乳53530制品批發(fā)商批發(fā)環(huán)節(jié)調味100919品集體食2012.223002991. 堂食品快餐食 品2002000使用冰100982粉絲腐100991竹即食果-H-71710疏2012.212012.(20)2012.19餐前小 吃 涼茶 即食蛋 制品 即

43、食熟 肉制品 即食水 產品 方便面 涼茶 小麥粉 茶葉 乳制品 及含乳 制品 生產領 域月餅 流通領 域月餅 餐飲領 域月餅 瓶裝飲 用水 固體飲 料 飲料食 品 飲料酒 淀粉及 淀粉制 品 冷凍飲 品 速凍米 面制品 乳制品 及含乳 制品 餐具洗 滌劑流通領 域暢銷 大米1321257484801010035323262333838061601393901009825050013212936058211811803533250491807826663379772109127270494904040077770食品和食品用2012.1875573817相關產品預包裝2012.1750491罐頭

44、谷物粉類制成80791品食糖60600蛋制品50482食用植86860物油乳制品及含乳50500制品福田區(qū)學校周37370邊食品餐飲茶2012.16葉2001991涼拌菜1531530現配奶茶2001964米面制 品917912沙縣小吃1691690商業(yè)中心餐飲48471原料2012.15市場4生產領2012.1411110域粽子流通領域粽子16151餐飲領域粽子88826餐飲酒2012.13類97943生食水1261242主口 產品沙津2001991壽司1001000餐飲原料飯盒120501205000煎炸用食物油100991保鮮膜學校食50446堂米油南山區(qū)商業(yè)中4004000心餐飲原料r

45、e 預包裝11011002012.12豆制品預包裝26224餅干食 品80737膨化食 品80782糖果40391濕紙巾白酒洋20200酒學校周37343邊小賣鋪食品63612葡萄酒乳制品27261及含乳 制品鹽焗食505002012.11品45414調味品100964油炸食 品89890小麥粉福田區(qū)商業(yè)區(qū)1001000專項執(zhí)法抽查494902012.（10）蜂產品46442調味品796811炒貨38371餐巾紙35350十貨45441乳制品 及含乳 制品乳制品505002012.9及含乳制品40400巧克力及制品50500預包裝果凍50500預包裝水產品40355餐飲乳2012.8制品及含乳

46、制40400品大米食用植物油100200100199012012.7散裝水殘品66633加油站便利店28315268蜜餞473215清洗產 品50500罐頭食 品交通樞20018614紐便利店散裝和200181192012.6預包裝食用油1531503腌臘肉原料79709菜肴98620966一次性用5預包裝 肉制品15014010糖果10090102012.42012.32012.22012.12011.192011.18炒貨食 品 乳制品 大米 油炸食 品 鮮榨果 蔬汁 即食果 蔬 佐餐飲 料 速凍米 面食品 水餃和 湯圓 預包裝 乳制品 及含乳 制品 米面制 品

47、培炒咖 啡及其 制品 龜苓膏 食用油 工業(yè)食 堂原料 商業(yè)中 心餐飲 單位原 料飲料單 位酒類 油條 集中式 消毒餐 具 預包裝 乳制品 及含乳 制品 果凍96942299299010710618281110099130030001009911131130121121015014916059150500257251628127652982944200197324186549510393939073649巧克力及其制57543品方便面68671食用植艾2092063物油即食豆_r制品31301腌品肉80755十貨堅79736果調味品15513124預包裝熟肉制81801品2011.172011.

48、16食品1559151148乳制品40400豆制品80800燒鹵熟18217111肉一次性包裝材604911料預包裝2011.15乳制品80及含乳80制品大米78762茶葉80800腐竹粉77前m絲776710餅干806812十貨食68662品蜂產品806218膨化食如,80755品小麥粉48462預包裝糕點面80746包糖果40400葡萄酒20200熟制鴨脖.生食水102594314 檸60591產品壽司60600西餐菜60600肴冰淇淋44431餐前小80719吃小麥粉47470大米50500即食現配奶茶50500沙拉44440涼拌菜63630茶葉80791涼茶53530預包裝13乳制品1

49、551550及含乳制品調味品60591預包裝即食豆987810制品蛋制品49472食糖60600涼茶60573谷物及淀粉制60600品罐頭食51501品飲料80791瓶裝飲48453用水固體飲50491料冷凍飲16151品即使水503812產品散裝糕4030102011.2011.2011.122011.112011(10)2011.92011.82011.7點面包 保鮮膜 月餅 預包裝 乳制品 及含乳 制品 兒童食 品 肉制品 十貨 烘烤食 品 餐飲單 位廚房 一次性 用品 餐飲單 位使用 十貨 餐飲單 位制作 油炸食 品 餐飲經 營單位 使用乳 制品 餐飲經 營單位 使用調 味品 市民日

50、常消費 的各種 原輔料 保鮮膜 鮮榨果 蔬汁 使用油 預包裝 乳制品 及含乳 制品31310280278239390100982797904847150464494545047350428353505049111951150453030010010002011983404002011.62011.52011.4蜂蜜 餐飲服 務單位 的即食 熟肉制 品 餐飲經 營單位 的食品 原料腌 臘肉 餐飲經 營單位 的食品 原料使 用油 餐飲經 營單位 使用糧 食和調 味品 餐飲經 營單位 使用水 產十制 品 預包裝 茶葉 預包裝 飲料酒 預包裝 糕點面 包餅干 預包裝 乳制品 及含乳 制品 十貨一次性

51、餐具 餐飲服 務單位 使用紙 品52484115833220164201914040010098291910969249585107877197841383785200175252011.32011.22011.12011.元旦前2010 豆制品2010 小麥粉2010 商業(yè)中心餐飲經 營單位 制作糕 點八、 餐飲經 營單位 制作面 類食品 加油站 便利店 雞精燕 窩飲品 米面制 品 十果堅 果 休閑食 品嬰幼兒 用奶瓶 奶嘴 預包裝 冷凍食 品生產環(huán) 節(jié)抽查 流通環(huán) 節(jié)抽查 餐飲消 費環(huán)節(jié) 抽查 生產環(huán) 節(jié)抽查 流通環(huán) 節(jié)抽查 餐飲消 費環(huán)節(jié) 抽查豆制品小麥粉 商業(yè)中 心餐飲29817812

52、01981663239535936343133823691350473649522274646020219844194109126212085421802085954694541551549916615588278573125050026825810餐飲食 品原料 2010 雙 節(jié)前2011 雙節(jié)全2010.32010.22010.32012.2單位原 料 應節(jié)食 品 生產環(huán) 節(jié)抽查 流通環(huán) 節(jié)抽查 餐飲服 務環(huán)節(jié) 應節(jié)食 品 生產環(huán) 節(jié)抽查 流通環(huán) 節(jié)抽查 餐飲服 務環(huán)節(jié) 大米 生產企 業(yè) 市場暢 銷 餐飲單 位 中秋佳 節(jié)將全 月餅 月餅生 產企業(yè) 餐飲單 位制作 市場暢 銷 醬油和 耗油

53、 乳飲料 和植物 蛋白飲 料 啤酒 預包裝 沙拉醬 預包裝159315326119519148458202555352132143613756119519148458202539636432333332178771747401811810377376112812711981980515102220217170111101414040382固體飲 料 食飲具 瓶裝飲 用水 洗滌消 毒劑 保鮮膜 西餐菜 肴 涼菜 罐頭食 品 現凍動 物性水 產品 預包裝 粗糧谷 物 預包裝 飲料 鮮榨果 蔬汁飲 料 水食水 產品 預包裝 水產品 預包裝 肉制品 預包裝 冷凍飲 品 預包裝 淀粉及 谷物制 品 洗滌

54、劑 產品 預包裝 涼茶 保鮮膜 一次性 飯盒固體飲 料 食飲具 瓶裝飲 用水 洗滌消 毒劑 保鮮膜 西餐菜 肴 涼菜 罐頭食 品 現凍動 物性水 產品 預包裝 粗糧谷 物 預包裝 飲料 鮮榨果 蔬汁飲 料 水食水 產品 預包裝 水產品 預包裝 肉制品 預包裝 冷凍飲 品 預包裝 淀粉及 谷物制 品 洗滌劑 產品 預包裝 涼茶 保鮮膜 一次性 飯盒1331312312921816232302251246518218116059116116015050080800200198260582342861391345505001651650474704848029290917615一次性筷子16315211餐飲使 用油2482426餐飲使用茶葉32275生產環(huán) 節(jié)粽子18162流通環(huán) 節(jié)粽子44440流通環(huán) 節(jié)咸蛋 餐飲服40400務環(huán)節(jié)49481粽子 生產環(huán) 節(jié)兒童1121075食品 流通環(huán) 節(jié)兒童食品2192127生產環(huán)節(jié)食品1483141568流通環(huán) 節(jié)食品3022948餐飲消費環(huán)節(jié)79790生產環(huán)節(jié)應品流通環(huán)節(jié)應節(jié)食品1391390餐飲消 費環(huán)節(jié) 食端午前2010 六一將全2010 春節(jié)