Add DurationDiD for absorbing individual outcomes

Author: igerber

diff_diff/duration_did_results.py

@@ -254,23 +254,24 @@ def to_dataframe(self, level: str = "event_study") -> pd.DataFrame:
         raise ValueError("level must be 'event_study', 'simple', 'survival', or 'diagnostics'")
 
     def summary(self, alpha: Optional[float] = None) -> str:
-        """Summarize stored whole-population absorption ATT and availability."""
+        """Summarize stored absorption ATT and availability, printing shared reasons once."""
         _require_fit_alpha(alpha, self.alpha, message=_SUMMARY_ALPHA_MESSAGE)
         units = (
             "hazard gap per normalized observation interval"
             if self.method == "common_dynamics"
             else "dimensionless hazard ratio"
         )
+        reasons = dict.fromkeys(
+            reason for values in self.inference_reasons.values() for reason in values
+        )
         return (
             f"DurationDiD ({self.method})\n"
             f"Coefficient ({units}): {self.coefficient:.6g}\n"
             f"Mean cumulative absorption ATT: {self.att:.6g}; SE: {self.se:.6g}\n"
             f"{_coverage_pct(self.alpha)}% CI: {self.conf_int}; p-value: {self.p_value:.6g}\n"
             f"Estimation: {self.estimation_status}; inference: {self.inference_status}\n"
             f"Pooled individual bootstrap: {self.n_bootstrap_valid}/{self.n_bootstrap} valid effect draws"
-            + "".join(
-                "\n" + reason for values in self.inference_reasons.values() for reason in values
-            )
+            + "".join("\n" + reason for reason in reasons)
             + "".join("\n" + warning for warning in self.support_warnings)
             + "\n"
             + self.pretrend_results.summary()

diff_diff/guides/llms-full.txt

@@ -115,6 +115,8 @@ Owns `effects`, `survival_curve`, `bootstrap_failures`, B-by-post `bootstrap_eff
 `inference_reasons` have pointwise/simultaneous/simple keys; available/unavailable,
 plus partial for pointwise. `support_warnings` and `pretrend_results` remain
 inspectable. Properties: `inference_method='bootstrap'`, uniform raw post `raw_att`.
+`summary()` prints distinct inference-unavailability reasons once in first appearance
+order; the full per-family `inference_reasons` mapping remains in the result and `to_dict()`.
 
 Native table column order (chronological RangeIndex):
 - effects (post only): period, event_time, att, se, t_stat, p_value, conf_int_lower,

docs/api/duration_did.rst

@@ -135,6 +135,10 @@ mappings as period/value records. Native tables have chronological RangeIndex.
   ``support_warnings`` and stored ``pretrend_results`` remain inspectable.
 - Properties: ``inference_method='bootstrap'`` and uniform raw post mean ``raw_att``.
 
+``summary()`` prints each distinct inference-unavailability reason once, in first
+appearance order. The full per-family ``inference_reasons`` mapping is preserved
+on the result and in ``to_dict()``.
+
 Native schemas, in column order:
 
 - ``effects`` (post only): period, event_time, att, se, t_stat, p_value,

tests/test_duration_did_results.py

@@ -93,6 +93,19 @@ def test_serialization_confidence_and_metadata(result):
     assert "dimensionless hazard ratio" in ph.summary()
 
 
+def test_summary_reports_shared_inference_failure_once():
+    data = duration_panel(control=(3, 2, 1, 1, 1, 1), treated=(3, 2, 1, 1, 1, 1), n=3, scale=1)
+    r = fit_duration(data, n_bootstrap=15, seed=2)
+    stored = r.to_dict()
+    reason = "one or more effect bootstrap draws failed"
+    assert all(reason in r.inference_reasons[family] for family in r.inference_status)
+    text = r.summary()
+    assert text.splitlines().count(reason) == 1
+    assert all(warning in text for warning in r.support_warnings)
+    assert r.pretrend_results.summary() in text
+    assert r.to_dict() == stored
+
+
 @pytest.mark.parametrize(
     "clock", [pd.date_range("2020", periods=6), pd.timedelta_range("0 days", periods=6)]
 )