rd2d
Estimate causal effects at a two-dimensional cutoff boundary using local polynomial regression discontinuity, with data-driven bandwidth selection and diagnostics.
Instructions
2D boundary regression discontinuity estimation. Assumptions: Conditional expectations of potential outcomes are continuous at the cutoff; Units cannot precisely manipulate the running variable around the cutoff (no sorting); For fuzzy designs: monotonicity of treatment take-up at the cutoff. Pre-conditions: A continuous running/forcing variable with a known cutoff that (sharply or fuzzily) assigns treatment; Enough observations in a neighbourhood of the cutoff to fit a local polynomial. Failure modes: Density of the running variable jumps at the cutoff (manipulation / sorting) -> Run a McCrary / density test (rdplotdensity); if manipulation is present the design is invalid near the cutoff; Estimate swings with the bandwidth -- results are not robust -> Report a bandwidth-sensitivity curve and use a data-driven MSE-optimal bandwidth. Alternatives: sp.rdrobust, sp.rdrandinf, sp.rdbwselect. Typical minimum N: 500.
Input Schema
| Name | Required | Description | Default |
|---|---|---|---|
| h | No | Manual bandwidth. If None, MSE-optimal bandwidth is selected. | |
| p | No | Polynomial order for point estimation (1 = local linear). | |
| y | Yes | Outcome variable name. | |
| x1 | Yes | First running variable (score dimension 1). | |
| x2 | Yes | Second running variable (score dimension 2). | |
| alpha | No | Significance level for confidence intervals. | |
| detail | No | Payload depth: 'minimal' (~150 tokens) for sub-step calls where only the point estimate is needed; 'standard' (~1K tokens) for diagnostics + coefficient table; 'agent' (~2K tokens, default) adds violations / next_steps / suggested_functions so the LLM can plan its next call without another round-trip. | agent |
| kernel | No | Kernel function: 'triangular', 'uniform', or 'epanechnikov'. | triangular |
| n_eval | No | Number of evaluation points when ``eval_points`` is None. Use 1 for a single pooled effect. | |
| approach | No | ``'distance'``: project onto signed distance to boundary, then apply univariate local polynomial RD. ``'location'``: fit bivariate local polynomial on each side of the boundary at evaluation points. | distance |
| boundary | No | Function ``f(x1) -> x2`` defining the boundary curve. If None, the boundary is the vertical line ``x1 = 0``. | |
| bwselect | No | Bandwidth selection method (used when ``h`` is None). | mserd |
| as_handle | No | If true, cache the fitted result on the server and return result_id + result_uri alongside the JSON payload so a subsequent tools/call can chain without re-running. | |
| data_path | Yes | Absolute path or URL to a data file. Supported: .csv / .tsv / .txt (delimited), .parquet / .pq, .feather / .arrow, .xlsx / .xls, .dta (Stata), .json / .jsonl. Schemes: file://, s3://, gs://, https://. | |
| result_id | No | Optional handle to a previously-fitted result (returned by an earlier call when as_handle=true). Tools that operate on a fitted object accept this in place of re-supplying data_path + columns. | |
| treatment | Yes | Binary treatment indicator (1 = treated, 0 = control). | |
| eval_points | No | Shape ``(k, 2)`` array of boundary evaluation points. If None, points are automatically selected along the boundary. | |
| data_columns | No | Optional column projection. Parquet/Feather/Stata loaders honour this for fast partial reads. | |
| data_sample_n | No | Optional uniform random subsample size (seed=0, deterministic) — useful on huge panels. |
Output Schema
| Name | Required | Description | Default |
|---|---|---|---|
No arguments | |||