Description Usage Arguments Details Value Author(s) References See Also Examples
Tools for fitting proportional hazards models to clustered bivariate recurrent events data. Nested frailties are modeled by their best linear unbiased predictors under an auxiliary Poisson model.
1 2 | ## S3 method for class 'formula'
bivrec(formula, data = parent.frame(), ...)
|
formula |
See details |
data |
See details |
... |
See details |
Supported parameters are:
1 2 3 4 5 |
a formula object, similar to coxph
. The response to the
left of the ~
should be a survival object generated by Surv2
. The
right side must contain an id(x)
term, where x
is a variable that takes a
unique value for each subject within a cluster. It may also contain a cluster(y)
term if y
is the variable that indicates cluster membership, and a strata(z)
term.
a data.frame
with columns corresponding to the terms
in the formula.
either an integer, a vector of integers, or a value between 0 and 1,
to determine the level of discretization. If it is an integer, K1
gives the
number of breakpoints in the baseline hazard for the first process.
If it is a vector of integers, K1
should have length equal to the
number of strata, and each value gives the number of breakpoints in the
baseline hazard for each stratum. If it is a number between 0 and 1, it gives
the ratio of the number of breakpoints relative to the maximum possible.
Defaults to 10.
analogous to K1
, for the second process.
a vector of strings giving the names of variables or interactions that should be excluded from the model for the first process.
analogous to excludevars1
for the second process.
an integer from 0 to 3 that determines the quantity of output
printed to the screen. Setting verbose=0
is completely silent.
logical, describing the at-risk function for the model. If
FALSE
, patients are assumed to be constantly at risk for both processes,
if TRUE
, they are only at-risk for one process at a time. Defaults to
FALSE
.
a string, determining the method used to estimate the dispersion
parameters. Possible values are "pearson", "marginal", "ohlsson"
, defaults
to "pearson"
.
a string describing the degree-of-freedom correction proposed
by Ma. It only applies when dispest="pearson"
.
Possible values are "single", "double", "none"
, defaults to "none"
a boolean determining whether standard errors should be computed.
Defaults to TRUE
.
logical determining whether to use the full sensitivity matrix
in computing standard errors, or only the covariate portion. Using fullS=TRUE
leads to more accurate results but increases computer time. Defaults to TRUE
.
a vector of strings, listing any dispersion parameters that should
be set to 0. Can contain any subset of "clust1", "clust2", "subj1", "subj2", "cov"
,
all else is ignored. Defaults to NULL
.
logical, determines whether the baseline hazard should be smoothed
at each iteration. Defaults to FALSE
.
an integer giving the maximum number of iterations permitted.
double, determines the value of the convergence criterion required at termination.
a list of initial values in the format used by the code internally.
An object of class bivrec
with the following components:
call |
the original call to the model-fitting function |
regression |
a list containing results from the regression fit in the last iteration. It has components
|
frailty |
a list containing results from the frailty estimation in the last iteration. It has components
|
dispersion |
a list containing results from the dispersion parameter estimation in the last iteration. It has components
|
hazard |
a list describing the baseline hazard. It has components
|
summaries |
a list of summary matrices. It has components
|
Emmanuel Sharef ess28@cornell.edu
E. Sharef and R. Strawderman. “A nested frailty model for clustered bivariate recurrent events”, in preparation.
1 2 3 4 5 6 7 |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.