penmodel_cmp {FamEvent} | R Documentation |
Fit a penetrance model for competing risks data
Description
Fits a competing risks model for family data with ascertainment correction and provides model parameter estimates.
Usage
penmodel_cmp(formula1, formula2, cluster = "famID", gvar = "mgene",
parms, cuts = NULL, data, design = "pop", base.dist = "Weibull",
frailty.dist = "none", agemin = NULL, robust = FALSE)
Arguments
formula1 |
A formula expression for event 1 as for other regression models. The response should be a survival object as returned by the |
formula2 |
A formula expression for event 2 as for other regression models. The response should be a survival object as returned by the |
cluster |
Name of cluster variable. Default is |
gvar |
Name of genetic variable. Default is |
parms |
list of Vectors of initial values for the parameters in each model including baseline parameters and regression coefficients and frailty parameters.
|
cuts |
Vector of cut points that define the intervals when |
data |
Data frame generated from |
design |
Study design of the family data. Possible choices are: |
base.dist |
Vector of two baseline hazard distributions to be fitted for competing events. Possible choices for each event are: |
frailty.dist |
Choice of frailty distribution to fit a shared frailty model for competing events. Possible choices are: |
agemin |
Minimum age of disease onset or minimum age. Default is |
robust |
Logical; if |
Details
The shared frailty comepting risks model is fitted to family data with specified baseline hazard distributions and frailty distribution
Event 1:
h1(t|X,Z) = h01(t - t0) Z1 exp(βs1 * xs + βg1 * xg),
Event 2:
h2(t|X,Z) = h02(t - t0) Z2 exp(βs2 * xs + βg2 * xg),
where h01(t) and h02(t) are the baseline hazard functions for event 1 and event 2, respectively, which can be specified by base.dist
. t0 is a minimum age of disease onset, Z1 and Z2 are frailties shared within families for each event and follow either a gamma, log-normal, correlateg gamma, or correlated log-normal distributions, which can be specified by frailty.dist
. xx and xg indicate male (1) or female (0) and carrier (1) or non-carrier (0) of a main gene of interest, respectively. Additional covariates can be added to formula1
for event 1 and formula2
for event 2 in the model.
Choice of frailty distributions for competing risk models
frailty.dist = "gamma"
shares the frailties within families generated from a gamma distribution independently for each competing event, where
Zj follows Gamma(kj, 1/kj).
frailty.dist = "lognormal"
shares the frailties within families generated from a log-normal distribution independently for each competing event, where
Zj follows log-normal distribution with mean 0 and variance 1/kj.
frailty.dist = "cgamma"
shares the frailties within families generated from a correlated gamma distribution to allow the frailties between two events to be correlated, where the correlated gamma frailties (Z1, Z2) are generated with three independent gamma frailties (Y0, Y1, Y2) as follows:
;
Z2 = k0/(k0 + k2) Y0 + Y2,
where Y0 from Gamma(k0, 1/k0), Y1 from Gamma(k1, 1/(k0 + k1)), Y2 from Gamma(k2, 1/(k0 + k2)).
frailty.dist = "clognormal"
shares the frailties within families generated from a correlated log-normal distribution where
log(Zj) follows a normal distribution with mean 0, variance 1/kj and correlation between two events k0.
depend
should specify the values of related frailty parameters: c(k1, k2)
with frailty.dist = "gamma"
or frailty.dist = "lognormal"
; c(k1, k2, k0)
for frailty.dist = "cgamma"
or frailty.dist = "clognormal"
.
More details about the competing risks model for family data arising from population-based study designs (design="pop", "pop+"
and their inference procedure based on the ascertainment-corrected likelihood approach can be found in Choi et al., 2021.
Note that the baseline parameters include lambda
and rho
, which represent the scale and shape parameters, respectively, and eta
, additional parameter to specify for "logBurr"
distribution. For the "lognormal"
baseline distribution, lambda
and rho
represent the location and scale parameters for the normally distributed logarithm, where lambda
can take any real values and rho
> 0. For the other baselinse distributions, lambda
> 0, rho
> 0, and eta
> 0. When a piecewise constant distribution is specified for the baseline hazards, base.dist="piecewise"
, baseparm
should specify the initial interval-constant values, one more than the cut points specified bycuts
.
Transformed baseline parameters are used for estimation; log transformation is applied to both scale and shape parameters () for
"Weibull"
, "loglogistic"
, "Gompertz"
and "gamma"
baselines, to () for
"logBurr"
and to the piecewise constant parameters for a piecewise
baseline hazard. For "lognormal"
baseline distribution, the log transformation is applied only to , not to
, which represents the location parameter for the normally distributed logarithm.
Calculations of penetrance estimates and their standard errors and 95% confidence intervals at given ages can be obtained by penetrance
function via Monte-Carlo simulations of the estimated penetrance model.
Value
Returns an object of class 'penmodel_cmp'
, including the following elements:
estimates |
Parameter estimates of transformed baseline parameters and regression coefficients. |
varcov |
Variance-covariance matrix of parameter estimates obtained from the inverse of Hessian matrix. |
varcov.robust |
Robust ‘sandwich’ variance-covariance matrix of parameter estimates when |
se |
Standard errors of parameter estimates obtained from the inverse of Hessian matrix. |
se.robust |
Robust ‘sandwich’ standard errors of parameter estimates when |
logLik |
Loglikelihood value for the fitted penetrance model. |
AIC |
Akaike information criterion (AIC) value of the model; AIC = 2*k - 2*logLik, where k is the number of parameters used in the model. |
Author(s)
Yun-Hee Choi
References
Choi, Y.-H., Jung, H., Buys, S., Daly, M., John, E.M., Hopper, J., Andrulis, I., Terry, M.B., Briollais, L. (2021) A Competing Risks Model with Binary Time Varying Covariates for Estimation of Breast Cancer Risks in BRCA1 Families, Statistical Methods in Medical Research 30 (9), 2165-2183. https://doi.org/10.1177/09622802211008945.
Choi, Y.-H., Briollais, L., He, W. and Kopciuk, K. (2021) FamEvent: An R Package for Generating and Modeling Time-to-Event Data in Family Designs, Journal of Statistical Software 97 (7), 1-30. doi:10.18637/jss.v097.i07
See Also
simfam_cmp
, penplot_cmp
, print.penmodel_cmp
,
summary.penmodel_cmp
, print.summary.penmodel_cmp
, plot.penmodel_cmp
Examples
# Competing risk family data simulated from population-based design
# using Weibull baseline hazards with gamma frailty distribution.
## Not run:
set.seed(4321)
fam1 <- simfam_cmp(N.fam = 200, design = "pop+", variation = "frailty",
base.dist = "Weibull", frailty.dist = "cgamma", depend=c(0.5, 1, 0.5),
allelefreq = 0.02, base.parms = list(c(0.01, 3), c(0.01, 3)),
vbeta = list(c(-1.13, 2.35), c(-1, 2)))
# Fitting shared gamma frailty Penetrance model for simulated competing risk data
fit1 <- penmodel_cmp(
formula1 = Surv(time, status==1) ~ gender + mgene,
formula2 = Surv(time, status==2) ~ gender + mgene,
cluster = "famID", gvar = "mgene", design = "pop+",
parms = list(c(0.01, 3, -1, 2), c(0.01, 3, -1, 2), c(0.5, 1)),
base.dist = "Weibull", frailty.dist = "gamma", data = fam1, robust = TRUE)
# Fitting shared correlated gamma frailty Penetrance model for simulated competing risk data
fit2 <- penmodel_cmp(
formula1 = Surv(time, status==1) ~ gender + mgene,
formula2 = Surv(time, status==2) ~ gender + mgene,
cluster = "famID", gvar = "mgene", design = "pop+",
parms = list(c(0.01, 3, -1, 2), c(0.01, 3, -1, 2), c(0.5, 1, 0.5)),
base.dist = "Weibull", frailty.dist = "cgamma", data = fam1, robust = TRUE)
# Summary of the model parameter estimates from the model fit
summary(fit1)
# Plot the lifetime penetrance curves with 95
# gender and mutation status groups along with their nonparametric penetrance curves
# based on data excluding probands.
plot(fit1, add.CIF = TRUE, conf.int = TRUE, MC = 100)
## End(Not run)