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

编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第44页 共44页 食品质量安全抽检数据分析摘要本文根据文章提出的不同问题，建立相应的数学模型，利用matlab软件进行求解，对食品质量进行评价和找规律以及合理抽检方法。针对问题一，对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对数据处理按季节分为17个子样本点的抽样值进行分析，得到各子样本食品安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。可以得到深圳市近三年来的食品质量情况应是明显提高。年份201120122013食品安全系数0.270.170.095等级级级级针对问题二，我们先通过MATLAB对原始数据进行检验，对残差向量进行分析，得到了残差向量分析图，剔除其中的异常点。运用MATLAB进行编程，得到各因素的偏回方和：（）0.33840.00300.36850.10090.14291.24180.14720.1963根据的大小可判断各因素对食品安全系数的影响程度结论：在食品质量影响因素中食品产地影响最大，食品加工次之，季节影响最小，抽查地点几乎无影响。 针对问题三，根据问题一、二数据结果以定量比较评估的方法分析了各类影响食品安全的因素及其可能造成的危害性的问题。改进后的食品抽检的办法以主要食品为准则层建立了层次分析（AHP）模型，对影响食品安检的危害性因素做出定量分析如问题二个影响因素大小。由问题一、二结果可以通过建立抽检模型，即改进后的规准型抽样检验模型，并以蒙特卡罗法对抽检的全过程进行模拟，得到相对误差逐渐趋向于0。关键词： 方差与回归分析 残差向量分析 评价指数 层次分析 一 问题重述 “民以食为天”，食品安全关系到千家万户的生活与健康。随着人们对生活质量的追求和安全意思的提高，食品安全已成为社会关注的热点，也是政府民生工程的一个主题。城市食品的来源越来越广泛，人们消费加工好的食品的比例也越来越高，因此除食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。深圳是食品抽检、监督最统一、最规范、最公开的城市之一。请下载2010年、2011年和2012年深圳市的食品抽检数据（注意蔬菜、鱼类、鸡鸭等抽检数据的获取），并根据这些资料来讨论：1. 如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势；2. 从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等；3. 能否改进食品抽检的办法，使之更科学更有效地反映食品质量状况且不过分增加监管成本（食品抽检是需要费用的），例如对于抽检结果稳定且抽检频次过高的食品领域该作怎样的调整？ 二、问题分析本问题是在监督最统一，最规范，最公开的城市之一深圳，通过对食品安全抽查数据，抽样，合格，不合格以及与近几年数据比较的数据比例进行研究，探讨使人们追求的食品质量安全与可持续科学发展观加快吻合的方案。 对于市场食品安全的分析，我们知道，和有关的变量有8个，研究与变量，,之间的定量关系的问题为多元回归问题。又因为许多多元非线性回归问题都可以化为多元线性回归问题，所以对于本问题我们建立了多元线性回归的数学模型。第一问：由于食品安全生成这一过程涉及到诸多环节，所以我们必须考虑食材的生产收获外，食品的运输、加工、包装、贮存、销售以及餐饮等每一个环节皆可能影响食品的质量与安全。另一方面，食品质量与安全又是一个专业性很强的问题，其标准的制定和抽样检测及评价都需要科学有效的方法。运用统计分析中的方差，分析各个影响因子对食品安全的影响程度。同时我们确定出哪些因素是主要的，哪些是次要的。通过表格可以看出，由于抽样的各种不同方法以及子样空间的不同导致的子样均差以及子样方差的不同，导致了食品安全系数的标准不同。首先，考虑抽样方法的影响。通过对抽样数据和合格数据的观察，我们发现在不同抽样方法以及不同子样本空间下，对食品安全系数有较大的影响。通过对整个表格的分析我们把抽样数据分了四个数据梯度,，分别为2013年的，2012年的，以及2011年和2010年的，然后列出四组数据，再应用方差分析模型进行分析，判断出评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势。其次，我们在来考虑诸多环节对食品安全系数的影响。我们把数据分为四个组进行分析，分别为取材环节，生产环节，流通环节，餐饮环节。然后我们同样应用方差分析模型来对数据进行处理，来判断各个环节对评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势。第二问：首先对附表1的数据进行处理，对残差向量进行，剔除其中的异常点。然后我们建立了多元线性回归的数学模型，并采用了最小二乘法来估计参数。把模型写成矩阵的形式，化简整理得其正规方程组，通过对正规方程组的求解，最后得到回归方程。对于各因素对食品安全的影响程度，由于利用偏回归平方和可以衡量每个变量在回归中所起的作用大小（即影响程度），我们对每个变量的偏回归平方和进行了计算，最后把影响程度由深到浅的各因素进行了排序。回归方程的显著性检验：事先我们并不能断定随机变量与一般变量，,之间是否确有线性关系。在求线性回归方程前线性回归模型只是一种假设，所以在求出线性回归方程之后，我们需要对其进行统计检验。将总的平方和分解为回归平方和和剩余平方和，运用F检验法来判定线性回归方程的显著性。回归系数的显著性检验：由于回归方程显著并不意味着每个自变量，,对因变量的影响都是重要的。而我们要找出影响食品安全的主要因素，即从回归方程中剔除那些次要的、可有可无的变量，这就需要我们对每个变量进行考察。显然，如果某个变量对的作用不显著，那么在多元线性回归模型中，它前面的系数就可以取值为零。因此，检验因子是否显著等价于检验假设 。最后再运用T检验法来辨别模型中哪些因子是显著的及从这些数据中能否找出某些规律性的东西：如食品产地与食品质量的关系；食品销售地点（即抽检地点）与食品质量的关系；季节因素等等。第三问：由于回归系数之间存在相关性，当从原回归方程中剔除一个变量时，其他变量，特别是与它密切相关的一些变量的回归系数就会受到影响，剔除一个变量后，这个变量对的影响很大部分转加到另一个变量对的影响上。所以，我们对回归系数进行一次检验后，只能剔除所有不显著因子中值最小的，然后重新建立新的回归方程，再对新的回归系数逐个进行检验，直到余下的回归系数都显著为止。 三、问题假设1、土壤对取材生长状况基本相同，从而对生成食品的安全系数的影响相同。2、假设周围自然环境因素除去土壤外，其他的对食品安全系数的生成没有影响。3、在食品生产过程中，材料都已经转化成了食品。也就是说产生是完全的，没有发生材料的剩余。4、同类食品是相同的，也就是说抽查的食品对生成食品的安全系数影响相同。5、餐饮服务等条件始终是不变的。6、抽查是在随机下进行，不存在任何的人为干扰。四、 参数及符号说明表一符号表示的意义单位备注第i个子样点与样本空间的差值个第i个空间的大小个第i个子样本第j种食品的合格%第i个样本空间第j子样本经管理后下一年的合格率%第i个样本空间中第j子样本总量个第i个子样本第j个食品合格率%管理系数1/月样本空间与子样本之间差距个PI食品安全质量综合评价指数1Ai准则层项目因子Bi项目层项目因子C.I.一致性指标R.I.平均随机一次性指标C.R.一致性比例指标为第i年需要处理的不合格食品个i=110为预测的i年的食品不合格数个i=110 表二符号说明影响食品质量安全的因素。（）各食品安全系数（）相互独立且服从同一正态分布的随机变量（）变量的偏回归平方和总平方和回归平方和剩余平方和待估计系数参数（）的回归值参数的最小二乘估计（）食材生长土壤，食材的生产收获，食品的运输、加工、包装、贮存、销售以及餐饮八个影响因素 五、模型的分析、建立与求解 问题一的模型问题一为：如何评价深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势 1. 问题的分析与模型的建立对深圳市这三年各主要食品领域微生物、重金属、添加剂含量等安全情况的变化趋势做出定量的综合评价，首先对数据处理按季节分为17个子样本点的抽检值进行分析，得到各子样本食品质量安全情况。运用层次分析法决定食品安全单位指标在综合评价中权重。1）构造各季节食品安全质量综合评价的阶梯层次结构在构造深圳市食品领域综合评价模型时，考虑地理位置，人员结构主要考虑老中青人员比例，经济结构等4项食品安全指标的综合效果。构造其综合评价的递阶层次结构如下图：食品安全综合评价2011年数据2010年数据2012年数据2013年数据局重金属含量微生物含量添加剂含量其他含量 图：食品综合评价递阶层次结构2）建立V-A层判断矩阵建立V-A层判断矩阵的依据是：考虑A层的四个因素对食品质量安全综合评价的重要性时，认为年均指标A1最重要，且A1比微生物、重金属、添加剂含量等值重要。而在微生物、重金属、添加剂含量等值中，而各类食品都用到添加剂，所以，添加剂比微生物、重金属等重要。 根据判断矩阵标度内容及其含义，可得V-A层判断矩阵A： A1A2 A3 A4 3）A-B层判断矩阵该层的矩阵元素确定如下： 将每条准则下各指标抽查值用分指数公式 （1.1）标准化，其中是在准则下第j个指标抽查值的均值（或最大值），是第j个抽查指标的最大允许不合格值。 对于每条准则下，令计算 （1.2） 设 （1.3）对于任意的 与相比，其判断矩阵元素由以下确定此外 （1.4） 按式（1.1）可得相应于表1.1中各指标抽查值的标准化值由式子（1.2）（1.3）（1.4）得A-B层的判断矩阵4）用方根法计算出判断矩阵A ， 的归一化特征相量: VA、VB(n)5）计算判断矩阵的最大特征值并进行一致性检验一致性指标 其中，n为判断矩阵的阶数计算一致性比例其中为平均一致性指标，当时判断矩阵具有好的一直性。6）食品质量的综合评价指数的计算2.模型的求解及结果以2013年第一至四期食品安全抽样检验产品抽查样本的32种食品的抽查值为例，求该抽查样本的食品质量安全综合评价指数。根据题目给出的32种食品抽样数据值，统计得到下面的平均值表：表1.12013年第一至四期食品安全抽样检验产品指标统计表表项目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.683330.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.12642111.0899B30.111110.553740.91751B4B(4):B1B2B3B417.724498.3325B10.1294611.27561B20.111110.7839311B30.12001111B44）用方根法计算出判断矩阵A ， 的归一化特征相量:=(，,)=0.461539 0.2307690.1538460.153846;=( ,，,)=0.732352 0.0984589 0.0889053 0.0802843=(，)=0.72337 0.104756 0.0994301 0.0724441=( ，,)=0.723233 0.108991 0.0928252 0.0749515=(,)=0.734622 0.0956177 0.0814922 0.0882682其中，的4个分量分别表示A层的4个因子(A1平均值、A2微生物平均值、A3、重金属平均值、A4、添加剂平均值)在综合评价中的各自的权重。VB(1)的4个分量分别表示只考虑年均值A1的条件下4个指标各自应占的权重；VB(2)、VB(3)、VB(4)类推。5）计算判断矩阵的最大特征值并进行一致性检验 一致性指标 其中，n为判断矩阵的阶数计算一致性比例其中为平均一致性指标，查表=0.9，故 当时判断矩阵具有好的一致性。6）市场食品质量安全的综合评价指数的计算=0.2610550.26根据深圳市食品安检局标准处推荐的评价标准（GB38382002食品安检质量标准以及食品质量分级）为参考依据，见下表1.3表1.3深圳市食品质量系数（PI）分级质量安全值级别食品质量评价<0.1良好0.10.25较好0.250.4一般0.40.55轻劣食品0.550.80重劣食品>0.80劣严重劣食品 这样就得到了2013年第一至四期食品安全抽样检验产品抽查样本的食品安全的优劣状况，并做出定量的综合评价。其它年份的食品安全情况同理可求。求解程序见附录1。表1.4为求得的各个抽查样本的食品质量安全状况和定量综合评价。其中C.R.为一致性指标(当C.R.<0.1时有好的一致性)详细的计算结果见附录2。抽查产品食品综合评价指数C.R.（食品质量安全系数）超市速食面制品0.2651150.00000000000.00000000000.00000008830.0000000000预包装调味0.2659770.00000000000.00000008830.00000008830.0000000000预包装即食小食品0.2846380.00000008830.00000000000.00000008830.0000000000一次性纸杯塑料杯0.2486960.00000000000.00000008830.00000000000.0000000883饼干及膨化食品0.258630.00000000000.00000000000.00000000000.0000000000速冻面制品0.2478460.00000008830.00000008830.00000000000.0000000000巧克力及其制品0.2610550.00000000000.00000000000.00000000000.0000000000乳制品及含乳制品0.3574770.00000000000.00000000000.00000000000.0000000000酒店奶制品0.2901320.00000008830.00000008830.00000000000.0000000000大米0.3788740.00000000000.00000008830.00000008830.0000000000食用油0.2380920.00000000000.00000008830.00000017660.0000000883广式烧卤熟肉制品0.3214150.00000000000.00000000000.00000000000.0000000000烟熏肉制品0.2946080.00000008830.00000000000.00000000000.0000000883预包装茶叶0.2751090.00000000000.00000000000.00000008830.0000000883大型企事业单位食堂0.8441380.00000000000.00000008830.00000000000.0000000000速冻汤圆食品0.2857430.00000008830.00000000000.00000000000.0000000000幼儿园食堂食品0.2705070000.00000008832013年第一至四期食品安全抽样检验产品抽查样本的食品安全情况的综合评价可由加权平均法来作出判断。 表1.4 2013年各抽查产品合格评价表 （1.5） 其中权值可以反映不同食品质量对整个深圳市食品质量的影响的不同作用。我们采用目标应用最多的因子实测值与标准值为双重判定依据的赋权方法超标指数法或指数赋权法，即根据各食品的质量分级指数来确定权重。其公式为： （i=1,2,n）式中为各食品质量安全分级的指数，为对应分级指数的最大允许值。将我们计算出的值代入式1.5计算得 故可以认为深圳市近三年来的食品质量情况应是明显提高。年份201120122013食品安全系数0.270.170.095等级问题二：“多元线性回归的数学模型”1、处理数据我们先通过MATLAB（程序见附录1）对原始数据进行检验，对残差向量进行分析，得到了残差向量分析图，剔除其中的异常点。2、设随机变量假如变量与另外8个变量，,的内在联系是线性的，它的第次试验数据是 =1,2，,8 （1）那么这一组数据可以假设有如下的结构式： （2）其中，是9个待估计参数，,是8个可以精确测量的一般变量，是38个相互独立且服从同一正态分布的随机变量，这就是多元线性回归的数学模型。 令， ， ， 那么多元线性回归的数学模型（2）可以写成矩阵形式 （3）其中是38维随机向量，它的分量是相互独立的。3、参数的最小二乘估计为了估计参数，我们采用最小二乘估计法。设分别是参数，的最小二乘估计，则回归方程为 （4）由最小二乘法知道，应使得全部观察值与回归值的偏差平方和达到最小，即使 (5)所以是的非负二次式，最小值一定存在。根据微积分学中的极值原理，应是下列正规方程组的解： (6)显然，正规方程组的系数矩阵是对称矩阵，用来表示，则，且其右端常数项矩阵亦可采用矩阵和来表示：。所以可以得到回归方程的回归系数： (7)4、由于利用偏回归平方和可以衡量每个变量在回归中所起的作用大小（即影响程度），设是p个变量所引起的回归平方和，是p-1个变量所引起的回归平方和（即除去），则偏回归平方和为：=-=-= （8）就是去掉变量后，回归平方和所减少的量。5.1.2 模型的求解1、数据筛选通过MATLAB（程序见附录1）作图如下：此时可见第八个点、第十四个点和第二十八个点是异常点，于是删除原始数据中第八行和第十四行和第二十八行数据。2、回归方程的求解 由附表1和所得的公式（7），运用MATLAB进行编程（程序见附录2），可得正规矩阵的系数矩阵为:回归系数为：， ， ， 回归方程为：3、偏回归平方和的比较运用MATLAB进行编程（程序见附录2），得到各因素的偏回方和：（）0.33840.00300.36850.10090.14291.24180.14720.1963根据的大小可判断各因素对食品安全系数的影响程度结论：在食品质量影响因素中食品产地影响最大，食品加工次之，季节影响最小，抽查地点几乎无影响。问题三：规准型抽样检验 “民以食为天，食以安为先”，食品安全和卫生问题一直是关系国计民生的重大问题。 根据问题一、二数据结果以定量比较评估的方法分析了各类影响食品安全的因素及其可能造成的危害性的问题。改进后的食品抽检的办法以主要食品为准则层建立了层次分析（AHP）模型，对影响食品安检的危害性因素做出定量分析如问题二个影响因素大小。由问题一、二结果可以通过建立抽检模型，即改进后的规准型抽样检验模型，并以蒙特卡罗法对抽检的全过程进行模拟，得到相对误差逐渐趋向于0。 就目前“营养强化.(查看更多)面粉”如何抽检问题，本题在将各检测项目与各类食品安全影响因素一一对应的基础上，建立逐步优化模型，最终得出一次全面检验的最小总成本。 依据问题一二结果通过参照计量统计检验程序给出了一套完备的连续多次跟踪抽检策略，建立更加切合实际的多目标规划的抽检策略。最终得到当检测数量达到100以上时，该部门进行抽检工作时可达到检测可靠性高、成本低、工时少的目的。故针对至少720个总批次的食品本题最终建议抽检125批次。 综上，本文较好地解决了当前食品安全问题中存在的矛盾关系。不仅提高了消费者的满意度，也兼顾了生产者的利益，具有较高的社会经济价值。6 模型的评价6.1模型的优点 本文对于各种因素对于食品质量安全的影响建立了多元线性回归模型，全面综合考虑了各个方面的因素，避免了单一因素分析的不准确性，得出了合理的数学模型。并且通过各因素的显著性分析，找到了影响食品质量安全的主要因素，较符合实际情况，模型可靠，并且模型相对简单，利于操作；该方法不仅适用于本题，也适用于其他方面的数据预测，有实际背景，可运用于实践，具有广泛适用性。6.2模型的缺点 本文忽略了除了所给因素之外的因素对食品安全的影响，与实际问题存在偏差。同时是在假设各因素相互独立的情况下对食品质量的影响进行分析，可能会导致误差7 模型的改进与推广应用统计推断先判断食品的安全性，又有中心极限定律假设食品的合格服从正态分布使其食品安全量化，并且量化成数据。比如超过量化后的数据为不合格，没有超过为合格。各个国家的食品安全等级不一样，反映深圳市也是如此。运用2010年、2011年和2012年深圳市的食品抽检数据统计量化后的数据分别作为各类产品安全等级。生活食品质量标准是指按特定生产方式生产，并经国家有关的专门机构认定，或许使用绿色食品标志的无污染、无公害、安全、优质、营养型的食品。在许多国家，绿色食品又有著许多相似的名称和叫法，诸如“生态食品”、“自然食品”、“蓝色天使食品”、“健康食品”、“有机农业食品”等。影响因素有原料生长土壤环境及食品的生产收获，此外食品的运输,加工，包装，储存，销售，以及餐饮等每个环节皆可能影响食品的质量与安全。运用统计分析中的方差，分析各个影响因子对食品安全运用统计分析中的方差，分析各个影响因子对食品安全的影响程度的影响程度。最终使人们心中无污染、无公害、安全、优质、营养型的食品与市场食品安全等级数据差极限为零。相反如果人们认识的食品安全数据比市场标准化低，或没有很好认识安全食品概念及本质。我们可以通过媒体，广告，食品安全讲座等。让人们认识和提高对安全食品实质。如果人们心目中的标准已经是市场安全等级，若是的话，就可以得出合理的结论：人们追求更加持续性，更加科学化的生活质量标准。正好和可持续科学发展观相吻合。模型中得到最优回归方程的方法是从包含全部变量的回归方程中逐次剔除不显著因子，这种方法是在不显著因子不多时采用，当不显著因子较多时，则工作量将会相当大，因为每剔除一个变量就得重新计算回归系数。鉴于以上问题，我们引入了逐步回归分析的方法，它的基本思想是将因子一个个引入，引入因子的条件是，该因子的偏回归平方和经检验时显著的。同时，每引入一个新因子后，要对老因子逐个检验，将偏回归平方和变为不显著的因子剔除。这种方法不需要计算偏相关系数，计算较简便，并且由于每步都作检验，因而保证了最后所得的方程中所有因子都是显著的。若回归方程是拟合好的，就可以进一步利用它来进行预报和控制。预报问题，用统计数学的语言来说就是一个区间估计问题。在建立气象预报、地震预报、自动控制等数学模型时，都可以用到本文的模型。 8 参考文献1 马新民，王逸迅. 概率与数理统计M. 北京：机械工业出版社，2010.2 刘卫国. MATLAB程序设计与应用M. 北京：高等教育出版社，2009.3 茆诗松. 回归分析及其试验设计M. 上海：华东师范大学出版社，1986.4陈巧玲,沈晓斌，长江水质污染预测的灰色模型，泉州师范学院学报(自然科学)，2005.115孙靖南 ,邹志红 ,任广平.模糊综合评价在天然水体水质评价中的应用研究6 .环境污染治理技术与设备 ,2005 ,6 (2) :45248.7李小燕 灰色神经网络预测模型的优化研究,378.html 8邓聚龙，灰理论基础M ，武汉:华中科技大学出版社,2002.9傅立，灰色系统理论及其运用M，北京:科学技术文献出版社, 1992.10刘荣珍,赵军,模糊评价模型在长江水质评价中的应用, 兰州交通大学学报,2007.1212C.Henry Edwards &David E.penney微分方程及边值问题，北京：清华大学出版社，2007附录：一数据统计名称抽取合格不合格2013.4超市速食面制品59518预包装调味品12010416预包装即食小食品544311一次性纸杯塑料杯50500饼干及膨化食品100928糕点42402速冻面制品48480年糕21201巧克力及其制品30300乳制品及含乳制品505002013.3酒店奶制品21210面制品877611盐焗食品795920年糕1073大米51510食用油49490广式烧卤熟肉制品1007822烟熏肉制品807642013.2预包装茶叶100991熟肉制品53521干货食品20018515速冻汤圆食品19190炒货1511510糖果37370乳制品及含乳制品505002013.1特色菜肴1221226甜品店食品1001000港式早餐20017228烧烤肉制品15013515幼儿园食堂食品1731730自制饮料2001991大型企事业单位食堂20019912012.26即食蔬菜2002000学校学生餐3003000娱乐场所食品30255水发水产34340潮式打冷食品60546自制米面34340油炸食品60600即食水果1001000快餐食品1501500城中村快餐食品1501500中央厨房半成品61601餐饮火锅汤底83830学校食堂面制品1008713西式菜肴2001982餐前小吃食品2001928餐饮领域盐焗类食品880工业区餐饮1001000餐饮领域乳及其乳制品27270一次性餐饮具69654腌腊肉制品55532餐饮领域茶叶20019912012.25饼干79736干果及坚果食品65641预包装肉制品65623动物性水产干制品、藻类干制品807010交通枢纽商铺预包装食品1981899糕点面包1009010预包装龟苓膏27270茶叶1201200焙烤咖啡及其制品27270饮料酒47470乳制品及含乳制品1461460火锅原料食品20016238罐头50482调味品100964塑料奶瓶橡胶奶嘴16160大米1501500食用植物油44440火锅底料产品27252腌腊肉制品252502012.24餐饮自制饮料2002000餐饮大米1001000餐饮食用油100982餐饮西式菜肴100973餐饮生食水产品1361333企事业单位食堂分装菜肴.饭盒1001000罗湖区.龙岗区餐饮肉丸181622012.23食用植物油80773腐竹及淀粉制品5553