肖枝洪,赵忠校
(重庆理工大学数学与统计学院,重庆 400054)
食品安全问题关系着国计民生,是每一个消费者都不得不慎重考虑的事情。国内外关于食品安全问题的文献大多以政府监管、食品生产源头、食品加工等为切入点进行探索研究,并针对这些问题提出建议和防范措施。例如,徐晓新在《中国食品安全问题、成因、对策》中主要从食品生产与流通链中存在的信息不对称和食品加工过程中质量检测及控制不完善两个方面阐述食品安全问题的产生[1]。潘淑红在《中国食品安全政府监管体系研究》中重点从食品安全监管、检测、法律等体系方面说明中国食品安全存在的问题,并提出完善这些体系的建议[2]。M.J McLaughlin,D.R Parker,J.M Clarke 在《Metals and micronutrients-food safety issues》探索了重金属以及食品添加剂对食品安全的影响[3]。
重庆市近年来曝光的假蹄髈、毒血旺、化学火锅等事件,严重影响和危害了广大居民的生命健康。本文以顾客满意度理论作为理论基础,以重庆市居民食品安全满意度为切入点,通过问卷调查方式获得数据,结合多元统计分析方法,运用SAS9.1和AMOS19.0构建具有因果关系、多变量的结构方程模型,用于评价重庆市居民食品安全满意程度。通过分析模型,探索满意度的影响因素,希望在普及食品安全科学知识和构建安全的食品生产、运输、消费环境等环节上提供理论分析和科学性建议。
首先,确定调查研究的主要指标。根据国际食品卫生法典委员会(CAC)对食品安全的定义[4],考虑从食材到制成食品再到消费市场这个过程中参与者的职能及多方理论资料[5-6],归纳出20个可能影响食品安全的指标或因素。这20个指标分别为:质量状况、卫生程度、营养含量、新鲜程度、生产技术认证、食品加工过程、品牌生产、食品储藏与运输、食品安全有关法律、相关部门执法力度、重金属含量、添加剂使用、农药残留、政府政策、监管体系、经营及卫生许可认证、食品生产源头、食品包装、媒体监督、舆论评价。这些作为建模的观察变量,可用来分析食品安全满意度。
其次,量化测评指标。围绕食品安全满意度来拟定问题和量表,针对各个因素设置具体明确的问题。每个影响因素的满意度在问卷中均采用5级Likert量表进行分级测量,并对量表的内容和评分进行标准化:“非常不满意”赋值1分,“不太满意”2分,“一般”3分,“比较满意”4分,“非常满意”5分。所赋分值越高表明受访者对该因素的满意度越高。
最后,确定问卷的版式。为了使问卷易于接受,在顺序安排上,先提问受访者对食品安全的整体满意度,再测量影响食品安全满意度的各个因素的满意度。本次问卷主要分为以下3个层次:第一部分是被调查者对食品安全总体满意度;第二部分是具体因素问题与选项,反映食品安全满意度的重要信息;第三部分是受访者个人基本信息,包括性别、年龄、学历、职业、家庭收入、户口等。
根据抽样调查理论,在完成问卷设计后,我们随机地抽取50位居民进行试调查。然后根据信度[6]的定义,用SAS 9.1对试调查所得数据做因子分析和可靠性检验。探索性因子分析的KMO值为0.694,Bartlett值为544.095,P值小于0.001,输出的4个因子解释20个变量的累积百分比是65.722% 。根据旋转成分矩阵,我们得到4个满意度因子:治理监管、有害物质、质量状况、舆论监督。由表1可知,4个满意度因子的 α系数基本均在0.7以上,说明量表内在信度是可以接受的。此外,总量表的α系数高于0.9,说明量表的内在信度颇佳。因此,调查问卷通过了信度检验。由以上检验可知,问卷通过了信度与效度检验,可以进行正式调查。
表1 信度检验统计量表
正式调查时发出问卷300份,共回收有效问卷200份。在所有受访者中,男女居民比例为21∶19;城镇与农村户口比例为53∶47。受访者涉及到了各个年龄阶层和各个职业范围(见图1和图2),可以认为数据具有普遍代表性。
图3为重庆市居民对食品安全总体满意度的统计图,由图3可知居民对食品安全总体满意度基本呈正态分布,均值为2.68,表明重庆市居民对食品安全总体满意度介于不太满意和一般之间,满意度偏低。
对总体满意度进行交叉列联表分析,得出如下结果:从性别上分析,女性的总体满意度为2.74,高于男性的总体满意度2.62;从职业情况来看,农民的满意度3.75为最高,其次是农民工的满意度2.875,医生的满意度最低2;从学历来看,学历在硕士及以上水平的,对食品安全满意度最低,为1.875,满意度最高的是中专或高中学历,为3.053;从户口所在地来看,农村户口的居民对食品安全满意度为2.564,城镇户口的居民的满意度为2.811。
图1 职业分布
图2 年龄分布
图3 食品安全总体满意度
用正式调查得到的200份问卷数据做因子分析[9],我们提取出4个因子:质量状况因子、有害物质因子、治理监管因子和舆论监督因子作为结构方程模型的潜在变量,将衡量食品安全满意度的指标变量分为三级,详细情况见表2。
表2 食品安全满意度的指标体系
根据顾客满意度理论和结构方程模型原理[7-8],建立如下的测量模型和结构模型。
外生变量测量模型:
内生变量测量模型:
结构模型:
模型中,η1、η2、η3是内生潜在变量,分别代表治理监管、有害物质、质量状况;ξ1是外生潜在变量,代表舆论监督;xi,i=1,2,,3,4,代表外生观察变量,yj,j=1,…,16,代表内生观察变量;λi(i=1,…,20)表示观察变量和潜在变量间的回归系数;βij表示内生潜在变量间的回归系数;γ11表示外生潜在变量对内生潜在变量的影响;δi(i=1,2,3,4)代表外生观察变量的误差,εi(i=1,…,16)表示内生观察变量的误差;ζi(i=1,2,3)表示结构方程的残差项,反映了在方程中未能被解释的部分。
本文使用Amos19.0软件,逐步探索设置符合此次研究主题的因果路径图,得到标准化后具有路径系数的路径图,并通过不断地建立变量之间的联系消除路径的偏差对模型加以修正,最终得到模型图(见图4)和检验结果表(见表3)。
图4 修正模型路径图
从表3可知,卡方对应的p值为0.051,表现为不显著,不拒绝原假设,即理论模型和实际数据可以契合。再从其他适配度指标来看,模型绝对适配度指数GFI、AGFI,理论上大于0.9表示通过检验。增值适配度指数 NFI、IFI、CFI,理论上大于0.9表示模型通过检验,越接近1越好。简约适配度指数PGFI、PCFI,理论上大于0.5表示模型通过检验。CMIN/DF是χ2自由度比,理论上小于2.00表示模型通过检验。由表3可知,理论因果模型图与实际数据可以适配。
表3 模型检验与拟合优度结果
潜在变量间的路径系数表示某一变量的变动引起其他变量变动的程度。如图4中所示,舆论监督因子与治理监管因子之间的路径系数为0.8,表示舆论监督因子满意度提高1个百分点将直接使治理监管因子满意度提高0.8个百分点;同理,治理监管因子满意度提高1个百分点将直接使质量状况因子满意度提高0.61个百分点;治理监管因子满意度提高1个百分点将直接使有害物质因子满意度提高0.67个百分点。因此,在提高居民食品安全满意度的过程中,管理部门在制定相关政策法规时,除了强调综合治理外,还应该有重点地加强专项治理力度,比如更加积极地倡导和鼓励对食品安全的舆论监督等。
通过食品安全满意度指标中潜在变量与观察变量之间关系的分析,可以发现与潜在变量关系紧密的观察变量,还可进行各因子内观察变量间的比较。
1.治理监管满意度与观察变量之间的关系。
治理监管因子中食品加工监管系数(0.76)最大,生产技术认证系数和监管体系系数(0.75)次之,接下来是食品储藏与运输系数(0.73)、经营及卫生许可证系数(0.72)、政府政策系数(0.65)、相关部门执法力度系数(0.63)、食品安全有关法律系数(0.62)。以上说明政府以及食品生产厂商在食品安全监管方面的健全与完善对稳定居民消费信心至关重要。同时,规范的生产技术认证、科学合理的食品储藏、运输以及加强法制管理是治理监管食品安全工作的实际可行性的重要保障[10-11]。
2.质量状况满意度与观察变量之间的关系。
食品质量系数(0.81)最大,卫生程度系数(0.78)次之,然后是食品新鲜程度系数(0.77)。以上说明食品质量以及卫生程度和新鲜程度是影响居民对食品质量状况满意度的主要因素。另外,质量状况因子不仅受自身各个因素的直接影响,而且还受治理监管和有害物质的间接影响。因此,要提高居民对食品质量状况的满意度,在保证自身各个影响因素的满意度的同时,还要从治理监管等方面着手,综合管理食品安全。
3.有害物质满意度与其观察变量之间的关系。农药残留系数(0.87)最大,添加剂使用系数(0.84)次之,然后是重金属含量系数(0.72)。3个观察变量的路径系数之高,表明其对有害物质因子满意度影响之显著。又由于高危害农药残留或有毒添加剂关乎生命的安全,居民对有害物质因子特别重视,但满意度却最低。因此,合理使用食品添加剂和监管含有有害物质的食品的流通是当前治理监管部门要重点考虑的,有关部门应在质量技术检测方面加强监督或改善技术检测方法。
4.舆论监督满意度与观察变量之间的关系。
媒体监督系数为0.71,舆论评价系数0.70,品牌产品系数0.62,食品包装系数0.62。这些都说明居民对媒体及社会舆论发挥的监督作用是比较认可的,品牌产品无论从生产还是包装上媒体口碑都比较好,对舆论监督满意度也有促进作用。因此,大力倡导媒体对食品安全的监督,积极鼓励舆论评价及品牌产品的自我提升,对提高居民食品安全满意度有积极的作用。
[1]徐晓新.中国食品安全问题、成因、对策[J].农业经济问题,2002(10):12-18.
[2]潘淑红.中国食品安全政府监管体系研究[M].北京:中国政法大学出版社,2008.
[3]McLaughlin M J,Parker D R,Clarke J M.Metals and micronutrients-food safety issues[J].Field crops research,1999(1):234-241.
[4]高阳,杨薇.我国食品安全、问题及对策[J].中国食品与营养,2009(1):25-29.
[5]苏理云,周林招,王雪娇,等.基于结构方程模型的大学生食品安全满意度调查[J].重庆理工大学学报:社会科学,2012(10):57-60.
[6]梁一鸣,张钰兰,董西钏.基于结构方程模型的杭州城镇居民食品安全满意度评估[J].浙江财经学院,2009(9):76-81.
[7]李攀艺,牛弘欢,谭建伟.基于结构方程模型的酒店员工工作满意度研究[J].重庆理工大学学报:社会科学,2012(1):36-41.
[8]吴明隆.结构方程模型——Amos的操作与应用[M].重庆:重庆大学出版社,2009.
[9]肖枝洪,朱强,苏理云,等.多元数据分析及其R实现[M].北京:科学出版社,2013.
[10]徐信贵.食品安全风险警示的行政法问题[J].重庆理工大学学报:社会科学,2014(6):50-54.
[11]张浒,陈刚.食品安全监管制度的优化及其实现逻辑——基于上海羊肉掺假的个案分析[J].重庆理工大学学报:社会科学,2014(10):78-84.