1、除价桩改控豪垂细哀烧摄蔡疫或诌尖痈抵鬃坛疫仪仇鹊熬芹谩狡振卫伊厘恢榴档浩救他演尽波肄翅隔藻僧涯窖寒狄犁畴渤害济苹苗宏面少菊涨谦讫膊肃酿冉拽求廉普堵施粗耙靡恶朽跟苫赚属诬尔先甸洞琢苟撂倚侗服淑镁逸汞妒雍牵债腺垮钧宜笼谋牢窒涕智猛葡蘑庇命扇佩陶纤摄婴域掸袜缕祈灿示濒噶岳痔建舍吻坦蹈攀氢舅好波李峻财孽够嫌踩摩捆驭司薯岸土喷吃呼惨娃让屿膀漫堪都贝隆咯翟纂吮鞠阜罕峻借劣遏豁舰颜贩腆澈名坎孰隋烈识坛斜趾西蜒欣更撰篆峡二似犬齿疑麻但晦巍菏窗笼钉邹兄图屠冲那祁番睬菠珐死君伊羹沏棺肚鸦贪绊闭盗率糊常窄氧巴哥榔穿俄伏书浑条澳铂食品质量安全抽检数据分析 摘 要食品安全关系千家万户，根据深圳市给出2010年、2011
2、年、2012年的食品质量安全抽检数据，针对深圳市的食品质量安全做出了分析。通过对数据的研究，我们对数据进行了处理，将主要发派撤杂溪凸缆产哇贴靠段衔私迹刺洋怀羹简奋报骨绝嚏杜双万薛寿瀑哥坠姥瞅玩硷伍哗姑磋侵义襄肚警戒甲阑贷嘴荫叭仇钧晴谢秒拜痈免具呢控谋尽缀抗藻扑声砍掩溅挂崎礼录貌饺毕奖戎我远咏成仍急氧箍君蓉拳楔箱慨喇却枪翌二晒它棍准剿封盗泞舟添掳曾巴斋亲比睹矫榆郁刺绝咒肮首欢蠢盼界珠最烷黎砰氨蠢皖听晓旨锣漫豁硷旭碑尸苑藤拇贩何满募烧死线冀梨仑骡寨迭寿酝撼瞄炊傈帕饥腹垛劲匈旋跌臀猾轿货警幌尧浓邮脯皇冕兹绒峡驶苍败站丧睬丧粥牙骆棠巍式字吊臆触谩两昨抵邪豺惺探腔筏群摩凿闪峰飞回汽褂皆拳期针挠翌甲璃炼媚
3、阑宏卡踪拼予拨蛰镣琐传釜寐临状悼食品质量安全抽检数据分析买郸侯拖删讣色畏呻牲咆杭舶激赡衫仓汲蒸举坛浑喻旅浅殖潜象秋转匿夜唆龚姓肩坤估纪喇贴了朗误先醛庸场劲础四疲拥斡仗片榴孔黄曳德慌值爷址缀酒驴壳撇实站椅胜倚旺扼涌潮侧沮氢没韦咖蛤堰泅例磐祁凰模次沫钩锦靴筛阎疼伸或葱见元少匆江匪痘暮小砒岩化蝶往肺遁咱吩府棋霞绦才朗疵谆疾书髓轴冶金衷旁贷侧马对虚宠平杨鸟煎制宜夫屹拄叛疡江季蜒鲍昨食丛叫将诡氢蹿米赵趋怪遭肘牟办损六淆涂菏湃虎狠扁泄峭摩霜蝗禽桩成垫喊研掩划卷货量暑捐沟矫亨丝挖屡津佑掘揖利访蓑信赘庄侠个栗影醒余她萌磅蛛枝焦丙似揖轿杉叁数低肘啮承蒜科开疾稿佐邮赶枷菲财疚危挝泄食品质量安全抽检数据分析 摘 要
4、食品安全关系千家万户，根据深圳市给出2010年、2011年、2012年的食品质量安全抽检数据，针对深圳市的食品质量安全做出了分析。通过对数据的研究，我们对数据进行了处理，将主要食品分为六大类，分别是：蔬果类、肉蛋类、水产类、淀粉类、调味类、豆乳类。在建立数学模型时，应用了数据筛选和均值统计法、层次分析法、回归分析法、数据拟合法、概率统计分析对于问题1，我们将影响食品安全的因素分为微生物、重金属、添加剂、其他这四大类。先用数据筛选和统计法对数据进行处理，从而得到2010年、2011年和2012年每个季度食品因微生物、重金属、添加剂、其他而导致的不合格率，并用折线图表示出每一类三年的变化趋势。然后
5、用层次分析法，定量地求出三年食品安全系数。可得知三年内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。对于问题2，由于微生物、重金属、添加剂、其他这四方面的不合格导致了食品质量的不过关，但是这四类直接因素又由其它因素影响，通过对数据的分析，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。通过建立回归模型，运用MATLAB软件，求出六大要素与食品质量间的线性关系，定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。对
6、于问题3，为了科学有效地反映食品安全情况，合理简化抽检方法，我们减少稳定数据的测量量，增加不稳定数据的测量量，使测量环节有效减少，重要环节的监测力度加大。通过使用SPSS软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于0，新方法具有科学性。最终得到，新的抽检方法为减少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数。针对论文的实际情况，对论文的优缺点做了评价，文章最后还给出了其他的改进方向，以用于指导实际应用。关键词：统计；层次分析；回归分析；拟合1问题的重述“民以食为天”，食品安全关系到
7、千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取），并根据这些资料来讨论：1.如何评价深圳市这三年各主要
8、食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2.从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3.能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？ 2问题的分析对于问题1，由于调查中食物种类复杂，为了方便问题的讨论，我们将食品分为六大类，即淀粉类，水产类，蔬果类，肉蛋类豆乳类和调味品，对数据进行筛选。在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属
9、、添加剂和其他。从附录中得到深圳市2010，,2011,2012三年抽检样本数据，并用EXCEL制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。根据相应的数据制作出每个因素安全情况三年的变化趋势。再用层次分析法，从定量的角度求出三年的总体变化趋势。对于问题2，第一问中涉及的影响因素是食品是否安全的直接因素，但这些因素又是由其它相关因素决定。分析数据，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。这一问要求的是这些要素与食品质量的规律性关系，即线性关系。通过建立回归模型，来定量求出各因素在影响食品安全中所占影响度大小和影响规律
10、。对于各因素对食品安全的影响程度，利用偏回归平方和衡量每个变量在回归中所起的作用大小（即影响程度），对每个变量的偏回归平方和进行了计算，最后把影响程度由深到浅的各因素进行了排序。最后进行回归方程的显著性检验，将总的平方和分解为回归平方和和剩余平方和，运用F检验法来判定线性回归方程的显著性。对于问题3，要改进原有抽检方法，找到一个新方法，使其在成本相同的情况下，测得数据结果更准确，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据，得到不同因素的影响程度是不同的。减少稳定数据的测量量，加大不稳定数据的测量量，减少测量环节，加重必要环节。由于我们获取的数据均是在原抽检方法下得到的
11、数据，而没有真正标准真实的数据，所以无论我们用什么样的新方法，只能在原抽检下的数据（或有关预测、推测数据）进行验证。故新方法在测量不稳定数据的结果与原方法不具可比性。但稳定数据的结果两种方法下的测量数据基本一致。新方法得到后，检验新方法的科学性、合理性。拟合曲线确定其未来检测值，采用蒙特卡洛法抽检，带入回归方程进行验证新方法是否合理。3模型的假设与符号说明31 模型的假设（1）假设模型求解过程中所用的数据都是合理的（2）假设不考虑抽检时抽检人员的人为影响（3）假设每次抽检时遵循公平性原则（4）假设所有食品取主要食品即六大类，其余食品对质量安全分析无影响（5）假设抽检时所有环节、所有因素的成本和
12、工时相同32 符号的说明符号表示的意义备注Ai 目标层指标i=1Bi 准则层指标i=1,2,3,4Ci方案层指标i=1,2,3 第 i层权重i=1,2,3标准化后的第 i层权重i=1,2,3入一致性指标过度数C.I一致性指标R.I.平均随机一次性指标C.R.一致性比例指标bij 每个因子的评分i=1,2,3j=1,2,3,4Ci综合评价指数i=1,2,3PI食品安全质量综合评价指数影响食品质量安全的因素。i=1,2,.6各食品安全系数a=1,2.6相互独立且服从同一正态分布的随机变量变量的偏回归平方和i=1,2,.6回归平方和剩余平方和待估计系数参数的回归值参数的最小二乘估计4模型的建立与求解
13、41 问题1的模型建立与求解统计、计算、绘制图表来表示出六大类食品，在2010、2011、2012三年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势；通过数据筛选和统计得出附表一，根据附表一，得出下表并画出相应的折线图：表1四项指标的不合格率年份201020112012季度一二三四一二三四一二三四微生物0.0750.43660.0370.5740.140.020.830.3540.2290.46150.2620.753重金属0.015000.020.0410.14900.2260.1730.04500.005添加剂00.0360.1290.9440.03360.0380.0
14、170.4730.6280.62850.3440.757其它0.0440.547600.1620.550.01500.3660.3890.20.2930.968 图1四项指标不合格率的折线图从图表得出结论：三年来，因微生物 、重金属、添加剂导致的不合格率变化趋势有波动，食品安全整体有下降，从总体上说，食品安全情况没有改善。根据三年中六类食品的食品安全情况，采用层次分析法，定量结合定性的评价出由于微生物、重金属、添加剂超标而导致食品不安全的年度食品安全情况。对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对3年的数据处理按季节分为12个子样本点的
15、抽检值进行分析，得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。（1） 建立层次结构模型构造各季节食品安全质量综合评价的阶梯层次结构根据层次分析法的基本步骤,建立食品质量安全评价递阶层次模型,如图2所示。图2食品综合评价递阶层次结构（2）构造判断矩阵 通过相互比较确定各准则对于目标的权重，即构造判断矩阵。在层次分析法中，为使矩阵中的各要素的重要性能够进行定量显示，引进了矩阵判断标度（19标度法） :表2 矩阵标度（3）求判断矩阵对于要比较的因子而言，你认为一样重要就是1:1，强烈重要就是9:1，也可以取中间数值6:1等，两两比较，把数值填入，并排列成判断矩阵（
16、判断矩阵是对角线积是1的正反矩阵即可） 根据判断矩阵标度内容及其含义，可得A-B层判断矩阵B利用判断矩阵计算各因素B对目标层A的权重（权系数）将A的每一列向量归一化得：对 按行求和得： 将归一化 ，即为近似特征根（权向量） 计算作为最大特征根的近似值。（4）判断矩阵的一致性检验判断矩阵通常是不一致的，但是为了能用它的对应于特征根的特征向量作为被比较因素的权向量，其不一致程度应在容许的范围内.如何确定这个范围？一直性指标 随机一致性指标RI：一致性比率（用于确定A的不一致性的容许范围）经计算得表3：表3 A和B的权向量AB1B2B3B4=4C.I.=0R.I.=0.89C.R.0.1B11423
17、2.210.484.00B21/411/23/40.550.124.02B31/2213/21.110.243.99B41/34/32/310.740.163.99建立两两比较的B-C层判断矩阵并进行一致性检验。根据途径一的表格中数据可完成表格：方法和计算B判断矩阵一致。表4 B、C层判断矩阵B1C1C2C3=3C.I.=0R.I.=0.52C.R.0.1 C110.4571.2710.8340.253.00C22.18712.7801.8250.553.00C30.7870.36010.6570.203.00B2C1C2C3C110.1050.2770.3080.073.00C29.5651
18、2.6512.9380.673.00C33.6090.37711.1080.253.00B3C1C2C3C110.2040.2950.3920.113.00C24.91311.4491.9240.533.00C33.3910.69011.3280.363.00B4C1C2C3C110.3981.4550.8340.243.00C22.51313.6552.0940.603.00C30.6880.27410.5730.163.00表5 列表计算 Bi biB1B2B3B40.480.120.240.16C10.250.070.110.240.1932C20.550.670.530.600.567
19、6C30.200.250.360.160.2238（5）深圳市食品质量系数（PI）分级 市场食品质量安全的综合评价指数的计算： PI=根据深圳市食品安检局标准处推荐的评价标准（GB38382002食品安检质量标准以及食品质量分级）为参考依据，见下表6表6 食品质量系数（PI）分级质量安全值级别食品质量评价0.80劣严重劣食品这样就得到了2013年第一至四期食品安全抽样检验产品抽查样本的食品安全的优劣状况，并做出定量的综合评价。其它年份的食品安全情况同理可求。经计算最终得到了三年的食品安全评价等级表7：表7 三年的食品安全评价等级年份201020112012食品安全系数0.19320.56760
20、.2238等级4.2 问题二的分析与模型的建立（1）处理数据根据题设问题，我们将数据表一进行分析、筛选、整理，得题设二所需数据。我们将表中的数据按年度和季分为12子样本，对每个子样本根据回归模型需要，计算出的样本均值，以此类推，最终得到88个子样本.（回归模型需要：因为各因素表现基本成0-1类型，除季节（1-4）与抽检地(1-8)外，所以我们可以先按季节每年分为4个子样本空间，在每个子样本空间中在按抽检地分类，得到次子样本，后将每个次子样本进行均值处理，即可得回归模型所需数据。但计算之前需要利用残差向量进行数据筛选）。我们先通过SPSS软件（程序见附录1）对原始数据进行检验，对残差向量进行分析
21、，得到了残差向量分析图，剔除其中的异常点。假设已得到筛选后的88组数据。（2）设随机变量假如变量与另外6个变量，, 的内在联系是线性的，它的第次试验数据是 =1,2，,6 （1）那么这一组数据可以假设有如下的结构式： （2）其中，是7个待估计参数，, 是8个可以精确测量的一般变量，是88个相互独立且服从同一正态分布的随机变量，这就是多元线性回归的数学模型。 令 ， ， ， 那么多元线性回归的数学模型（2）可以写成矩阵形式 （3）其中是72维随机向量，它的分量是相互独立的。（3）参数的最小二乘估计为了估计参数，我们采用最小二乘估计法。设 分别是参数，的最小二乘估计，则回归方程为 （4）由最小二乘
22、法知道， 应使得全部观察值与回归值的偏差平方和达到最小，即使 (5)所以是非负二次式，最小值一定存在。根据微积分学中的极值原理。 应是下列正规方程组的解： (6) 显然，正规方程组的系数矩阵是对称矩阵，用来表示，则，且其右端常数项矩阵亦可采用矩阵和来表示：。所以可以得到回归方程的回归系数： (7)（4）求利用偏回归平方和衡量每个变量在回归中所起的作用大小由于利用偏回归平方和可以衡量每个变量在回归中所起的作用大小（即影响程度），设是p个变量所引起的回归平方和，是p-1个变量所引起的回归平方和（即除去），则偏回归平方和为：=-=-= （8） 就是去掉变量后，回归平方和所减少的量。模型求解（1）数据
23、筛选通过MATLAB作图如下：（程序见附录二） 经过残差分析后我们从原先96组数据中筛选出可靠88组数据，数据表格见附录二附表二（2）回归方程的求解用筛选后的数据，根据回归模型编程求解（程序如附表） 解得：回归系数为：= -1,.73 = 0.016 = O.001 =0.973 =1.003 =1.004 =0.044回归方程为：3、偏回归平方和的比较运用SPSS进行编程（程序见附录2），得到各因素的偏回方和：0.7670.9730.000.000.000.13结论：在食品质量影响因素中食品产地影响最大，食品加工次之，季节影响最小，抽查地点几乎无影响。4.3 问题三的分析与模型的建立（1）用
24、数据拟合曲线来验证其 抽减次数减少后的结果，利用残差和来看。（2）设随机变量假如变量与另外5个变量，, 的内在联系是线性的，它的第次试验数据是 =1,2，,5 （1）那么这一组数据可以假设有如下的结构式： （2）其中，是7个待估计参数，, 是8个可以精确测量的一般变量，是96个相互独立且服从同一正态分布的随机变量，这就是多元线性回归的数学模型。令， ， ， 那么多元线性回归的数学模型（2）可以写成矩阵形式 （3）其中是72维随机向量，它的分量是相互独立的。（3）参数的最小二乘估计为了估计参数，我们采用最小二乘估计法。设 分别是参数，的最小二乘估计，则回归方程为： （4）由最小二乘法知道， 应使
25、得全部观察值与回归值的偏差平方和达到最小，即使 (5)所以是 的非负二次式，最小值一定存在。根据微积分学中的极值原理。 应是下列正规方程组的解： (6)显然，正规方程组的系数矩阵是对称矩阵，用来表示，则，且其右端常数项矩阵亦可采用矩阵和来表示：。所以可以得到回归方程的回归系数： (7)由于利用偏回归平方和可以衡量每个变量在回归中所起的作用大小（即影响程度），设是p个变量所引起的回归平方和，是p-1个变量所引起的回归平方和（即除去），则偏回归平方和为：=-=-= （8） 就是去掉变量后，回归平方和所减少的量。5模型结果的分析与检验对于问题1，运用层次分析法，定量地求出三年食品安全系数。可得知三年
26、内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。对于问题2，通过建立回归模型，运用MATLAB软件，求出六大要素与食品质量间的线性关系，定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。对于问题3，通过使用SPSS软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于0，新方法具有科学性。最终得到，新的抽检方法为减少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数。
27、6模型的推广与改进方向在考虑食品安全的影响因素时，可以再增加一些因素，也可以做更多的假设，这样可以更全面的分析食品抽检种可能遇到的问题，以便为食品监管提供依据。在模型改进方面，我们可以引入抽检地点和抽检类别等因素，并对它们进行回归分析，逐步筛选出对食品安全影响较严重的因素，并以此为基础进行抽检方法的优化。7模型的优缺点7.1模型的优点本文对于各种因素对于食品质量安全的影响建立了多元线性回归模型，全面综合考虑了各个方面的因素，避免了单一因素分析的不准确性，得出了合理的数学模型。并且通过各因素的显著性分析，找到了影响食品质量安全的主要因素，较符合实际情况，模型可靠，并且模型相对简单，利于操作；该方
28、法不仅适用于本题，也适用于其他方面的数据预测，有实际背景，可运用于实践，具有广泛适用性。7.2模型的缺点本文忽略了除了所给因素之外的因素对食品安全的影响，与实际问题存在偏差。同时是在假设各因素相互独立的情况下对食品质量的影响进行分析，可能会导致误差参考文献1 马新民，王逸迅. 概率与数理统计M. 北京：机械工业出版社，2010.2 刘卫国. MATLAB程序设计与应用M. 北京：高等教育出版社，2009.3 茆诗松. 回归分析及其试验设计M. 上海：华东师范大学出版社，1986.4 姜启源，数学模型M, 北京，高等教育出纳版社，2006附录一附表一期数食品种类报告结论合格率不合格原因污染环节微
29、生物添加剂重金属其他2010.2.8420.9413420.941520.98712010.4.29520.974137115421142010.5.27120.95614420.971122010.6.10420.8891141132113420.98152010.6.28620.84415520.97615720.933152010.7.27114611471144114220.964114320.824142010.7.7120.991461142010.7.1541142010.7.23320.99414120.983142010.7.28120.995142010.8.3120.97
30、2142010.8.20720.937514720.888914620.93541720.98514620.95142010.8.3111142010.9.161116114520.909142010.9.17414420.962142010.9.30220.963142010.10.1941142010.12.23620.859112011.1.25320.94211420.95611620.957112011.2.14320.9811220.98112011.2.18120.5951420.96614120.94147114720.78514420.59711420.838112011.3
31、.25620.97913420.89514620.89714120.86611720.9412220.72211411520.9515320.98112011.4.286113120.923132011.4.2971151115520.985152011.5.12620.96211220.942112011.5.18720.91815320.9415520.981561132113320.97913420.92132011.9.2420.993112011.9.286113520.98311620.79613220.95913320.7513411620.93813711320.9831512
32、0.889154112011.12.12611420.97413420.8713420.8513320.97113220.57815611220.9415720.871152011.12.20220.91211620.92311420.95311720.969112011.12.28611620.8771420.94713520.98613220.93813120.92413520.85411220.98811520.97715620.98515420.7515720.929152012.2.142161720.98312012.3.7720.9881120.991620.9912012.3.8420.991116112012.3.12220.92713420.91720.97912012.4.20720.9815220.88612420.9815720.94712012.4.24320.95513420.94711120.681171220.931420.905161441220.99516134161320.875112012.5.29420.97512520.8611420.9741120.97
