diff-diff 3.8.0__tar.gz → 3.9.1__tar.gz

{diff_diff-3.8.0 → diff_diff-3.9.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: diff-diff
-Version: 3.8.0
+Version: 3.9.1
 Classifier: Development Status :: 5 - Production/Stable
 Classifier: Intended Audience :: Science/Research
 Classifier: Operating System :: OS Independent
@@ -20,9 +20,9 @@ Requires-Dist: scipy>=1.10
 Requires-Dist: pytest>=7.0 ; extra == 'dev'
 Requires-Dist: pytest-xdist>=3.0 ; extra == 'dev'
 Requires-Dist: pytest-cov>=4.0 ; extra == 'dev'
-Requires-Dist: black==26.3.1 ; extra == 'dev'
-Requires-Dist: ruff==0.15.13 ; extra == 'dev'
-Requires-Dist: mypy==2.1.0 ; extra == 'dev'
+Requires-Dist: black==26.5.1 ; extra == 'dev'
+Requires-Dist: ruff==0.16.2 ; extra == 'dev'
+Requires-Dist: mypy==2.3.0 ; extra == 'dev'
 Requires-Dist: maturin>=1.4,<2.0 ; extra == 'dev'
 Requires-Dist: matplotlib>=3.5 ; extra == 'dev'
 Requires-Dist: nbmake>=1.5 ; extra == 'dev'
@@ -33,6 +33,7 @@ Requires-Dist: sphinxext-opengraph>=0.9 ; extra == 'docs'
 Requires-Dist: sphinx-sitemap>=2.5 ; extra == 'docs'
 Requires-Dist: nbsphinx>=0.9 ; extra == 'docs'
 Requires-Dist: myst-parser>=2.0 ; extra == 'docs'
+Requires-Dist: sphinx-design>=0.6.1 ; extra == 'docs'
 Requires-Dist: matplotlib>=3.5 ; extra == 'docs'
 Requires-Dist: ipython>=8.0 ; extra == 'docs'
 Requires-Dist: plotly>=5.0 ; extra == 'plotly'
@@ -95,7 +96,7 @@ data = pd.DataFrame({
 })
 
 did = DifferenceInDifferences()
-results = did.fit(data, outcome='outcome', treatment='treated', time='post')
+results = did.fit(data, outcome='outcome', treatment='treated', post='post')
 print(results)              # DiDResults(ATT=3.0000, SE=1.7321, p=0.1583)
 results.print_summary()     # full statsmodels-style table
 ```
@@ -143,9 +144,9 @@ For rigorous DiD analysis, follow these 8 steps. Skipping diagnostic steps produ
 2. **State identification assumptions** - which parallel trends variant (unconditional, conditional, PT-GT-Nev, PT-GT-NYT), no-anticipation, overlap.
 3. **Test parallel trends** - simple 2x2: `check_parallel_trends()`, `equivalence_test_trends()`; staggered: inspect CS event-study pre-period coefficients (generic PT tests are invalid for staggered designs). Insignificant pre-trends do NOT prove PT holds.
 4. **Choose estimator** - staggered adoption -> CS/SA/BJS (NOT plain TWFE); few treated units -> SDiD; factor confounding -> TROP; simple 2x2 -> DiD. Run `BaconDecomposition` to diagnose TWFE bias.
-5. **Estimate** - `estimator.fit(data, ...)`. Always print the cluster count first and choose inference method based on the result (cluster-robust if >= 50 clusters, wild bootstrap if fewer).
-6. **Sensitivity analysis** - `compute_honest_did(results)` for bounds under PT violations (MultiPeriodDiD, CS, or dCDH), `run_all_placebo_tests()` for 2x2 falsification, specification comparisons for staggered designs.
-7. **Heterogeneity** - CS: `aggregate='group'`/`'event_study'`; SA: `results.event_study_effects` / `to_dataframe(level='cohort')`; subgroup re-estimation.
+5. **Estimate** - `estimator.fit(data, ...)`. Always print the cluster count first and choose inference method based on the result (cluster-robust if >= 50 clusters, wild bootstrap if fewer - for DifferenceInDifferences pass `cluster=`; TwoWayFixedEffects auto-clusters at unit level).
+6. **Sensitivity analysis** - `compute_honest_did(results)` for bounds under PT violations (MultiPeriodDiD, CS, or dCDH natively; the TwoWayFixedEffects `event_study=True` surface and a StackedDiD `results.aggregate('event_study')` container also admit - Stacked needs `kappa_pre >= 2`), `run_all_placebo_tests()` for 2x2 falsification, specification comparisons for staggered designs.
+7. **Heterogeneity** - CS: `results.aggregate('group')`/`'event_study'` (post-fit, no refit); SA: `results.event_study_effects` / `to_dataframe(level='cohort')`; Stacked: `results.aggregate('event_study')`/`'simple'` post-fit views (surface always computed since 3.9); EDiD: `results.aggregate(...)` post-fit from retained EIFs (3.9); ImputationDiD/TwoStageDiD: `results.aggregate(...)` post-fit from panel-backed kits (3.9); ContinuousDiD: `results.aggregate('dose'/'simple'/'event_study')` post-fit (3.9; dose/simple are views, event_study recomputes); subgroup re-estimation.
 8. **Robustness** - compare 2-3 estimators (CS vs SA vs BJS), report with and without covariates (shows whether conditioning drives identification), present pre-trends and sensitivity bounds.
 
 Full guide: `diff_diff.get_llm_guide("practitioner")`.
@@ -154,7 +155,7 @@ Full guide: `diff_diff.get_llm_guide("practitioner")`.
 
 - [DifferenceInDifferences](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - basic 2x2 DiD with robust/cluster-robust SEs, wild bootstrap, formula interface, and fixed effects
 - [TwoWayFixedEffects](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - panel data DiD with unit and time fixed effects via within-transformation or dummies
-- [MultiPeriodDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - event study design with period-specific treatment effects for dynamic analysis
+- [MultiPeriodDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - event study design with period-specific treatment effects for dynamic analysis (deprecated 3.9 - use TwoWayFixedEffects `event_study=True`)
 - [CallawaySantAnna](https://diff-diff.readthedocs.io/en/stable/api/staggered.html) - Callaway & Sant'Anna (2021) group-time ATT estimator for staggered adoption
 - [ChaisemartinDHaultfoeuille](https://diff-diff.readthedocs.io/en/stable/api/chaisemartin_dhaultfoeuille.html) - de Chaisemartin & D'Haultfœuille (2020/2022) for **reversible (non-absorbing) treatments** with multi-horizon event study, normalized effects, cost-benefit delta, sup-t bands, and dynamic placebos. The most general option for treatments that switch on AND off (see also `LPDiD`/`TROP` `non_absorbing`). Alias `DCDH`.
 - [SunAbraham](https://diff-diff.readthedocs.io/en/stable/api/staggered.html) - Sun & Abraham (2021) interaction-weighted estimator for heterogeneity-robust event studies
@@ -163,26 +164,29 @@ Full guide: `diff_diff.get_llm_guide("practitioner")`.
 - [SpilloverDiD](https://diff-diff.readthedocs.io/en/stable/api/spillover.html) - Butts (2021) ring-indicator spillover-aware DiD identifying direct effect on treated + per-ring spillover on near-control units; handles non-staggered and staggered timing; supports survey-design variance under `survey_design=` for HC1 / CR1 (Wave E.1 Binder TSL) and Conley (Wave E.2 panel-aware stratified-Conley sandwich on per-period PSU totals; extended in Wave E.2 follow-up to `conley_lag_cutoff > 0` via panel-block composition with within-PSU serial Bartlett HAC — `lag>0` requires an effective PSU via explicit `survey_design.psu` or injected `cluster=<col>`); `SurveyDesign.subpopulation()` preserves full-design `n_psu` / `df_survey` via zero-padded scores (Wave E.3, R `svyrecvar(subset())` form)
 - [SyntheticDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - Synthetic DiD combining standard DiD and synthetic control for few treated units
 - [SyntheticControl](https://diff-diff.readthedocs.io/en/stable/api/synthetic_control.html) - Abadie, Diamond & Hainmueller (2010) classic synthetic control for a single treated unit (donor-weight counterfactual, nested/cv/inverse-variance/custom V; in-space placebo permutation inference via `in_space_placebo()`, plus ADH-2015 `leave_one_out()` + `in_time_placebo()` robustness, Firpo-Possebom (2018) test-inversion confidence sets, and Chernozhukov-Wüthrich-Zhu (2021) conformal inference)
-- [TripleDifference](https://diff-diff.readthedocs.io/en/stable/api/triple_diff.html) - triple difference (DDD) estimator for designs requiring two criteria for treatment eligibility
+- [TripleDifference](https://diff-diff.readthedocs.io/en/stable/api/triple_diff.html) - triple difference (DDD) estimator for designs requiring two criteria for treatment eligibility; serves both the 2x2x2 and the staggered-adoption design from one signature (`fit(..., first_treat=)` selects the staggered engine)
 - [ContinuousDiD](https://diff-diff.readthedocs.io/en/stable/api/continuous_did.html) - Callaway, Goodman-Bacon & Sant'Anna (2024) continuous treatment DiD with dose-response curves
 - [HeterogeneousAdoptionDiD](https://diff-diff.readthedocs.io/en/stable/api/had.html) - de Chaisemartin, Ciccia, D'Haultfœuille & Knau (2026) for designs where **no unit remains untreated**; local-linear estimator at the dose support boundary returning Weighted Average Slope (WAS) on Design 1' (`d̲ = 0` / QUG) or `WAS_{d̲}` on Design 1 (`d̲ > 0`, continuous-near-d̲ or mass-point), with a multi-period event-study extension (last-treatment cohort, pointwise CIs). **Panel-only** in this release - repeated cross-sections rejected by the validator. Alias `HAD`.
-- [RegressionDiscontinuity](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2014) sharp, fuzzy, AND covariate-adjusted regression discontinuity with robust bias-corrected inference and rdrobust-parity bandwidth selection (all 10 selectors, mass-point handling; fuzzy via `treatment_col=` with a first-stage block and weak-identification warning; covariates via `covariates=` - CCFT 2019, same estimand, covariate-aware bandwidths). Canonical `att` is the bias-corrected estimate with a coherent robust CI (rdrobust's printed headline is `att_conventional`). Alias `RDD`.
+- [RegressionDiscontinuity](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2014) sharp, fuzzy, AND covariate-adjusted regression discontinuity with robust bias-corrected inference and rdrobust-parity bandwidth selection (all 10 selectors, mass-point handling; fuzzy via `takeup=` with a first-stage block and weak-identification warning; covariates via `covariates=` - CCFT 2019, same estimand, covariate-aware bandwidths). Canonical `att` is the bias-corrected estimate with a coherent robust CI (rdrobust's printed headline is `att_conventional`). Alias `RDD`.
 - [StackedDiD](https://diff-diff.readthedocs.io/en/stable/api/stacked_did.html) - Wing, Freedman & Hollingsworth (2024) stacked DiD with Q-weights and sub-experiments; optional covariate balancing (Ustyuzhanin 2026)
 - [EfficientDiD](https://diff-diff.readthedocs.io/en/stable/api/efficient_did.html) - Chen, Sant'Anna & Xie (2025) efficient DiD with optimal weighting for tighter SEs
 - [TROP](https://diff-diff.readthedocs.io/en/stable/api/trop.html) - Triply Robust Panel estimator (Athey et al. 2025) with nuclear norm factor adjustment
-- [StaggeredTripleDifference](https://diff-diff.readthedocs.io/en/stable/api/staggered.html#staggeredtripledifference) - Ortiz-Villavicencio & Sant'Anna (2025) staggered DDD with group-time ATT
+- [StaggeredTripleDifference](https://diff-diff.readthedocs.io/en/stable/api/staggered.html#staggeredtripledifference) - Ortiz-Villavicencio & Sant'Anna (2025) staggered DDD with group-time ATT (deprecated 3.9 - use `TripleDifference` with `first_treat=`)
 - [WooldridgeDiD](https://diff-diff.readthedocs.io/en/stable/api/wooldridge_etwfe.html) - Wooldridge (2023, 2025) ETWFE: saturated OLS, logit/Poisson QMLE (ASF-based ATT). Alias `ETWFE`.
 - [LPDiD](https://diff-diff.readthedocs.io/en/stable/api/lpdid.html) - Dube, Girardi, Jorda & Taylor (2025) Local Projections DiD: per-horizon long-difference event study on clean controls (no negative weighting), variance- or equally-weighted ATT, for absorbing or non-absorbing (reversible) treatment
-- [ChangesInChanges](https://diff-diff.readthedocs.io/en/stable/api/changes_in_changes.html) - Athey & Imbens (2006) nonlinear/distributional DiD for the 2x2 design: full counterfactual distribution and quantile treatment effects via CDF transformation, plus the QDiD comparison estimator; bootstrap inference; R qte parity. Alias `CiC`
+- [ChangesInChanges](https://diff-diff.readthedocs.io/en/stable/api/changes_in_changes.html) - Athey & Imbens (2006) nonlinear/distributional DiD for the 2x2 design: full counterfactual distribution and quantile treatment effects via CDF transformation, plus the QDiD comparison estimator via `method="qdid"`; bootstrap inference; R qte parity. Alias `CiC`
 - [BaconDecomposition](https://diff-diff.readthedocs.io/en/stable/api/bacon.html) - Goodman-Bacon (2021) decomposition for diagnosing TWFE bias in staggered settings
 
 ## Diagnostics & Sensitivity
 
+- [RD Plots](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2015) optimal data-driven RD plots (`RDPlot`): all 8 rdrobust `binselect` bin selectors, implied-scale/WIMSE-weight reporting, optional matplotlib rendering
+- [Manipulation Testing](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Cattaneo, Jansson & Ma (2020) density-discontinuity test (`RDDensityTest`): rddensity 3.0 parity, robust bias-corrected inference, unrestricted/restricted models, mass-point adjustment
 - [Parallel Trends Testing](https://diff-diff.readthedocs.io/en/stable/api/diagnostics.html) - simple and Wasserstein-robust parallel trends tests, equivalence testing (TOST)
 - [Placebo Tests](https://diff-diff.readthedocs.io/en/stable/api/diagnostics.html) - placebo timing, group, permutation, leave-one-out
 - [Honest DiD](https://diff-diff.readthedocs.io/en/stable/api/honest_did.html) - Rambachan & Roth (2023) sensitivity analysis: robust CI under PT violations, breakdown values
 - [Pre-Trends Power Analysis](https://diff-diff.readthedocs.io/en/stable/api/pretrends.html) - Roth (2022) minimum detectable violation and power curves
 - [Power Analysis](https://diff-diff.readthedocs.io/en/stable/api/power.html) - analytical and simulation-based MDE, sample size, power curves for study design
+- [MMM Calibration Export](https://diff-diff.readthedocs.io/en/stable/api/mmm.html) - convert experiment results into MMM calibration inputs: PyMC-Marketing lift-test frames and Google Meridian lognormal ROI priors
 - Conley spatial HAC SE (`vcov_type="conley"`) on cross-sectional `LinearRegression` / `compute_robust_vcov` plus panel `DifferenceInDifferences` / `MultiPeriodDiD` / `TwoWayFixedEffects` (with `conley_lag_cutoff` for within-unit Bartlett temporal HAC) - Conley (1999) spatial-correlation-aware SEs with parity vs R `conleyreg` on cross-sectional + panel fixtures, optional combined spatial + cluster product kernel via explicit `cluster=`, auto-activating sparse k-d-tree fast path for `n > 5_000`
 
 ## Survey Support

{diff_diff-3.8.0 → diff_diff-3.9.1}/README.md

@@ -41,7 +41,7 @@ data = pd.DataFrame({
 })
 
 did = DifferenceInDifferences()
-results = did.fit(data, outcome='outcome', treatment='treated', time='post')
+results = did.fit(data, outcome='outcome', treatment='treated', post='post')
 print(results)              # DiDResults(ATT=3.0000, SE=1.7321, p=0.1583)
 results.print_summary()     # full statsmodels-style table
 ```
@@ -89,9 +89,9 @@ For rigorous DiD analysis, follow these 8 steps. Skipping diagnostic steps produ
 2. **State identification assumptions** - which parallel trends variant (unconditional, conditional, PT-GT-Nev, PT-GT-NYT), no-anticipation, overlap.
 3. **Test parallel trends** - simple 2x2: `check_parallel_trends()`, `equivalence_test_trends()`; staggered: inspect CS event-study pre-period coefficients (generic PT tests are invalid for staggered designs). Insignificant pre-trends do NOT prove PT holds.
 4. **Choose estimator** - staggered adoption -> CS/SA/BJS (NOT plain TWFE); few treated units -> SDiD; factor confounding -> TROP; simple 2x2 -> DiD. Run `BaconDecomposition` to diagnose TWFE bias.
-5. **Estimate** - `estimator.fit(data, ...)`. Always print the cluster count first and choose inference method based on the result (cluster-robust if >= 50 clusters, wild bootstrap if fewer).
-6. **Sensitivity analysis** - `compute_honest_did(results)` for bounds under PT violations (MultiPeriodDiD, CS, or dCDH), `run_all_placebo_tests()` for 2x2 falsification, specification comparisons for staggered designs.
-7. **Heterogeneity** - CS: `aggregate='group'`/`'event_study'`; SA: `results.event_study_effects` / `to_dataframe(level='cohort')`; subgroup re-estimation.
+5. **Estimate** - `estimator.fit(data, ...)`. Always print the cluster count first and choose inference method based on the result (cluster-robust if >= 50 clusters, wild bootstrap if fewer - for DifferenceInDifferences pass `cluster=`; TwoWayFixedEffects auto-clusters at unit level).
+6. **Sensitivity analysis** - `compute_honest_did(results)` for bounds under PT violations (MultiPeriodDiD, CS, or dCDH natively; the TwoWayFixedEffects `event_study=True` surface and a StackedDiD `results.aggregate('event_study')` container also admit - Stacked needs `kappa_pre >= 2`), `run_all_placebo_tests()` for 2x2 falsification, specification comparisons for staggered designs.
+7. **Heterogeneity** - CS: `results.aggregate('group')`/`'event_study'` (post-fit, no refit); SA: `results.event_study_effects` / `to_dataframe(level='cohort')`; Stacked: `results.aggregate('event_study')`/`'simple'` post-fit views (surface always computed since 3.9); EDiD: `results.aggregate(...)` post-fit from retained EIFs (3.9); ImputationDiD/TwoStageDiD: `results.aggregate(...)` post-fit from panel-backed kits (3.9); ContinuousDiD: `results.aggregate('dose'/'simple'/'event_study')` post-fit (3.9; dose/simple are views, event_study recomputes); subgroup re-estimation.
 8. **Robustness** - compare 2-3 estimators (CS vs SA vs BJS), report with and without covariates (shows whether conditioning drives identification), present pre-trends and sensitivity bounds.
 
 Full guide: `diff_diff.get_llm_guide("practitioner")`.
@@ -100,7 +100,7 @@ Full guide: `diff_diff.get_llm_guide("practitioner")`.
 
 - [DifferenceInDifferences](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - basic 2x2 DiD with robust/cluster-robust SEs, wild bootstrap, formula interface, and fixed effects
 - [TwoWayFixedEffects](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - panel data DiD with unit and time fixed effects via within-transformation or dummies
-- [MultiPeriodDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - event study design with period-specific treatment effects for dynamic analysis
+- [MultiPeriodDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - event study design with period-specific treatment effects for dynamic analysis (deprecated 3.9 - use TwoWayFixedEffects `event_study=True`)
 - [CallawaySantAnna](https://diff-diff.readthedocs.io/en/stable/api/staggered.html) - Callaway & Sant'Anna (2021) group-time ATT estimator for staggered adoption
 - [ChaisemartinDHaultfoeuille](https://diff-diff.readthedocs.io/en/stable/api/chaisemartin_dhaultfoeuille.html) - de Chaisemartin & D'Haultfœuille (2020/2022) for **reversible (non-absorbing) treatments** with multi-horizon event study, normalized effects, cost-benefit delta, sup-t bands, and dynamic placebos. The most general option for treatments that switch on AND off (see also `LPDiD`/`TROP` `non_absorbing`). Alias `DCDH`.
 - [SunAbraham](https://diff-diff.readthedocs.io/en/stable/api/staggered.html) - Sun & Abraham (2021) interaction-weighted estimator for heterogeneity-robust event studies
@@ -109,26 +109,29 @@ Full guide: `diff_diff.get_llm_guide("practitioner")`.
 - [SpilloverDiD](https://diff-diff.readthedocs.io/en/stable/api/spillover.html) - Butts (2021) ring-indicator spillover-aware DiD identifying direct effect on treated + per-ring spillover on near-control units; handles non-staggered and staggered timing; supports survey-design variance under `survey_design=` for HC1 / CR1 (Wave E.1 Binder TSL) and Conley (Wave E.2 panel-aware stratified-Conley sandwich on per-period PSU totals; extended in Wave E.2 follow-up to `conley_lag_cutoff > 0` via panel-block composition with within-PSU serial Bartlett HAC — `lag>0` requires an effective PSU via explicit `survey_design.psu` or injected `cluster=<col>`); `SurveyDesign.subpopulation()` preserves full-design `n_psu` / `df_survey` via zero-padded scores (Wave E.3, R `svyrecvar(subset())` form)
 - [SyntheticDiD](https://diff-diff.readthedocs.io/en/stable/api/estimators.html) - Synthetic DiD combining standard DiD and synthetic control for few treated units
 - [SyntheticControl](https://diff-diff.readthedocs.io/en/stable/api/synthetic_control.html) - Abadie, Diamond & Hainmueller (2010) classic synthetic control for a single treated unit (donor-weight counterfactual, nested/cv/inverse-variance/custom V; in-space placebo permutation inference via `in_space_placebo()`, plus ADH-2015 `leave_one_out()` + `in_time_placebo()` robustness, Firpo-Possebom (2018) test-inversion confidence sets, and Chernozhukov-Wüthrich-Zhu (2021) conformal inference)
-- [TripleDifference](https://diff-diff.readthedocs.io/en/stable/api/triple_diff.html) - triple difference (DDD) estimator for designs requiring two criteria for treatment eligibility
+- [TripleDifference](https://diff-diff.readthedocs.io/en/stable/api/triple_diff.html) - triple difference (DDD) estimator for designs requiring two criteria for treatment eligibility; serves both the 2x2x2 and the staggered-adoption design from one signature (`fit(..., first_treat=)` selects the staggered engine)
 - [ContinuousDiD](https://diff-diff.readthedocs.io/en/stable/api/continuous_did.html) - Callaway, Goodman-Bacon & Sant'Anna (2024) continuous treatment DiD with dose-response curves
 - [HeterogeneousAdoptionDiD](https://diff-diff.readthedocs.io/en/stable/api/had.html) - de Chaisemartin, Ciccia, D'Haultfœuille & Knau (2026) for designs where **no unit remains untreated**; local-linear estimator at the dose support boundary returning Weighted Average Slope (WAS) on Design 1' (`d̲ = 0` / QUG) or `WAS_{d̲}` on Design 1 (`d̲ > 0`, continuous-near-d̲ or mass-point), with a multi-period event-study extension (last-treatment cohort, pointwise CIs). **Panel-only** in this release - repeated cross-sections rejected by the validator. Alias `HAD`.
-- [RegressionDiscontinuity](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2014) sharp, fuzzy, AND covariate-adjusted regression discontinuity with robust bias-corrected inference and rdrobust-parity bandwidth selection (all 10 selectors, mass-point handling; fuzzy via `treatment_col=` with a first-stage block and weak-identification warning; covariates via `covariates=` - CCFT 2019, same estimand, covariate-aware bandwidths). Canonical `att` is the bias-corrected estimate with a coherent robust CI (rdrobust's printed headline is `att_conventional`). Alias `RDD`.
+- [RegressionDiscontinuity](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2014) sharp, fuzzy, AND covariate-adjusted regression discontinuity with robust bias-corrected inference and rdrobust-parity bandwidth selection (all 10 selectors, mass-point handling; fuzzy via `takeup=` with a first-stage block and weak-identification warning; covariates via `covariates=` - CCFT 2019, same estimand, covariate-aware bandwidths). Canonical `att` is the bias-corrected estimate with a coherent robust CI (rdrobust's printed headline is `att_conventional`). Alias `RDD`.
 - [StackedDiD](https://diff-diff.readthedocs.io/en/stable/api/stacked_did.html) - Wing, Freedman & Hollingsworth (2024) stacked DiD with Q-weights and sub-experiments; optional covariate balancing (Ustyuzhanin 2026)
 - [EfficientDiD](https://diff-diff.readthedocs.io/en/stable/api/efficient_did.html) - Chen, Sant'Anna & Xie (2025) efficient DiD with optimal weighting for tighter SEs
 - [TROP](https://diff-diff.readthedocs.io/en/stable/api/trop.html) - Triply Robust Panel estimator (Athey et al. 2025) with nuclear norm factor adjustment
-- [StaggeredTripleDifference](https://diff-diff.readthedocs.io/en/stable/api/staggered.html#staggeredtripledifference) - Ortiz-Villavicencio & Sant'Anna (2025) staggered DDD with group-time ATT
+- [StaggeredTripleDifference](https://diff-diff.readthedocs.io/en/stable/api/staggered.html#staggeredtripledifference) - Ortiz-Villavicencio & Sant'Anna (2025) staggered DDD with group-time ATT (deprecated 3.9 - use `TripleDifference` with `first_treat=`)
 - [WooldridgeDiD](https://diff-diff.readthedocs.io/en/stable/api/wooldridge_etwfe.html) - Wooldridge (2023, 2025) ETWFE: saturated OLS, logit/Poisson QMLE (ASF-based ATT). Alias `ETWFE`.
 - [LPDiD](https://diff-diff.readthedocs.io/en/stable/api/lpdid.html) - Dube, Girardi, Jorda & Taylor (2025) Local Projections DiD: per-horizon long-difference event study on clean controls (no negative weighting), variance- or equally-weighted ATT, for absorbing or non-absorbing (reversible) treatment
-- [ChangesInChanges](https://diff-diff.readthedocs.io/en/stable/api/changes_in_changes.html) - Athey & Imbens (2006) nonlinear/distributional DiD for the 2x2 design: full counterfactual distribution and quantile treatment effects via CDF transformation, plus the QDiD comparison estimator; bootstrap inference; R qte parity. Alias `CiC`
+- [ChangesInChanges](https://diff-diff.readthedocs.io/en/stable/api/changes_in_changes.html) - Athey & Imbens (2006) nonlinear/distributional DiD for the 2x2 design: full counterfactual distribution and quantile treatment effects via CDF transformation, plus the QDiD comparison estimator via `method="qdid"`; bootstrap inference; R qte parity. Alias `CiC`
 - [BaconDecomposition](https://diff-diff.readthedocs.io/en/stable/api/bacon.html) - Goodman-Bacon (2021) decomposition for diagnosing TWFE bias in staggered settings
 
 ## Diagnostics & Sensitivity
 
+- [RD Plots](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Calonico, Cattaneo & Titiunik (2015) optimal data-driven RD plots (`RDPlot`): all 8 rdrobust `binselect` bin selectors, implied-scale/WIMSE-weight reporting, optional matplotlib rendering
+- [Manipulation Testing](https://diff-diff.readthedocs.io/en/stable/api/regression_discontinuity.html) - Cattaneo, Jansson & Ma (2020) density-discontinuity test (`RDDensityTest`): rddensity 3.0 parity, robust bias-corrected inference, unrestricted/restricted models, mass-point adjustment
 - [Parallel Trends Testing](https://diff-diff.readthedocs.io/en/stable/api/diagnostics.html) - simple and Wasserstein-robust parallel trends tests, equivalence testing (TOST)
 - [Placebo Tests](https://diff-diff.readthedocs.io/en/stable/api/diagnostics.html) - placebo timing, group, permutation, leave-one-out
 - [Honest DiD](https://diff-diff.readthedocs.io/en/stable/api/honest_did.html) - Rambachan & Roth (2023) sensitivity analysis: robust CI under PT violations, breakdown values
 - [Pre-Trends Power Analysis](https://diff-diff.readthedocs.io/en/stable/api/pretrends.html) - Roth (2022) minimum detectable violation and power curves
 - [Power Analysis](https://diff-diff.readthedocs.io/en/stable/api/power.html) - analytical and simulation-based MDE, sample size, power curves for study design
+- [MMM Calibration Export](https://diff-diff.readthedocs.io/en/stable/api/mmm.html) - convert experiment results into MMM calibration inputs: PyMC-Marketing lift-test frames and Google Meridian lognormal ROI priors
 - Conley spatial HAC SE (`vcov_type="conley"`) on cross-sectional `LinearRegression` / `compute_robust_vcov` plus panel `DifferenceInDifferences` / `MultiPeriodDiD` / `TwoWayFixedEffects` (with `conley_lag_cutoff` for within-unit Bartlett temporal HAC) - Conley (1999) spatial-correlation-aware SEs with parity vs R `conleyreg` on cross-sectional + panel fixtures, optional combined spatial + cluster product kernel via explicit `cluster=`, auto-activating sparse k-d-tree fast path for `n > 5_000`
 
 ## Survey Support

{diff_diff-3.8.0 → diff_diff-3.9.1}/diff_diff/__init__.py

@@ -20,6 +20,9 @@ This library provides sklearn-like estimators for causal inference using
 the difference-in-differences methodology.
 """
 
+import warnings as _warnings
+from typing import Any as _Any
+
 # Import backend detection from dedicated module (avoids circular imports)
 from diff_diff._backend import (
     HAS_RUST_BACKEND,
@@ -30,6 +33,9 @@ from diff_diff._backend import (
 )
 from diff_diff._guides_api import get_llm_guide
 from diff_diff.agent_workflow import agent_workflow
+from diff_diff.aggregation import (
+    AggregationResult,
+)
 from diff_diff.bacon import (
     BaconDecomposition,
     BaconDecompositionResults,
@@ -154,6 +160,12 @@ from diff_diff.local_linear import (
 )
 from diff_diff.lpdid import LPDiD
 from diff_diff.lpdid_results import LPDiDResults
+from diff_diff.mmm import (
+    MeridianROIPrior,
+    meridian_calibration_mask,
+    to_meridian_roi_prior,
+    to_pymc_marketing_lift_test,
+)
 from diff_diff.power import (
     PowerAnalysis,
     PowerResults,
@@ -212,6 +224,14 @@ from diff_diff.rdd import (
     RegressionDiscontinuity,
     RegressionDiscontinuityResults,
 )
+from diff_diff.rddensity import (
+    RDDensityTest,
+    RDDensityTestResult,
+)
+from diff_diff.rdplot import (
+    RDPlot,
+    RDPlotResult,
+)
 from diff_diff.results import (
     DiDResults,
     MultiPeriodDiDResults,
@@ -219,6 +239,11 @@ from diff_diff.results import (
     SpilloverDiDResults,  # re-export
     SyntheticDiDResults,
 )
+from diff_diff.results_base import (
+    BaseResults,
+    Diagnostic,
+    EventStudyResults,
+)
 from diff_diff.spillover import (
     SpilloverDiD,
 )
@@ -301,13 +326,10 @@ TWFE = TwoWayFixedEffects
 EventStudy = MultiPeriodDiD
 SDiD = SyntheticDiD
 CS = CallawaySantAnna
-CDiD = ContinuousDiD
 SA = SunAbraham
 BJS = ImputationDiD
-Gardner = TwoStageDiD
 DDD = TripleDifference
 SDDD = StaggeredTripleDifference
-Stacked = StackedDiD
 Bacon = BaconDecomposition
 EDiD = EfficientDiD
 ETWFE = WooldridgeDiD
@@ -315,8 +337,36 @@ DCDH = ChaisemartinDHaultfoeuille
 HAD = HeterogeneousAdoptionDiD
 CiC = ChangesInChanges
 RDD = RegressionDiscontinuity
+SCM = SyntheticControl
+
+# Alias diet (rows M-132..M-134, mechanism M-135): CDiD / Gardner /
+# Stacked are deprecated in 3.9 and removed in 4.0. They deliberately
+# do NOT live in module globals — dir()/vars() no longer list them —
+# but stay importable (and in __all__) through 3.9, served by the PEP
+# 562 module __getattr__ below, which emits the FutureWarning naming
+# the surviving class.
+_DEPRECATED_ALIASES = {
+    "CDiD": "ContinuousDiD",
+    "Gardner": "TwoStageDiD",
+    "Stacked": "StackedDiD",
+}
+
+
+def __getattr__(name: str) -> _Any:
+    """PEP 562 warning shim for the dieted aliases (row M-135)."""
+    target = _DEPRECATED_ALIASES.get(name)
+    if target is not None:
+        _warnings.warn(
+            f"diff_diff.{name} is deprecated and will be removed in 4.0; "
+            f"use diff_diff.{target}.",
+            FutureWarning,
+            stacklevel=2,
+        )
+        return globals()[target]
+    raise AttributeError(f"module 'diff_diff' has no attribute {name!r}")
+
 
-__version__ = "3.8.0"
+__version__ = "3.9.1"
 __all__ = [
     # Estimators
     "DifferenceInDifferences",
@@ -336,7 +386,9 @@ __all__ = [
     "TROP",
     "SyntheticControl",
     "StackedDiD",
-    # Estimator aliases (short names)
+    # Estimator aliases (short names). CDiD / Gardner / Stacked are
+    # deprecated (M-132..M-134): still importable through 3.9 via the
+    # module __getattr__ (M-135), gone from module globals/dir().
     "DiD",
     "TWFE",
     "EventStudy",
@@ -350,6 +402,7 @@ __all__ = [
     "Gardner",
     "DDD",
     "SDDD",
+    "SCM",
     "Stacked",
     "Bacon",
     # Bacon Decomposition
@@ -519,6 +572,12 @@ __all__ = [
     "RegressionDiscontinuity",
     "RegressionDiscontinuityResults",
     "RDD",
+    # RDPlot (data-driven RD plots, rdplot parity)
+    "RDPlot",
+    "RDPlotResult",
+    # RDDensityTest (manipulation testing, rddensity parity)
+    "RDDensityTest",
+    "RDDensityTestResult",
     # HeterogeneousAdoptionDiD pre-test diagnostics (Phase 3)
     "qug_test",
     "stute_test",
@@ -529,7 +588,7 @@ __all__ = [
     "YatchewTestResults",
     "HADPretestReport",
     # HAD joint pre-tests (Phase 3 follow-up) — multi-period event-study
-    # workflow dispatch via did_had_pretest_workflow(aggregate="event_study")
+    # workflow dispatch selected from the panel shape (M-139)
     "stute_joint_pretest",
     "joint_pretrends_test",
     "joint_homogeneity_test",
@@ -559,8 +618,18 @@ __all__ = [
     "Alert",
     "OutcomeShape",
     "TreatmentDoseShape",
+    # MMM calibration export (interop)
+    "to_pymc_marketing_lift_test",
+    "to_meridian_roi_prior",
+    "meridian_calibration_mask",
+    "MeridianROIPrior",
     # LLM guide accessor
     "get_llm_guide",
+    # Results-contract foundations (4.0 program Phase 2, spec sections 3.5/5)
+    "BaseResults",
+    "Diagnostic",
+    "EventStudyResults",
+    "AggregationResult",
 ]
 
 # Agent-facing entrypoints surface first in dir(diff_diff). LLM agents

diff_diff-3.9.1/diff_diff/_base.py

@@ -0,0 +1,169 @@
+"""Shared estimator parameter surface (4.0 program, Phase 2 PR 2(c)-i).
+
+This module is a deliberate LEAF: it imports only the standard library, never
+other ``diff_diff`` modules, so every estimator module can import it without
+cycles (the ``results_base.py`` precedent).
+
+One public symbol:
+
+``BaseEstimator``
+    Mixin providing the sklearn-compatible ``get_params`` / ``set_params``
+    pair for every estimator class, replacing the 25 hand-rolled copies
+    (see ``docs/v4-design.md`` section 7, "Constructor hygiene").
+
+Design contract (the locked transactional rule):
+
+- Parameter names are introspected from the concrete class's ``__init__``
+  signature, so the parameter surface can never drift from the constructor.
+- ``get_params`` returns ``{name: getattr(self, alias_or_name)}`` for every
+  constructor parameter. ``deep`` is accepted for sklearn compatibility
+  (``sklearn.base.clone`` calls ``get_params(deep=False)``) and ignored:
+  no diff-diff estimator nests another estimator.
+- ``set_params`` is TRANSACTIONAL via probe re-init: it validates the merged
+  configuration by constructing a throwaway ``type(self)(**merged)`` - which
+  raises before ``self`` is touched - then adopts the probe's parameter
+  attributes plus each class's declared derived-config attributes. Fitted
+  state (attributes ``fit()`` sets) is never touched, because only declared
+  attributes are copied. The behavioral contract, stated post-normalization:
+
+        est.set_params(**p)  ==  config state of
+        type(est)(**{**est.get_params(), **cls._normalize_set_params(p)})
+
+  Validation is therefore exactly ``__init__``'s validation, eagerly, and
+  can never drift from it.
+
+Per-class accommodation hooks (declarative; default no-op):
+
+``_PARAM_ATTR_ALIASES``
+    Maps a constructor parameter name to the attribute that stores its RAW
+    value when the two differ (e.g. ``DifferenceInDifferences`` stores the
+    raw ``vcov_type`` argument under ``_vcov_type_arg`` and the resolved
+    value under ``vcov_type``).
+``_DERIVED_CONFIG_ATTRS``
+    Extra attributes ``__init__`` computes FROM the parameters that must be
+    re-adopted from the probe after a successful ``set_params`` (e.g.
+    ``_vcov_type_explicit``).
+``_normalize_set_params``
+    Classmethod rewriting the user-supplied ``params`` dict before the
+    merge (e.g. the ``robust``-alone alias re-derivation on
+    ``DifferenceInDifferences``).
+"""
+
+import inspect
+from typing import Any, ClassVar, Dict, Mapping, Tuple, TypeVar
+
+__all__ = ["BaseEstimator"]
+
+TSelf = TypeVar("TSelf", bound="BaseEstimator")
+
+# Signature kinds accepted as constructor parameters. VAR_POSITIONAL /
+# VAR_KEYWORD are rejected loudly: an estimator constructor forwarding
+# through *args/**kwargs has no introspectable parameter surface.
+_ACCEPTED_KINDS = (
+    inspect.Parameter.POSITIONAL_OR_KEYWORD,
+    inspect.Parameter.KEYWORD_ONLY,
+)
+
+
+class BaseEstimator:
+    """Mixin providing sklearn-compatible ``get_params`` / ``set_params``.
+
+    See the module docstring for the full contract. Subclasses only declare
+    the accommodation hooks when their constructor stores state under
+    different names than its parameters; most estimators are plain drop-ins.
+    """
+
+    _PARAM_ATTR_ALIASES: ClassVar[Mapping[str, str]] = {}
+    _DERIVED_CONFIG_ATTRS: ClassVar[Tuple[str, ...]] = ()
+
+    @classmethod
+    def _param_names(cls) -> Tuple[str, ...]:
+        """Constructor parameter names, introspected once per class.
+
+        Cached via ``cls.__dict__`` (never a plain class attribute, which a
+        subclass with its own ``__init__`` would inherit; never an instance
+        attribute, which would break ``vars(est)`` stability across fit).
+        """
+        cached = cls.__dict__.get("_param_names_cache")
+        if cached is not None:
+            return cached
+        names = []
+        for name, param in inspect.signature(cls.__init__).parameters.items():
+            if name == "self":
+                continue
+            if param.kind not in _ACCEPTED_KINDS:
+                raise TypeError(
+                    f"{cls.__name__}.__init__ uses *args/**kwargs "
+                    f"(parameter {name!r}); its parameter surface cannot be "
+                    "introspected for get_params/set_params."
+                )
+            names.append(name)
+        result: Tuple[str, ...] = tuple(names)
+        setattr(cls, "_param_names_cache", result)
+        return result
+
+    @classmethod
+    def _normalize_set_params(cls, params: Dict[str, Any]) -> Dict[str, Any]:
+        """Rewrite the user-supplied params before the merge (default no-op)."""
+        return params
+
+    def get_params(self, deep: bool = True) -> Dict[str, Any]:
+        """
+        Get estimator parameters (sklearn-compatible).
+
+        Parameters
+        ----------
+        deep : bool, default True
+            Accepted for sklearn compatibility (``sklearn.base.clone``
+            calls ``get_params(deep=False)``) and ignored: no diff-diff
+            estimator nests another estimator.
+
+        Returns
+        -------
+        Dict[str, Any]
+            Estimator parameters suitable for passing to ``__init__``.
+            Keys follow the ``__init__`` signature order. Where a class
+            stores a parameter's raw value under a different attribute
+            (``_PARAM_ATTR_ALIASES``), the raw value is returned so that
+            ``type(est)(**est.get_params())`` reconstructs the same
+            configuration.
+        """
+        del deep
+        aliases = self._PARAM_ATTR_ALIASES
+        return {name: getattr(self, aliases.get(name, name)) for name in self._param_names()}
+
+    def set_params(self: TSelf, **params: Any) -> TSelf:
+        """
+        Set estimator parameters (sklearn-compatible, transactional).
+
+        Validates the merged configuration by constructing a throwaway
+        instance of ``type(self)`` - so a rejected call raises the same
+        error ``__init__`` would and leaves this estimator unchanged - then
+        adopts the probe's parameter attributes and the class's declared
+        derived-config attributes. Fitted state is never touched.
+
+        Parameters
+        ----------
+        **params
+            Estimator parameters. Unknown names raise ``ValueError``
+            before any validation or mutation.
+
+        Returns
+        -------
+        self
+        """
+        names = self._param_names()
+        for key in params:
+            if key not in names:
+                raise ValueError(f"Unknown parameter: {key}")
+        cls = type(self)
+        normalized = cls._normalize_set_params(dict(params))
+        merged = {**self.get_params(), **normalized}
+        probe = cls(**merged)
+        aliases = self._PARAM_ATTR_ALIASES
+        for name in names:
+            attr = aliases.get(name, name)
+            setattr(self, attr, getattr(probe, attr))
+        for attr in self._DERIVED_CONFIG_ATTRS:
+            setattr(self, attr, getattr(probe, attr))
+        return self