Announcement

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

  • Likelihood ratio test for exact poisson regression

    Good afternoon,

    I am trying to assess model fit for multivariate regression analysis using the liklihood ratio test on exact poisson regression models (using exact poisson due to a small sample size).

    See code below:
    expoisson Outcome Exposure1 Expousre2 Exposure3 Exposure4, irr
    est store m1
    expoisson Outcome Exposure1 Expousre2 Exposure3 Exposure4 Exposure5 , irr
    est store m2
    lrtest m1 m2

    I get the following error messages:
    (m1 does not contain matrix e(V); rank = 0 assumed)
    m1 does not contain scalar e(ll)

    Is there an alternative method for exact models?

    Many thanks

  • #2
    exact models, as used here, are not maximum likelihood models and thus a likelihood ratio test is not possible

    Comment


    • #3
      Thanks Richard. So would I have to use another model or are there other tests for models such as these? I wondered if the AIC (estat ic) could be used here. Many thanks
      Last edited by Mona Dave; 02 Aug 2022, 10:30.

      Comment


      • #4
        AIC, etc., are also based on likelihoods and thus us also not appropriate

        the issue is in what sense do you mean "model fit"? i.e., what exactly are you worried about? remember that, in a sense, this is an exact estimator

        Comment


        • #5
          The test you want to do is equivalent to the test that the coefficient of exposure 5 equals 0 in model 2.
          ​​
          ---------------------------------
          Maarten L. Buis
          University of Konstanz
          Department of history and sociology
          box 40
          78457 Konstanz
          Germany
          http://www.maartenbuis.nl
          ---------------------------------

          Comment


          • #6
            In reply to Richard, I'm trying to do forward stepwise regression. What I've been taught is to use a test based on liklihood to arrive at the final model in multivariate analysis. However, I've not used exact regression before which is where I've come a bit stuck with methodology. That's useful to know that the other tests based on liklihood won't work here.
            Last edited by Mona Dave; 02 Aug 2022, 13:15.

            Comment


            • #7
              what is your N? you appear to have either 4 or 5 predictors/covariates (or whatever you want to call them) and my worry right now is that your N is too small for this many predictors

              Comment


              • #8
                N=23

                Comment

                Working...
                X