mets (version 1.3.4)

BinAugmentCifstrata: Augmentation for Binomial regression based on stratified NPMLE Cif (Aalen-Johansen)

Description

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)$$.

Usage

BinAugmentCifstrata(
  formula,
  data = data,
  cause = 1,
  cens.code = 0,
  km = TRUE,
  time = NULL,
  weights = NULL,
  offset = NULL,
  ...
)

Arguments

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.

Author

Thomas Scheike

Details

Standard errors computed under assumption of correct $$G_c(s)$$ model.

Examples

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)