species<-factor(c(rep('s',50),rep('v',50)))
irsv<-rbind(iris[,,1],iris[,,2])
irsvdf<-data.frame(species,irsv)
irsvdf
irsv.discrim<-discrim(species,data=irsvdf)
irsv.discrim
ircoef<-coef(irsv.discrim)
ircoef1<-ircoef$linear.coefficients[,1]
ircoef2<-ircoef$linear.coefficients[,2]
ira<-ircoef1-ircoef2
这是以鸢尾花数据为例的两总体Fisher判别

请求高手指教！！！！

主要是到了discrim就运行不出来了