Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • logistic regression standard errors

    Hi All,

    When I run the logit command I get the following result

    outcome | Coef. Std. Err. z P>|z| [95% Conf. Interval]
    --------------------------------+----------------------------------------------------------------
    indep_var | 3.147702 1.095392 2.87 0.004 1.000772 5.294632


    When I run the same command with the or option I get the following output

    outcome | Odds Ratio Std. Err. z P>|z| [95% Conf. Interval]
    --------------------------------+----------------------------------------------------------------
    indep_var| 23.2825 25.50347 2.87 0.004 2.720382 199.2642


    I understand how the confidence interval is calculated from the log scale results - for example -

    exp(3.147702 + 1.96*1.095392) approx.= exp(5.294632) = 199.2642


    How does the numerical scale Std Err of 25.50347 relate to the log scale Std Err of 1.095392 and how would it be used to calculate the confidence interval on the numerical scale?

    Thanks in advance,

    Don


  • #2
    See https://www.stata.com/support/faqs/s...cs/delta-rule/
    -------------------------------------------
    Richard Williams, Notre Dame Dept of Sociology
    StataNow Version: 19.5 MP (2 processor)

    EMAIL: [email protected]
    WWW: https://www3.nd.edu/~rwilliam

    Comment


    • #3
      Thanks Richard.

      Using the link you gave, I can see that, in my case Variance[exp(b)] = Variance(b)*exp(b) = 1.095392*23.2825 = 25.50347.

      I don't quite understand how exp(b)*Variance(b) relates to the Taylor expansion for the variance

      Variance[exp(b)] = Variance(b)*[exp'(b)]2 , where exp' means the derivative, d(exp(b))/d(b).

      I understand that the derivative of the exponential is itself so I don't understand why we multiply by exp(b) instead of [exp'(b)]2 ?

      Thanks again in advance,

      Don

      Comment


      • #4
        Sorry, I can't help you. This is one of those things I take on blind faith.
        -------------------------------------------
        Richard Williams, Notre Dame Dept of Sociology
        StataNow Version: 19.5 MP (2 processor)

        EMAIL: [email protected]
        WWW: https://www3.nd.edu/~rwilliam

        Comment


        • #5
          Hi Don
          I think you are forgetting the following:
          exp'(b)=exp(b)
          and indeed:
          Variance[exp(b)] = Variance(b)*[exp(b)]2
          But what the "table" shows is the standard error (square root of the above expressions)
          STDDev(exp(b))=STDDEV(b)*[exp(b)]

          HTH
          Last edited by FernandoRios; 03 Oct 2019, 21:23.

          Comment


          • #6
            Of course! Thanks very much Fernando and thanks again Richard for taking the time to consider this.


            Don

            Comment

            Working...
            X