Usage
idm(formula01, formula02, formula12, data, maxiter=200, eps=c(5,5,3),
n.knots=c(7,7,7), irec=0,kappa0 = c(1000000,500000,20000), igraph=1,
hazard="Weib", print.iter=FALSE, subset=NULL, na.action=na.omit)
Arguments
formula01
a formula object for transition '0 --> 1' (from the initial state to the transient state), with the response on the left of a $\texttildelow$
operator, and the terms on the right. The response must be a
survival object or Hist object as returned by the 'S
formula02
a formula object for transition '0 --> 2' (from the initial state to the absorbant state), with the response on the left of a $\texttildelow$
operator, and the terms on the right. The response must be a
survival object or Hist object as returned by the 'S
formula12
a formula object for transition '1 --> 2' (from the transient state to the absorbant state). The response on the left of the $\texttildelow$
operator is not required. If missing formula12, variables of transition '1--->2' are the same as the ones of trans
data
a data frame in which to interpret the variables named in formula01, formula02 and formula12.
maxiter
maximum number of iterations. The default is 200.
eps
a vector of length 3 for the convergence criteria (criterion for parameters, criterion for likelihood, criterion for second derivatives). The default is 'c(5,5,3)' and corresponds to criteria equals to $10^{-5}$, $10^{-5}$ and $10^{-3}$.
n.knots
a vector of length 3 for the number of knots for the splines to use to approximate the intensities for transitions 0 --> 1, 0 --> 2 and 1 --> 2 respectively. Argument for the penalized likelihood approach.
The default is c(7,7,7).
irec
binary variable equals to 1 when search (by approximated cross validation) of the smoothing parameters kappa and 0 otherwise. Argument for the penalized likelihood approach. The default is 0.
kappa0
a vector of length 3. If irec=0, smoothing parameters for the transition 0 --> 1, 0 --> 2 and 1 --> 2. If irec=1, initial values of the smoothing parameters for the cross validation search. Argument for the penalized likelihood approach.
igraph
a binary variable equals to 1 when you want to output graphs of the transition intensities, 0 otherwise. The default is 1.
hazard
type of estimation method: "Splines" for a penalized likelihood approach with approximation of the transition intensities by M-splines, "Weib" for a parametric approach with a Weibull distribution on the transition intensities. Default is "Weib".
print.iter
boolean parameter. Equals to TRUE to print the likelihood during the iteration process, FALSE otherwise. Default is FALSE. This option is not running on Windows.
subset
expression indicating the subset of the rows of data to be used in the fit. All observaation are included by default.
na.action
how NAs are treated. The default is first, any na.action attribute of data, second a na.action setting of options, and third 'na.fail' if that is unset. The 'factory-fresh' default is na.omit. Another possible value is NULL.