Retail Pharmacies and Drug Diversion

All papers and comparison plots · Results JSON · Quarto source

2 of 2 registered specifications have estimates. Every registered specification appears below.

Janssen and Zhang (2023), AER. https://doi.org/10.1257/aer.20210357

Data and estimation

Outcomes are divided by their unweighted source-sample standard deviation (ddof=1), before exclusions; they are not centered. Treatment adoption is the first observed treated period. Models share an eligible sample, fixed effects, controls and weights. Initially treated units, missing −1 baselines and unsupported untreated-outcome predictions are removed. For switching treatments, untreated rows after first adoption are removed. Exact counts and settings are in the JSON.

Static TWFE uses the treatment indicator. Dynamic TWFE omits event time −1. BJS/did2s fits untreated outcomes, then averages imputation residuals. Sun–Abraham estimates cohort-by-event-time effects. The three dynamic estimates share post-period aggregation weights. Event-study tails are unbinned. These are common-sample reruns, not Peter’s original estimator-specific sample and tail-binning rules.

The homogeneity test restricts identified cohort deviations at positive event times to zero; leads and event time 0 remain unrestricted. It is a homoskedastic nested-model F-test. Estimate intervals are separately cluster-robust, with the did2s two-stage correction. Large panels use algebraically equivalent sparse regressions; backend and validation results are recorded in JSON.

Adoption timing

Treatment adoption timing

S15-P002: seeded sample of 100 of 82,873 eligible units. Purple untreated; yellow treated; white unobserved. Calendar periods follow the configured time grid. Each specification uses its own full eligible sample.

Specification Status F_pvalue Static TWFE Dynamic TWFE BJS / did2s Sun–Abraham
S15-P001 estimated 1 -0.2053 -0.1536 -0.1719 -0.1462
S15-P002 estimated 1 -0.1894 -0.07098 -0.1375 -0.06488

Specifications

Open a specification for its estimates, event studies and homogeneity test. Plots display −10 through 10; JSON contains the full paths. Large-panel intervals are calculated over the displayed window.

S15-P001 · RPD-T3-C4-ALL-OPIOIDS

Specification JSON · Configuration

RPD-T3-C4-ALL-OPIOIDS

Source rows: 5,055,761; fitted rows: 5,054,615; units: 82,873; clusters: 15,493; treated cohorts: 73. Source outcome SD: 541.10911.

method estimate se lower upper
Static TWFE -0.20531 0.031204 -0.26648 -0.14415
Dynamic TWFE -0.15362 0.034359 -0.22096 -0.086269
BJS / did2s -0.17192 0.027282 -0.2254 -0.11844
Sun–Abraham -0.14615 0.029035 -0.20306 -0.089238

Post-period weighting: observations. Homogeneity F(2749, 4965650) = 0.4001; p = 1.

TWFE, imputation, saturated event studies and cohort deviations
S15-P002 · RPD-T3-C9-OXYCONTIN

Specification JSON · Configuration

RPD-T3-C9-OXYCONTIN

Source rows: 5,055,761; fitted rows: 5,054,615; units: 82,873; clusters: 15,493; treated cohorts: 73. Source outcome SD: 75.909813.

method estimate se lower upper
Static TWFE -0.18942 0.054416 -0.29608 -0.082755
Dynamic TWFE -0.070978 0.032999 -0.13566 -0.006297
BJS / did2s -0.1375 0.036135 -0.20833 -0.066671
Sun–Abraham -0.064879 0.029923 -0.12353 -0.0062269

Post-period weighting: observations. Homogeneity F(2749, 4965650) = 0.3805; p = 1.

TWFE, imputation, saturated event studies and cohort deviations

Reproduce

Code, configurations and environments · Generic Quarto template · Editable Jot template. Source datasets remain local. Run python templates/full_archive/run_archive.py, then python templates/full_archive/build_collection.py.