葡萄酒评价问题2012a研究

112：葡萄和葡萄酒的理化指标（2）问题一1281所在xls选中“计算结果”表单。鼠标单击【视图】->【】->【控件工具CommandButton1VBAPrivateSubEndVBA代码了。VBAPrivateSubCommandButton1_Click()Dimi,j,kAsIntegerDimInfo(30)AsDimx(30,10),All(30)10Total=Cells(5,1)Fork=1Topos=3+14*(k-1) '获得第kInfo(k)Sheets("第一组红葡萄酒品尝评分").Cells(pos,1)Forj=1To1010s=ForItem=1TossSheets("第一组红葡萄酒品尝评分").Cells(pos1Item,2NextItem'计算第jx(k,j)=s'第k个样品10个品酒员的得分NextjNextkFork=1ToCells(6k,1) '显示第一组红葡萄酒第kForj=1ToCells(6+k,2+j)=x(k,'显示第一组红葡萄酒第k10NextjNextkCells(37,1)Fork=1ToTotalpos=3+14*(k-1)Info(k)=Sheets("第二组红葡萄酒品尝评分").Cells(pos,Forj1To10s=0ForItem=1TossSheets("第二组红葡萄酒品尝评分").Cells(pos1Item,2NextItemx(k,j)=sNextjNextFork=1ToTotalCells(38+k,1)=Info(k)Forj=1ToCells(38+k,2+j)=x(k,j)NextjNextkEndPrivateSubCommandButton2_Click()Dimi,j,kAsIntegerDimInfo(30)AsDimx(30,10),All(30)10Total=Cells(68,1)Fork=1Topos=4+13*(k-1)'获得第kInfo(k)=Sheets("第一组白葡萄酒品尝评分").Cells(pos,Forj=1To 10s=ForItem=1TossSheets("第一组白葡萄酒品尝评分").Cells(posItem,3Next '计算第jx(k,j)=s'第k个样品10个品酒员的得分NextjNextFork=1ToCells(69k,1)Info(k)'显示第一组白葡萄酒第kForj=1ToCells(69+k,2+j)=x(k,'显示第一组白葡萄酒第k10NextjNextkCells(100,1)="第二组白葡萄酒得分"Fork=1Topos=4+12*(k-Info(k)=Sheets("第二组白葡萄酒品尝评分").Cells(pos,Forj=1Tos=ForItem=1TossSheets("第二组白葡萄酒品尝评分").Cells(pos1Item,4NextItemx(k,j)=sNextjNextFork=1ToTotalCells(101+k,1)=Info(k)Forj=1ToCells(101+k,2+j)=x(k,j)NextjNextkEnd图 xls表单中的样14品酒品酒员品酒品酒品酒品酒员品酒品酒品酒品酒品酒品酒品酒12员员员员6员员9员号号号号号号23456789品酒品酒品酒品酒品酒品酒品品酒品品酒品酒员号员号3员号5员号员号酒员号员号23456789品酒品酒品酒品酒品酒品酒品品酒品品酒品酒员号员号3员号5员号员号酒员号员号23456789品酒品酒品酒品品酒品酒品酒品品酒品品酒品酒员号员号3员号5员号员号酒员号员号23456789A：酒样品；B：组别； m SST(xijki1j1kSSA是酒样品因素，SSB是组别因素，SSC是品酒员因素，SSEA1:SAS8TheGLMDependentVariable:SumMeanFPr>Corrected R-CoeffRootValueTypeIIIMeanFPr>19F=19.46,p<0.0001，说明对红葡萄酒，两组人员A2: Mean F Pr>1CorrectedR-CoeffRootValueAnovaMeanFPr>1 A：酒样品；B：组别； 白葡萄酒SAS8结果TheGLMDependentVariable:SumMeanFPr>CorrectedR-CoeffRootValueTypeIIIMeanFPr>19F=22.27,p<0.0001，说明对白葡萄酒，两组人员 TheANOVAProcedureDependentVariable:ValueSumMeanFPr>1 CorrectedR-CoeffRootValueAnovaMeanFPr1A：酒样品；B：m m

x)2

i1SSASSBSSESum Mean FPr> Corrected CoeffVar RootMSE ValueMean TypeIII Mean FPr> F1=11.39,p1<0.0001，说明红F2=7.48，p2<0.0001，说明人员之间有差异。SumMeanFPr>CorrectedR- Coeff Root Value 9 9F2=15.45，p2<0.0001，说明人员之间有差异。F1F2FFFSum Mean F Pr> Corrected CoeffVar RootMSE ValueMean TypeIIIMeanFPr>9SumMeanFPr>Corrected CoeffVar RootMSE ValueMean TypeIII Mean FPr> 对白葡萄酒，第二组人员评价指标F2=7.4199<第一组人员评价指标F1=8.2426，SAS操作,Solutions->ysis->yst.Statistics->ANOVA->FactorialANOVA.问题二：按酿酒葡萄的理化指标与葡萄酒的质量评价对酿酒葡萄进行分级。我们先对酿酒葡萄的理化指进行主分分析于本题给的数关于酿酒萄的理指标数繁多而很多理化标反映因素是似的或是互相含，例如果颜色主通过花苷反映所以我们以剔除部分不要的理指标我们过主成分析法可提取出干个互不相关的综指标以替原指标按照各标的贡率进行得到各葡萄的化指标评价，最后葡萄酒质量评进行 运算后以到酿酒葡的综合分。

,x28，评价对象有n=28个第ijxijxijxijxijxxijs

，(i1, ,m

j1, ,n 1其中x

，s

(x

)2xsjnn nn

xi=xixii1 Rr)mn

nnxki ,(i,j1,2, ,mn式中rii

rijrji

rij是第ij计算相关系数矩阵R的特征值12m0及对应的特征向量u1,u2 ,um其中

(u,u

)Tm 1 2 y1u11x1u21x2un1y2u12x1u22x2un2y x x

2m

nm对于红、白葡萄，我们提取的综合指标如下表所示（110号葡萄的数据，：4.2-19个主成分（部分123456789123-----------4-----------5678-----------9p(p≤m)计算特征值j(j1, bj

jmmp

,(j1 m)yj 为主成分y1,y2, ,yp 表格方差bjp由上表可看出，当取到第9p已经达到了86.263%>85%)9y1y2,,y9928个9个主成分进行综合分析。表格方差bjp由上表可看出，当取到第11p已经达到了87.225%85%)11z1z2,,z11112811个主成分进行综合分析。ppZ=bjyj jZ进行了归一化处理，得到如下数4.2-5Z12345678904.2-6Z1234567890.S1Sa1S1a2S2,其中a1a20.5。表格123456789表格S2SS2S12-34567894.2-90.7-0.5-0.2-0.3-0-0-4.2-100.75-0.7-0.65-0.6-0.4-0.4-0-0-分级结果与葡萄的质量进行对比，可以看出我们的分级比较好地符合了评酒员从身经验发得出较为符酒人喜好评价结果这些属感官指这种客因素和酒员评的感官标种因素而使得级结果加科SAS软0.01水平上双侧显著相关的指标作为联系紧密的标准。表16DPPH表17表18DPPH表19DPPH表20表21DPPHDPPH表22表23表240.05水平（双侧）上有部分显著相关：25DPPH（0.05DPPH,XqT,XqT 设随机向量xX,X ,XT，y

xyx

12y

22pq，11pp22q T是两组变量之间的协方差阵。且当 12与21UaTxaXaX a 11 1

xy及给定的条件下，选取ab使U1与V1 covU1,V1

covaTx,bTyvaraTxvarbTy后进行线性回归分析。所有数据处理在统计软件SAS8中，都先采用逐步回归，1生成用于SAS导入的程序a2012_2,m，SAS数据文件a2012_44R20.96RMSE=0.12709表

概率Pr--x13表示白藜芦醇含量，x14代表黄糖醇含量，x19代表可滴定酸含量，x27代表5.4.1R20.9026RMSE=0.19014z0.65067y60.43483表

概率PrDPPH半抑制体积（IV50