1

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

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

Stepwise_ML_NLOGIT5_results_7.01.13.docx

Author: Abt Associates Inc.

Created: 7/1/13

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

based on the 316(b) stated preference survey data with stepwise fish saved. 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(999) $

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

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

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

|-> REJECT; region = 2 $

|-> REJECT; region = 5 $

|-> 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 421 9 9 9.0 |

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

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> REJECT; w_var = -999 $

|-> CREATE ; fish_sq = fish_sav*fish_sav $

|-> CREATE ; If(fish_sav = 5) fish_5=1; (Else) fish_5=0 $

|-> CREATE ; If(fish_sav = 50) fish_50=1; (Else) fish_50=0 $

|-> CREATE ; If(fish_sav = 95) fish_95=1; (Else) fish_95=0 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = com_fish, fish_pop, fish_5, fish_50, fish_95, aquatic,

cost

; RPL

; Fcn = com_fish(N), fish_pop(N), aquatic(N)

; Parameters

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -1090.81361

Estimation based on N = 1078, K = 7

Inf.Cr.AIC = 2195.6 AIC/N = 2.037

Model estimated: Jul 01, 2013, 18:24:18

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

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$
3

Response data are given as ind. choices

Number of obs.= 1078, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

COM_FISH| .05254** .02319 2.27 .0235 .00708 .09800

FISH_POP| .01001 .03654 .27 .7841 -.06160 .08162

AQUATIC| .09350** .03697 2.53 .0114 .02104 .16595

FISH_5| .02297 .18708 .12 .9023 -.34370 .38965

FISH_50| .88462*** .17833 4.96 .0000 .53510 1.23414

FISH_95| 1.04920*** .18399 5.70 .0000 .68859 1.40981

COST| -.01320*** .00236 -5.60 .0000 -.01782 -.00858

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -968.74364

Restricted log likelihood -1175.88015

Chi squared [ 10 d.f.] 414.27303

Significance level .00000

McFadden Pseudo R-squared .1761544

Estimation based on N = 1078, K = 10

Inf.Cr.AIC = 1957.5 AIC/N = 1.816

Model estimated: Jul 01, 2013, 18:39:36

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -1090.8136 .1119******

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 372 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 1078, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

|Random parameters in utility functions

COM_FISH| .03229 .04034 .80 .4234 -.04677 .11135

FISH_POP| .12996 .08591 1.51 .1303 -.03842 .29834

AQUATIC| .21497** .09542 2.25 .0243 .02795 .40199

|Nonrandom parameters in utility functions

FISH_5| -.24201 .27171 -.89 .3731 -.77456 .29053

FISH_50| 1.09973*** .28561 3.85 .0001 .53994 1.65952

FISH_95| 1.39847*** .30151 4.64 .0000 .80752 1.98942

COST| -.01173*** .00332 -3.53 .0004 -.01824 -.00521

|Distns. of RPs. Std.Devs or limits of triangular

NsCOM_FI| .21162** .09479 2.23 .0256 .02582 .39741

NsFISH_P| .92959*** .14490 6.42 .0000 .64559 1.21358
4

NsAQUATI| .91559*** .14981 6.11 .0000 .62196 1.20922

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

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

-----------------------------------------------------------------------------
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(123) $

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

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

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

|-> REJECT; region = 1 $

|-> REJECT; region = 5 $

|-> 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 506 9 9 9.0 |

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

|-> CREATE; num_r=z*9 $

|-> REJECT; num_r = 0 $

|-> REJECT; choice = -999 $

|-> REJECT; response = -999 $

|-> REJECT; w_var = -999 $

|-> CREATE ; If(fish_sav = 25) fish_25=1; (Else) fish_25=0 $

|-> CREATE ; If(fish_sav = 55) fish_55=1; (Else) fish_55=0 $

|-> CREATE ; If(fish_sav = 90) fish_90=1; (Else) fish_90=0 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = com_fish, fish_pop, fish_25, fish_55, fish_90, aquatic,

cost

; RPL

; Fcn = com_fish(N), fish_pop(N), aquatic(N)

; Parameters

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -1365.48362

Estimation based on N = 1321, K = 7

Inf.Cr.AIC = 2745.0 AIC/N = 2.078

Model estimated: Jul 01, 2013, 18:26:35

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

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$
6

Response data are given as ind. choices

Number of obs.= 1321, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

COM_FISH| .03486 .02147 1.62 .1045 -.00722 .07695

FISH_POP| -.02291 .03346 -.68 .4936 -.08850 .04268

AQUATIC| .10604*** .03406 3.11 .0018 .03930 .17279

FISH_25| .32491** .16357 1.99 .0470 .00432 .64549

FISH_55| .64920*** .16639 3.90 .0001 .32308 .97532

FISH_90| 1.17608*** .16954 6.94 .0000 .84379 1.50836

COST| -.02008*** .00218 -9.22 .0000 -.02435 -.01581

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -1206.67087

Restricted log likelihood -1441.96170

Chi squared [ 10 d.f.] 470.58166

Significance level .00000

McFadden Pseudo R-squared .1631741

Estimation based on N = 1321, K = 10

Inf.Cr.AIC = 2433.3 AIC/N = 1.842

Model estimated: Jul 01, 2013, 18:42:48

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -1365.4836 .1163******

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 458 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 1321, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

|Random parameters in utility functions

COM_FISH| .12915*** .04140 3.12 .0018 .04801 .21029

FISH_POP| .04549 .07020 .65 .5170 -.09210 .18309

AQUATIC| .30748*** .07066 4.35 .0000 .16899 .44597

|Nonrandom parameters in utility functions

FISH_25| .31780 .21999 1.44 .1486 -.11337 .74897

FISH_55| .68790*** .23475 2.93 .0034 .22781 1.14800

FISH_90| 1.36035*** .25588 5.32 .0000 .85882 1.86187

COST| -.02828*** .00341 -8.29 .0000 -.03496 -.02159

|Distns. of RPs. Std.Devs or limits of triangular

NsCOM_FI| .40460*** .07425 5.45 .0000 .25906 .55013

NsFISH_P| .71690*** .12040 5.95 .0000 .48092 .95288
7

NsAQUATI| .67675*** .11427 5.92 .0000 .45279 .90072

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

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

-----------------------------------------------------------------------------
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(123) $

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

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

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

|-> REJECT; region = 1 $

|-> REJECT; region = 5 $

|-> REJECT; region = 2 $

|-> REJECT; region = 4 $

|-> 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 $

|-> REJECT; w_var = -999 $

|-> CREATE ; If(fish_sav = 55) fish_55=1; (Else) fish_55=0 $

|-> CREATE ; If(fish_sav = 75) fish_75=1; (Else) fish_75=0 $

|-> CREATE ; If(fish_sav = 95) fish_95=1; (Else) fish_95=0 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = com_fish, fish_pop, fish_55, fish_75, fish_95, aquatic,

cost

; RPL

; Fcn = com_fish(N), fish_pop(N), aquatic(N)

; Parameters

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -2281.35196

Estimation based on N = 2180, K = 7

Inf.Cr.AIC = 4576.7 AIC/N = 2.099

Model estimated: Jul 01, 2013, 18:28:59

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

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

Response data are given as ind. choices
9

Number of obs.= 2180, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

COM_FISH| .00760 .01647 .46 .6445 -.02468 .03989

FISH_POP| .00185 .02574 .07 .9428 -.04860 .05230

AQUATIC| .07116*** .02591 2.75 .0060 .02039 .12194

FISH_55| .57841*** .12481 4.63 .0000 .33379 .82303

FISH_75| .90525*** .12754 7.10 .0000 .65527 1.15522

FISH_95| 1.13624*** .13112 8.67 .0000 .87925 1.39322

COST| -.01938*** .00167 -11.58 .0000 -.02266 -.01610

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -2058.51630

Restricted log likelihood -2380.09089

Chi squared [ 10 d.f.] 643.14917

Significance level .00000

McFadden Pseudo R-squared .1351102

Estimation based on N = 2180, K = 10

Inf.Cr.AIC = 4137.0 AIC/N = 1.898

Model estimated: Jul 01, 2013, 19:14:17

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -2281.3520 .0977******

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

COM_FISH| -.03519 .03234 -1.09 .2765 -.09858 .02820

FISH_POP| -.01603 .06383 -.25 .8017 -.14113 .10906

AQUATIC| .26696*** .04418 6.04 .0000 .18037 .35356

|Nonrandom parameters in utility functions

FISH_55| .59433*** .18494 3.21 .0013 .23186 .95681

FISH_75| 1.10540*** .19512 5.67 .0000 .72296 1.48783

FISH_95| 1.38457*** .20258 6.83 .0000 .98752 1.78162

COST| -.02593*** .00262 -9.90 .0000 -.03107 -.02079

|Distns. of RPs. Std.Devs or limits of triangular

NsCOM_FI| .37020*** .05488 6.75 .0000 .26264 .47776

NsFISH_P| 1.02313*** .10819 9.46 .0000 .81108 1.23518

NsAQUATI| .49541*** .08636 5.74 .0000 .32614 .66469

--------+--------------------------------------------------------------------
10

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

-----------------------------------------------------------------------------
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(123) $

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

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

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

|-> REJECT; region = 1 $

|-> REJECT; region = 5 $

|-> REJECT; region = 2 $

|-> REJECT; region = 3 $

|-> 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 $

|-> REJECT; w_var = -999 $

|-> CREATE ; If(fish_sav = 2) fish_2=1; (Else) fish_2=0 $

|-> CREATE ; If(fish_sav = 50) fish_50=1; (Else) fish_50=0 $

|-> CREATE ; If(fish_sav = 95) fish_95=1; (Else) fish_95=0 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = com_fish, fish_pop, fish_50, fish_95, aquatic, cost

; RPL

; Fcn = com_fish(N), fish_pop(N), aquatic(N)

; Parameters

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -682.11512

Estimation based on N = 720, K = 6

Inf.Cr.AIC = 1376.2 AIC/N = 1.911

Model estimated: Jul 09, 2013, 18:01:02

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

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

Response data are given as ind. choices

Number of obs.= 720, skipped 0 obs
12

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

| Standard Prob. 95% Confidence

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

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

COM_FISH| .03070 .02802 1.10 .2733 -.02423 .08563

FISH_POP| .05548 .04219 1.31 .1886 -.02722 .13818

AQUATIC| .15072*** .04282 3.52 .0004 .06679 .23465

FISH_50| .80059*** .14191 5.64 .0000 .52245 1.07873

FISH_95| 1.51057*** .14324 10.55 .0000 1.22982 1.79132

COST| -.01804*** .00300 -6.02 .0000 -.02391 -.01217

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -602.34800

Restricted log likelihood -779.04062

Chi squared [ 9 d.f.] 353.38523

Significance level .00000

McFadden Pseudo R-squared .2268080

Estimation based on N = 720, K = 9

Inf.Cr.AIC = 1222.7 AIC/N = 1.698

Model estimated: Jul 09, 2013, 18:08:41

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -682.1151 .1169******

Response data are given as ind. choices

Replications for simulated probs. = 300

Halton sequences used for simulations

RPL model with panel has 252 groups

Variable number of obs./group =NUM_R

BHHH estimator used for asymp. variance

Number of obs.= 720, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

|Random parameters in utility functions

COM_FISH| .04594 .06607 .70 .4869 -.08356 .17544

FISH_POP| .16759 .12439 1.35 .1779 -.07620 .41138

AQUATIC| .48854*** .08097 6.03 .0000 .32984 .64724

|Nonrandom parameters in utility functions

FISH_50| .92662*** .22718 4.08 .0000 .48136 1.37187

FISH_95| 1.82906*** .25662 7.13 .0000 1.32610 2.33202

COST| -.02383*** .00490 -4.86 .0000 -.03344 -.01422

|Distns. of RPs. Std.Devs or limits of triangular

NsCOM_FI| .45357*** .12460 3.64 .0003 .20935 .69778

NsFISH_P| 1.22225*** .17914 6.82 .0000 .87114 1.57336

NsAQUATI| .26818 .19488 1.38 .1688 -.11378 .65013

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

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

-----------------------------------------------------------------------------
13

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

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(123) $

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

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

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

|-> REJECT; region = 1 $

|-> REJECT; region = 4 $

|-> REJECT; region = 2 $

|-> REJECT; region = 3 $

|-> 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 $

|-> REJECT; w_var = -999 $

|-> CREATE ; If(fish_sav = 25) fish_25=1; (Else) fish_25=0 $

|-> CREATE ; If(fish_sav = 55) fish_55=1; (Else) fish_55=0 $

|-> CREATE ; If(fish_sav = 95) fish_95=1; (Else) fish_95=0 $

|-> NLOGIT ; Lhs = Choice

; Choices = A,B,Neither

; Rhs = com_fish, fish_pop, fish_25, fish_55, fish_95, aquatic,

cost

; RPL

; Fcn = com_fish(N), fish_pop(N), aquatic(N)

; Parameters

; Pds = num_r

; Pts = 300

; Maxit = 1000

; Halton

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -814.46699

Estimation based on N = 789, K = 7

Inf.Cr.AIC = 1642.9 AIC/N = 2.082

Model estimated: Jul 01, 2013, 18:33:27

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

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

Response data are given as ind. choices
14

Number of obs.= 789, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

COM_FISH| .05297* .02781 1.90 .0568 -.00154 .10748

FISH_POP| .01231 .04301 .29 .7747 -.07198 .09660

AQUATIC| .01304 .04334 .30 .7635 -.07191 .09798

FISH_25| .39314* .20589 1.91 .0562 -.01040 .79668

FISH_55| .66466*** .21173 3.14 .0017 .24968 1.07965

FISH_95| 1.32261*** .21608 6.12 .0000 .89910 1.74612

COST| -.01766*** .00277 -6.37 .0000 -.02310 -.01223

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -717.39916

Restricted log likelihood -864.47725

Chi squared [ 10 d.f.] 294.15618

Significance level .00000

McFadden Pseudo R-squared .1701353

Estimation based on N = 789, K = 10

Inf.Cr.AIC = 1454.8 AIC/N = 1.844

Model estimated: Jul 01, 2013, 18:39:00

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -814.4670 .1192******

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

Number of obs.= 789, skipped 0 obs

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

| Standard Prob. 95% Confidence

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

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

|Random parameters in utility functions

COM_FISH| .20843*** .07404 2.82 .0049 .06332 .35355

FISH_POP| .04687 .09498 .49 .6217 -.13929 .23303

AQUATIC| .03409 .12791 .27 .7898 -.21660 .28478

|Nonrandom parameters in utility functions

FISH_25| .32266 .47315 .68 .4953 -.60469 1.25001

FISH_55| .79403 .54173 1.47 .1427 -.26774 1.85581

FISH_95| 1.71456** .69948 2.45 .0142 .34359 3.08552

COST| -.02624*** .00799 -3.29 .0010 -.04189 -.01059

|Distns. of RPs. Std.Devs or limits of triangular

NsCOM_FI| .25697* .13935 1.84 .0652 -.01615 .53008

NsFISH_P| 1.01673** .41534 2.45 .0144 .20268 1.83078

NsAQUATI| .79238* .48003 1.65 .0988 -.14846 1.73323

--------+--------------------------------------------------------------------
15

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

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