Augmentation for Binomial regression based on stratified NPMLE Cif (Aalen-Johansen)
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.
Augmentation term only computed for standard FG model, since strata is used to specify working models for CIF's.
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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.