This PR adds a script that assesses vaccine effectiveness using a sequential trials approach.
model_stcox.R fits the models. In the interests of time, it uses the full one-row-per-person-per-day person-time dataset (created by the data_stset.R script) to determine who is eligible for inclusion and to get "baseline" covariates for each trial (may have done differently if writing from scratch).
report_stcox.R creates a plot of the estimated effects, and is very similar to this script in the comparative effectiveness repo.
I haven't yet written anything than combines estimates from across different brands/outcomes into one plot.
This PR adds a script that assesses vaccine effectiveness using a sequential trials approach.
model_stcox.Rfits the models. In the interests of time, it uses the full one-row-per-person-per-day person-time dataset (created by thedata_stset.Rscript) to determine who is eligible for inclusion and to get "baseline" covariates for each trial (may have done differently if writing from scratch).report_stcox.Rcreates a plot of the estimated effects, and is very similar to this script in the comparative effectiveness repo.I haven't yet written anything than combines estimates from across different brands/outcomes into one plot.