1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
|
R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-pc-linux-gnu (32-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(survival)
Loading required package: splines
> #
> # A test of nesting. It makes sure the model.frame is built correctly
> #
> tfun <- function(fit, mydata) {
+ survfit(fit, newdata=mydata)
+ }
>
> myfit <- coxph(Surv(time, status) ~ age + factor(sex), lung)
>
> temp1 <- tfun(myfit, lung[1:5,])
> temp2 <- survfit(myfit, lung[1:5,])
> indx <- match('call', names(temp1)) #the call components won't match
>
> all.equal(unclass(temp1)[-indx], unclass(temp2)[-indx])
[1] TRUE
>
>
> proc.time()
user system elapsed
0.196 0.032 0.225
|
