frail.fit {eha} | R Documentation |
This function is called from mlreg
. Can be called directly
by a user who knows what (s)he is doing!
frail.fit(X, Y, rs, strats, offset, init, max.survs, frailty, control)
X |
Design matrix |
Y |
Survival object |
rs |
Risk set pointers |
strats |
Stratum variable |
offset |
Offset variable |
init |
Initial regression parameter values |
max.survs |
Sampling of risk sets? |
frailty |
The frailty grouping variable |
control |
Values for the optimization |
Frailty is analyzed in a marginal framework, with numerical integration. This is very slow and somewhat unstable.
See the code!
Not very much tested. Can fail, and is very slow.
Göran Broström
fam.size <- rep(20, 5) n <- sum(fam.size) fam <- rep(1:length(fam.size), fam.size) beta <- 1 sigma <- 1 frail <- rep(rnorm(length(fam.size), 0, sigma), fam.size) x <- runif(n) exit <- rexp(n, exp(x * beta + frail)) event <- rep(1, n) enter <- numeric(n) mlreg(Surv(enter, exit, event) ~ x, frailty = frail)