View source: R/binomial.regression.R

binregATE | R Documentation |

Under the standard causal assumptions we can estimate the average treatment effect E(Y(1) - Y(0)). We need Consistency, ignorability ( Y(1), Y(0) indep A given X), and positivity.

binregATE( formula, data, cause = 1, time = NULL, beta = NULL, treat.model = ~+1, cens.model = ~+1, offset = NULL, weights = NULL, cens.weights = NULL, se = TRUE, kaplan.meier = TRUE, cens.code = 0, no.opt = FALSE, method = "nr", augmentation = NULL, ... )

`formula` |
formula with outcome (see |

`data` |
data frame |

`cause` |
cause of interest |

`time` |
time of interest |

`beta` |
starting values |

`treat.model` |
logistic treatment model given covariates |

`cens.model` |
only stratified cox model without covariates |

`offset` |
offsets for partial likelihood |

`weights` |
for score equations |

`cens.weights` |
censoring weights |

`se` |
to compute se's with IPCW adjustment, otherwise assumes that IPCW weights are known |

`kaplan.meier` |
uses Kaplan-Meier for IPCW in contrast to exp(-Baseline) |

`cens.code` |
gives censoring code |

`no.opt` |
to not optimize |

`method` |
for optimization |

`augmentation` |
to augment binomial regression |

`...` |
Additional arguments to lower level funtions |

The first covariate in the specification of the competing risks regression model must be the treatment effect that is a factor. If the factor has more than two levels then it uses the mlogit for propensity score modelling. If there are no censorings this is performing ordinary logistic regression modelling.

This is then model using a logistic regresssion using the standard binary double robust estimating equations that are then IPCW censoring adjusted using binomial regression.

Rather than binomial regression we also consider a IPCW weighted version of standard logistic regression logitIPCWATE.

The original version of the program with only binary treatment binregATEbin take binary-numeric as input for the treatment, and also computes the ATT and ATC, average treatment effect on the treated (ATT), E(Y(1) - Y(0) | A=1), and non-treated, respectively. Experimental version.

Thomas Scheike

data(bmt) dfactor(bmt) <- ~. brs <- binregATE(Event(time,cause)~tcell.f+platelet+age,bmt,time=50,cause=1, treat.model=tcell.f~platelet+age) summary(brs) brsi <- binregATE(Event(time,cause)~tcell.f+tcell.f*platelet+tcell.f*age,bmt,time=50,cause=1, treat.model=tcell.f~platelet+age) summary(brsi)

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.