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)