[R] sbrier (Brier score) and coxph

From: Stephen Henderson <s.henderson_at_ucl.ac.uk>
Date: Thu 30 Jun 2005 - 06:51:11 EST


Hello
I've decided to try and distill an earlier rather ill focused question to try and elicit a response. Any help is greatly appreciated. Why does mod.cox not work with sbrier whilst mod.km does? Can I make it work?

> data(DLBCL)
> DLBCL.surv<-Surv(DLBCL$time,DLBCL$cens)
>
> mod.km<-survfit(DLBCL.surv)
> mod.cox<-survfit(coxph(DLBCL.surv~IPI, data=DLBCL))
>
>
> sbrier(DLBCL.surv, mod.km)

integrated Brier score

             0.2076454
attr(,"time")
[1] 1.3 129.9
> sbrier(DLBCL.surv, mod.cox)

Error in switch(ptype, survfit = { : switch: EXPR must return a length 1 vector

Thanks in advance
Stephen Henderson



This email and any files transmitted with it are confidentia...{{dropped}}

R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jun 30 07:06:30 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:33:07 EST