coxph
and
survReg
models.
frailty.gaussian(x, sparse=(ngroup>5), theta, df, eps=1e-05, method=c("reml", "aic", "df", "fixed"), ...)
"fixed"
method requires that
theta
be given explicitly (and is assumed
when
theta
is specified), the
"df"
method is used when degrees of freedom
is explicitly set.
The
"reml"
method uses McGilchrist's restricted maximum likelihood equation
to choose a value of the variance.
This is an iterative method that searches for the zero of an equation.
The
"aic"
method chooses the variance (and hence the degrees of freedom) based
on Akaike's information criteria.
trace=T
is the most
usual.
(The set of available arguments depends on the specific control function).
For the AIC method the
"caic=T"
argument may be used to choose the corrected
AIC criteria.
"coxph.penal"
.
When using frailty terms with a
survReg
model,
only
method="aic"
should be used.
The
method
argument can be set
in the call to
frailty
in the formula.
The
method
argument will be passed down to the specific
frailty function via the
...
argument.
fit1 <- coxph(Surv(time, status) ~ age + sex + frailty(inst, dist='gauss', df=4), data=lung, na.action=na.exclude) round(fit1$frail,3)