[R] new object

[ronggui]

> sum.out<-summary(fit)
> names(sum.out)
 [1] "surv"     "time"     "n.risk"   "n.event"  "conf.int" "std.err"
 [7] "lower"    "upper"    "strata"   "call"
> sum.out$n.risk
 [1] 11 10  8  7  5  4  2 12 10  8  6  5  4  3  2  1
> sum.out$n.event
 [1] 1 1 1 1 1 1 1 2 2 1 1 1 1 1 1 1
> cbind(sum.out$n.risk,sum.out$n.event,sum.out$strata)
      [,1] [,2] [,3]
 [1,]   11    1    1
 [2,]   10    1    1
 [3,]    8    1    1
 [4,]    7    1    1
 [5,]    5    1    1
 [6,]    4    1    1
 [7,]    2    1    1
 [8,]   12    2    2
 [9,]   10    2    2
[10,]    8    1    2
[11,]    6    1    2
[12,]    5    1    2
[13,]    4    1    2
[14,]    3    1    2
[15,]    2    1    2
[16,]    1    1    2


2006/7/11, Mauricio Cardeal <mcardeal at ufba.br>:
> 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
>


-- 
»ÆÈÙ¹ó
Department of Sociology
Fudan University