《食品质量安全分析》由会员分享，可在线阅读，更多相关《食品质量安全分析（26页珍藏版）》请在金锄头文库上搜索。
1、 食品质量安全抽检数据分析食品质量安全抽检数据分析 摘 要 食品安全关系千家万户，根据深圳市给出 2010 年、2011 年、2012 年的食品质量安全抽检数据， 针对深圳市的食品质量安全做出了分析。 通过对数据的研究，我们对数据进行了处理，将主要食品分为六大类，分别是：蔬果类、肉蛋类、水产类、淀粉类、调味类、豆乳类。在建立数学模型时，应用了数据筛选和均值统计法、层次分析法、回归分析法、数据拟合法、概率统计分析 对于问题 1，我们将影响食品安全的因素分为微生物、重金属、添加剂、其他这四大类。先用数据筛选和统计法对数据进行处理，从而得到 2010 年、2011年和 2012 年每个季度食品因微生
2、物、重金属、添加剂、其他而导致的不合格率，并用折线图表示出每一类三年的变化趋势。 然后用层次分析法，定量地求出三年食品安全系数。可得知三年内食品安全质量有波动，总体上说，安全情况有下降趋势，食品安全有待提高。 对于问题 2，由于微生物、重金属、添加剂、其他这四方面的不合格导致了食品质量的不过关， 但是这四类直接因素又由其它因素影响， 通过对数据的分析，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。通过建立回归模型，运用 MATLAB 软件，求出六大要素与食品质量间的线性关系， 定量求出各因素在影响食品安全中所占影响度大小和影响规律。最终得
3、到，在食品质量影响因素中食品产地影响最大，食品加工环节影响次之，季节影响最小，抽查地点几乎无影响。 对于问题 3，为了科学有效地反映食品安全情况，合理简化抽检方法，我们减少稳定数据的测量量，增加不稳定数据的测量量，使测量环节有效减少，重要环节的监测力度加大。通过使用 SPSS 软件，求解出新选出的各项因素与食品安全的线性关系，并通过回归曲线进行检验，论证了新的检验方法求得出的数据与原检验方法下的数据误差趋于 0，新方法具有科学性。最终得到，新的抽检方法为减少调料类、豆乳类、蔬果类、淀粉类食品抽检次数，增加水产类、肉蛋类抽检次数。 针对论文的实际情况， 对论文的优缺点做了评价，文章最后还给出了其
4、他的改进方向，以用于指导实际应用。 关键词关键词：统计；层次分析；回归分析；拟合 语句不通顺数据筛选和均值 统计能和层次分 析、回归分析会 馆分析等并列 为“法”吗这个法大了点，起码 包括回归分析法这一句话说得好有分析吗 数据合适才可以用此方法在问题三中，应有体现1问题的重述 “民以食为天” ，食品安全关系到千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高， 食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安
5、全。 另一方面，食品质量与安全又是一个专业性很强的问题， 其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。请下载 2010年、2011 年和 2012 年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取） ，并根据这些资料来讨论： 1.如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势； 2.从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等； 3.能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（
6、食品抽检是需要费用的） ，例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？ 2问题的分析 对于问题 1，由于调查中食物种类复杂，为了方便问题的讨论，我们将食品分为六大类，即淀粉类，水产类，蔬果类，肉蛋类豆乳类和调味品，对数据进行筛选。在考虑食品安全情况的影响因素时，根据数据分析，我们主要分析四大主要因素，分别为微生物、重金属、添加剂和其他。从附录中得到深圳市2010，,2011,2012 三年抽检样本数据，并用 EXCEL 制表进行统计和整理，求出每个因素导致的每个季度的食品不合格率。 根据相应的数据制作出每个因素安全情况三年的变化趋势。 再用层次分析法， 从定量的角度求出三年的总
7、体变化趋势。 对于问题 2，第一问中涉及的影响因素是食品是否安全的直接因素，但这些因素又是由其它相关因素决定。 分析数据，我们将影响食品安全的这些因素归结为检测地点、食品产地、抽查季节、生产环节、流通环节、餐饮环节六个要素。这一问要求的是这些要素与食品质量的规律性关系， 即线性关系。通过建立回归模型，来定量求出各因素在影响食品安全中所占影响度大小和影响规律。对于各因素对食品安全的影响程度，利用偏回归平方和衡量每个变量在回归中所起的作用大小（即影响程度） ，对每个变量的偏回归平方和进行了计算，最后把影响程度由深到浅的各因素进行了排序。 最后进行回归方程的显著性检验，将总的平方和分解为回归平方和和
8、剩余平方和，运用 F 检验法来判定线性回归方程的显著性。 对于问题 3，要改进原有抽检方法，找到一个新方法，使其在成本相同的情况下，测得数据结果更准确， ，在测得数据结果和老方法接近下，成本更低。根据第一问和第二问得到的数据， 得到不同因素的影响程度是不同的。减少稳定数据的测量量，加大不稳定数据的测量量，减少测量环节，加重必要环节。由于我们获取的数据均是在原抽检方法下得到的数据， 而没有真正标准真实的数据，所以无论我们用什么样的新方法， 只能在原抽检下的数据 （或有关预测、 推测数据）进行验证。故新方法在测量不稳定数据的结果与原方法不具可比性。但稳定数据的结果两种方法下的测量数据基本一致。新方
9、法得到后，检验新方法的科学性、合理性。拟合曲线确定其未来检测值，采用蒙特卡洛法抽检，带入回归方程进行验证新方法是否合理。 3模型的假设与符号说明 31 模型的假设 （1）假设模型求解过程中所用的数据都是合理的 （2）假设不考虑抽检时抽检人员的人为影响 （3）假设每次抽检时遵循公平性原则 （4）假设所有食品取主要食品即六大类，其余食品对质量安全分析无影响 （5）假设抽检时所有环节、所有因素的成本和工时相同 32 符号的说明 符号 表示的意义 备注 Ai 目标层指标 i=1 Bi 准则层指标 i=1,2,3,4 Ci 方案层指标 i=1,2,3 第 i 层权重 i=1,2,3 标准化后的第 i 层
10、权重 i=1,2,3 入 一致性指标过度数 C.I 一致性指标 R.I. 平均随机一次性指标 C.R. 一致性比例指标 bij 每个因子的评分 i=1,2,3 j=1,2,3,4 Ci 综合评价指数 i=1,2,3 PI 食品安全质量综合评价指数 影响食品质量安全的因素。 i=1,2,.6 各食品安全系数 a=1,2.6 相互独立且服从同一正态分布的随机变量 变量的偏回归平方和 i=1,2,.6 回归平方和 剩余平方和 待估计系数参数 的回归值 参数的最小二乘估计 4模型的建立与求解 41 问题 1 的模型建立与求解 统计、计算、绘制图表来表示出六大类食品，在 2010、2011、2012 三
11、年中各季度下抽检样品中不合格食品里，微生物、重金属、添加剂的变化趋势；通过数据筛选和统计得出附表一，根据附表一，得出下表并画出相应的折线图： 表表 1 1 四项四项指标的不合格率指标的不合格率 年 份 2010 2011 2012 季 度 一 二 三 四 一 二 三 四 一 二 三 四 微 生 物 0.0 75 0.4 366 0.0 37 0.5 74 0.1 4 0.0 2 0.8 3 0.3 54 0.2 29 0.4 615 0.2 62 0.7 53 重 金 属 0.0 15 0 0 0.0 2 0.0 41 0.1 49 0 0.2 26 0.1 73 0.0 45 0 0.0 0
12、5 添 加 剂 0 0.0 36 0.1 29 0.9 44 0.0 336 0.0 38 0.0 17 0.4 73 0.6 28 0.6 285 0.3 44 0.7 57 其 它 0.0 44 0.5 476 0 0.1 62 0.5 5 0.0 15 0 0.3 66 0.3 89 0.2 0.2 93 0.9 68 图图 1 1 四项指标不合格率的折线图四项指标不合格率的折线图 从图表得出结论：三年来，因微生物 、重金属、添加剂导致的不合格率变化趋势有波动，食品安全整体有下降，从总体上说，食品安全情况没有改善。 根据三年中六类食品的食品安全情况， 采用层次分析法，定量结合定性的评价出
13、由于微生物、重金属、添加剂超标而导致食品不安全的年度食品安全情况。 对深圳市这三年各主要食品领域微生物、 重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对 3 年的数据处理按季节分为 12 个子样本点的抽检值进行分析， 得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。 （1 1） 建立层次结构模型建立层次结构模型 构造各季节食品安全质量综合评价的阶梯层次结构 根据层次分析法的基本步骤,建立食品质量安全评价递阶层次模型,如图 2所示。 图图 2 2 食品综合评食品综合评价递阶层次结构价递阶层次结构 （2 2）构造判断矩阵构造判断矩阵 通过相互比较
14、确定各准则对于目标的权重， 即构造判断矩阵。在层次分析法中，为使矩阵中的各要素的重要性能够进行定量显示，引进了矩阵判断标度（19 标度法） : 表表 2 2 矩阵标度矩阵标度 （3 3）求判断矩阵求判断矩阵 对于要比较的因子而言，你认为一样重要就是 1:1，强烈重要就是 9:1，也可以取中间数值 6:1 等，两两比较，把数值填入，并排列成判断矩阵（判断矩阵是对角线积是 1 的正反矩阵即可） 根据判断矩阵标度内容及其含义，可得 A-B 层判断矩阵 B 利用判断矩阵计算各因素 B 对目标层 A 的权重（权系数） 将 A 的每一列向量归一化得： 对 按行求和得： 将归一化 ，即为近似特 征根（权向量
15、） 计算作为最大特征根的近似值。 （4 4）判断矩阵的一致性检验判断矩阵的一致性检验 判断矩阵通常是不一致的， 但是为了能用它的对应于特征根的特征向量作为被比较因素的权向量，其不一致程度应在容许的范围内.如何确定这个范围？ 一直性指标一直性指标 1nnCICI=0 时A一致； CI 越大，A的不一致性程度 越严重。随机一致性指标 RI： n 1 2 3 4 5 6 7 8 9 10 11RI 0 0 0.58 0.90 1.12 1.24 1.32 1.41 1.45 1.49 1.51一致性比率一致性比率（用于确定用于确定 A A 的不一致性的容许范围的不一致性的容许范围） RICICR当CR0.80 劣 严重劣食品 这样就得到了 2013 年第一至四期食品安全抽样检验产品抽查样本的食品安全的优劣状况，并做出定量的综合评价。其它年份的食品安全情况同理可求。 经计算最终得到了三年的食品安全评价等级表 7： 表表 7 7 三年的食品安全评价等级三年的食品安全评价等级 年份 2010 2011 2012 食品安全系数 0.1932 0.5676 0.2238 等级 4.2 问题二的分析与模型的建立问题二的分析与模型的建立
