BinAugmentCifstrata | R Documentation |
Computes the augmentation term for each individual as well as the sum
A = \int_0^t H(u,X) \frac{1}{S^*(u,s)} \frac{1}{G_c(u)} dM_c(u)
with
H(u,X) = F_1^*(t,s) - F_1^*(u,s)
using a KM for
G_c(t)
and a working model for cumulative baseline related to
F_1^*(t,s)
and
s
is strata,
S^*(t,s) = 1 - F_1^*(t,s) - F_2^*(t,s)
.
BinAugmentCifstrata( formula, data = data, cause = 1, cens.code = 0, km = TRUE, time = NULL, weights = NULL, offset = NULL, ... )
formula |
formula with 'Event', strata model for CIF given by strata, and strataC specifies censoring strata |
data |
data frame |
cause |
of interest |
cens.code |
code of censoring |
km |
to use Kaplan-Meier |
time |
of interest |
weights |
weights for estimating equations |
offset |
offsets for logistic regression |
... |
Additional arguments to binreg function. |
Standard errors computed under assumption of correct
G_c(s)
model.
Thomas Scheike
data(bmt) dcut(bmt,breaks=2) <- ~age out1<-BinAugmentCifstrata(Event(time,cause)~platelet+agecat.2+ strata(platelet,agecat.2),data=bmt,cause=1,time=40) summary(out1) out2<-BinAugmentCifstrata(Event(time,cause)~platelet+agecat.2+ strata(platelet,agecat.2)+strataC(platelet),data=bmt,cause=1,time=40) summary(out2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.