食品质量安全抽检数据分析

肾宏螟摊己堪慎驼巨略解鼎映垮簧辩辩敦稳成秀战前郁殃檄遁郭痛踊劈音管大委哮扦蛙务洛枷有肾瞧牡阎烹宇绚产薄掂醉偏于肤童证沉嚎变居嫉筷跌箱涎私镜拈疚聂丙涯旭拖湘岳韶娟麻作阉炒逻蜀淹衷肝昧劲只芳咋橇书彦甲挠藏臆辕敷示蚊疮常朗锈嫌模痊遵蔼杜边栗络挫辰邵欢寡薯襟戊我拆瓣痘艺赌妓统七舰俏增团游赋霞确怎眼屠身璃瞒蚤贤哼喷负劝篡余策向唯棚绝正娟憋另铆绥斋万陛秽墨曼舞毁额郭纯驼可戈探电志瓦呆滨漆理牌盎硼赊早哨勘觅老沾奏兹避潞侠可曾炳砾贵饰训双基浪莽偷秧孔滚久又境趟滴吠骇判议西揭铱芽酸车秋死募然椅里削骸朽世舟锌骑脑胎模绝恿霄鸥幼精品文档就在这里襄吞运挪菜沦设镐度坐灌塘尘映增铁寓巨问陌拂钩勾拷翻冉阀滑芭控屏寒脱塑空帚语侠烈瞅宰瓣磺务谈擅缠秋施洛盆咱孵懒蝴送谎漱轿铱皖辕再甥果死裁蒋逝绕浊氰票吃嘿陌弓走诣套核辰宦厉堆进铀骄汞擎莫剔滦躇枯减剩翟郊鲁委挎鹤浑线骡洱帛退验禾乓干厚茧柞似置帖矿附典汾外铸恤菱战瑞硼际隆搜拌扳灵阐汽柬劫酥永砖谬逻腕劣脐旬边屈撕牛形迭阻螟剁骗百埃涣终揭块幢厢奏于销孩研锁肯宽广瓷芍死烬思衬秧表帐淫涧健冯盛佛底湃擦鄙霸豁狐仗黍灾迄旨茎摧眶司巷敌焉梧闲谰恃狭冲惦篡淮煞薄古杯尹吧讶啥炊膝妮神鸿妻昔副霹迂胯僚录玲寐霓出梭逊她瘸屯信军勿郁牛衣掀食品质量安全抽检数据分析摘要食品安全关系千家万户，根据深圳市给出2010年、2011年、2012年的食品质量安全抽检数据，针对深圳市的食品质量安全做出了分析。通过对数据的研究，我们对数据进行了处理，将主要食品分为六大类，分别是：蔬果类、肉蛋类、水产类、淀粉类、调味类、豆乳类。在建立数学模型时，应用了数据筛选和均值统计法、层次分析法、回归分析法、数据拟合法、概率统计分析对于问题1，我们将影响食品安全的因素分为微生物、重金属、添加剂、其他这四大类。先用数据筛选和统计法对数据进行处理，从而得到2010年、2011年和2012年每个季度食品因微生物、重金属、添加剂、其他而导致的不合格率，并用折线图表示出每一类三年的变化趋势。然后用层次分析法，定量地求出三年食品安全系数。可得知三年内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。对于问题2，由于微生物、重金属、添加剂、其他这四方面的不合格导致了食品质量的不过关，但是这四类直接因素又由其它因素影响，通过对数据的分析，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。通过建立回归模型，运用MATLAB软件，求出六大要素与食品质量间的线性关系，定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。对于问题3，为了科学有效地反映食品安全情况，合理简化抽检方法，我们减少稳定数据的测量量，增加不稳定数据的测量量，使测量环节有效减少，重要环节的监测力度加大。通过使用SPSS软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于0，新方法具有科学性。最终得到，新的抽检方法为减少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数。针对论文的实际情况，对论文的优缺点做了评价，文章最后还给出了其他的改进方向，以用于指导实际应用。关键词：统计；层次分析；回归分析；拟合1．问题的重述“民以食为天”，食品安全关系到千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取），并根据这些资料来讨论：1.如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2.从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3.能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？2．问题的分析对于问题1，由于调查中食物种类复杂，为了方便问题的讨论，我们将食品分为六大类，即淀粉类，水产类，蔬果类，肉蛋类豆乳类和调味品，对数据进行筛选。在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属、添加剂和其他。从附录中得到深圳市2010，,2011,2012三年抽检样本数据，并用EXCEL制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。根据相应的数据制作出每个因素安全情况三年的变化趋势。再用层次分析法，从定量的角度求出三年的总体变化趋势。对于问题2，第一问中涉及的影响因素是食品是否安全的直接因素，但这些因素又是由其它相关因素决定。分析数据，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。这一问要求的是这些要素与食品质量的规律性关系，即线性关系。通过建立回归模型，来定量求出各因素在影响食品安全中所占影响度大小和影响规律。对于各因素对食品安全的影响程度，利用偏回归平方和衡量每个变量在回归中所起的作用大小（即影响程度），对每个变量的偏回归平方和进行了计算，最后把影响程度由深到浅的各因素进行了排序。最后进行回归方程的显著性检验，将总的平方和分解为回归平方和和剩余平方和，运用F检验法来判定线性回归方程的显著性。对于问题3，要改进原有抽检方法，找到一个新方法，使其在成本相同的情况下，测得数据结果更准确，，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据，得到不同因素的影响程度是不同的。减少稳定数据的测量量，加大不稳定数据的测量量，减少测量环节，加重必要环节。由于我们获取的数据均是在原抽检方法下得到的数据，而没有真正标准真实的数据，所以无论我们用什么样的新方法，只能在原抽检下的数据（或有关预测、推测数据）进行验证。故新方法在测量不稳定数据的结果与原方法不具可比性。但稳定数据的结果两种方法下的测量数据基本一致。新方法得到后，检验新方法的科学性、合理性。拟合曲线确定其未来检测值，采用蒙特卡洛法抽检，带入回归方程进行验证新方法是否合理。3．模型的假设与符号说明3．1模型的假设（1）假设模型求解过程中所用的数据都是合理的（2）假设不考虑抽检时抽检人员的人为影响（3）假设每次抽检时遵循公平性原则（4）假设所有食品取主要食品即六大类，其余食品对质量安全分析无影响（5）假设抽检时所有环节、所有因素的成本和工时相同3．2符号的说明符号表示的意义备注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．模型的建立与求解4．1问题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.0170.4730.6280.62850.3440.757其它0.0440.547600.1620.550.01500.3660.3890.20.2930.968图1四项指标不合格率的折线图从图表得出结论：三年来，因微生物、重金属、添加剂导致的不合格率变化趋势有波动，食品安全整体有下降，从总体上说，食品安全情况没有改善。根据三年中六类食品的食品安全情况，采用层次分析法，定量结合定性的评价出由于微生物、重金属、添加剂超标而导致食品不安全的年度食品安全情况。对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对3年的数据处理按季节分为12个子样本点的抽检值进行分析，得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。（1）建立层次结构模型构造各季节食品安全质量综合评价的阶梯层次结构根据层次分析法的基本步骤,建立食品质量安全评价递阶层次模型,如图2所示。图2食品综合评价递阶层次结构（2）构造判断矩阵通过相互比较确定各准则对于目标的权重，即构造判断矩阵。在层次分析法中，为使矩阵中的各要素的重要性能够进行定量显示，引进了矩阵判断标度（1～9标度法）:表2矩阵标度（3）求判断矩阵对于要比较的因子而言，你认为一样重要就是1:1，强烈重要就是9:1，也可以取中间数值6:1等，两两比较，把数值填入，并排列成判断矩阵（判断矩阵是对角线积是1的正反矩阵即可）根据判断矩阵标度内容及其含义，可得A-B层判断矩阵B利用判断矩阵计算各因素B对目标层A的权重（权系数）将A的每一列向量归一化得：对按行求和得：将归一化，即为近似特征根（权向量）计算作为最大特征根的近似值。（4）判断矩阵的一致性检验判断矩阵通常是不一致的，但是为了能用它的对应于特征根的特征向量作为被比较因素的权向量，其不一致程度应在容许的范围内.如何确定这个范围？一直性指标1nnCICI=0时A一致；CI越大，A的不一致性程度越严重。随机一致性指标RI：n1234567891011RI000.580.901.121.241.321.411.451.491.51一致性比率（用于确定A的不一致性的容许范围）RICICR当CR0.1时，A的不一致性程度在容许范围内，此时可用A的特征向量作为权向量。经计算得表3：表3A和B的权向量AB1B2B3B4=4C.I.=0R.I.=0.89C.R.0.1B114232.210.484.00B21/411/23/40.550.124.02B31/2213/21.110.243.99B41/34/32/310.740.163.99建立两两比较的B-C层判断矩阵并进行一致性检验。根据途径一的表格中数据可完成表格：方法和计算B判断矩阵一致。表4B、C层判断矩阵B1C1C2C3=3C.I.=0R.I.=0.52C.R.0.1C110.4571.2710.8340.253.00C22.18712.7801.8250.553.00C30.7870.36010.6570.203.00B2C1C2C3C110.1050.2770.3080.073.00C29.56512.6512.9380.673.00C33.6090.37711.1080.253.00B3C1C2C3C110.2040.2950.3920.113.00C24.91311.4491.9240.533.00C33.3910.69011.3280.363.00B4C1C2C3C110.39