2
$\begingroup$

Some major commercial statistical packages (e.g., SPSS) do not report a CI for the intercept term in logistic regression. [Based on answer below R does provide CI for intercept]

Why might confidence intervals for the intercept term not be included by default?

UPDATE:

Based on feedback, Confidence Intervals for odds associated with intercept term are reported in some stat. packages. And, obviously, they can be easily computed manually, knowing the standard error. Therefore, reformulated question:

What is the interpretation of intercept's CIs (for odds) in logistic regression?

Improve this question
$\endgroup$
5
  • 1
    $\begingroup$ Although I can find plenty of annotated outputs of SPSS binary logistic regression on the Web, none of them directly provide confidence intervals, but they all report standard errors (from which CIs are easily computed), including for the constant term. Do you think you could document the behavior you are asking about by means of a reproducible example? $\endgroup$
    – whuber ♦
    Commented May 14, 2014 at 19:01
  • $\begingroup$ SPSS is good for behavioral sciences, where they normally don't care about those C.I that much ;) $\endgroup$
    – Stat
    Commented May 14, 2014 at 19:06
  • $\begingroup$ I agree with both above comments. But still I found it a bit suspicious - please see my comments to gung's answer below $\endgroup$
    – Oleg Shirokikh
    Commented May 14, 2014 at 19:07
  • $\begingroup$ Use Generalized Linear Model, you'll get 95% CI that way. But I agree with another answer that 95%CI in terms of Odds Ratio for the intercept is pretty hard to interpret. $\endgroup$
    – Penguin_Knight
    Commented May 14, 2014 at 19:24
  • $\begingroup$ Is there an option to run the model without the intercept in SPSS? (I know there is for linear regression, not sure about logistic). If so, create a variable in SPSS which is equal to 1 for everyone, put that in as a predictor, and ask to remove the intercept. The new variable will be the intercept and will have CIs. $\endgroup$
    – Jeremy Miles
    Commented May 14, 2014 at 20:09

2 Answers 2

Reset to default
4
$\begingroup$

Well, that's not correct. You can use confint funtion in both S-plus and R to obtain C.I. for the estimated parameters. I will give an example for logistic regression in R since I don't have the S-plus:

> mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv")
> mylogit <- glm(admit ~ gre + gpa + rank, data = mydata, family = "binomial")
> confint(mylogit)
Waiting for profiling to be done...
                    2.5 %       97.5 %
(Intercept) -5.7109591680 -1.260314066
gre          0.0001715446  0.004461385
gpa          0.1415710585  1.428341503
rank        -0.8149612229 -0.315479733
>
Improve this answer
$\endgroup$
1
  • $\begingroup$ Interesting. I haven't checked R.. E.g. SPSS doesn't provide CIs for intercept. Please see updated question $\endgroup$
    – Oleg Shirokikh
    Commented May 14, 2014 at 18:48
2
$\begingroup$

I don't have access to SPSS, but I strongly suspect that it can output a confidence interval for the intercept in logistic regression. You probably have to know how to use the underlying syntax (e.g., PASTE) to call for it.

As to why it isn't output by default, you'd really have to contact the company. But I would guess that they believe people aren't very interested in seeing the CI for the intercept and want to minimize the volume of statistical output.

Improve this answer
$\endgroup$
4
  • $\begingroup$ Thanks. I've checked documentation and all UI options and I'm pretty sure that there is no option to display CI for intercepts. They provide a table with coefficients and statistics including CIs for predictors - but leave space for CI for intercept empty. For instance, ats.ucla.edu/stat/spss/output/logistic.htm - last table. $\endgroup$
    – Oleg Shirokikh
    Commented May 14, 2014 at 19:05
  • $\begingroup$ Therefore I became curious - maybe this can be justified by some reasoning (though I don't see any) $\endgroup$
    – Oleg Shirokikh
    Commented May 14, 2014 at 19:06
  • $\begingroup$ I don't know SPSS's syntax, but I believe it can be done the UCLA page notwithstanding. You can also compute it yourself (remember that the SEs will be on the scale of the linear predictor) fairly easily. $\endgroup$
    – gung - Reinstate Monica
    Commented May 14, 2014 at 19:11
  • 5
    $\begingroup$ At the related page ats.ucla.edu/stat/stata/faq/oratio.htm (about 2/3 the way down) the SPSS output includes a CI for the constant by default. This is because no odds ratio calculations were requested. I'm not surprised that CIs are not provided for an odds ratio for a constant: after all, what would that even mean? Thus it would behoove us all to distinguish carefully between output in terms of the raw parameter estimates and output in terms of odds ratios. $\endgroup$
    – whuber ♦
    Commented May 14, 2014 at 19:23

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.