[R] FW: Pec Function in R - syntax

K F Pearce k.f.pearce at newcastle.ac.uk
Tue Aug 25 13:13:27 CEST 2009


 
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



More information about the R-help mailing list