tests/surv-Ex.R

library("mlt")
library("survival")
set.seed(29)

y <- rexp(10)
d <- rep(c(TRUE, FALSE), length = length(y))
### right
(s <- Surv(y, d, type = "right"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

### right censored, left truncated
(s <- Surv(y, y + 1, d))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

### left censored
(s <- Surv(y, d, type = "left"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

### interval
dd <- rep(0:3, length = length(y))
(s <- Surv(y, y + 1, dd, type = "interval"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

### interval2
(s <- Surv(y, y + 1, type = "interval2"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

(s <- Surv(y, ifelse(d, Inf, y + 1), type = "interval2"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

(s <- Surv(y, ifelse(d, NA, y + 1), type = "interval2"))
attr(s, "type")
cbind(s, R(s))
stopifnot(all.equal(s, as.Surv(R(s))))

### this is not the same as R(s) because of censoring!
as.numeric(R(s))
as.double(R(s))

R(list(runif(10), rnorm(10)))

Try the mlt package in your browser

Any scripts or data that you put into this service are public.

mlt documentation built on Dec. 1, 2023, 7:16 p.m.