葡萄酒评价解题思路.docx

葡萄酒评价解题思路
第一问
先对数据标号，标号方式参考如下：
红葡萄酒
澄清度
色调
香气纯正度
香气浓度
香气质量
口感纯正度
口感浓度
持久性
口感质量
平衡/整体评价

白葡萄酒
澄清度
色调
香气纯正度
香气浓度
香气质量
口感纯正度
口感浓度
持久性
口感质量
平衡/整体评价

还可以更细得标号，用表示第组第位品酒师对第种红葡萄酒第项指标的评分，用表示第组第位品酒师对第种白葡萄酒第项指标的评分，比如第1组第0位品酒师对第21种红葡萄酒第7项指标（口感浓度）的评分就是，第2组第5位品酒师对第8种白葡萄酒第2项指标（色调）的评分就是。
若是评价“哪一组结果更可信”的话，应该考虑一下系统误差和偶然误差，系统误差小的结果比系统误差大的结果可信，偶然误差小（数据比较集中）的结果比偶然误差大的结果可信。
比如说，第号红酒澄清度的实际分值（带有主观性，不过根据大数定理，无穷多个品酒师的评分的数学期望就是实际分值）为4分，那么同样是对第号红酒澄清度的打分，第一组打了7个4分，2个5分和1个3分，第二组打了6个4分，1个5分和3个3分，那么第一组的评分的数学期望就是4.1分，第二组的评分的数学期望就是3.8分，第二组的系统误差更大，第一组更可信；再比如说第一组打了5个4分，3个5分和2个3分，第二组打了7个4分，2个5分和1个3分，那么两组的评分的数学期望都是4.1分，不过第一组的数据比较分散，偶然误差比较大，第二组更可信。
不同的品酒师的个人感受不可能完全一样，评分标准掌握尺度也有差异，因此难免有主观误差（系统误差的一类），不过如果品酒师是随机分配到两组的话，多数情况下可以认为不同的人的系统误差相互抵消（下文说不能相互抵消的话怎么办），因此重点考虑偶然误差。
可以通过统计学中的理论（需要用到分布）得出同样的置信水平（可以设）下每一个统计量（同一组人对同一种酒同一项指标的评分）的置信区间，然后求出置信区间跨度（置信上限与置信下限的查，设为和，与和对应，比如就是第二组的品酒师对第七种红葡萄酒口感浓度的评分的置信区间跨度），跨度小的偶然误差小。
可以求出同一组品酒师评价同样的指标的置信区间跨度的平均值（但目前不知道应不应该分红葡萄酒和白葡萄酒），即、，然后比较和的大小，以及和的大小，小者偶然误差小，更可靠，这样就可以得出两组品酒师对同一个品种（一共两个品种）的酒的同一项指标（一共十个指标）的评分哪个更可靠了。
如果计算不方便的话，可以不算置信区间跨度，而算标准差，这样和就是第组品酒师对第种酒第项指标的评分的标准差了，而平均标准差越小的数据越可靠。
当然，系统误差也可能存在，如果两组品酒师对同一个品牌的酒的同一项指标的平均打分（即和，或和）差距比较大，说明其中一组存在较大的系统误差，或者两组都存在较大的系统误差，此时怎么办？有两种办法，一种是在比较和的时候删去这个品牌的酒的该指标数据，另一种是比较分析哪种更有可能有系统误差（但如何分析还没想好）。
此外，品酒师水平参差不齐，评价尺度也不一样，不过对于主观打分评判来说，都宽松和都严格是公平的，但宽严尺度不一（不一定是“黑哨”）就有问题了，比如三种酒的“客观”评分应该是70、80、90，这样给它们分别打75、85、95分和分别打65、75、85分对它们来说是公平的，但都打80分就不公平了。
因此，需要检验一下品酒师的水平。我的方法是，把同一个人对同一项指标的打分以品牌为单位一一列出来（、……、、……），然后求出二十个人对各个品牌这项指标打的分数的平均值，也一一列出来，再求二者的相关系数，如果相关系数太小，说明这位品酒师对这项指标的打分存在不公平的情况，数据有问题（但我没想好如何删数据）。
至于显著差异，不好意思，我没带书，不记得怎么求显著差异了。
第二问
第一问已经求出哪组品酒师对同一个品种的酒的同一项指标（红葡萄酒指标0、红葡萄酒指标1、红葡萄酒指标2……白葡萄酒指标9、白葡萄酒指标0、白葡萄酒指标1、白葡萄酒指标2……白葡萄酒指标9）的评分更可靠了，这样可以得到一个相对客观的打分（个人建议，如果两组数据存在显著差异的话，用更可靠的那组；如果不存在显著差异的话，用两组平均值），把同一个品牌的酒的十项指标得分累加起来，就是酒的总分。
分级的时候，总分越高的酒对应的葡萄的等级越高。
第三问
首先，附件2的数据该求平均值的一律求出来平均值。
对葡萄的各类理化指标分别编号、……，对葡萄酒的各类理化指标分别编号、……，然后进行拟合数据
令

……

则对于每一个，都有28组数据可以用来拟合，从而得到葡萄酒理化指标与葡萄理化指标的关系。
第四问
还是拟合，只不过这次酒的评分是因变量，葡萄的理化参数、葡萄酒的理化参数都是自变量。
根据我的经验，应该是同一个理化指标在一定范围内与得分正相关（不一定是线性相关），但达到“峰值”之后就负相关了。
这里还存在一个问题：可能有的理化指标对有利于提高葡萄酒指标甲（色调、香气浓度等）的得分，却不利于提高指标乙（色调、香气浓度等）的得分，因此需要分指标讨论；但这样一来计算量将难以估量。
我目前能想到的就是这些了，我的思路需要的计算量非常大（虽然是计算机处理），希望高人能简化一下。

判定哪组更可信的相关程序
function y= xlsread（’D \cumcm2012problems\A’）
global a10 a20 a30 total k;
syms k a10;
x1=dsolve('Dx1=-0.75*x1','x1(0)=a10');
t=1;a20=subs(x1);
%************************************************************
x21=dsolve('Dx21=-(0.75+0.32*k)*x21','x21(0)=a20');
t=2/3;a21=subs(x21);
x22=dsolve('Dx22=-0.75*x22','x22(2/3)=a21');
t=1;a30=subs(x22);
x31=dsolve('Dx31=-(0.75+k)*x31','x31(0)=a30');
t=2/3;a31=subs(x31);
x32=dsolve('Dx32=-0.75*x32','x32(2/3)=a31');
nn=1.109*10^11*(0.5*a21+a31);
eq1=a10-nn*1.22*10^11/(1.22*10^11+nn);
S=solve(eq1,a10);a10=S(2);
syms t;
t3=subs(subs(int(0.32*k*x21,t,0,2/3)));
t4=subs(subs(int(k*a31,t,0,2/3)));
total=17.86*t3+22.99*t4;
k=x;
y=subs(-total);
%**************************
function y= problem1(x)
global a10 a20 a30 total;
[k,mtotal]=fminbnd('buyu',16,18);
ezplot(total,0,25)
xlabel('组别')
ylabel('系列')
title(' ')
format long;
k
total=-mtotal
a10=eval(a10)
a20=eval(a20)
a30=eval(a30)
format short
clear
大家可根据实际情况做相应修改！