拟合统计量
回归方程检验
回归系数检验
回归诊断
图12-回归分析2
结果分析:
从结果可以分析可以分析出,回归模型中各参数都通过了显著性检验,而且模型整体的拟合效果也不错,所以我们所建立的回归方程是合理的。从回归方程中可以看出,b5-酒石酸、b10-DPPH自由基1/IC50、b20-可滴定酸、b29=果皮颜色a*(+红;-绿)、b30-果皮颜色b*(+黄;-蓝)对葡萄酒质量的影响比较大。
5.4.2 论证能否用葡萄和葡萄酒的理化指标来评价葡萄酒的质量:
查找了相关文献,葡萄酒的质量可以用理化指标和感官指标进行评价,所以我们将酿酒葡萄和葡萄酒的感官指标当做变量引入方程,从而去验证能否只用葡萄和葡萄酒的理化指标来评价葡萄酒的质量,如果添加了感官指标后的拟合度比没有添加的拟合度高,各参数都通过显著性检验,则说明感官指标对葡萄酒的质量的评价也有影响。 建立回归模型:
z??0??1x1??2x2???nxn??n?1a1??n?2a2???2nan?? (6) 其中an为感官指标
运用SAS软件,利用附录的理性指标和感官指标对葡萄酒质量解出回归方程结果:
16
图13-回归分析3
上面是红葡萄和红葡萄酒的理性指标和感官指标对葡萄酒质量解出回归方程结果,回归方程的拟合度为98.62%,拟合度明显比没加上感官指标时78.89%高出很多,而且模型中的各参数都通过了显著性检验,所以上面的结果论证了不能用葡萄和葡萄酒的理化指标来评价葡萄酒的质量。
六、 模型评价与检验
6.1 模型的优点:
本文从题意出发,充分运用SPSS、SAS、MATLAB等程序并绘制了较多图表,在解题过程中采用了多种算法,包括配对样本T检验、聚类分析、因子分析、典型相关分析等法。使得模型叙述较为流畅。 6.2 模型的缺点:
由于时间紧迫,数据量巨大,我们在计算的过程中忽略了二级指标进行分析,这不免使得模型相对来说比较粗糙,同时由于数据量较大的关系使得计算过程中难免存在误差。 6.3 模型的检验
经检验证实模型均基本符合模型建立初衷。
七、 参考文献
【1】 李华等,运用多元统计分析确定葡萄感官特性的描述符,中国食品学报第七
卷第四期,2007.8
【2】 刘曙东等,葡萄酒感官评价结果的统计分析方法研究,中国食品学报第六卷
第二期,2006.4
【3】 李运等,统计分析在葡萄酒质量评价中的应用,酿酒科技2009年第四期,
2009
【4】 王英臣,基于主成分分析和人工神经网络的酒类辨识,河北工程大学学报第
26卷第2期,2009.6
17
八、 附录
附录1
正态性检验的Matlab程序:
x1=[62.7 80.3 80.4 68.6 73.3 72.2 71.5 72.3 81.5 74.2 70.1 53.9 74.6 73 58.7 74.9 79.3 60.1 78.6 79.2 77.1 77.2 85.6 78 69.2 73.8 73];
x2=[68.1 74 74.6 71.2 72.1 66.3 65.3 66 78.2 68.8 61.6 68.3 68.8 72.6 65.7 69.9 74.5 65.4 72.6 75.8 72.2 71.6 77.1 71.5 68.2 72 71.5];
x3=[82 74.2 79.4 71 68.4 77.5 71.4 72.9 74.3 72.3 73.3 65.9 72 72.4 74 78.8 73.1 72.2 77.8 76.4 71 75.9 73.3 77.1 81.3 64.8 81.3];
x4=[77.9 75.8 76.9 81.5 75.5 74.2 72.3 80.4 79.8 71.4 72.4 73.9 77.1 78.4 67.3 80.3 76.7 76.4 76.6 79.2 79.4 77.4 76.1 79.5 74.3 77 79.6]; subplot(2,2,1);h=normplot(x1); title('第一组红酒的正态性检验'); subplot(2,2,2);h=normplot(x2); title('第二组红酒的正态性检验'); subplot(2,2,3);h=normplot(x3); title('第一组白酒的正态性检验'); subplot(2,2,4);h=normplot(x4); title('第二组白酒的正态性检验')
18
附录2
各组葡萄酒平常评分表: 酒样品1 酒样品2 酒样品3 酒样品4 酒样品5 酒样品6 酒样品7 酒样品8 酒样品9 酒样品10 酒样品11 酒样品12 酒样品13 酒样品14 酒样品15 酒样品16 酒样品17 酒样品18 酒样品19 酒样品20 酒样品21 酒样品22 酒样品23 酒样品24 酒样品25 酒样品26 酒样品27
评酒员1 51 71 80 52 74 72 63 64 77 67 73 54 69 70 69 72 70 63 76 78 73 73 83 70 60 73 70 评酒员2 66 81 85 64 74 69 70 76 78 82 60 42 84 77 50 80 79 65 84 84 90 83 85 85 78 80 77 第一组红葡萄酒品尝评分 评酒评酒评酒评酒员3 员4 员5 员6 49 54 77 61 86 74 91 80 89 76 69 89 65 66 58 82 72 62 84 63 71 61 82 69 76 64 59 84 65 65 76 72 76 82 85 90 83 68 75 73 72 63 63 71 40 55 53 60 79 59 73 77 70 70 80 59 50 58 51 50 80 71 69 71 91 68 97 82 49 55 52 57 84 66 68 87 76 68 82 79 96 71 69 60 72 68 93 72 86 80 95 93 90 68 90 84 81 62 70 67 71 61 78 71 63 64 80 76 评酒员7 72 83 73 76 68 69 72 69 76 75 70 47 77 76 56 80 69 62 80 76 79 75 81 70 64 72 73 评酒员8 61 79 83 63 84 64 59 85 92 68 66 61 76 76 60 74 80 58 78 76 73 77 91 75 62 76 67 评酒员9 74 85 84 83 81 81 84 75 80 76 90 58 75 76 67 78 81 70 82 86 86 79 84 78 81 79 85 评酒员10 62 73 76 77 71 84 84 76 79 75 73 69 77 76 76 74 76 68 81 81 74 80 78 70 67 77 75 19
酒样品1 酒样品2 酒样品3 酒样品4 酒样品5 酒样品6 酒样品7 酒样品8 酒样品9 酒样品10 酒样品11 酒样品12 酒样品13 酒样品14 酒样品15 酒样品16 酒样品17 酒样品18 酒样品19 酒样品20 酒样品21 酒样品22 酒样品23 酒样品24 酒样品25 酒样品26 酒样品27 酒样28 评酒评酒员1 员2 85 80 78 47 85 67 75 77 84 47 61 45 84 81 75 46 79 69 75 42 79 46 64 42 82 42 78 48 74 48 69 49 81 54 86 44 75 66 80 68 84 49 65 48 71 66 82 56 86 80 75 66 58 40 66 75 第一组白葡萄酒品尝评分 评酒评酒评酒评酒评酒员3 员4 员5 员6 员7 88 61 76 93 83 86 54 79 91 85 89 75 78 75 136 80 65 77 83 88 77 60 79 62 74 83 65 78 56 80 83 66 74 80 80 81 54 81 59 73 81 60 70 55 73 86 60 87 75 83 85 60 74 71 86 75 52 67 62 77 83 49 66 65 76 84 67 79 64 78 87 71 81 61 79 86 65 70 91 87 90 70 78 71 87 83 71 72 71 85 83 68 73 64 80 82 71 83 81 84 85 59 76 86 83 90 58 72 77 76 80 69 80 82 78 79 73 67 59 68 82 69 74 67 77 82 75 93 91 81 79 67 59 55 66 89 69 88 87 85 评酒员8 80 68 79 78 74 67 68 77 81 73 62 56 62 68 67 62 74 64 63 62 70 70 71 78 78 76 74 76 评酒员9 95 73 90 85 79 65 77 85 76 91 88 68 65 81 74 84 92 74 73 87 88 80 87 86 77 90 73 88 评酒员10 79 81 79 86 74 84 82 83 85 71 72 70 69 73 82 77 91 81 77 80 84 74 75 85 81 84 77 90
20