《数据分析》实验报告三

＞> d=（2＊ｐ^２+3*p-1)*(1/(ｎ1-1)+1／（n2-1)-1/(ｎ-k)）／（6*(p＋1)*（k-1)）;
>> p1=n１／n;p2=n2／n;
>> ｍ１=ｍeaｎ(G1);m2=meaｎ(G2);
＞＞ s1=cov(G1）;ｓ2=cｏv(G2);
＞> s＝((ｎ1－1)＊ｓ1＋（n２-1)*s2)／（ｎ-k);
Ｐ＾=(N1+N2）/(m+n)
误判率的回代估计易于计算。但是，p^是由建立判别函数的数据反过来用作评估准则的数据而得到的。所以有偏，往往比真实误判率小。当训练样本容量较大时,p^可以作为真实误判率的一种估计。
４)交叉误判率估计是每次剔除一个样品,利用m+n-1个训练样本建立判别准则，再利用建立的准则对删除的样本进行判别。对每个样品做如上分析,以其误判的比例作为误判率，步骤；
N1１=Ｎ1１+１;
enｄ
end
end
eｎd
尝试引用非结构体数组的字段（我没辙了,实在不知道哪错了)
五、实验结果及分析
此题用贝叶斯判别法分析效果明显
说明:此部分的内容和格式可根据实验课程的具体需要、要求自行设计和确定相关栏目。
3．对第一类与第二类的先验概率分别为7／１5、8/15，计算误判概率。
一、实验目的
1．熟练掌握MATLAＢ软件进行距离判别与贝叶斯判别的方法与步骤。
2．掌握判别分析的回代误判率与交叉误判率的编程。
３．掌握贝叶斯判别的误判率的计算。
二、实验原理
1)在MATLＡＢ中,进行数据的判别分析命令为classｉｆy,其调用格式为:
＞> N１１＝0;Ｎ２２=0；
>> ｆor k＝１:n1
A=G1([1:ｋ-1,k+1:n1],:）；
N1＝ｌenｇth(A(:,１));
M１＝mean(A,1);s11=cｏv(A）;
S１=（(N1-1)*ｓ1１+（n2－1)*s2）/(N1+n2-ｋ);
ｐo1＝N1/（n-1)；ｐｏ2=ｎ２/(n-1)；
for j=1：2
if ｗ（ｊ)＝=mａx（w）
diｓｐ([＇待判样品属于第',nuｍ２sｔr（j），＇类污染'］);
enｄ
enｄ
end
待判样品属于第2类污染
待判样品属于第2类污染
待判样品属于第2类污染
待判样品属于第2类污染
3、计算回代误判率
>> n１1=0；n22＝0;
>> foｒ ｉ=1:n1
ｗ1(ｉ,１)=m1*inｖ（s)*Ｇ1（i,:)＇－1/２*m1*inv（s）*ｍ1'+lｏg(p1);
ｗ2(i，2）=m2*ｉｎv（s)*Ｇ2(i,：)＇-1/2*m2*ｉnv(s)*m２＇+ｌｏg(p2);
fｏr ｊ=1:2
if ｗ２(ｉ,ｊ）==maｘ(w２(i，:))＆j~＝2
n２2＝ｎ２2+1;
end
end
enｄ
>> ｐｏo=(ｎ１１＋n２2)/(n１+n2)
ｐoｏ =
0.1000
４、计算交叉误判率
w１(ｉ，２)＝m2＊inv(s)*G１(i,:)＇-１/２*ｍ2*iｎv(ｓ）*m2'+ｌog(p２);
ｆor j=1:2
if w1(ｉ,j)=＝mａx（w1(i,:))&j~＝1
ｎ11＝n11+１;
ｅnd
end
ｅnd
>> ｆor ｉ=1:n2
w2(ｉ,1）=ｍ１*iｎv(ｓ）*Ｇ2(i,:)'－1／2＊ｍ1*ｉnv(s)＊ｍ１'+loｇ(p1);
>> Ｍ＝（ｎ-k)＊ｌog（det（s))-（(n1-１）*log(ｄet(s1))+（n２-１)*log(ｄet(s２)）)；
>> T＝（1-d)＊Ｍ
T =
-4４.8２３7 + 0.９２88i
>> C＝ｃhｉ2iｎv(０.９5,f)
C =
３2.67０6
>> iｆ T<ｃhi2inv(０.９5,f)
disp('两组数据协方差相等');
elｓe
ｄisp（＇两组数据协方差不相等'）;
ｅｎd
两组数据协方差相等
2、构造判别函数，得出判别结果
＞> foｒ i=１:４
ｗ(1)＝ｍ１*ｉｎv（s)*x（i,:)'-1／2*ｍ1*iｎｖ(ｓ)*ｍ1'+log(ｐ1);
ｗ(２)=m2*inv(s）＊x(i,：）＇-1/2*m1*inv(ｓ)*ｍ２'＋ｌoｇ(p2）;
ｃlass=clａssifｙ(ｓample,training,gｒouｐ’tｙpe’)
将ｓａmplｅ数据的每一行指定到训练集tｒaｉnｉng的一个类中。Samｐle和ｔraｉｎiｎｇ必须具有相同的列数。grｏup向量包含从1到组数的正整数，它指明训练营集中的每一行属于哪一类。group和trａiniｎg必须具有相同的行数。’ｔype’是可选项，选’ｌineａｒ’表示总体为多元正态总体，选’ｑuａdｒａtｉc’与’mahalanobiｓ’。该函数返回class,它是一个与sampｌe具有相同行数的向量。Clasｓ的每一个元素指定sａmpｌe中对应元素的分类。通过计算sample和trainｉng中每一行的马氏距离,classify函数决定ｓａｍple中的每一行属于哪一个分类。
课程名称
数据分析方法
课程编号
实验地点
系统建模与仿真实验室SＬ11０
实验时间
校外指导教师
无
校内指导教师
实验名称
实验3 距离判别与贝叶斯判别分析
评阅人签字
成绩
实验数据与内容
我国山区某大型化工厂，在厂区及邻近地区挑选有代表性的１5个大气取样点，每日４次同时抽取大气样品，测定其中含有的6种气体的浓度,前后共４天，每个取样点每种气体实测16次,计算每个取样点每种气体的平均浓度,数据见表4－8。气体数据对应的污染地区分类见表4-8中最后一列。现有两个取自该地区的4个气体样本,气体指标见表4-8中后4行,试解决以下问题:
fｏr ｉ＝1：n1
w1(i,1）=M1*inv(Ｓ1)*Ｇ1(ｉ,:)'-1/2*M１*inv(S1）＊M1'＋loｇ(ｐo1);
w1(ｉ，2）=ｍ2*inv（S１)*G1(i，:)'-1/2＊m2*inｖ（Ｓ１)*m2＇＋loｇ(po2）;
fｏr j=1:2
if w1(ｉ.ｊ)==maｘ（Ｗ１(i,:))＆ｊ~=1
0.０740 0．0830 ０.105０ 0.1９0０ ０.０2０0 1.0000］
>> G１=A（[1:4 ７：８ １5],：）;
>＞ G2=A([５：6 ９:４ １5］,:);
>＞ ｎ1=size(G１,1）;
＞> n2=size(G２,1）;
>＞ n=ｎ１+n２;
>> k=2；
>> p＝6;
>> f＝p＊(ｐ＋1)＊(ｋ-1)/2;
从总体G1的训练样本开始,剔除其中一个样品,剩余的ｍ-1个样品与G2中的全部样品建立判别函数;
用建立的判别函数对剔除的样品进行判别;
重复以上步骤,直到Ｇ1中的全部样本依次被删除又进行判别，其误判的样品个数记为N１*；
对Ｇ2的样品重复以上步骤,直到G２中的全部样本依次被删除又进行判别，其误判的样品个数记为N2*。
0．０５2０ 0.０840 0.0２10 0.037０ ０.007１ ０．０220]
>> x=［0.052 0.084 0.0２１ 0．０３7 0.0071 ０.022
0.０41０ 0．0550 0.1１０0 ０.1１０0 0.0210 0.0073
0.030０ 0.１120 ０．07２0 ０.1６00 0．05６0 0.0210
于是交叉误判率估计为：
p^*=(Ｎ1*+Ｎ2*)/(m+ｎ)
5)贝叶斯判别的有效性可以通过平均误判率来确定。判别准则的误判率在一定程度上依赖于所考虑的各总体间的差异程度。各总体间差异越大,就越有可能建立有效的判别准则。如果各总体间差异很小,做判别分析的意义不大。
三、实验步骤ቤተ መጻሕፍቲ ባይዱ
输入数据，判别两类总体的协方差阵是否相等，用马氏距离判断判别污染类别，计算回代误判率与交叉误判率,贝叶斯判别污染分类。
１．判别两类总体的协方差矩阵是否相等，然后用马氏距离差别这4个未知气体样本的污染类别，并计算回代误判率与交叉误判率；若两类总体服从正态分布,第一类与第二类的先验概率分别为7/１5、8/1５,利用贝叶斯判别样本的污染分类。
2.先验概率为多少时,距离判别与贝时斯判别相同?调整先验概率对判别结果的影响是什么?
2)贝叶斯判别方法步骤
第1步,验证两个总体服从二元正态分布;第2步,检验两个总体的协方差矩阵相等；估计两个总体的先验概率p１、p2;利用MATLAＢ软件计算。
3）回代误判率
设G1,G2为两个总体,x１,x2…和y1,y２…是分别来自G１,Ｇ2的训练样本，以全体训练样本作为ｍ+n个新样品,逐个代入已建立的判别准则中判别其归属,这个过程称为回判。回判结果中若属于G１的样品被误判为属于G2的个数为N1个，属于G２的样品被误判为属于G1的个数为N2个，则误判估计为:
四、实验过程原始记录(数据、图表、计算等)
1、输入矩阵,计算协方差矩阵是否相等
>> Ａ＝［0.0560 ０.08４０ 0.0310 0.038０ ０．0081 0．0２２0
０.０４00 0.0５5０ 0.1000 0.1100 0.022０ 0.007３
……
0．0６90 ０.0870 ０.02７0 0.0500 0．089０ 0.021０