idmTPreg (version 1.1)

TPreg: Regression Modeling of Transition Probabilities in a Progressive Illness-Death Model

Description

Fits a semi-parametric regression model to estimate the effects on transition probabilities in a (possibly non-Markov) progressive illness-death model for a sequence of time.

Usage

TPreg(formula, data, link, s = 0, t = NULL, R = 199, by = NULL, trans, ncores = NULL)

Arguments

formula

an object of class formula which specifies the covariates. For example formula = ~ age + sex.

data

a data.frame of iddata class or a data.frame in which other than covariates five variables; id, Zt, delta1, Tt, delta are included.

link

a link function for binomial family which are logit, probit and cauchit, (corresponding to logistic, normal and Cauchy CDFs respectively).

s

the current time for the transition probabilities; default is zero which reports the occupation probabilities.

t

the Future time for the transition probabilities; default is NULL which is the largest uncensored sojourn time in the initial state.

R

the number of bootstrap replicates. Default is 199.

by

number: increment of the sequence from time s to time t.The default is NULL which is \(\lfloor (\max({ Zt}) - \min({Zt}))/q_{0.01}({ Zt }) \rfloor\) , where \(q_{0.01}(.)\) is the sample quantile corresponding to 0.01 probability and \(\lfloor x\rfloor\) gives the largest integer less than or equal to \(x\). A binomial regression at every byth time between s and t is performed. by=1 reports all binomial regression results for each jump time, corresponding to the specified transition(s), between s and t. By increasing by we skip some times. In order to save the time, for a relatively large dataset a relatively big by is recommended.

trans

the possible transition(s) for a progressive illness-death model. For trans argument there are five options: "11", "12", "13", "23", and "all".

ncores

the number of cores to use for parallel execution. Default is the number of CPU cores on the current host.

Value

TPreg returns an object of class TPreg. An object of class TPreg is a list containing at least the following components:

co

the list of:

• ‘transition’ the transition,

• ‘time’ the jump times,

• ‘coefficients’ the estimated effects ,

• ‘SD’ standard errors,

• ‘LWL’ lower confidence limits,

• ‘UPL’ upper confidence limits,

• ‘p.value’ p-values.

call

the matched call.

transition

the transition, this is equal to the transition in co unless for trans="all"

s

the current time for the transition probability.

t

the future time for the transition probability.

n.misobs

the number of missing observations.

In addition, trans="all" will have four lists: co11, co12, co13, and co23 instead of co and will give the information for all possible transitions in the progressive illness death model

References

Azarang, L. Scheike, TH. and de Una-Alvarez, J. (2017) Direct modeling of regression effects for transition probabilities in the progressive illness-death model, Statistics in Medicine 36, \(1964-1976\).

See Also

print.TPreg, summary.TPreg, and plot.TPreg as generic functions.

Examples

# NOT RUN {
data(colonTPreg)

co11 <- TPreg( ~ Age + Nodes + treatment, colonTPreg, link = "logit", s = 50, R = 19, t = 200,
trans = "11", ncores = 1)
co11
summary(co11)
# }