Thanks to Dieter Menne for suggesting cph in the Design package. I had been
looking at cph, but I can't seem to find an analog to the "expected" argument
in the survival package. Am I missing something? I am using R 2.6.1 and Windows
XP. Cheers,Reid> To: [EMAIL PROTECTED]> From: [EMAIL PROTECTED]> Date: Sat, 7
Jun 2008 08:15:48 +0000> Subject: Re: [R] expected risk from coxph (survival)>
> Reid Tingley <r_tingley <at> hotmail.com> writes:> > When I try to to obtain
the expected risk for a new dataset using coxph in the> survival package I get
an error.> > Using the example from ?coxph:> > # Example rewritten by DM;
please do not use HTML mail> library(survival)> test1 <- list(time= c(4,
3,1,1,2,2,3),> status=c(1,NA,1,0,1,1,0),> x= c(0, 2,1,1,1,0,0),> sex=
c(0,0,0,0,1,1,1))> cox<-coxph( Surv(time, status) ~ x + strata(sex), test1)
#stratified model> new<-list(time= c(5, 1,1,2,2,4,3),>
status=c(1,NA,1,0,0,1,1),> x= c(0, 2,1,1,1,0,0),> sex= c(0,0,0,0,1,1,1))> !
predict(cox,new,type="expected")> > > # }> # else if (type == "expected") {> #
if (missing(newdata)) > # pred <- y[, ncol(y)] - object$residuals> # else
stop("Method not yet finished")> Looks like this is "by design"; see the code
above.> > You might try to use cph and predict.Design from Frank Harrell's
Design package> instead.> > Dieter> >
______________________________________________> R-help@r-project.org mailing
list> https://stat.ethz.ch/mailman/listinfo/r-help> PLEASE do read the posting
guide http://www.R-project.org/posting-guide.html> and provide commented,
minimal, self-contained, reproducible code.
_________________________________________________________________
If you like crossword puzzles, then you'll love Flexicon, a game which
comb[[elided Hotmail spam]]
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
