[R] Obtaining partial output from a function that does not run to completion.

Eric rmailbox at justemail.net
Thu Dec 31 20:31:39 CET 2009


To heck with print(), use R's debugging capabilities instead:

trace ( minBMI, browser )

Be sure to ?trace and ?browser so you can figure out how to 
interactively debug.

Eric


On 12/31/09 6:29 AM, John Sorkin wrote:
> I have written a function that contains runs
> lm()
> vif()  and
> glm()
>
> When the glm() blows up with an error message, I don't get the output from either the lm() or vf() even thought neither lm() nor vif() have any problems . How can I force the function to print sequential results rather than wait for the entire function to complete before listing the functhion's output?
> Thanks,
> John
>
>
>
> minBMI<-function(SS,SimData)
> {
> SampleData<-sample(1:SS,size=SS,replace=TRUE)
> fitBMIEpiRevlm<-lm(AAMTCARE~BMIEpiRevAdjc+BMIEpiRevAdjcSq+SEX+jPHI+jMEDICAID+H_AGE+jMARSTAT+factor(jEDUCATION)+factor(jsmokercat)+factor(jrace)+log(INCOME_C+1),data=SimData[SampleData,],x=TRUE)
> print(summary(fitBMIEpiRevlm))
> print(vif(fitBMIEpiRevlm))
> fitBMIEpiRev<- glm(AAMTCARE~BMIEpiRevAdjc+BMIEpiRevAdjcSq+SEX+jPHI+jMEDICAID+H_AGE+jMARSTAT+factor(jEDUCATION)+factor(jsmokercat)+factor(jrace)+log(INCOME_C+1),data=SimData[SampleData,],family=Gamma(link="log"))
>
> print(summary(fitBMIEpiRev))
> }
> minBMI(SS,SimData)
>
>
> John David Sorkin M.D., Ph.D.
> Chief, Biostatistics and Informatics
> University of Maryland School of Medicine Division of Gerontology
> Baltimore VA Medical Center
> 10 North Greene Street
> GRECC (BT/18/GR)
> Baltimore, MD 21201-1524
> (Phone) 410-605-7119
> (Fax) 410-605-7913 (Please call phone number above prior to faxing)
>
> Confidentiality Statement:
> This email message, including any attachments, is for th...{{dropped:6}}
>
> ______________________________________________
> R-help at 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.
>




More information about the R-help mailing list