Hi everyone, I ran into an error message when I use the cov.* (rob, mcd, mve) function in the package lqs. The data matrix is 761 rows by 661 columns. I get the following message: NA/NaN/Inf in foreign function call (arg 1) Can anyone provide some insight on why this occurs, and better yet, how to get around it? R 1.2.0 compiled from source on Linux Mandrake 7.1, running on Compaq SP750 with 2 P3 Xeons. TIA! Regards, Andy Liaw -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
On Wed, 24 Jan 2001, Liaw, Andy wrote:> Hi everyone, > > I ran into an error message when I use the cov.* (rob, mcd, mve) function in > the package lqs. The data matrix is 761 rows by 661 columns. I get the > following message: > > NA/NaN/Inf in foreign function call (arg 1) > > Can anyone provide some insight on why this occurs, and better yet, how to > get around it?Not without more information. But that is probably in z <- .C("mve_fitlots", as.double(x), as.integer(n), as.integer(p), in which case it probably means some column has IQR 0. Generally, you need to help yourself on these things. 1) traceback() tells you where the error is. 2) debug(cov.rob) and re-run the call, and step though a line at a time In particular, look at divisor. -- Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
> From: Prof Brian Ripley [mailto:ripley at stats.ox.ac.uk] > On Wed, 24 Jan 2001, Liaw, Andy wrote: > > I ran into an error message when I use the cov.* (rob, mcd, > mve) function in > > the package lqs. The data matrix is 761 rows by 661 > columns. I get the > > following message: > > > > NA/NaN/Inf in foreign function call (arg 1) > > > > Can anyone provide some insight on why this occurs, and > better yet, how to > > get around it? > > Not without more information. But that is probably in > > z <- .C("mve_fitlots", as.double(x), as.integer(n), > as.integer(p), > > in which case it probably means some column has IQR 0. > > Generally, you need to help yourself on these things. > > 1) traceback() tells you where the error is. > 2) debug(cov.rob) and re-run the call, and step though a line > at a time > In particular, look at divisor.Thanks a lot for the pointer. There are five columns that had IQR=0! BTW, traceback() didn''t give any more information than what I provided. It simply said the call bombed inside cov.rob. I''ll have to give debug() a try... Regards, Andy -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._