先读取iris数据,再用程序包MASS(记得要在使用MASS前下载好该程序包)中的线性函数lda()作判别分析:
data(iris) #读入数据
iris #展示数据
attach(iris) #用变量名绑定对应数据
library(MASS) #加载MASS程序包
ld=lda(Species~Sepal.Length+Sepal.Width+Petal.Len gth+Petal.Width) #进行Fisher判别
ld #输出Fisher判别结果
分析上述数据,对150个原始数据的回代判别中,只有3个错误,误判率为2%。
采集到3朵新的鸢尾花,其测量值分别为(5.1, 3.5, 1.5, 0.25)、(5.9, 2.8, 4.3, 1.3)和(6.6, 2.9, 5.6, 2.1),对它们的品种进行判别: