用R进行数据判别时出现了问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

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

求大神解释原因啊- -

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：数据判别

兄弟，你的代码可以注释下，再往上传吗？看不懂的。

Solution is attached. Little comments for your coding:

1. Always comment your code, as mentioned up;
2.  is.vector() this type is already logical value with TRUE or FALSE, no need to specify, i.e.,
is.vector() ==TRUE ;
3. check input values for left to right, i.e., TrnX1, TrnX2, TstX