Hi folks, I know, this is a fairly common question and I am really disappointed that I could not find a solution. I am trying to calculate Mahanalobis distances in a data frame, where I have several hundreds groups and several hundreds of variables.
Whatever I do, however I subset it I get the "system is computationally singular: reciprocal condition number" error. I know what it means and I know what should be the problem, but there is no way this is a singular matrix. I have uploaded the input file to my ftp: http://mkk.szie.hu/dep/talt/lv/CentInpDuplNoHeader.txt It is a tab delimited txt file with no headers. I tried the StatMatch Mahanalobis function and also this function: mahal_dist <-function (data, nclass, nvariable) { dist <- matrix(0, nclass, nclass) n=0 w <- cov(data) print(w) for(i in 1:nclass) { for(c in 1:nclass){ diffl <- vector(length = nvariable) for(l in 1:nvariable){ diffl[l]=abs(data[i,l]-data[c,l]) } ### matrixes print(diffl) dist[i,c]= (t(diffl))%*%(solve(w))%*%(diffl) } n=n+1 print(n) } return(dist) sqrt_dist <- sqrt(dist) print(sqrt_dist) } I have a deadline for this project (not a homework:)), and I could always use this codes, so I thought I will be able to quit the calculations short, but now I am just lost. I would really appreciate any help. Thanks for any help -- View this message in context: http://r.789695.n4.nabble.com/system-is-computationally-singular-reciprocal-condition-number-tp4647472.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
