先读取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判别结果

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

分析上述数据,对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),对它们的品种进行判别:

Fisher判别示例:鸢尾花(iris)数据(R)-LMLPHP

04-30 12:00