多元统计分析---判别分析

library(car) iris<-as.matrix(IRIS[,3:6]) scatterplotMatrix(~SL+SW+PL+PW|IRIS,data=IRIS,smooth=FALSE,reg.line=FALSE,ellipse=TRUE,levels=0.95,by.groups=TRUE,diagonal="none")

#线性判别分析 iris.lda<-lda(iris,factor(IRIS$IRIS)) #iris.lda是建立的判别准则 z1<-predict(iris.lda,dim=2)$x #基于2个判别准则预测得分 eqscplot(z1,type="n",xlab="first linear discriminant",ylab="second linear discriminant",main="IRIS data") #用了下面的text后才有图像 text(z1,labels=as.character(iris.Sp)) #给z1按照iris.Sp贴标签，按顺序排的前50个就是s,50-100是c，100-150是v z2<-predict(iris.lda,dim=2)$class #分类结果

#计算错判数 table(z2,IRIS$IRIS)