[R] Running cox models

Geoff Russell geoffrey.russell at gmail.com
Sun Sep 3 12:50:38 CEST 2006


Hi,

I'm reading van Belle et al "Biostatistics" and trying to run a cox test using
a dataset from:

http://faculty.washington.edu/~heagerty/Books/Biostatistics/chapter16.html

(Primary Biliary Cirrhosis data link at top of the page),

I'm using the following code:

--------------- start of code
library(survival)
liver <- scan("liver2.txt",list(age=0,albumin=0,alkphos=0,ascites=0,bili=0,
        cholest=0,edema=0,edmadj=0,hepmeg=0,obstime=0,platelet=0,protime=0,
        sex=0,sgot=0,spiders=0,stage=0,status=0,treatmnt=0,
        triglyc=0,urinecu=0))
fit<-coxph(Surv(obstime,status)~bili+edmadj+albumin+protime+age,data=liver)
summary(fit)
----------------- End of code

but the answer is rather different from that in the book (p.688 - for
anyone with the book).

The book refers to EDTRT, but the dataset has EDMADJ and EDMEMA, also
the book talks about 312 patients and the dataset has 418 lines.

Has anybody else used this dataset?

Cheers,
Geoff Russell



More information about the R-help mailing list