Hi all, and thanks for the responses.
As per Bill's suggestion, I tried:
x <- try(flowClust(tFrame, K=30, B=1000, varNames=c('CD4', 'CD8','KI67',
'CD45RO', 'CD28', 'CD57', 'CCR5', 'CD19','CD27', 'CCR7','CD127')),
silent=TRUE)
and it caught the error.
Indeed,
x <- try(res30 <- flowClust(tFrame, K=30, B=1000, varNames=c('CD4',
'CD8','KI67', 'CD45RO', 'CD28', 'CD57', 'CCR5', 'CD19','CD27',
'CCR7','CD127')), silent=TRUE)
caught the error. (What I really meant by "res30 = flowClust(...)" was
"res30 <- flowClust(...)". I've fallen into the probably bad habit of
using '=' in place of '<-' ) It had completely slipped my mind that '='
has a dual meaning, and that using it there might result in that part of
my expression being evaluated as the passing of an argument.
Anyway, the point is that it works, and the error is caught when I use '<-'.
Thanks again for all your help (and an interesting peek into the nuts
and bolts of R).
-Kieran
William Dunlap wrote:
Note that Kieren's example labelled the first
argument to try() with an improper label res30=,
not expr= (or is that a mailer turning something
into '30='?). If it really is an improper argument
tag then this could be showing a buglet in reporting
on wrongly named arguments:
> invisible(rm(x,y))
> x<-try(silent=TRUE, badTag=stop("Oops"))
Error in try(silent = TRUE, badTag = stop("Oops")) : Oops
> x
Error: object "x" not found
> y<-try(silent=TRUE, expr=stop("Oops"))
> y
[1] "Error in try(silent = TRUE, expr = stop(\"Oops\")) : Oops\n"
attr(,"class")
[1] "try-error"
In the first example I would expect an error message like
unused argument(s) (badTag = stop("Oops"))
but it is appropriate that try() would abort if it
is called in a bad way. Perhaps it is trying to make that
error message and that triggered the evaluation of the argument,
as in
> grep(mypattern=stop("Oops"), "wxyz")
Error in grep(mypattern = stop("Oops"), "wxyz") : Oops
where one might expect an error message regarding the wrongly
named argument, as in:
> grep(mypattern="x", "wxyz")
Error in grep(mypattern = "x", "wxyz") :
unused argument(s) (mypattern = "x")
Bill Dunlap
TIBCO Software Inc - Spotfire Division
wdunlap tibco.com
-Original Message-
From: r-devel-boun...@r-project.org
[mailto:r-devel-boun...@r-project.org] On Behalf Of Dirk Eddelbuettel
Sent: Wednesday, April 15, 2009 7:14 PM
To: Kieran O'Neill
Cc: r-devel@r-project.org
Subject: Re: [Rd] How can I catch errors thrown from c via
the Rcpperror() function?
Kieran,
On 15 April 2009 at 18:03, Kieran O'Neill wrote:
| I am using the flowClust package from BioConductor, which
is largely
| implemented in c. For some of my data, the package
occasionally (and
| quite stochastically) encounters a particular condition
which halts its
| operation. At this point, it calls the error() function
defined by Rcpp,
| and halts.
|
| What I would like to be able to do is to catch the error
thrown, and
| retry the operation a few times before giving up.
|
| However, when I wrap the call to flowClust in try() or
tryCatch(), the
| error seems to completely bypass them:
|
| Examples:
|
| 1. This is a trivial example just to test the try() function, and
| correctly assigns the error to the variable x:
|
| > x <- try(stop(simpleError('blah')))
| Error : blah
| > x
| [1] "Error : blah\n"
| attr(,"class")
| [1] "try-error"
|
| 2. This is an example using flowClust (using real data, set up to
| guarantee that the error is thrown):
|
| > x <- try(res30 = flowClust(tFrame, K=30, B=1000,
varNames=c('CD4',
| 'CD8','KI67', 'CD45RO', 'CD28', 'CD57', 'CCR5', 'CD19',
'CD27', 'CCR7',
| 'CD127')))
| Error in flowClust(tFrame, K = 30, B = 1000, varNames =
c("CD4", "CD8", :
|
| The covariance matrix is near singular!
| Try running the program with a different initial
configuration or less
| clusters
| > x
| Error: object "x" not found
|
|
| The c code throwing the error is as follows (from flowClust.c):
|
| if(status!=0)
|{
|error("\n The covariance matrix is near singular! \n
Try running
| the program with a different initial configuration or less clusters
| \n"); }
|
|
| I looked up
