frail.fit {eha}R Documentation

Fits a frailty proportional hazards model

Description

This function is called from mlreg. Can be called directly by a user who knows what (s)he is doing!

Usage

frail.fit(X, Y, rs, strats, offset, init, max.survs, frailty, control)

Arguments

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

Details

Frailty is analyzed in a marginal framework, with numerical integration. This is very slow and somewhat unstable.

Value

See the code!

Warning

Not very much tested. Can fail, and is very slow.

Author(s)

Göran Broström

See Also

mlreg

Examples

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)

[Package eha version 0.95-4 Index]