Hello everyone - sorry, I forgot to add the syntax for my program in the last
message:
My syntax is as follows:
>library(Hmisc)
> library(survival)
> library(pec)
> AL<-read.table("PredErr2.txt",header=TRUE)
>
Models<-list("Kaplan.Meier"=survfit(Surv(OVS,dead)~1,data=AL),"EBMT+SN
>
Pgroup"=coxph(Surv(OVS,dead)~GroupMod,method="efron",data=AL),"EBMTgro
> up"=coxph(Surv(OVS,dead)~GroupGrat,
>
method="efron",data=AL),"EBMToriginal"=coxph(Surv(OVS,dead)~Gratwohl,
> method="efron",data=AL),"EBMT+SNPoriginal">
coxph(Surv(OVS,dead)~Gratwohl+PIL1rAnyC+PIL4AnyT,
> method="efron",data=AL))
> PredError.632plus<-pec(object=Models,formula=Surv(OVS,dead)~1,data=AL,
>
exact=TRUE,cens.model="marginal",replan="boot632plus",B=100,verbose=TR
> UE)
Thankyou so much.
Kim
Dr Kim Pearce CStat
Industrial Statistics Research Unit (ISRU) School of Mathematics and Statistics
Herschel Building University of Newcastle Newcastle upon Tyne United Kingdom
NE1 7RU
Tel. 0044 (0)191 222 6244 (direct)
Fax. 0044 (0)191 222 8020
Email: K.F.Pearce at ncl.ac.uk