某企生产的产品,其造型、性能和价位及所属级别数据如下表所示:
某企业产品的造型、性能、价位、级别等指标
题目来自《多元统计分析-基于R》课后习题
下面分别用Fisher判别法和Bayes判别法进行判别分析。
先进性Fisher判别,R程序如下:
#Fisher判别
ex5.3<-read.csv("ex5.3.csv",header = T)
ex5.3
library(MASS) #加载MASS程序包,以便使用其中lda函数
ld=lda(G~x1+x2+x3,data=ex5.3)
ld
运行结果如下:
程序输出了lda()所用的函数,各组的先验概率,各组的均值向量,第一及第二判别函数的系数
,两个判别2函数对判别的贡献大小。
下面我们用predict()函数对原始数据进行回判,将lda()函数的分类与原数据分类进行对比,考察误差的大小,R程序如下: