Kaplan-Meier with robust standard errors Robust variance is default variance and obtained from the predict call
km(formula, data = data, ...)Arguments
Examples
library(mets)
data(sTRACE)
sTRACE$cluster <- sample(1:100,500,replace=TRUE)
out1 <- km(Surv(time,status==9)~strata(vf,chf),data=sTRACE)
out2 <- km(Surv(time,status==9)~strata(vf,chf)+cluster(cluster),data=sTRACE)
summary(out1,times=1:3)
#> $pred
#> [,1] [,2] [,3]
#> [1,] 0.9216911 0.8606809 0.8171022
#> [2,] 0.7180838 0.6043700 0.5032911
#> [3,] 0.6930406 0.6930406 0.6930406
#> [4,] 0.5320165 0.4044111 0.4044111
#>
#> $se.pred
#> [,1] [,2] [,3]
#> [1,] 0.01768107 0.02279610 0.02545049
#> [2,] 0.02911535 0.03163815 0.03234100
#> [3,] 0.16275118 0.16275118 0.16275118
#> [4,] 0.09970269 0.09739478 0.09739478
#>
#> $lower
#> [,1] [,2] [,3]
#> [1,] 0.8876802 0.8171413 0.7687123
#> [2,] 0.6632273 0.5454355 0.4437332
#> [3,] 0.4373867 0.4373867 0.4373867
#> [4,] 0.3684728 0.2522477 0.2522477
#>
#> $upper
#> [,1] [,2] [,3]
#> [1,] 0.9570051 0.9065405 0.8685383
#> [2,] 0.7774775 0.6696724 0.5708430
#> [3,] 1.0000000 1.0000000 1.0000000
#> [4,] 0.7681477 0.6483640 0.6483640
#>
#> $times
#> [1] 1 2 3
#>
#> attr(,"class")
#> [1] "summarypredictrecreg"
summary(out2,times=1:3)
#> $pred
#> [,1] [,2] [,3]
#> [1,] 0.9216911 0.8606809 0.8171022
#> [2,] 0.7180838 0.6043700 0.5032911
#> [3,] 0.6930406 0.6930406 0.6930406
#> [4,] 0.5320165 0.4044111 0.4044111
#>
#> $se.pred
#> [,1] [,2] [,3]
#> [1,] 0.01570548 0.02197520 0.02391423
#> [2,] 0.02841423 0.03050283 0.03315955
#> [3,] 0.16275118 0.16275118 0.16275118
#> [4,] 0.09520495 0.09294011 0.09294011
#>
#> $lower
#> [,1] [,2] [,3]
#> [1,] 0.8914173 0.8186703 0.7715502
#> [2,] 0.6644977 0.5474474 0.4423209
#> [3,] 0.4373867 0.4373867 0.4373867
#> [4,] 0.3746292 0.2577528 0.2577528
#>
#> $upper
#> [,1] [,2] [,3]
#> [1,] 0.9529931 0.9048474 0.8653437
#> [2,] 0.7759911 0.6672113 0.5726656
#> [3,] 1.0000000 1.0000000 1.0000000
#> [4,] 0.7555245 0.6345163 0.6345163
#>
#> $times
#> [1] 1 2 3
#>
#> attr(,"class")
#> [1] "summarypredictrecreg"
par(mfrow=c(1,2))
plot(out1,se=TRUE)
plot(out2,se=TRUE)
