[R] Defining reference category for a cph model summary inside of a "for" loop

From: Wells, Brian <WELLSB_at_ccf.org>
Date: Fri, 28 Mar 2008 13:48:25 -0400

I have the following code.    

> f <- cph(formula = Surv(TimeToDeath, Dead == "Yes")
~1,data=single.dat, x=T, y=T, surv=T)

> for(i in c('A', 'B', 'C', 'D', 'E', 'F')){

> f <-update(f,as.formula(paste('Surv(TimeToDeath, Dead ==

> print(summary(f, paste(i,"=1st Quartile", sep='')))

There is no error message generated in R, but R ignores the reference category defined with paste in the summary function for the cph model.  

The output uses the "1st Quartile" as the reference category to calculate hazards for some of the variables defined by i, but not all of them.    

Any help would be greatly appreciated.  


Brian J. Wells, MD, MS

Research Associate

Quantitative Health Sciences

Cleveland Clinic    

P Please consider the environment before printing this e-mail

Cleveland Clinic is ranked one of the top hospitals in America by U.S. News & World Report (2007). Visit us online at http://www.clevelandclinic.org for a complete listing of our services, staff and locations.

Confidentiality Note: This message is intended for use only by the individual or entity to which it is addressed and may contain information that is privileged, confidential, and exempt from disclosure under applicable law. If the reader of this message is not the intended recipient or the employee or agent responsible for delivering the message to the intended recipient, you are hereby notified that any dissemination, distribution or copying of this communication is strictly prohibited. If you have received this communication in error, please contact the sender immediately and destroy the material in its entirety, whether electronic or hard copy. Thank you.

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.
Received on Fri 28 Mar 2008 - 17:52:31 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sat 29 Mar 2008 - 01:30:25 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive