[R] How to flag those iterations which yield a warning?

[Ravi Varadhan]

Nope - that does not work.  The value of last.warning is not reset after the initial NULL.

Ravi.

____________________________________________________________________

Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine and Gerontology
School of Medicine
Johns Hopkins University

Ph. (410) 502-2619
email: rvaradhan at jhmi.edu


----- Original Message -----
From: Ravi Varadhan <rvaradhan at jhmi.edu>
Date: Friday, February 18, 2011 11:35 am
Subject: Re: [R] How to flag those iterations which yield a warning?
To: Ravi Varadhan <rvaradhan at jhmi.edu>
Cc: r-help at r-project.org


> May be I could do:
>  
>  for (i in 1:nsim) {
>  
>  last.warning <- NULL
>  
>  # do model fitting
>  
>  if(!is.null(last.warning)) # discard simulation result
>  
>  }
>  
>  I think this might work.  Any other ideas?
>  
>  Ravi.
>  
>  ____________________________________________________________________
>  
>  Ravi Varadhan, Ph.D.
>  Assistant Professor,
>  Division of Geriatric Medicine and Gerontology
>  School of Medicine
>  Johns Hopkins University
>  
>  Ph. (410) 502-2619
>  email: rvaradhan at jhmi.edu
>  
>  
>  ----- Original Message -----
>  From: Ravi Varadhan <rvaradhan at jhmi.edu>
>  Date: Friday, February 18, 2011 11:29 am
>  Subject: [R] How to flag those iterations which yield a warning?
>  To: r-help at r-project.org
>  
>  
>  > Hi,
>  >  
>  >  I am running a simulation study with the survival::coxph.  Some of 
> 
>  > the simulations result in problematic fits due to flat partial 
>  > likelihood. So, you get the warning message:
>  >  
>  >  Warning message:
>  >  In fitter(X, Y, strats, offset, init, control, weights = weights,  
> 
>  > ... :
>  >    Loglik converged before variable  2 ; beta may be infinite. 
>  >  
>  >  How can I keep track of the simulations which yield any kind of 
>  > warning message, so that I can appropriately handle them in the 
>  > post-processing of simulation results?
>  >  
>  >  I know that `try' can be used to capture errors, but I would like 
> to 
>  > capture and detect warnings.
>  >  
>  >  Any pointers would be appreciated.
>  >  
>  >  Thanks,
>  >  Ravi.
>  >  ____________________________________________________________________
>  >  
>  >  Ravi Varadhan, Ph.D.
>  >  Assistant Professor,
>  >  Division of Geriatric Medicine and Gerontology
>  >  School of Medicine
>  >  Johns Hopkins University
>  >  
>  >  Ph. (410) 502-2619
>  >  email: rvaradhan at jhmi.edu
>  >  
>  >  ______________________________________________
>  >  R-help at r-project.org mailing list
>  >  
>  >  PLEASE do read the posting guide 
>  >  and provide commented, minimal, self-contained, reproducible code. 
>