Dealing with Separation in Logistic Regression Models
Carlisle Rainey Assistant Professor Texas A&M University crainey@tamu.edu
paper, data, and code at crain.co/research
Dealing with Separation in Logistic Regression Models Carlisle - - PowerPoint PPT Presentation
Dealing with Separation in Logistic Regression Models Carlisle Rainey Assistant Professor Texas A&M University crainey@tamu.edu paper, data, and code at crain.co/research The prior matters a lot, so choose a good one. The prior matters
Dealing with Separation in Logistic Regression Models
Carlisle Rainey Assistant Professor Texas A&M University crainey@tamu.edu
paper, data, and code at crain.co/research
Variable Coefficient Confidence Interval Democratic Governor
[-126,979.03; 126,926.33] % Uninsured (Std.) 0.92 [-3.46; 5.30] % Favorable to ACA 0.01 [-0.17; 0.18] GOP Legislature 2.43 [-0.47; 5.33] Fiscal Health 0.00 [-0.02; 0.02] Medicaid Multiplier
[-2.45; 1.80] % Non-white 0.05 [-0.12; 0.21] % Metropolitan
[-0.17; 0.02] Constant 2.58 [-7.02; 12.18]
Variable Coefficient Confidence Interval Democratic Governor
[-126,979.03; 126,926.33] % Uninsured (Std.) 0.92 [-3.46; 5.30] % Favorable to ACA 0.01 [-0.17; 0.18] GOP Legislature 2.43 [-0.47; 5.33] Fiscal Health 0.00 [-0.02; 0.02] Medicaid Multiplier
[-2.45; 1.80] % Non-white 0.05 [-0.12; 0.21] % Metropolitan
[-0.17; 0.02] Constant 2.58 [-7.02; 12.18]
Variable Coefficient Confidence Interval Democratic Governor
[-126,979.03; 126,926.33] % Uninsured (Std.) 0.92 [-3.46; 5.30] % Favorable to ACA 0.01 [-0.17; 0.18] GOP Legislature 2.43 [-0.47; 5.33] Fiscal Health 0.00 [-0.02; 0.02] Medicaid Multiplier
[-2.45; 1.80] % Non-white 0.05 [-0.12; 0.21] % Metropolitan
[-0.17; 0.02] Constant 2.58 [-7.02; 12.18]
This is a failure of maximum likelihood.
For
the posterior distribution of βs is proportional to the prior distribution for βs, so that p(βs|y) ∝ p(βs|σ).
For
the posterior distribution of βs is proportional to the prior distribution for βs, so that p(βs|y) ∝ p(βs|σ).
for Choosing a Good Prior
Pr(yi) = Λ(βc + βssi + β1xi1 + ... + βkxik)
Transforming the Prior Distribution
We Already Know Few Things
β1 ≈ ˆ βmle
1
β2 ≈ ˆ βmle
2
. . . βk ≈ ˆ βmle
k
p∗(β|βs < 0, β−s = ˆ βmle
−s ),
where ˆ βmle
s
= −∞
(on GitHub)
and choose a good one.
range of prior distributions.