[Limdep Nlogit List] MNL in WTP space
William Greene
wgreene at stern.nyu.edu
Fri May 7 00:38:08 AEST 2021
Thao. When there are no random coefficients in the model, and tau=gamma=0,
then the resulting model is a simple MNL. Specifying this as an GMXL model
makes
the program act peculiarly as the GMXL is written specifically for the case
when those
are nonzero. It's hard to predict how it will behave. The right way to
fit that model is
as an MNL. But, it is also important that the MNL model in WTP space is a
one to one
transformation of the same model in preference space. That means that you
get the
identical answers in the two cases. This is a theoretical property of the
MLE known as
"invariance." The MLE is invariant to 1:1 transformations of the parameter
space.
/Bill Greene
On Thu, May 6, 2021 at 6:05 AM Thao Thai via Limdep <
limdep at mailman.sydney.edu.au> wrote:
> 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
>
> -----------------------------------------------------------------------------
> _______________________________________________
> Limdep site list
> Limdep at mailman.sydney.edu.au
> https://protect-au.mimecast.com/s/_P-bClx1Nji2AD49zTGHxbu?domain=limdep.itls.usyd.edu.au
>
>
--
William Greene
Department of Economics, emeritus
Stern School of Business, New York University
44 West 4 St.
New York, NY, 10012
URL: https://protect-au.mimecast.com/s/8pTZCmO5glujAG40RFOi6g5?domain=people.stern.nyu.edu
Email: wgreene at stern.nyu.edu
Editor in Chief: Journal of Productivity Analysis
<https://protect-au.mimecast.com/s/d2uTCnx1jni7KQBy6UN015w?domain=springer.com>
<https://protect-au.mimecast.com/s/d2uTCnx1jni7KQBy6UN015w?domain=springer.com>
Editor in Chief: Foundations and Trends in Econometrics
Associate Editor: Economics Letters
Associate Editor: Journal of Business and Economic Statistics
More information about the Limdep
mailing list