linearmodels Skill
linearmodels: panel data, IV/GMM, system regression, and asset pricing models in Python. Covers PanelOLS (FE/RE), BetweenOLS, FirstDifferenceOLS, Fama-MacBeth, IV2SLS/LIML/GMM, SUR, IV3SLS, and Driscoll-Kraay SEs. Use for random effects estimation, between or first-difference panel models, system estimation (SUR, 3SLS), LIML/GMM instrumental variables, Fama-MacBeth regressions, or Driscoll-Kraay standard errors. Complements pyfixest (high-dimensional FE + DiD) and statsmodels (GLM + time series).
Comprehensive skill for panel data estimation, instrumental variables, system regression, and asset pricing with linearmodels (Kevin Sheppard). Use decision trees below to find the right guidance, then load detailed references.
What is linearmodels?
linearmodels extends statsmodels with specialized model classes for structured data:
- Panel data: PanelOLS (fixed effects), RandomEffects, BetweenOLS, FirstDifferenceOLS, PooledOLS, FamaMacBeth
- Instrumental variables: IV2SLS, IVLIML (k-class), IVGMM, IVGMMCUE (continuously updating), AbsorbingLS
- System estimation: SUR (Seemingly Unrelated Regression), IV3SLS, IVSystemGMM
- Asset pricing: LinearFactorModel, LinearFactorModelGMM, TradedFactorModel
- Rich inference: Driscoll-Kraay, clustered (1- and 2-way), HAC kernels (Bartlett, Parzen, Quadratic Spectral)
- Dual API: Formula-based (via formulaic) and array-based interfaces
How to Use This Skill
Reference File Structure
| File |
Purpose |
When to Read |
quickstart.md |
Installation, MultiIndex setup, formula vs array API, first model |
Starting with linearmodels |
panel-models.md |
PanelOLS, RandomEffects, BetweenOLS, FD, Pooled, FamaMacBeth |
Panel data estimation |
iv-models.md |
IV2SLS, IVLIML, IVGMM, IVGMMCUE, AbsorbingLS |
IV / GMM estimation |
system-models.md |
SUR, IV3SLS, IVSystemGMM, cross-equation constraints |
System estimation |
asset-pricing.md |
LinearFactorModel, TradedFactorModel, GMM estimation |
Asset pricing tests |
covariance-inference.md |
All SE types, Driscoll-Kraay, clustering, GMM weights |
Choosing standard errors |
gotchas.md |
MultiIndex requirement, pyfixest/statsmodels boundary, limits |
Debugging issues |
Reading Order
- New to linearmodels? Start with
quickstart.md then panel-models.md
- Need IV/GMM? Read
quickstart.md then iv-models.md
- System estimation (SUR/3SLS)? Read
quickstart.md then system-models.md
- Asset pricing? Read
quickstart.md then asset-pricing.md
- Choosing SEs? Read
covariance-inference.md
- Coming from pyfixest? Read
quickstart.md then gotchas.md
Related Skills
| Skill |
Relationship |
pyfixest |
Preferred for high-dimensional FE, FE + IV, DiD, fast demeaning, publication tables. Use linearmodels when pyfixest cannot do what you need (RE, system models, LIML/GMM, Fama-MacBeth) |
statsmodels |
Foundation library. Use statsmodels for GLM, time series, diagnostics. linearmodels extends statsmodels for panel/IV/system models |
svy |
Survey-weighted regression with complex survey designs. linearmodels supports weights for population/precision weighting in panel models, but this is NOT equivalent to design-based survey inference — it does not handle stratification, clustering as a design feature, or replicate weights. If your data comes from a complex probability survey, use svy |
data-scientist |
Methodology guidance — load for "why and when" behind model choices |
polars |
Data preparation before estimation; convert to pandas with .to_pandas() before passing to linearmodels |
Quick Decision Trees
"I need a panel model"
What panel estimation method?
├─ Fixed effects (within estimator)
│ ├─ 1-2 way FE, no IV → linearmodels PanelOLS or pyfixest feols
│ ├─ 3+ way FE → pyfixest (linearmodels max 2-way in PanelOLS)
│ ├─ FE + IV combined → pyfixest (linearmodels has no Panel IV)
│ └─ FE + DiD → pyfixest (linearmodels has no DiD)
├─ Random effects (GLS) → linearmodels RandomEffects
│ └─ → ./references/panel-models.md
├─ FE vs RE comparison → linearmodels (run both, compare)
│ └─ → ./references/panel-models.md
├─ Between estimator → linearmodels BetweenOLS
│ └─ → ./references/panel-models.md
├─ First difference → linearmodels FirstDifferenceOLS
│ └─ → ./references/panel-models.md
├─ Pooled OLS (panel-aware SEs) → linearmodels PooledOLS
│ └─ → ./references/panel-models.md
└─ Fama-MacBeth → linearmodels FamaMacBeth
└─ → ./references/panel-models.md
"I need IV / GMM estimation"
What IV method?
├─ 2SLS (standard IV)
│ ├─ With fixed effects → pyfixest (linearmodels has no Panel IV)
│ └─ Without FE → linearmodels IV2SLS or pyfixest
│ └─ → ./references/iv-models.md
├─ LIML / k-class (better finite-sample) → linearmodels IVLIML
│ └─ → ./references/iv-models.md
├─ GMM-IV (efficient, overidentified) → linearmodels IVGMM
│ └─ → ./references/iv-models.md
├─ Continuously updating GMM → linearmodels IVGMMCUE
│ └─ → ./references/iv-models.md
└─ High-dimensional absorbed FE (OLS) → linearmodels AbsorbingLS
└─ → ./references/iv-models.md
"I need system estimation"
System of equations?
├─ Multiple equations, correlated errors → SUR
│ └─ → ./references/system-models.md
├─ Multiple equations + endogenous variables → IV3SLS
│ └─ → ./references/system-models.md
├─ System GMM → IVSystemGMM
│ └─ → ./references/system-models.md
├─ Cross-equation parameter restrictions → LinearConstraint
│ └─ → ./references/system-models.md
└─ Not sure which → Start with SUR
└─ → ./references/system-models.md
"Something isn't working"
Having issues?
├─ TypeError about DataFrame index → ./references/gotchas.md
├─ Need FE + IV in one model → ./references/gotchas.md
├─ Need 3+ way fixed effects → ./references/gotchas.md
├─ Constant term confusion → ./references/gotchas.md
├─ Formula parsing errors → ./references/gotchas.md
├─ Want to compare with pyfixest → ./references/gotchas.md
└─ SUR performance issues → ./references/gotchas.md
File-First Execution in Research Workflows
Important: In data research pipelines (see CLAUDE.md), linearmodels estimation is executed through script files, not interactively. This ensures auditability and reproducibility.
The pattern:
- Write model code to
scripts/stage8_analysis/{step}_{task-name}.py
- Execute via Bash with automatic output capture wrapper script
- Validation results get automatically embedded in scripts as comments
- If failed, create versioned copy for fixes
Closely read agent_reference/SCRIPT_EXECUTION_REFERENCE.md for the mandatory file-first execution protocol covering complete code file writing, output capture, and file versioning rules.
The examples below show linearmodels syntax. In research workflows, wrap them in scripts following the file-first pattern.
Quick Reference
Essential Imports
from linearmodels.panel import PanelOLS, RandomEffects, BetweenOLS
from linearmodels.panel import FirstDifferenceOLS, PooledOLS, FamaMacBeth
from linearmodels.iv import IV2SLS, IVLIML, IVGMM, IVGMMCUE, AbsorbingLS
from linearmodels.system import SUR, IV3SLS, IVSystemGMM
from linearmodels.panel import compare # Panel model comparison tables
Data Setup (Critical — MultiIndex Required for Panel Models)
import pandas as pd
# Panel data MUST have a MultiIndex with (entity, time)
df = df.set_index(["entity_id", "year"])
# Verify the index
print(f"Index names: {df.index.names}")
print(f"Index levels: {df.index.nlevels}")
Core Operations
| Operation |
Code |
| Panel FE (formula) |
PanelOLS.from_formula("y ~ x1 + x2 + EntityEffects", data=df).fit() |
| Panel FE (array) |
PanelOLS(df.y, df[["x1","x2"]], entity_effects=True).fit() |
| Two-way FE |
PanelOLS.from_formula("y ~ x1 + EntityEffects + TimeEffects", data=df).fit() |
| Random effects |
RandomEffects.from_formula("y ~ 1 + x1 + x2", data=df).fit() |
| Between OLS |
BetweenOLS.from_formula("y ~ 1 + x1 + x2", data=df).fit() |
| First difference |
FirstDifferenceOLS.from_formula("y ~ x1 + x2", data=df).fit() |
| Fama-MacBeth |
FamaMacBeth.from_formula("y ~ 1 + x1 + x2", data=df).fit() |
| IV / 2SLS |
IV2SLS.from_formula("y ~ 1 + exog + [endog ~ inst]", data=df).fit() |
| LIML |
IVLIML.from_formula("y ~ 1 + exog + [endog ~ inst]", data=df).fit() |
| Clustered SEs |
mod.fit(cov_type="clustered", cluster_entity=True) |
| Driscoll-Kraay |
mod.fit(cov_type="kernel", kernel="bartlett", bandwidth=5) |
| Summary |
results.summary |
| Model comparison |
compare({"FE": fe_res, "RE": re_res}) |
Formula Syntax
# Panel FE keywords (appear in formula, not after |)
"y ~ x1 + x2 + EntityEffects" # Entity FE
"y ~ x1 + x2 + EntityEffects + TimeEffects" # Two-way FE
"y ~ x1 + x2 + TimeEffects" # Time FE only
# IV bracket notation
"y ~ 1 + exog + [endog ~ instrument1 + instrument2]"
# Suppress intercept
"y ~ x1 + x2 - 1"
Topic Index
| Topic |
Reference File |
| Installation |
./references/quickstart.md |
| MultiIndex data setup |
./references/quickstart.md |
| Formula vs array API |
./references/quickstart.md |
| First model |
./references/quickstart.md |
| Syntax comparison (pyfixest, statsmodels) |
./references/quickstart.md |
| PanelOLS (entity/time effects) |
./references/panel-models.md |
| RandomEffects |
./references/panel-models.md |
| BetweenOLS |
./references/panel-models.md |
| FirstDifferenceOLS |
./references/panel-models.md |
| PooledOLS |
./references/panel-models.md |
| FamaMacBeth |
./references/panel-models.md |
| FE vs RE decision |
./references/panel-models.md |
| Variance decomposition |
./references/panel-models.md |
| Weighted panel estimation |
./references/panel-models.md |
| R-squared types (within, between, overall) |
./references/panel-models.md |
| IV2SLS |
./references/iv-models.md |
| IVLIML and k-class estimators |
./references/iv-models.md |
| IVGMM (1-step, 2-step, iterative) |
./references/iv-models.md |
| IVGMMCUE |
./references/iv-models.md |
| AbsorbingLS (high-dim FE OLS) |
./references/iv-models.md |
| First-stage diagnostics |
./references/iv-models.md |
| Overidentification tests |
./references/iv-models.md |
| SUR (Seemingly Unrelated Regression) |
./references/system-models.md |
| IV3SLS |
./references/system-models.md |
| IVSystemGMM |
./references/system-models.md |
| Cross-equation constraints |
./references/system-models.md |
| LinearFactorModel |
./references/asset-pricing.md |
| TradedFactorModel |
./references/asset-pricing.md |
| Factor model GMM |
./references/asset-pricing.md |
| Driscoll-Kraay SEs |
./references/covariance-inference.md |
| Clustered SEs (entity, time, both) |
./references/covariance-inference.md |
| HAC / kernel covariance |
./references/covariance-inference.md |
| GMM weight matrices |
./references/covariance-inference.md |
| Debiased inference |
./references/covariance-inference.md |
| MultiIndex requirement |
./references/gotchas.md |
| Maximum 2-way FE limit |
./references/gotchas.md |
| No Panel IV |
./references/gotchas.md |
| pyfixest vs linearmodels boundary |
./references/gotchas.md |
| statsmodels vs linearmodels boundary |
./references/gotchas.md |
| Constant term handling |
./references/gotchas.md |
Citation
When this library is used as a primary analytical tool, include in the report's
Software & Tools references:
Sheppard, K. linearmodels: Econometric models for panel data, IV/GMM, and system regression [Computer software]. https://bashtage.github.io/linearmodels/
Cite when: linearmodels is used for panel estimation (RE, between), IV/GMM, Fama-MacBeth, or system regression (SUR, 3SLS).
Do not cite when: Only imported but no estimation performed.
For method-specific citations (e.g., individual estimators or techniques),
consult the reference files in this skill and agent_reference/CITATION_REFERENCE.md.
1---2name: linearmodels3description: Panel data, IV/GMM, system regression. PanelOLS (FE/RE), BetweenOLS, Fama-MacBeth, IV2SLS/LIML/GMM, SUR, 3SLS, Driscoll-Kraay SEs. Use for RE/between, system estimation, or GMM. Complements pyfixest (FE + DiD) and statsmodels (GLM + time series).4---56# linearmodels Skill78linearmodels: panel data, IV/GMM, system regression, and asset pricing models in Python. Covers PanelOLS (FE/RE), BetweenOLS, FirstDifferenceOLS, Fama-MacBeth, IV2SLS/LIML/GMM, SUR, IV3SLS, and Driscoll-Kraay SEs. Use for random effects estimation, between or first-difference panel models, system estimation (SUR, 3SLS), LIML/GMM instrumental variables, Fama-MacBeth regressions, or Driscoll-Kraay standard errors. Complements pyfixest (high-dimensional FE + DiD) and statsmodels (GLM + time series).910Comprehensive skill for panel data estimation, instrumental variables, system regression, and asset pricing with linearmodels (Kevin Sheppard). Use decision trees below to find the right guidance, then load detailed references.1112## What is linearmodels?1314linearmodels extends statsmodels with specialized model classes for structured data:15- **Panel data**: PanelOLS (fixed effects), RandomEffects, BetweenOLS, FirstDifferenceOLS, PooledOLS, FamaMacBeth16- **Instrumental variables**: IV2SLS, IVLIML (k-class), IVGMM, IVGMMCUE (continuously updating), AbsorbingLS17- **System estimation**: SUR (Seemingly Unrelated Regression), IV3SLS, IVSystemGMM18- **Asset pricing**: LinearFactorModel, LinearFactorModelGMM, TradedFactorModel19- **Rich inference**: Driscoll-Kraay, clustered (1- and 2-way), HAC kernels (Bartlett, Parzen, Quadratic Spectral)20- **Dual API**: Formula-based (via formulaic) and array-based interfaces2122## How to Use This Skill2324### Reference File Structure2526| File | Purpose | When to Read |27|------|---------|--------------|28| `quickstart.md` | Installation, MultiIndex setup, formula vs array API, first model | Starting with linearmodels |29| `panel-models.md` | PanelOLS, RandomEffects, BetweenOLS, FD, Pooled, FamaMacBeth | Panel data estimation |30| `iv-models.md` | IV2SLS, IVLIML, IVGMM, IVGMMCUE, AbsorbingLS | IV / GMM estimation |31| `system-models.md` | SUR, IV3SLS, IVSystemGMM, cross-equation constraints | System estimation |32| `asset-pricing.md` | LinearFactorModel, TradedFactorModel, GMM estimation | Asset pricing tests |33| `covariance-inference.md` | All SE types, Driscoll-Kraay, clustering, GMM weights | Choosing standard errors |34| `gotchas.md` | MultiIndex requirement, pyfixest/statsmodels boundary, limits | Debugging issues |3536### Reading Order37381. **New to linearmodels?** Start with `quickstart.md` then `panel-models.md`392. **Need IV/GMM?** Read `quickstart.md` then `iv-models.md`403. **System estimation (SUR/3SLS)?** Read `quickstart.md` then `system-models.md`414. **Asset pricing?** Read `quickstart.md` then `asset-pricing.md`425. **Choosing SEs?** Read `covariance-inference.md`436. **Coming from pyfixest?** Read `quickstart.md` then `gotchas.md`4445## Related Skills4647| Skill | Relationship |48|-------|-------------|49| `pyfixest` | Preferred for high-dimensional FE, FE + IV, DiD, fast demeaning, publication tables. Use linearmodels when pyfixest cannot do what you need (RE, system models, LIML/GMM, Fama-MacBeth) |50| `statsmodels` | Foundation library. Use statsmodels for GLM, time series, diagnostics. linearmodels extends statsmodels for panel/IV/system models |51| `svy` | Survey-weighted regression with complex survey designs. linearmodels supports `weights` for population/precision weighting in panel models, but this is NOT equivalent to design-based survey inference — it does not handle stratification, clustering as a design feature, or replicate weights. If your data comes from a complex probability survey, use `svy` |52| `data-scientist` | Methodology guidance — load for "why and when" behind model choices |53| `polars` | Data preparation before estimation; convert to pandas with `.to_pandas()` before passing to linearmodels |5455## Quick Decision Trees5657### "I need a panel model"5859```60What panel estimation method?61├─ Fixed effects (within estimator)62│ ├─ 1-2 way FE, no IV → linearmodels PanelOLS or pyfixest feols63│ ├─ 3+ way FE → pyfixest (linearmodels max 2-way in PanelOLS)64│ ├─ FE + IV combined → pyfixest (linearmodels has no Panel IV)65│ └─ FE + DiD → pyfixest (linearmodels has no DiD)66├─ Random effects (GLS) → linearmodels RandomEffects67│ └─ → ./references/panel-models.md68├─ FE vs RE comparison → linearmodels (run both, compare)69│ └─ → ./references/panel-models.md70├─ Between estimator → linearmodels BetweenOLS71│ └─ → ./references/panel-models.md72├─ First difference → linearmodels FirstDifferenceOLS73│ └─ → ./references/panel-models.md74├─ Pooled OLS (panel-aware SEs) → linearmodels PooledOLS75│ └─ → ./references/panel-models.md76└─ Fama-MacBeth → linearmodels FamaMacBeth77 └─ → ./references/panel-models.md78```7980### "I need IV / GMM estimation"8182```83What IV method?84├─ 2SLS (standard IV)85│ ├─ With fixed effects → pyfixest (linearmodels has no Panel IV)86│ └─ Without FE → linearmodels IV2SLS or pyfixest87│ └─ → ./references/iv-models.md88├─ LIML / k-class (better finite-sample) → linearmodels IVLIML89│ └─ → ./references/iv-models.md90├─ GMM-IV (efficient, overidentified) → linearmodels IVGMM91│ └─ → ./references/iv-models.md92├─ Continuously updating GMM → linearmodels IVGMMCUE93│ └─ → ./references/iv-models.md94└─ High-dimensional absorbed FE (OLS) → linearmodels AbsorbingLS95 └─ → ./references/iv-models.md96```9798### "I need system estimation"99100```101System of equations?102├─ Multiple equations, correlated errors → SUR103│ └─ → ./references/system-models.md104├─ Multiple equations + endogenous variables → IV3SLS105│ └─ → ./references/system-models.md106├─ System GMM → IVSystemGMM107│ └─ → ./references/system-models.md108├─ Cross-equation parameter restrictions → LinearConstraint109│ └─ → ./references/system-models.md110└─ Not sure which → Start with SUR111 └─ → ./references/system-models.md112```113114### "Something isn't working"115116```117Having issues?118├─ TypeError about DataFrame index → ./references/gotchas.md119├─ Need FE + IV in one model → ./references/gotchas.md120├─ Need 3+ way fixed effects → ./references/gotchas.md121├─ Constant term confusion → ./references/gotchas.md122├─ Formula parsing errors → ./references/gotchas.md123├─ Want to compare with pyfixest → ./references/gotchas.md124└─ SUR performance issues → ./references/gotchas.md125```126127## File-First Execution in Research Workflows128129**Important:** In data research pipelines (see `CLAUDE.md`), linearmodels estimation is executed through **script files**, not interactively. This ensures auditability and reproducibility.130131**The pattern:**1321. Write model code to `scripts/stage8_analysis/{step}_{task-name}.py`1332. Execute via Bash with automatic output capture wrapper script1343. Validation results get automatically embedded in scripts as comments1354. If failed, create versioned copy for fixes136137Closely read `agent_reference/SCRIPT_EXECUTION_REFERENCE.md` for the mandatory file-first execution protocol covering complete code file writing, output capture, and file versioning rules.138139The examples below show linearmodels syntax. In research workflows, wrap them in scripts following the file-first pattern.140141---142143## Quick Reference144145### Essential Imports146147```python148from linearmodels.panel import PanelOLS, RandomEffects, BetweenOLS149from linearmodels.panel import FirstDifferenceOLS, PooledOLS, FamaMacBeth150from linearmodels.iv import IV2SLS, IVLIML, IVGMM, IVGMMCUE, AbsorbingLS151from linearmodels.system import SUR, IV3SLS, IVSystemGMM152from linearmodels.panel import compare # Panel model comparison tables153```154155### Data Setup (Critical — MultiIndex Required for Panel Models)156157```python158import pandas as pd159160# Panel data MUST have a MultiIndex with (entity, time)161df = df.set_index(["entity_id", "year"])162163# Verify the index164print(f"Index names: {df.index.names}")165print(f"Index levels: {df.index.nlevels}")166```167168### Core Operations169170| Operation | Code |171|-----------|------|172| Panel FE (formula) | `PanelOLS.from_formula("y ~ x1 + x2 + EntityEffects", data=df).fit()` |173| Panel FE (array) | `PanelOLS(df.y, df[["x1","x2"]], entity_effects=True).fit()` |174| Two-way FE | `PanelOLS.from_formula("y ~ x1 + EntityEffects + TimeEffects", data=df).fit()` |175| Random effects | `RandomEffects.from_formula("y ~ 1 + x1 + x2", data=df).fit()` |176| Between OLS | `BetweenOLS.from_formula("y ~ 1 + x1 + x2", data=df).fit()` |177| First difference | `FirstDifferenceOLS.from_formula("y ~ x1 + x2", data=df).fit()` |178| Fama-MacBeth | `FamaMacBeth.from_formula("y ~ 1 + x1 + x2", data=df).fit()` |179| IV / 2SLS | `IV2SLS.from_formula("y ~ 1 + exog + [endog ~ inst]", data=df).fit()` |180| LIML | `IVLIML.from_formula("y ~ 1 + exog + [endog ~ inst]", data=df).fit()` |181| Clustered SEs | `mod.fit(cov_type="clustered", cluster_entity=True)` |182| Driscoll-Kraay | `mod.fit(cov_type="kernel", kernel="bartlett", bandwidth=5)` |183| Summary | `results.summary` |184| Model comparison | `compare({"FE": fe_res, "RE": re_res})` |185186### Formula Syntax187188```python189# Panel FE keywords (appear in formula, not after |)190"y ~ x1 + x2 + EntityEffects" # Entity FE191"y ~ x1 + x2 + EntityEffects + TimeEffects" # Two-way FE192"y ~ x1 + x2 + TimeEffects" # Time FE only193194# IV bracket notation195"y ~ 1 + exog + [endog ~ instrument1 + instrument2]"196197# Suppress intercept198"y ~ x1 + x2 - 1"199```200201## Topic Index202203| Topic | Reference File |204|-------|---------------|205| Installation | `./references/quickstart.md` |206| MultiIndex data setup | `./references/quickstart.md` |207| Formula vs array API | `./references/quickstart.md` |208| First model | `./references/quickstart.md` |209| Syntax comparison (pyfixest, statsmodels) | `./references/quickstart.md` |210| PanelOLS (entity/time effects) | `./references/panel-models.md` |211| RandomEffects | `./references/panel-models.md` |212| BetweenOLS | `./references/panel-models.md` |213| FirstDifferenceOLS | `./references/panel-models.md` |214| PooledOLS | `./references/panel-models.md` |215| FamaMacBeth | `./references/panel-models.md` |216| FE vs RE decision | `./references/panel-models.md` |217| Variance decomposition | `./references/panel-models.md` |218| Weighted panel estimation | `./references/panel-models.md` |219| R-squared types (within, between, overall) | `./references/panel-models.md` |220| IV2SLS | `./references/iv-models.md` |221| IVLIML and k-class estimators | `./references/iv-models.md` |222| IVGMM (1-step, 2-step, iterative) | `./references/iv-models.md` |223| IVGMMCUE | `./references/iv-models.md` |224| AbsorbingLS (high-dim FE OLS) | `./references/iv-models.md` |225| First-stage diagnostics | `./references/iv-models.md` |226| Overidentification tests | `./references/iv-models.md` |227| SUR (Seemingly Unrelated Regression) | `./references/system-models.md` |228| IV3SLS | `./references/system-models.md` |229| IVSystemGMM | `./references/system-models.md` |230| Cross-equation constraints | `./references/system-models.md` |231| LinearFactorModel | `./references/asset-pricing.md` |232| TradedFactorModel | `./references/asset-pricing.md` |233| Factor model GMM | `./references/asset-pricing.md` |234| Driscoll-Kraay SEs | `./references/covariance-inference.md` |235| Clustered SEs (entity, time, both) | `./references/covariance-inference.md` |236| HAC / kernel covariance | `./references/covariance-inference.md` |237| GMM weight matrices | `./references/covariance-inference.md` |238| Debiased inference | `./references/covariance-inference.md` |239| MultiIndex requirement | `./references/gotchas.md` |240| Maximum 2-way FE limit | `./references/gotchas.md` |241| No Panel IV | `./references/gotchas.md` |242| pyfixest vs linearmodels boundary | `./references/gotchas.md` |243| statsmodels vs linearmodels boundary | `./references/gotchas.md` |244| Constant term handling | `./references/gotchas.md` |245246## Citation247248When this library is used as a primary analytical tool, include in the report's249Software & Tools references:250251> Sheppard, K. linearmodels: Econometric models for panel data, IV/GMM, and system regression [Computer software]. https://bashtage.github.io/linearmodels/252253**Cite when:** linearmodels is used for panel estimation (RE, between), IV/GMM, Fama-MacBeth, or system regression (SUR, 3SLS).254**Do not cite when:** Only imported but no estimation performed.255256For method-specific citations (e.g., individual estimators or techniques),257consult the reference files in this skill and `agent_reference/CITATION_REFERENCE.md`.