Erik Iverson
2010-Sep-13  18:30 UTC
[R] using survexp and ratetable with coxph object that includes a factor term
Hello,
I'm attempting to use the ratetable argument to
survexp in the survival package.  I use
the example from the ?survexp help page below,
and then slightly modify it to produce an error.
library(survival)
data(pbc)
#fit a model without any factors
pfit1 <- coxph(Surv(time, status > 0) ~ trt + log(bili) +
          log(protime) + age + platelet, data=pbc)
#this works as expected
survexp(~ trt, ratetable = pfit1, data = pbc)
#now I fit a similar model, simply adding the 'sex' term
#note sex is now the only factor in the model
pfit2 <- coxph(Surv(time, status > 0) ~ trt + log(bili) +
          log(protime) + age + platelet + sex, data = pbc)
#however, I now get an error with survexp
survexp(~ trt, ratetable = pfit2, data = pbc)
Error in ratetable(trt = trt, bili = bili, protime = protime, age = age,
: object 'sexf' not found
Does anyone have any suggestions?
Thank you,
Erik Iverson
sessionInfo()
R version 2.11.1 (2010-05-31)
x86_64-unknown-linux-gnu
locale:
  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
  [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8
  [5] LC_MONETARY=C              LC_MESSAGES=en_US.UTF-8
  [7] LC_PAPER=en_US.UTF-8       LC_NAME=C
  [9] LC_ADDRESS=C               LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] grDevices datasets  splines   graphics  utils     stats     methods
[8] base
other attached packages:
[1] Hmisc_3.7-0     survival_2.35-8
loaded via a namespace (and not attached):
[1] cluster_1.12.3 grid_2.11.1    lattice_0.18-8 tools_2.11.1
