1

*******************************************************************************************************

*******************************************************************************************************

Lin_ML_NLOGIT5_results_3.27.13.docx

Author: Abt Associates Inc.

Created: 3/27/2013

Purpose: This file presents the NLOGIT5 code and output for linear mixed logit

models estimated based on the 316(b) stated preference survey data. This file includes

results for all survey regions including the national survey version.

*******************************************************************************************************

*******************************************************************************************************
2

*******************************************************************************************************

Northeast Region

*******************************************************************************************************

|-> RESET $

|-> LOAD ; File="H:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits &

SP Study)\Task 4 - Stated Pref Survey\07 Peer Review\Post-peer review anal

Project file contained 20817 observations.

|-> CALC ; Ran(12345) $

|-> CREATE ; If(option=0) Neither=1; (Else) Neither=0 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> REJECT; region = 2 $

|-> REJECT; region = 3 $

|-> REJECT; region = 4 $

|-> REJECT; region = 5 $

|-> CREATE; newind=Ndx(id,0) $

|-> CREATE; z = GroupMean(choice, str=newind) $

+-----------------------------------------------------------------+

| Variable = ____________ Variable Groups Max Min Average |

| Z Group means CHOICE 421 9 9 9.0 |

+-----------------------------------------------------------------+

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = neither, com_fish, fish_pop, fish_sav, aquatic, cost

; RPL

; Fcn = neither(N), com_fish(N), fish_pop(N), fish_sav(N),

aquatic(N)

; Parameters

; Correlation

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

$

Normal exit: 4 iterations. Status=0, F= 1149.010

-----------------------------------------------------------------------------

Start values obtained using MNL model

Dependent variable Choice

Log likelihood function -1149.00961

Estimation based on N = 1127, K = 6

Inf.Cr.AIC = 2310.0 AIC/N = 2.050

Model estimated: Mar 27, 2013, 13:11:06

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

Constants only -1212.1328 .0521 .0432

Response data are given as ind. choices

Number of obs.= 1127, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

NEITHER| -.04758 .17615 -.27 .7871 -.39282 .29766

COM_FISH| .05078** .02287 2.22 .0264 .00596 .09561
3

FISH_POP| -.01563 .03606 -.43 .6646 -.08632 .05505

FISH_SAV| .01261*** .00126 10.00 .0000 .01014 .01509

AQUATIC| .09012** .03627 2.48 .0130 .01904 .16121

COST| -.01328*** .00230 -5.77 .0000 -.01779 -.00877

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

Line search at iteration 13 does not improve fn. Exiting optimization.

-----------------------------------------------------------------------------

Random Parameters Logit Model

Dependent variable CHOICE

Log likelihood function -946.89472

Restricted log likelihood -1238.13605

Chi squared [ 21 d.f.] 582.48265

Significance level .00000

McFadden Pseudo R-squared .2352256

Estimation based on N = 1127, K = 21

Inf.Cr.AIC = 1935.8 AIC/N = 1.718

Model estimated: Mar 27, 2013, 13:30:50

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

No coefficients -1238.1360 .2352 .2280

Constants only -1212.1328 .2188 .2115

At start values -1149.0096 .1759 .1682

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 391 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 1127, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

|Random parameters in utility functions

NEITHER| -.56919 .42972 -1.32 .1853 -1.41144 .27305

COM_FISH| .13145* .06765 1.94 .0520 -.00115 .26404

FISH_POP| .13602 .10508 1.29 .1955 -.06993 .34197

FISH_SAV| .02852*** .00537 5.31 .0000 .01799 .03905

AQUATIC| .17102 .10955 1.56 .1185 -.04368 .38573

|Nonrandom parameters in utility functions

COST| -.02677*** .00513 -5.21 .0000 -.03683 -.01670

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .37257 1.01103 .37 .7125 -1.60901 2.35415

NsCOM_FI| .33999*** .12781 2.66 .0078 .08948 .59050

NsFISH_P| .13862 .41291 .34 .7371 -.67067 .94791

NsFISH_S| .00905 .02755 .33 .7425 -.04494 .06305

NsAQUATI| .01195 .79388 .02 .9880 -1.54403 1.56793

|Below diagonal values in L matrix. V = L*Lt

COM_:NEI| -.27252 .17195 -1.58 .1130 -.60953 .06449

FISH:NEI| -.45360* .23161 -1.96 .0502 -.90754 .00034

FISH:COM| .19175 .23733 .81 .4191 -.27341 .65691

FIS0:NEI| -.03597** .01504 -2.39 .0167 -.06545 -.00650

FIS0:COM| .03279* .01948 1.68 .0922 -.00538 .07097

FISH:FIS| -.02933 .03232 -.91 .3642 -.09267 .03401
4

AQUA:NEI| -.56692* .29401 -1.93 .0538 -1.14316 .00932

AQUA:COM| -.16970 .27289 -.62 .5340 -.70456 .36516

AQUA:FIS| -.01169 .60928 -.02 .9847 -1.20586 1.18249

AQU0:FIS| -.18190 .55693 -.33 .7440 -1.27347 .90966

|Standard deviations of parameter distributions

sdNEITHE| .37257 1.01103 .37 .7125 -1.60901 2.35415

sdCOM_FI| .43573*** .12740 3.42 .0006 .18604 .68543

sdFISH_P| .51160** .21996 2.33 .0200 .08049 .94271

sdFISH_S| .05755** .02241 2.57 .0102 .01363 .10146

sdAQUATI| .61933*** .22791 2.72 .0066 .17263 1.06602

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

--------+--------------------------------------------

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

--------+--------------------------------------------

NEITHER| 1.00000 -.62544 -.88663 -.62514 -.91538

COM_FISH| -.62544 1.00000 .84697 .83564 .35871

FISH_POP| -.88663 .84697 1.00000 .62977 .70379

FISH_SAV| -.62514 .83564 .62977 1.00000 .37951

AQUATIC| -.91538 .35871 .70379 .37951 1.00000
5

*******************************************************************************************************

Southeast Region

*******************************************************************************************************

|-> RESET $

|-> LOAD ; File="H:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits &

SP Study)\Task 4 - Stated Pref Survey\07 Peer Review\Post-peer review anal

Project file contained 20817 observations.

|-> CALC ; Ran(12345) $

|-> CREATE ; If(option=0) Neither=1; (Else) Neither=0 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> REJECT; region = 1 $

|-> REJECT; region = 3 $

|-> REJECT; region = 4 $

|-> REJECT; region = 5 $

|-> CREATE; newind=Ndx(id,0) $

|-> CREATE; z = GroupMean(choice, str=newind) $

+-----------------------------------------------------------------+

| Variable = ____________ Variable Groups Max Min Average |

| Z Group means CHOICE 506 9 9 9.0 |

+-----------------------------------------------------------------+

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = neither, com_fish, fish_pop, fish_sav, aquatic, cost

; RPL

; Fcn = neither(N), com_fish(N), fish_pop(N), fish_sav(N),

aquatic(N)

; Parameters

; Correlation

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

$

Normal exit: 4 iterations. Status=0, F= 1444.082

-----------------------------------------------------------------------------

Start values obtained using MNL model

Dependent variable Choice

Log likelihood function -1444.08230

Estimation based on N = 1389, K = 6

Inf.Cr.AIC = 2900.2 AIC/N = 2.088

Model estimated: Mar 27, 2013, 13:11:44

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

Constants only -1519.7592 .0498 .0426

Response data are given as ind. choices

Number of obs.= 1389, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

NEITHER| .17407 .17028 1.02 .3067 -.15968 .50782

COM_FISH| .02977 .02095 1.42 .1552 -.01128 .07083
6

FISH_POP| .01207 .03268 .37 .7118 -.05197 .07612

FISH_SAV| .01387*** .00158 8.78 .0000 .01078 .01697

AQUATIC| .10495*** .03309 3.17 .0015 .04010 .16980

COST| -.02084*** .00213 -9.78 .0000 -.02502 -.01666

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

Line search at iteration 16 does not improve fn. Exiting optimization.

-----------------------------------------------------------------------------

Random Parameters Logit Model

Dependent variable CHOICE

Log likelihood function -1155.01295

Restricted log likelihood -1525.97247

Chi squared [ 21 d.f.] 741.91904

Significance level .00000

McFadden Pseudo R-squared .2430971

Estimation based on N = 1389, K = 21

Inf.Cr.AIC = 2352.0 AIC/N = 1.693

Model estimated: Mar 27, 2013, 13:36:40

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

No coefficients -1525.9725 .2431 .2373

Constants only -1519.7592 .2400 .2342

At start values -1444.0823 .2002 .1941

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 483 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 1389, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

|Random parameters in utility functions

NEITHER| -.19920 .30806 -.65 .5179 -.80299 .40459

COM_FISH| .12662*** .04660 2.72 .0066 .03528 .21796

FISH_POP| .14670** .06364 2.31 .0212 .02196 .27143

FISH_SAV| .02183*** .00462 4.73 .0000 .01278 .03089

AQUATIC| .24591*** .07741 3.18 .0015 .09419 .39762

|Nonrandom parameters in utility functions

COST| -.03868*** .00417 -9.27 .0000 -.04686 -.03050

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .00256 4.10600 .00 .9995 -8.04506 8.05018

NsCOM_FI| .26720*** .08472 3.15 .0016 .10115 .43325

NsFISH_P| .00053 .56151 .00 .9992 -1.10000 1.10106

NsFISH_S| .01474 .01730 .85 .3943 -.01917 .04864

NsAQUATI| .10110 .97909 .10 .9178 -1.81787 2.02008

|Below diagonal values in L matrix. V = L*Lt

COM_:NEI| .01697 .53383 .03 .9746 -1.02931 1.06325

FISH:NEI| -.02296 .76357 -.03 .9760 -1.51953 1.47361

FISH:COM| -.04206 .13328 -.32 .7523 -.30328 .21916

FIS0:NEI| -.00023 .05320 .00 .9965 -.10450 .10403

FIS0:COM| -.05434*** .01073 -5.07 .0000 -.07536 -.03332

FISH:FIS| .00870 .03932 .22 .8250 -.06837 .08576
7

AQUA:NEI| -.03153 1.02189 -.03 .9754 -2.03441 1.97134

AQUA:COM| -.08306 .16507 -.50 .6148 -.40660 .24048

AQUA:FIS| .00531 .81221 .01 .9948 -1.58659 1.59722

AQU0:FIS| -.43624 .30040 -1.45 .1465 -1.02501 .15254

|Standard deviations of parameter distributions

sdNEITHE| .00256 4.10600 .00 .9995 -8.04506 8.05018

sdCOM_FI| .26774*** .09727 2.75 .0059 .07709 .45838

sdFISH_P| .04792 .39545 .12 .9035 -.72714 .82298

sdFISH_S| .05697*** .01339 4.25 .0000 .03072 .08322

sdAQUATI| .45656*** .12343 3.70 .0002 .21463 .69849

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

--------+--------------------------------------------

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

--------+--------------------------------------------

NEITHER| 1.00000 .06338 -.47907 -.00409 -.06906

COM_FISH| .06338 1.00000 -.90630 -.95217 -.18594

FISH_POP| -.47907 -.90630 1.00000 .84082 .19289

FISH_SAV| -.00409 -.95217 .84082 1.00000 -.07155

AQUATIC| -.06906 -.18594 .19289 -.07155 1.00000
8

*******************************************************************************************************

Inland Region

*******************************************************************************************************

|-> RESET $

|-> LOAD ; File="H:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits &

SP Study)\Task 4 - Stated Pref Survey\07 Peer Review\Post-peer review anal

Project file contained 20817 observations.

|-> CALC ; Ran(12345) $

|-> CREATE ; If(option=0) Neither=1; (Else) Neither=0 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> REJECT; region = 1 $

|-> REJECT; region = 2 $

|-> REJECT; region = 4 $

|-> REJECT; region = 5 $

|-> CREATE; newind=Ndx(id,0) $

|-> CREATE; z = GroupMean(choice, str=newind) $

+-----------------------------------------------------------------+

| Variable = ____________ Variable Groups Max Min Average |

| Z Group means CHOICE 787 9 9 9.0 |

+-----------------------------------------------------------------+

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = neither, com_fish, fish_pop, fish_sav, aquatic, cost

; RPL

; Fcn = neither(N), com_fish(N), fish_pop(N), fish_sav(N),

aquatic(N)

; Parameters

; Correlation

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

$

Normal exit: 4 iterations. Status=0, F= 2312.653

-----------------------------------------------------------------------------

Start values obtained using MNL model

Dependent variable Choice

Log likelihood function -2312.65340

Estimation based on N = 2180, K = 6

Inf.Cr.AIC = 4637.3 AIC/N = 2.127

Model estimated: Mar 27, 2013, 13:12:11

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

Constants only -2375.8608 .0266 .0219

Response data are given as ind. choices

Number of obs.= 2180, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

NEITHER| .23330 .18451 1.26 .2061 -.12832 .59493

COM_FISH| .02490 .01634 1.52 .1277 -.00714 .05693
9

FISH_POP| -.01017 .02575 -.39 .6929 -.06065 .04031

FISH_SAV| .01242*** .00201 6.17 .0000 .00847 .01637

AQUATIC| .08714*** .02585 3.37 .0007 .03648 .13779

COST| -.01791*** .00166 -10.76 .0000 -.02118 -.01465

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

Line search at iteration 50 does not improve fn. Exiting optimization.

-----------------------------------------------------------------------------

Random Parameters Logit Model

Dependent variable CHOICE

Log likelihood function -1799.00579

Restricted log likelihood -2394.97479

Chi squared [ 21 d.f.] 1191.93800

Significance level .00000

McFadden Pseudo R-squared .2488415

Estimation based on N = 2180, K = 21

Inf.Cr.AIC = 3640.0 AIC/N = 1.670

Model estimated: Mar 27, 2013, 13:43:03

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

No coefficients -2394.9748 .2488 .2452

Constants only -2375.8608 .2428 .2391

At start values -2312.6534 .2221 .2183

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 756 groups

Variable number of obs./group =NUM_R

Number of obs.= 2180, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

|Random parameters in utility functions

NEITHER| -3.56478*** .58123 -6.13 .0000 -4.70396 -2.42559

COM_FISH| .06082** .02946 2.06 .0390 .00308 .11857

FISH_POP| .02175 .04405 .49 .6215 -.06458 .10808

FISH_SAV| .01483*** .00439 3.38 .0007 .00622 .02344

AQUATIC| .04630 .05036 .92 .3579 -.05240 .14500

|Nonrandom parameters in utility functions

COST| -.03356*** .00298 -11.24 .0000 -.03940 -.02771

|Diagonal values in Cholesky matrix, L.

NsNEITHE| 5.91389*** .97758 6.05 .0000 3.99787 7.82991

NsCOM_FI| .12211* .06541 1.87 .0619 -.00610 .25032

NsFISH_P| .04180 .07226 .58 .5629 -.09982 .18342

NsFISH_S| .00651 .00994 .65 .5127 -.01297 .02598

NsAQUATI| .01888 .12100 .16 .8760 -.21828 .25603

|Below diagonal values in L matrix. V = L*Lt

COM_:NEI| .00368 .05809 .06 .9495 -.11018 .11754

FISH:NEI| .06056 .07893 .77 .4430 -.09415 .21527

FISH:COM| -.00861 .09349 -.09 .9266 -.19185 .17462

FIS0:NEI| -.00280 .00815 -.34 .7314 -.01877 .01318

FIS0:COM| -.01221 .01058 -1.15 .2483 -.03295 .00852

FISH:FIS| -.04627*** .00657 -7.04 .0000 -.05915 -.03339

AQUA:NEI| -.30583*** .10392 -2.94 .0033 -.50951 -.10216
10

AQUA:COM| -.31353*** .10801 -2.90 .0037 -.52522 -.10183

AQUA:FIS| -.08142 .09646 -.84 .3986 -.27048 .10763

AQU0:FIS| .00598 .15266 .04 .9688 -.29324 .30520

|Standard deviations of parameter distributions

sdNEITHE| 5.91389*** .97758 6.05 .0000 3.99787 7.82991

sdCOM_FI| .12217* .06487 1.88 .0596 -.00497 .24931

sdFISH_P| .07409 .07769 .95 .3402 -.07817 .22635

sdFISH_S| .04838*** .00615 7.87 .0000 .03632 .06043

sdAQUATI| .44593*** .08826 5.05 .0000 .27295 .61891

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

--------+--------------------------------------------

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

--------+--------------------------------------------

NEITHER| 1.00000 .03014 .81738 -.05782 -.68583

COM_FISH| .03014 1.00000 -.09156 -.25410 -.72343

FISH_POP| .81738 -.09156 1.00000 -.55759 -.58188

FISH_SAV| -.05782 -.25410 -.55759 1.00000 .39361

AQUATIC| -.68583 -.72343 -.58188 .39361 1.00000
11

*******************************************************************************************************

Pacific Region

*******************************************************************************************************

|-> RESET $

|-> LOAD ; File="H:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits &

SP Study)\Task 4 - Stated Pref Survey\07 Peer Review\Post-peer review anal

Project file contained 20817 observations.

|-> CALC ; Ran(12345) $

|-> CREATE ; If(option=0) Neither=1; (Else) Neither=0 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> REJECT; region = 1 $

|-> REJECT; region = 2 $

|-> REJECT; region = 3 $

|-> REJECT; region = 5 $

|-> CREATE; newind=Ndx(id,0) $

|-> CREATE; z = GroupMean(choice, str=newind) $

+-----------------------------------------------------------------+

| Variable = ____________ Variable Groups Max Min Average |

| Z Group means CHOICE 311 9 9 9.0 |

+-----------------------------------------------------------------+

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = neither, com_fish, fish_pop, fish_sav, aquatic, cost

; RPL

; Fcn = neither(N), com_fish(N), fish_pop(N), fish_sav(N),

aquatic(N)

; Parameters

; Correlation

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

$

Normal exit: 5 iterations. Status=0, F= 785.6061

-----------------------------------------------------------------------------

Start values obtained using MNL model

Dependent variable Choice

Log likelihood function -785.60610

Estimation based on N = 812, K = 6

Inf.Cr.AIC = 1583.2 AIC/N = 1.950

Model estimated: Mar 27, 2013, 13:12:33

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

Constants only -867.2195 .0941 .0822

Response data are given as ind. choices

Number of obs.= 812, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

NEITHER| .25642 .21039 1.22 .2229 -.15595 .66878

COM_FISH| .03013 .02788 1.08 .2799 -.02452 .08478
12

FISH_POP| .06417 .04367 1.47 .1417 -.02142 .14977

FISH_SAV| .01730*** .00152 11.40 .0000 .01432 .02027

AQUATIC| .14388*** .04453 3.23 .0012 .05660 .23116

COST| -.01651*** .00289 -5.70 .0000 -.02218 -.01083

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

Line search at iteration 18 does not improve fn. Exiting optimization.

-----------------------------------------------------------------------------

Random Parameters Logit Model

Dependent variable CHOICE

Log likelihood function -602.70286

Restricted log likelihood -892.07318

Chi squared [ 21 d.f.] 578.74063

Significance level .00000

McFadden Pseudo R-squared .3243796

Estimation based on N = 812, K = 21

Inf.Cr.AIC = 1247.4 AIC/N = 1.536

Model estimated: Mar 27, 2013, 13:36:35

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

No coefficients -892.0732 .3244 .3155

Constants only -867.2195 .3050 .2959

At start values -785.6061 .2328 .2228

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 285 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 812, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

|Random parameters in utility functions

NEITHER| .18220 .45022 .40 .6857 -.70020 1.06461

COM_FISH| .15056** .07282 2.07 .0387 .00784 .29327

FISH_POP| .21176* .11078 1.91 .0559 -.00536 .42888

FISH_SAV| .04389*** .00698 6.29 .0000 .03022 .05756

AQUATIC| .34171*** .11007 3.10 .0019 .12597 .55744

|Nonrandom parameters in utility functions

COST| -.02756*** .00518 -5.32 .0000 -.03771 -.01741

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .09866 1.02078 .10 .9230 -1.90205 2.09936

NsCOM_FI| .05148 .54996 .09 .9254 -1.02641 1.12937

NsFISH_P| .03433 .89818 .04 .9695 -1.72606 1.79472

NsFISH_S| .03811 .03435 1.11 .2672 -.02921 .10543

NsAQUATI| .03108 1.29590 .02 .9809 -2.50885 2.57100

|Below diagonal values in L matrix. V = L*Lt

COM_:NEI| .29657 .22843 1.30 .1942 -.15114 .74428

FISH:NEI| .21922 .33024 .66 .5068 -.42803 .86647

FISH:COM| -.03956 1.30135 -.03 .9758 -2.59016 2.51105

FIS0:NEI| .03428* .01985 1.73 .0843 -.00463 .07319

FIS0:COM| .00758 .03421 .22 .8246 -.05946 .07462

FISH:FIS| .01381 .03904 .35 .7236 -.06271 .09033
13

AQUA:NEI| .28183 .33992 .83 .4070 -.38439 .94806

AQUA:COM| -.04331 1.06543 -.04 .9676 -2.13152 2.04490

AQUA:FIS| -.04100 .91836 -.04 .9644 -1.84095 1.75895

AQU0:FIS| -.21204 .43295 -.49 .6243 -1.06061 .63652

|Standard deviations of parameter distributions

sdNEITHE| .09866 1.02078 .10 .9230 -1.90205 2.09936

sdCOM_FI| .30101 .30183 1.00 .3186 -.29057 .89258

sdFISH_P| .22539 .41340 .55 .5856 -.58486 1.03564

sdFISH_S| .05362 .04330 1.24 .2156 -.03124 .13849

sdAQUATI| .35905 .57884 .62 .5351 -.77546 1.49355

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

--------+--------------------------------------------

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

--------+--------------------------------------------

NEITHER| 1.00000 .98527 .97263 .63920 .78495

COM_FISH| .98527 1.00000 .92828 .65396 .75275

FISH_POP| .97263 .92828 1.00000 .63610 .76724

FISH_SAV| .63920 .65396 .63610 1.00000 .03554

AQUATIC| .78495 .75275 .76724 .03554 1.00000
14

*******************************************************************************************************

National Survey Version

*******************************************************************************************************

|-> RESET $

|-> LOAD ; File="H:\ERD\TYCHE Contract\Assignments\WA 2-15 (316b Benefits &

SP Study)\Task 4 - Stated Pref Survey\07 Peer Review\Post-peer review anal

Project file contained 20817 observations.

|-> CALC ; Ran(12345) $

|-> CREATE ; If(option=0) Neither=1; (Else) Neither=0 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> REJECT; region = 1 $

|-> REJECT; region = 2 $

|-> REJECT; region = 3 $

|-> REJECT; region = 4 $

|-> CREATE; newind=Ndx(id,0) $

|-> CREATE; z = GroupMean(choice, str=newind) $

+-----------------------------------------------------------------+

| Variable = ____________ Variable Groups Max Min Average |

| Z Group means CHOICE 288 9 9 9.0 |

+-----------------------------------------------------------------+

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = neither, com_fish, fish_pop, fish_sav, aquatic, cost

; RPL

; Fcn = neither(N), com_fish(N), fish_pop(N), fish_sav(N),

aquatic(N)

; Parameters

; Correlation

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

$

Normal exit: 4 iterations. Status=0, F= 823.1674

-----------------------------------------------------------------------------

Start values obtained using MNL model

Dependent variable Choice

Log likelihood function -823.16744

Estimation based on N = 789, K = 6

Inf.Cr.AIC = 1658.3 AIC/N = 2.102

Model estimated: Mar 27, 2013, 13:13:04

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

Constants only -858.1651 .0408 .0278

Response data are given as ind. choices

Number of obs.= 789, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

NEITHER| .01976 .21987 .09 .9284 -.41117 .45069

COM_FISH| .05102* .02732 1.87 .0618 -.00253 .10456
15

FISH_POP| .01815 .04282 .42 .6717 -.06577 .10206

FISH_SAV| .01288*** .00194 6.64 .0000 .00908 .01668

AQUATIC| .00398 .04297 .09 .9262 -.08024 .08821

COST| -.01597*** .00274 -5.82 .0000 -.02134 -.01059

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

Line search at iteration 18 does not improve fn. Exiting optimization.

-----------------------------------------------------------------------------

Random Parameters Logit Model

Dependent variable CHOICE

Log likelihood function -660.03931

Restricted log likelihood -866.80510

Chi squared [ 21 d.f.] 413.53157

Significance level .00000

McFadden Pseudo R-squared .2385378

Estimation based on N = 789, K = 21

Inf.Cr.AIC = 1362.1 AIC/N = 1.726

Model estimated: Mar 27, 2013, 13:24:01

R2=1-LogL/LogL* Log-L fncn R-sqrd R2Adj

No coefficients -866.8051 .2385 .2283

Constants only -858.1651 .2309 .2205

At start values -823.1674 .1982 .1874

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 272 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 789, skipped 0 obs

--------+--------------------------------------------------------------------

| Standard Prob. 95% Confidence

CHOICE| Coefficient Error z |z|>Z* Interval

--------+--------------------------------------------------------------------

|Random parameters in utility functions

NEITHER| -2.09328*** .45875 -4.56 .0000 -2.99242 -1.19414

COM_FISH| .08925 .06077 1.47 .1419 -.02986 .20836

FISH_POP| .02838 .10735 .26 .7915 -.18202 .23879

FISH_SAV| .02241*** .00585 3.83 .0001 .01093 .03388

AQUATIC| -.12105 .11085 -1.09 .2748 -.33831 .09621

|Nonrandom parameters in utility functions

COST| -.03739*** .00520 -7.19 .0000 -.04758 -.02719

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .18221 2.51713 .07 .9423 -4.75128 5.11570

NsCOM_FI| .25789 .23968 1.08 .2819 -.21186 .72765

NsFISH_P| .42440 .70213 .60 .5456 -.95176 1.80055

NsFISH_S| .04638*** .01622 2.86 .0043 .01458 .07818

NsAQUATI| .09169 1.05598 .09 .9308 -1.97800 2.16138

|Below diagonal values in L matrix. V = L*Lt

COM_:NEI| .08478 .67488 .13 .9000 -1.23796 1.40752

FISH:NEI| .24647 1.59495 .15 .8772 -2.87958 3.37252

FISH:COM| -.23771 .41405 -.57 .5659 -1.04922 .57381

FIS0:NEI| .00720 .05298 .14 .8919 -.09665 .11104

FIS0:COM| -.02694 .02251 -1.20 .2315 -.07107 .01719

FISH:FIS| -.01058 .03803 -.28 .7808 -.08513 .06396
16

AQUA:NEI| .17645 1.11942 .16 .8748 -2.01757 2.37047

AQUA:COM| .01383 .45501 .03 .9757 -.87798 .90564

AQUA:FIS| -.08445 .69582 -.12 .9034 -1.44823 1.27934

AQU0:FIS| -.78306** .32141 -2.44 .0148 -1.41302 -.15311

|Standard deviations of parameter distributions

sdNEITHE| .18221 2.51713 .07 .9423 -4.75128 5.11570

sdCOM_FI| .27147 .41960 .65 .5176 -.55092 1.09387

sdFISH_P| .54531 1.33621 .41 .6832 -2.07362 3.16424

sdFISH_S| .05514*** .01773 3.11 .0019 .02040 .08989

sdAQUATI| .81244 .53981 1.51 .1323 -.24557 1.87045

--------+--------------------------------------------------------------------

Note: ***, **, * ==> Significance at 1%, 5%, 10% level.

-----------------------------------------------------------------------------

--------+--------------------------------------------

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

--------+--------------------------------------------

NEITHER| 1.00000 .31229 .45198 .13052 .21719

COM_FISH| .31229 1.00000 -.27296 -.42336 .08400

FISH_POP| .45198 -.27296 1.00000 .12261 .00985

FISH_SAV| .13052 -.42336 .12261 1.00000 -.77071

AQUATIC| .21719 .08400 .00985 -.77071 1.00000