[R] Catch errors

[Stephen Tucker]

?try

or

?tryCatch
http://www.maths.lth.se/help/R/ExceptionHandlingInR/

for example...

tryCatch(lme(Y ~ X1*X2, random = ~1|subj, Model[i]),
         error=function(err) return(0))

(you can do something with 'err' or just return 0 as above)

--- Gang Chen <gangchen at mail.nih.gov> wrote:

> I run a linear mixed-effects model in a loop
> 
> for (i in 1:N) {
>     fit.lme <- lme(Y ~ X1*X2, random = ~1|subj, Model[i]);
> }
> 
> As the data in some iterations are all (or most) 0's, leading to the  
> following error message from lme:
> 
> Error in chol((value + t(value))/2) : the leading minor of order 1 is  
> not positive definite
> 
> What is a good way to catch the error without spilling on the screen  
> so that I can properly stuff the corresponding output with some  
> artificial value such as 0?
> 
> Thanks,
> Gang
> 
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>