adjsurv | R Documentation |
Function for computing direct adjusted survival estimates
from a model fitted with the mexhaz
. It can be used to obtain
direct adjusted survival estimates for one or two populations. In the
latter case, survival difference estimates are also
computed. Corresponding variance estimates are based on the Delta Method
(based on the assumption of multivariate normality of the model
parameter estimates). When the model includes a random effect, two types
of predictions can be made: (i) marginal predictions (obtained by
integration over the random effect distribution) or (ii) conditional
predictions either for a particular cluster (using the corresponding
shrinkage estimate) or for the value 0 of the random effect.
adjsurv(object, time.pts, data, data.0 = NULL, weights = NULL, clust.name = NULL, marginal = TRUE, level = 0.95)
object |
an object of class |
time.pts |
a vector of numerical values representing the time points at which predictions are requested. Time values greater than the maximum follow-up time on which the model estimation was based are discarded. |
data |
a |
data.0 |
an optional |
clust.name |
name of the variable in |
weights |
optional argument specifying the weights to be
associated with each row of |
marginal |
logical value controlling the type of predictions
returned by the function when the model includes a random
intercept. When |
level |
a number in (0,1) specifying the level of confidence for
computing the confidence intervals of the hazard and the
survival. By default, |
An object of class resMexhaz
that can be used by the function
plot.resMexhaz
to produce graphics
of the direct adjusted survival curve. It contains the following
elements:
results |
a |
type |
type of results returned by the function. The value is
used by |
multiobs |
value used by
|
ci.method |
method used to compute confidence limits. Currently set
to |
level |
level of confidence used to compute confidence limits. |
Hadrien Charvat, Aurelien Belot
Charvat H, Remontet L, Bossard N, Roche L, Dejardin O, Rachet B, Launoy G, Belot A; CENSUR Working Survival Group. A multilevel excess hazard model to estimate net survival on hierarchical data allowing for non-linear and non-proportional effects of covariates. Stat Med 2016;35:3066-3084 (doi: 10.1002/sim.6881)
Booth JG, Hobert JP. Standard errors of prediction in generalized linear mixed models. J Am Stat Assoc 1998;93:262-272 (doi: 10.2307/2669622).
plot.resMexhaz
, lines.resMexhaz
data(simdatn1) ## Fit of a fixed-effect hazard model, with the baseline hazard ## described by a linear B-spline with two knots at 1 and 5 year and with ## effects of age (agecr), deprivation index (depindex) and sex (IsexH) Mod_bs1_2 <- mexhaz(formula=Surv(time=timesurv, event=vstat)~agecr+depindex+IsexH, data=simdatn1, base="exp.bs", degree=1, knots=c(1,5), verbose=0) ## Direct adjusted survival for the simdatn1 population DAS_Modbs1_2 <- adjsurv(Mod_bs1_2, time.pts=seq(1,10), data=simdatn1)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.