[R] ReL plot(cox.zph())

[Pu, Aihua]

Hello Terry,

Thanks very much for your prompt reply.
My data set is a big one containing 9233 unique patients records (one
patient per row).
Also I checked my x variable and none of them holds constant value. I
had no problem
in generating coxph and cox.zph object. The error message only came
after the plot statement. 

I would greatly appreciate if you can give me some advices on this. 

Thank you very much,

Aihua

PS:  My R code is shown below and the data set is in the attachment in
case you want to take a look.

revasFit <- coxph(Surv(myrevas2cath_bccr, revas_b365) ~
sex_num+age5564+age6574+age75
+ v3 +vlmd
+vless50+vmiss+ef3050+ef30+efmiss+dm+renal+HPD+HTN+PVD+CVA+CHF+copd
+liver+ cancer, data=mydata)
zph.revasFit <- cox.zph(revasFit,transform='log')
plot(zph.revasFit)

 
 
 




-----Original Message-----
From: Terry Therneau [mailto:therneau at mayo.edu]
Sent: Thursday, September 27, 2007 6:04 AM
To: Pu, Aihua
Cc: R-help at stat.math.ethz.ch
Subject: ReL plot(cox.zph())

You report an error message:

> plot(zph.revasFit[1])
Error in plot.window(xlim, ylim, log, asp, ...) :
        need finite 'ylim' values

I have never seen this error before, and I cannot guess what causes it.
You
need to provide more information, and likely a small data set that
produces
the problem.  Perhaps you have an x variable that is a constant?

        Terry Therneau
       


-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: mydata.txt
Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070927/00bfc322/attachment-0001.txt