[R] Coxph with factors

From: Kylie-Anne Richards <kar_at_itga.com.au>
Date: Fri 15 Jul 2005 - 23:39:37 EST


My apologies for lack of clarity.

Sample data:

      time.sec po vo done f.pom
      55 0 4402 1 -0.25
      378 -0.5 50000 1 -0.25
      51 0 1000 1 -0.25
      43 0 71581 1 -0.25

etc

>        DATA<-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo)
>        final<-coxph(Surv(time.sec,done)~f.pom+vo+po,data=DATA)
                      coef exp(coef) se(coef)       z       p
f.pom-2.5  1.30e+00      3.65 1.66e-01   7.822 5.2e-15
f.pom-2    1.11e+00      3.02 1.68e-01   6.604 4.0e-11
f.pom-1.5  1.37e+00      3.95 1.59e-01   8.656 0.0e+00
f.pom-1    1.76e+00      5.79 1.56e-01  11.286 0.0e+00
f.pom-0.5  2.10e+00      8.14 1.52e-01  13.757 0.0e+00
f.pom0     4.96e+00    142.10 1.95e-01  25.433 0.0e+00
f.pom0.5   4.69e+00    108.71 1.53e-01  30.573 0.0e+00
f.pom1     4.66e+00    105.50 1.58e-01  29.575 0.0e+00
f.pom1.5   4.66e+00    105.62 1.60e-01  29.051 0.0e+00
f.pom2     4.61e+00    100.77 1.60e-01  28.792 0.0e+00
vo           -6.93e-05      1.00 6.57e-06 -10.548 0.0e+00
po           2.65e-03      1.00 1.29e-03   2.047 4.1e-02

>final.surv<-survfit((final), individual=T,conf.type="log-log")
>print(final.surv)

  time n.risk n.event survival std.err lower 95% CI upper 95% CI

     0  13922    3914 8.10e-01 3.82e-03     8.02e-01      0.81712
     1   9726      91 7.98e-01 3.97e-03     7.90e-01      0.80547
     2   9457      46 7.92e-01 4.04e-03     7.84e-01      0.79949
     3   9297      36 7.87e-01 4.10e-03     7.79e-01      0.79477
etc

FIRST Q: The default uses the mean of 'vo' and mean of 'po', but what is it using for the factors?? Is it the sum of the coef of the factors divided by the number of factors??

SECOND Q: For a model with covariates I would normally specify: final.surv<-survfit((final), data.frame(po=0,vo=0,pom=0,individual=T)) to get the baseline survival prob.; what would I specify for a model with a factor, i.e., 'f.pom' ??

E.g. one of the many variations I have tried: > final.surv<-survfit((final),
> data.frame(po=0,vo=0,f.pom=0,individual=T)) Error in model.frame(formula, rownames, variables, varnames, extras, extranames, :

        variable lengths differ
In addition: Warning message:
'newdata' had 1 rows but variable(s) found have 13922 rows Execution halted

I appreciate the time anyone can spare to help me with my problem, and apologies for it's simplistic nature.

Many thanks,

Kylie-Anne

Adaikalavan Ramasamy ramasamy at cancer.org.uk Fri Jul 15 09:58:19 CEST 2005

  a.. Previous message: [R] Coxph with factors
  b.. Next message: [R] Pearson dispersion statistic
  c.. Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]

--------------------------------------------------------------------------------

Yes, and please show us a reproducible example or small section of thedata as well as the error output. On Thu, 2005-07-14 at 15:36 +0000, Dieter Menne wrote:> Kylie-Anne Richards <kar <at> itga.com.au> writes:> > > > > I am fitting a coxph model with factors. I am running into problems when> > using 'survfit'. I am unsure how R is treating the factors when I fit, say:> > > DATA<-data.frame(time.sec,done,f.pom=factor(f.pom),po,vo)> > > final<-coxph(Surv(time.sec,done)~f.pom*vo+po,data=DATA)> > > final.surv<-survfit((final), individual=T,conf.type="log-log")> > > print(final.surv)> ....> > Better chances to get a reply when you tell us what problems you are running > into.> > Dieter Menne



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 Fri Jul 15 23:48:37 2005

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