- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
5-9 设在某地区抽取了
14 块岩石标本, 其中 7 块含矿, 7 块不含矿。对每块岩石
测定了 Cu,Ag ,Bi 三种化学成分的含量,得到的数据如表 1。
表 1
岩石化学成分的含量数据
类型
序号
Cu
Ag
Bi
类型
序号
Cu
Ag
Bi
1
2.58
0.90
0.95
8
2.25
1.98
1.06
2
2.90
1.23
1.00
9
2.16
1.80
1.06
含
3
3.55
1.15
1.00
不
10
2.33
1.74
1.10
4
2.35
1.15
0.79
含
11
1.96
1.48
1.04
矿
矿
5
3.54
1.85
0.79
12
1.94
1.40
1.00
6
2.70
2.23
1.30
13
3.00
1.30
1.00
7
2.70
1.70
0.48
14
2.78
1.70
1.48
(1)假定两类样本服从正态分布, 使用广义平方距离判别法进行判别归类 (先验概率取为相等,并假定两类样本的协方差阵相等) ;
(2)今得一块标本,并测得其 Cu,Ag ,Bi 的含量分别为 2.95,2.15 和 1.54,试判断该标本是含矿还是不含矿?
问题求解
使用广义平方距离判别法对样本进行判别归类用 SAS 软件中的 DISCRIM 过程进行判别归类。 SAS 程序及结果如下。
data
d59;
input
group x1-x3@@;
cards
;
1 2.58
0.9
0.95
1 2.9 1.23
1
1 3.55
1.15
1
1 2.35
1.15
0.79
1 3.54
1.85
0.79
1 2.7 2.23
1.3
1 2.7 1.7 0.48
2 2.25
1.98
1.06
2 2.16
1.8
1.06
2 2.33
1.74
1.1
2 1.96
1.48
1.04
2 1.94
1.4
1
231.31
2 2.78
1.7
1.48
;
proc
print
data =d59;
run
;
proc
discrim
data =d59 pool =yes distance
list ;
class
group;
var
x1-x3;
run
;
由输出结果可知,两总体间的广义平方距离为 D2=3.19774。还可知两个三
元总体均值相等的检验结果: D=3.19774,F=3.10891,p=0.07560.10,故在显著
性水平
=0.10 时量总体的均值向量有显著差异,即认为讨论这两个三元总体的
判别问题是有意义的。
线性判别函数为:
Y1
31.1105
13.7895Cu
8.2120Ag
11.3311Bi,
Y2
28.7375
10.3139Cu
8.9904Ag
16.8578Bi.
判别结果为含矿的 6 号样本错判为不含矿;不含矿的 13 号样本错判为含矿。
对给定样本判别归类
将 Cu,Ag ,Bi 的含量数值 2.95、2.15、1.54 分别代入线性判别函数得:
Y1 44.6742,Y2 46.97888。
贝叶斯判别的解
D *
D1* ,
, D k
*
为
Dt*
X | Yt ( X )
Yj ( X ), j
t, j
1,
,k (t
1,
, k )
,
由于 Y
1
44.67422 Y2
46.97888,因此待判的样品判为不含矿。
5-10 已知某研究对象分为三类,每个样品考察 4 项指标,各类的观测样品数分
别为 7,4,6;类外还有 3
个待判样品(所有观测数据见表
2)。假定样本均来自正
态总体。
X1
表 2
判别分类的数据
样品号
X2
X3
X4
类别号
1
6.0
-11.5
19.0
90.0
1
2
-11.0
-18.5
25.0
-36.0
3
3
90.2
-17.0
17.0
3.0
2
4
-4.0
-15.0
13.0
54.0
1
5
0.0
-14.0
20.0
35.0
2
6
0.5
-11.5
19.0
37.0
3
7
-10.0
-19.0
21.0
-42.0
3
8
0.0
-23.0
5.0
-35.0
1
9
20.0
-22.0
8.0
-20.0
3
10
-100.0
-21.4
7.0
-15.0
1
11
-100.0
-21.5
15.0
-40.0
2
12
13.0
-17.2
18.0
2.0
2
13
-5.0
-18.5
15.0
18.0
1
14
10.0
-18.0
14.0
50.0
1
15
-8.0
-14.0
16.0
56.0
1
16
0.6
-13.0
26.0
21.0
3
17
-40.0
-20.0
22.0
文档评论(0)