Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

re: re:st: Hausman test


From   Ricardo Ovaldia <[email protected]>
To   [email protected]
Subject   re: re:st: Hausman test
Date   Tue, 3 Jan 2012 11:17:43 -0800 (PST)

Thank you again  Kit.

Yes I used sigmamore. I also get the message if I use -sigmaless- or nothing at all.

Here is the output, is the Hausman test valid?:

. xtset fam
       panel variable:  fam (unbalanced)

. 
.         xtreg chol l1 l2 l3  i.sex i.marker, fe

Fixed-effects (within) regression               Number of obs      =       416
Group variable: fam                             Number of groups   =         9

R-sq:  within  = 0.1368                         Obs per group: min =        27
       between = 0.7881                                        avg =      46.2
       overall = 0.5725                                        max =        73

                                                F(9,398)           =      7.01
corr(u_i, Xb)  = 0.3937                         Prob > F           =    0.0000

------------------------------------------------------------------------------
        chol |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          l1 |   4.064901   1.649737     2.46   0.014     .8216118    7.308189
          l2 |  -.1133366   .0582016    -1.95   0.052    -.2277575    .0010844
          l3 |   .0012481   .0006428     1.94   0.053    -.0000155    .0025118
       2.sex |   3.799992   .8166467     4.65   0.000     2.194512    5.405472
             |
      marker |
          2  |  -.9134657   1.173933    -0.78   0.437    -3.221351    1.394419
          3  |  -1.346317    1.20109    -1.12   0.263    -3.707591    1.014957
          4  |   .5833598   2.476976     0.24   0.814    -4.286232    5.452952
          5  |     .53107   1.352041     0.39   0.695    -2.126964    3.189104
          6  |   2.271173   3.366584     0.67   0.500    -4.347337    8.889682
             |
       _cons |   9.321154   14.70398     0.63   0.526    -19.58602    38.22833
-------------+----------------------------------------------------------------
     sigma_u |  5.0827273
     sigma_e |  8.1937869
         rho |  .27786908   (fraction of variance due to u_i)
------------------------------------------------------------------------------
F test that all u_i=0:     F(8, 398) =     9.45              Prob > F = 0.0000

.         est store fixed

.         xtreg chol l1 l2 l3  i.sex i.marker  

Random-effects GLS regression                   Number of obs      =       416
Group variable: fam                             Number of groups   =         9

R-sq:  within  = 0.1185                         Obs per group: min =        27
       between = 0.8693                                        avg =      46.2
       overall = 0.6005                                        max =        73

                                                Wald chi2(9)       =    610.15
corr(u_i, X)   = 0 (assumed)                    Prob > chi2        =    0.0000

------------------------------------------------------------------------------
        chol |      Coef.   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
          l1 |   6.341376   1.317354     4.81   0.000      3.75941    8.923342
          l2 |   -.170823   .0489787    -3.49   0.000    -.2668195   -.0748265
          l3 |   .0016893    .000569     2.97   0.003      .000574    .0028046
       2.sex |    3.77881   .8776635     4.31   0.000     2.058621    5.498999
             |
      marker |
          2  |  -1.579155   1.217692    -1.30   0.195    -3.965788    .8074783
          3  |  -1.128749    1.26016    -0.90   0.370    -3.598617     1.34112
          4  |   1.295496   2.350917     0.55   0.582    -3.312217    5.903209
          5  |  -.8996998   1.409237    -0.64   0.523    -3.661753    1.862353
          6  |   11.06828   3.289858     3.36   0.001     4.620282    17.51629
             |
       _cons |  -16.71301    11.0348    -1.51   0.130    -38.34082    4.914797
-------------+----------------------------------------------------------------
     sigma_u |          0
     sigma_e |  8.1937869
         rho |          0   (fraction of variance due to u_i)
------------------------------------------------------------------------------

.         noi hausman fixed ., sigmamore

Note: the rank of the differenced variance matrix (8) does not equal the number of coefficients being tested (9); be sure this is what you expect, or there may be problems computing
        the test.  Examine the output of your estimators for anything unexpected and possibly consider scaling your variables so that the coefficients are on a similar scale.

                 ---- Coefficients ----
             |      (b)          (B)            (b-B)     sqrt(diag(V_b-V_B))
             |     fixed          .          Difference          S.E.
-------------+----------------------------------------------------------------
          l1 |    4.064901     6.341376       -2.276476         1.19984
          l2 |   -.1133366     -.170823        .0574864         .039406
          l3 |    .0012481     .0016893       -.0004412        .0003977
       2.sex |    3.799992      3.77881        .0211814        .0878687
  2bn.marker |   -.9134657    -1.579155        .6656893         .353452
    3.marker |   -1.346317    -1.128749       -.2175685        .3081298
    4.marker |    .5833598     1.295496       -.7121362        1.276993
    5.marker |      .53107    -.8996998         1.43077        .3828838
    6.marker |    2.271173     11.06828       -8.797112        1.548825
------------------------------------------------------------------------------
                           b = consistent under Ho and Ha; obtained from xtreg
            B = inconsistent under Ha, efficient under Ho; obtained from xtreg

    Test:  Ho:  difference in coefficients not systematic

                  chi2(8) = (b-B)'[(V_b-V_B)^(-1)](b-B)
                          =       64.84
                Prob>chi2 =      0.0000




Ricardo Ovaldia, MS
Statistician 
Oklahoma City, OK


--- On Tue, 1/3/12, Christopher Baum <[email protected]> wrote:

> From: Christopher Baum <[email protected]>
> Subject: re: re:st: Hausman test
> To: "[email protected]" <[email protected]>
> Date: Tuesday, January 3, 2012, 12:39 PM
> <>
> "the rank of the differenced variance matrix (8) does not
> equal the number of coefficients being tested (9)..."
> 
> The two -xtreg- model are the same. Only difference is the
> -fe- option in the first model. I get estimates for all the
> covariates in both models so I am not sure why I am getting
> this message. Any ideas?
> 
> 
> Without seeing the output from both models and from the
> test, hard to tell. But note that xtreg, re can include
> time-invariant variables (such as race or gender)
> while xtreg,fe will drop them. That's probably not the
> case, as you say you get coefficients for all variables, but
> there has to be some reason that the 
> differenced VCE is rank-deficient. You are using the
> sigmamore or sigmaless option on hausman, per the help file
> suggestion, right?
> 
> KIt
> 
> 
> Kit Baum   |   Boston College
> Economics & DIW
> Berlin   |   http://ideas.repec.org/e/pba1.html
>                
>              An
> Introduction to Stata Programming 
> |   http://www.stata-press.com/books/isp.html
>   An Introduction to Modern Econometrics Using
> Stata  |   http://www.stata-press.com/books/imeus.html
> 
> 
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
> 

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index