[Limdep Nlogit List] MNL in WTP space

Thao Thai Thao.T.Thai at monash.edu
Thu May 6 20:03:59 AEST 2021


Dear Nlogit users,

I am estimating a MNL in WTP-space using the GMXL command in Nlogit

(1) If I specify all of my 20 coefficients as constant, the model does not
work. I specified ;tau = [0] as I think it would result in an MNL without
scale.

;fcn=  hos(c), pri(c), ind(c), gov(c), non(c), rlh1(c), rlh2(c), fl(c),
cr1(c), cr2(c), lo1(c), sa(*c), rlc1(c), rlc2(c), lo2(c), rlp1(c), rli1(c),
rli2(c), rlg1(c), rln1(c)
*;tau = [0]*

(2) If I specify some of my variables as non-random (which means I do not
specify them in ;fcn), the model works. However, the non-random variables
have coefficients in preference space (which means their coefficients are
the same as those in MNL model in preference space)

  ;fcn=  rlh1(c), rlh2(c), fl(c), cr1(c), cr2(c), lo1(c), sa(*c), rlc1(c),
rlc2(c), lo2(c), rlp1(c), rli1(c), rli2(c), rlg1(c), rln1(c)
*;tau = [0]*

(3) When I used the same specification as in (2) but with ;gamma = [0],
Nlogit reports "Error    805: Initial iterations cannot improve
function.Status=3"
 ;fcn=  rlh1(c), rlh2(c), fl(c), cr1(c), cr2(c), lo1(c), sa(*c), rlc1(c),
rlc2(c), lo2(c), rlp1(c), rli1(c), rli2(c), rlg1(c), rln1(c)
*;gamma = [0]*

My questions are:
1. Why does Nlogit produce coefficients in preference space for non-random
variables? Are there any ways to get an MNL in WTP space for all
coefficients in Nlogit ?
2. Please advise why the model with ;gamma = [0] can not run?
My detailed syntax is below.

Any suggestions you can share would be much appreciated. Thank you!
Best regards,
Thao

*(1) This did not work: *

|-> CALC ; Ran(29531928) $
|-> sample;all$
|-> reject;ALLCHO#30$
|-> reject;FCPC=1$
|-> GMXlogit;userp
    ;lhs = cho, cset, alts
    ;choices = a,b,c,d,e,f
    ;checkdata
    ;fcn= hos(c), pri(c), ind(c), gov(c), non(c), rlh1(c), rlh2(c), fl(c),
cr1(c), cr2(c), lo1(c), sa(*c), rlc1(c), rlc2(c), lo2(c), rlp1(c), rli1(c),
rli2(c), rlg1(c), rln1(c)
    ;pts=100
    ;pds=3
    ;halton
    ;maxit=300
    ;gmx
    ;par
    ;tau= [0]
    ;model:
    U(a) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(b) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(c) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(d) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(e) =hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(f) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    $
+----------------------------------------------------------+
| Inspecting the data set before estimation.               |
| These errors mark observations which will be skipped.    |
| Row Individual = 1st row then group number of data block |
+----------------------------------------------------------+
No bad observations were found in the sample

Iterative procedure has converged
Normal exit:   5 iterations. Status=0, F=    .3903337D+04

-----------------------------------------------------------------------------
Start values obtained using MNL model
Dependent variable               Choice
Log likelihood function     -3903.33652
Estimation based on N =   2370, K =  20
Inf.Cr.AIC  =   7846.7 AIC/N =    3.311
---------------------------------------
            Log likelihood R-sqrd R2Adj
ASCs  only  model must be fit separately
               Use NLOGIT ;...;RHS=ONE$
Note: R-sqrd = 1 - logL/Logl(constants)
Warning:  Model does not contain a full
set of ASCs. R-sqrd is problematic. Use
model setup with ;RHS=one to get LogL0.
---------------------------------------
Response data are given as ind. choices
Number of obs.=  2370, skipped    0 obs
--------+--------------------------------------------------------------------
        |                  Standard            Prob.      95% Confidence
     CHO|  Coefficient       Error       z    |z|>Z*         Interval
--------+--------------------------------------------------------------------
     HOS|     .11397         .14452      .79  .4303     -.16929    .39722
     PRI|     .34043***      .12631     2.70  .0070      .09286    .58800
     IND|    -.94307***      .16462    -5.73  .0000    -1.26571   -.62042
     GOV|     .08533         .13264      .64  .5200     -.17464    .34529
     NON|    -.08990         .14032     -.64  .5217     -.36492    .18512
    RLH1|     .08068         .15904      .51  .6119     -.23102    .39239
    RLH2|     .19219         .13233     1.45  .1464     -.06717    .45155
      FL|     .20004***      .05424     3.69  .0002      .09373    .30634
     CR1|     .43243***      .05697     7.59  .0000      .32077    .54410
     CR2|     .15373*        .08971     1.71  .0866     -.02210    .32957
     LO1|    -.56437***      .05770    -9.78  .0000     -.67746   -.45129
      SA|     .01263***      .00080    15.86  .0000      .01107    .01419
    RLC1|     .39288**       .16670     2.36  .0184      .06614    .71961
    RLC2|     .30506**       .15166     2.01  .0443      .00781    .60231
     LO2|    -.86848***      .08495   -10.22  .0000    -1.03499   -.70198
    RLP1|    -.01482         .12190     -.12  .9033     -.25373    .22410
    RLI1|     .65116***      .14499     4.49  .0000      .36699    .93534
    RLI2|     .77620***      .15085     5.15  .0000      .48054   1.07187
    RLG1|    -.34006***      .12687    -2.68  .0074     -.58873   -.09140
    RLN1|    -.08453         .13534     -.62  .5322     -.34980    .18073
--------+--------------------------------------------------------------------
***, **, * ==>  Significance at 1%, 5%, 10% level.
Model was estimated on May 06, 2021 at 03:04:54 PM
-----------------------------------------------------------------------------

Error    805: Initial iterations cannot improve function.Status=3
Function value was   .4246469942D+04 at entry -----------
                     .2835448956D+05 at exit  -----------
Error   1025: Failed to fit model. See earlier diagnostic.
*(2) This syntax worked and its results *

|-> CALC ; Ran(29531928) $
|-> sample;all$
|-> reject;ALLCHO#30$
|-> reject;FCPC=1$
|-> GMXlogit;userp
    ;lhs = cho, cset, alts
    ;choices = a,b,c,d,e,f
    ;checkdata
    ;fcn= rlh1(c), rlh2(c), fl(c), cr1(c), cr2(c), lo1(c), sa(*c), rlc1(c),
rlc2(c), lo2(c), rlp1(c), rli1(c), rli2(c), rlg1(c), rln1(c)
    ;pts=100
    ;pds=3
    ;halton
    ;maxit=300
    ;gmx
    ;par
    ;tau= [0]
    ;model:
    U(a) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(b) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(c) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(d) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(e) =hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    /
    U(f) = hos*HOS + pri*PRI + ind*IND + gov*GOV + non*NON
    + rlh1    * RL_H1  + rlh2  * RL_H2
    + fl      * FL_H
    + cr1     * CR_H1  + cr2   * CR_H2
    + lo1     * LO_H1
    + sa      * SA_H
    + rlc1    * RL_C1  + rlc2  * RL_C2
    + fl      * FL_C
    + cr1     * CR_C1  + cr2   * CR_C2
    + lo1    * LO_C1   + lo2   * LO_C2
    + sa      * SA_C
    + rlp1    * RL_P1
    + fl      * FL_P
    + cr1     * CR_P1  + cr2   * CR_P2
    + lo1     * LO_P1  + lo2   * LO_P2
    + sa      * SA_P
    + rli1    * RL_I1  + rli2  * RL_I2
    + fl      * FL_I
    + cr1     * CR_I1
    + lo1     * LO_I1
    + sa      * SA_I
    + rlg1    * RL_G1
    + fl      * FL_G
    + cr1     * CR_G1
    + lo1    * LO_G1
    + sa      * SA_G
    + rln1    * RL_N1
    + fl      * FL_N
    + cr1     * CR_N1
    + lo1     * LO_N1  + lo2   * LO_N2
    + sa      * SA_N
    $
+----------------------------------------------------------+
| Inspecting the data set before estimation.               |
| These errors mark observations which will be skipped.    |
| Row Individual = 1st row then group number of data block |
+----------------------------------------------------------+
No bad observations were found in the sample

Iterative procedure has converged
Normal exit:   5 iterations. Status=0, F=    .3903337D+04

-----------------------------------------------------------------------------
Start values obtained using MNL model
Dependent variable               Choice
Log likelihood function     -3903.33652
Estimation based on N =   2370, K =  20
Inf.Cr.AIC  =   7846.7 AIC/N =    3.311
---------------------------------------
            Log likelihood R-sqrd R2Adj
ASCs  only  model must be fit separately
               Use NLOGIT ;...;RHS=ONE$
Note: R-sqrd = 1 - logL/Logl(constants)
Warning:  Model does not contain a full
set of ASCs. R-sqrd is problematic. Use
model setup with ;RHS=one to get LogL0.
---------------------------------------
Response data are given as ind. choices
Number of obs.=  2370, skipped    0 obs
--------+--------------------------------------------------------------------
        |                  Standard            Prob.      95% Confidence
     CHO|  Coefficient       Error       z    |z|>Z*         Interval
--------+--------------------------------------------------------------------
    RLH1|     .08068         .15904      .51  .6119     -.23102    .39239
    RLH2|     .19219         .13233     1.45  .1464     -.06717    .45155
      FL|     .20004***      .05424     3.69  .0002      .09373    .30634
     CR1|     .43243***      .05697     7.59  .0000      .32077    .54410
     CR2|     .15373*        .08971     1.71  .0866     -.02210    .32957
     LO1|    -.56437***      .05770    -9.78  .0000     -.67746   -.45129
      SA|     .01263***      .00080    15.86  .0000      .01107    .01419
    RLC1|     .39288**       .16670     2.36  .0184      .06614    .71961
    RLC2|     .30506**       .15166     2.01  .0443      .00781    .60231
     LO2|    -.86848***      .08495   -10.22  .0000    -1.03499   -.70198
    RLP1|    -.01482         .12190     -.12  .9033     -.25373    .22410
    RLI1|     .65116***      .14499     4.49  .0000      .36699    .93534
    RLI2|     .77620***      .15085     5.15  .0000      .48054   1.07187
    RLG1|    -.34006***      .12687    -2.68  .0074     -.58873   -.09140
    RLN1|    -.08453         .13534     -.62  .5322     -.34980    .18073
     HOS|     .11397         .14452      .79  .4303     -.16929    .39722
     PRI|     .34043***      .12631     2.70  .0070      .09286    .58800
     IND|    -.94307***      .16462    -5.73  .0000    -1.26571   -.62042
     GOV|     .08533         .13264      .64  .5200     -.17464    .34529
     NON|    -.08990         .14032     -.64  .5217     -.36492    .18512
--------+--------------------------------------------------------------------
***, **, * ==>  Significance at 1%, 5%, 10% level.
Model was estimated on May 06, 2021 at 02:26:45 PM
-----------------------------------------------------------------------------

Iterative procedure has converged
Normal exit:  48 iterations. Status=0, F=    .3903337D+04

-----------------------------------------------------------------------------
Generalized Mixed (RP) Logit Model
Dependent variable                  CHO
Log likelihood function     -3903.33652
Restricted log likelihood   -4246.46994
Chi squared [ 20](P= .000)    686.26684
Significance level               .00000
McFadden Pseudo R-squared      .0808044
Estimation based on N =   2370, K =  20
Inf.Cr.AIC  =   7846.7 AIC/N =    3.311
---------------------------------------
            Log likelihood R-sqrd R2Adj
No coefficients -4246.4699  .0808 .0793
Constants only can be computed directly
               Use NLOGIT ;...;RHS=ONE$
At start values **********  .6987 .6982
Note: R-sqrd = 1 - logL/Logl(constants)
Warning:  Model does not contain a full
set of ASCs. R-sqrd is problematic. Use
model setup with ;RHS=one to get LogL0.
---------------------------------------
Response data are given as ind. choices
Replications for simulated probs. = 100
Used Halton sequences in simulations.
RPL model with panel has     790 groups
Fixed number of obsrvs./group=        3
Number of obs.=  2370, skipped    0 obs
--------+--------------------------------------------------------------------
        |                  Standard            Prob.      95% Confidence
     CHO|  Coefficient       Error       z    |z|>Z*         Interval
--------+--------------------------------------------------------------------
        |Random parameters in utility
functions..............................
    RLH1|    6.38639       12.59107      .51  .6120   -18.29166  31.06444
    RLH2|    15.2124       10.50848     1.45  .1477     -5.3839   35.8086
      FL|    15.8337***     4.24480     3.73  .0002      7.5141   24.1534
     CR1|    34.2284***     4.99186     6.86  .0000     24.4446   44.0123
     CR2|    12.1685*       7.13125     1.71  .0879     -1.8085   26.1455
     LO1|   -44.6720***     5.35949    -8.34  .0000    -55.1764  -34.1676
      SA|        1.0    .....(Fixed Parameter).....
    RLC1|    31.0976**     13.94301     2.23  .0257      3.7698   58.4254
    RLC2|    24.1465*      12.37301     1.95  .0510      -.1042   48.3971
     LO2|   -68.7436***     8.40285    -8.18  .0000    -85.2129  -52.2744
    RLP1|   -1.17281        9.63099     -.12  .9031   -20.04920  17.70357
    RLI1|    51.5419***    11.86105     4.35  .0000     28.2946   74.7891
    RLI2|    61.4391***    11.65835     5.27  .0000     38.5891   84.2890
    RLG1|   -26.9173***     9.86829    -2.73  .0064    -46.2588   -7.5758
    RLN1|   -6.69115       10.70810     -.62  .5321   -27.67865  14.29634
        |Nonrandom parameters in utility
functions...........................
     HOS|     .11397         .14452      .79  .4303     -.16929    .39722
     PRI|     .34043***      .12631     2.70  .0070      .09286    .58800
     IND|    -.94307***      .16462    -5.73  .0000    -1.26571   -.62042
     GOV|     .08533         .13264      .64  .5200     -.17464    .34529
     NON|    -.08990         .14032     -.64  .5217     -.36492    .18512
        |Distns. of RPs. Std.Devs or limits of
triangular....................
  CsRLH1|        0.0    .....(Fixed Parameter).....
  CsRLH2|        0.0    .....(Fixed Parameter).....
    CsFL|        0.0    .....(Fixed Parameter).....
   CsCR1|        0.0    .....(Fixed Parameter).....
   CsCR2|        0.0    .....(Fixed Parameter).....
   CsLO1|        0.0    .....(Fixed Parameter).....
    CsSA|        0.0    .....(Fixed Parameter).....
  CsRLC1|        0.0    .....(Fixed Parameter).....
  CsRLC2|        0.0    .....(Fixed Parameter).....
   CsLO2|        0.0    .....(Fixed Parameter).....
  CsRLP1|        0.0    .....(Fixed Parameter).....
  CsRLI1|        0.0    .....(Fixed Parameter).....
  CsRLI2|        0.0    .....(Fixed Parameter).....
  CsRLG1|        0.0    .....(Fixed Parameter).....
  CsRLN1|        0.0    .....(Fixed Parameter).....
        |Variance parameter tau in GMX scale
parameter.......................
TauScale|        0.0    .....(Fixed Parameter).....
        |Weighting parameter gamma in GMX
model..............................
GammaMXL|        0.0    .....(Fixed Parameter).....
        |Coefficient on SA       in preference space
form....................
Beta0WTP|     .01263***      .00080    15.86  .0000      .01107    .01419
S_b0_WTP|        0.0    .....(Fixed Parameter).....
        |  Sample Mean    Sample
Std.Dev.....................................
Sigma(i)|        1.0***   .2153D-05 ********  .0000  .10000D+01  .10000D+01
--------+--------------------------------------------------------------------
nnnnn.D-xx or D+xx => multiply by 10 to -xx or +xx.
***, **, * ==>  Significance at 1%, 5%, 10% level.
Fixed parameter ... is constrained to equal the value or
had a nonpositive st.error because of an earlier problem.
Model was estimated on May 06, 2021 at 02:57:53 PM
-----------------------------------------------------------------------------


More information about the Limdep mailing list