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)

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.