[R] new object

Jacques VESLOT jacques.veslot at good.ibl.fr
Tue Jul 11 14:28:11 CEST 2006


 > names(summary(fit))
  [1] "surv"     "time"     "n.risk"   "n.event"  "conf.int" "std.err"
  [7] "lower"    "upper"    "strata"   "call"
 > summary(fit)$n.risk
  [1] 11 10  8  7  5  4  2 12 10  8  6  5  4  3  2  1

-------------------------------------------------------------------
Jacques VESLOT

CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex

Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31

http://www-good.ibl.fr
-------------------------------------------------------------------


Mauricio Cardeal a écrit :
> Hi !
> 
> Please, how can I extract n.event and n.risk as a new object from 
> example below ?  Thanks in advance.
> Mauricio
> 
> require(survival)
> fit <- survfit(Surv(time, status) ~ x, data=aml)
> 
> summary(fit)
> 
> Call: survfit(formula = Surv(time, status) ~ x, data = aml)
> 
>                 x=Maintained
>  time n.risk n.event survival std.err lower 95% CI upper 95% CI
>     9     11       1    0.909  0.0867       0.7541        1.000
>    13     10       1    0.818  0.1163       0.6192        1.000
>    18      8       1    0.716  0.1397       0.4884        1.000
>    23      7       1    0.614  0.1526       0.3769        0.999
>    31      5       1    0.491  0.1642       0.2549        0.946
>    34      4       1    0.368  0.1627       0.1549        0.875
>    48      2       1    0.184  0.1535       0.0359        0.944
> 
>                 x=Nonmaintained
>  time n.risk n.event survival std.err lower 95% CI upper 95% CI
>     5     12       2   0.8333  0.1076       0.6470        1.000
>     8     10       2   0.6667  0.1361       0.4468        0.995
>    12      8       1   0.5833  0.1423       0.3616        0.941
>    23      6       1   0.4861  0.1481       0.2675        0.883
>    27      5       1   0.3889  0.1470       0.1854        0.816
>    30      4       1   0.2917  0.1387       0.1148        0.741
>    33      3       1   0.1944  0.1219       0.0569        0.664
>    43      2       1   0.0972  0.0919       0.0153        0.620
>    45      1       1   0.0000      NA           NA           NA
> 
> 	[[alternative HTML version deleted]]
> 
> ______________________________________________
> R-help at 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
>



More information about the R-help mailing list