DA1<-function
+ 　   (TrnX1, TrnX2, TstX = NULL, var.equal = FALSE){if (is.null(TstX) == TRUE) TstX<-rbind(TrnX1,TrnX2)
+ 　   if (is.vector(TstX) == TRUE)  TstX<-t(as.matrix(TstX)) else if (is.matrix(TstX) != TRUE)
+ 　      TstX<-as.matrix(TstX)
+ 　   if (is.matrix(TrnX1) != TRUE) TrnX1<-as.matrix(TrnX1)
+ 　   if (is.matrix(TrnX2) != TRUE) TrnX2<-as.matrix(TrnX2); nx<-nrow(TstX)
+ 　   blong<-matrix(rep(0, nx), nrow=1, byrow=TRUE, dimnames=list("blong", 1:nx))
+ 　   mu1<-colMeans(TrnX1); mu2<-colMeans(TrnX2)
+ 　   if (var.equal == TRUE  || var.equal == T){S<-var(rbind(TrnX1,TrnX2))
+ 　      w<-mahalanobis(TstX, mu2, S)-mahalanobis(TstX, mu1, S)} else{S1<-var(TrnX1); S2<-var(TrnX2)
+ 　      w<-mahalanobis(TstX, mu2, S2)-mahalanobis(TstX, mu1, S1)}
+ 　   for (i in 1:nx){if (w>0) blong<-1 else blong<-2}; blong}
> DA1(classX1, classX2, var.equal=TRUE) #样本协方差相同
      1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47
blong 1 1 1 1 1 1 1 1 1  1  1  1  1  1  1  1  1  1  1  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2
>  
>  
> DA1(classX1, classX2)#样本协方差不同
Error in solve.default(cov, ...) :
  system is computationally singular: reciprocal condition number = 7.66864e-17

求大神解释原因啊- -