1

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

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

Wgt_Lin_ML_NLOGIT5_results_11.14.12.docx

Author: Robert Johnston, Clark University

Created: 11/14/2012

Purpose: This file presents the NLOGIT5 code and output for weighted 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="C:\Rob Johnston\Applied Research and Development\Wetland\SP

and BT 2009\Data Analysis\all_reg_model_data_w_wgts_071112.lpj"$

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

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

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -1097.15198

Estimation based on N = 1078, K = 6

Inf.Cr.AIC = 2206.3 AIC/N = 2.047

Model estimated: Nov 14, 2012, 11:30:14

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.= 1078, skipped 0 obs

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

| Standard Prob. 95% Confidence
3

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

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

NEITHER| -.13870 .18107 -.77 .4437 -.49359 .21619

COM_FISH| .04855** .02308 2.10 .0354 .00331 .09380

FISH_POP| .00189 .03633 .05 .9585 -.06931 .07309

FISH_SAV| .01110*** .00127 8.72 .0000 .00861 .01359

AQUATIC| .08278** .03665 2.26 .0239 .01095 .15460

COST| -.01284*** .00234 -5.48 .0000 -.01743 -.00825

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -906.58721

Restricted log likelihood -1175.88015

Chi squared [ 21 d.f.] 538.58588

Significance level .00000

McFadden Pseudo R-squared .2290139

Estimation based on N = 1078, K = 21

Inf.Cr.AIC = 1855.2 AIC/N = 1.721

Model estimated: Nov 14, 2012, 11:37:33

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -1097.1520 .1737******

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

NEITHER| -.52259 .37095 -1.41 .1589 -1.24963 .20445

COM_FISH| .21316*** .05180 4.12 .0000 .11163 .31468

FISH_POP| .06502 .09369 .69 .4877 -.11862 .24866

FISH_SAV| .02900*** .00464 6.26 .0000 .01992 .03809

AQUATIC| .20135* .10306 1.95 .0507 -.00065 .40335

|Nonrandom parameters in utility functions

COST| -.02106*** .00453 -4.65 .0000 -.02995 -.01218

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .19465 1.01704 .19 .8482 -1.79872 2.18801

NsCOM_FI| .02436 .18777 .13 .8968 -.34366 .39237

NsFISH_P| .51866** .20836 2.49 .0128 .11029 .92703

NsFISH_S| .01194 .03594 .33 .7397 -.05850 .08237

NsAQUATI| .16464 .78656 .21 .8342 -1.37700 1.70628

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

COM_:NEI| -.13415 .16318 -.82 .4110 -.45399 .18568

FISH:NEI| -.26629 .30726 -.87 .3861 -.86851 .33592
4

FISH:COM| -.08380 .30845 -.27 .7859 -.68835 .52076

FIS0:NEI| -.01506 .02085 -.72 .4702 -.05593 .02581

FIS0:COM| .00782 .01851 .42 .6726 -.02846 .04410

FISH:FIS| -.03852** .01714 -2.25 .0246 -.07211 -.00493

AQUA:NEI| -.34493 .55957 -.62 .5376 -1.44168 .75181

AQUA:COM| -.49055 .38042 -1.29 .1972 -1.23616 .25507

AQUA:FIS| -.17135 .35684 -.48 .6311 -.87074 .52804

AQU0:FIS| -.07816 .86316 -.09 .9279 -1.76992 1.61360

|Standard deviations of parameter distributions

sdNEITHE| .19465 1.01704 .19 .8482 -1.79872 2.18801

sdCOM_FI| .13635 .16831 .81 .4179 -.19353 .46622

sdFISH_P| .58902*** .18116 3.25 .0011 .23395 .94409

sdFISH_S| .04375*** .01269 3.45 .0006 .01887 .06863

sdAQUATI| .64976 .44299 1.47 .1424 -.21848 1.51800

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

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

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

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

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

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

NEITHER| 1.00000 -.98392 -.45210 -.34422 -.53086

COM_FISH| -.98392 1.00000 .41941 .37062 .38746

FISH_POP| -.45210 .41941 1.00000 -.64504 .11520

FISH_SAV| -.34422 .37062 -.64504 1.00000 .24709

AQUATIC| -.53086 .38746 .11520 .24709 1.00000
5

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

Southeast Region

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

|-> RESET $

|-> LOAD ; File="D:\Rob Johnston\Applied Research and Development\Wetland\SP

and BT 2009\Data Analysis\all_reg_model_data_w_wgts_071112.lpj"$

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

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

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -1365.77901

Estimation based on N = 1321, K = 6

Inf.Cr.AIC = 2743.6 AIC/N = 2.077

Model estimated: Nov 14, 2012, 11:37:35

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.= 1321, skipped 0 obs

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

| Standard Prob. 95% Confidence
6

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

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

NEITHER| .03059 .17514 .17 .8614 -.31268 .37385

COM_FISH| .03566* .02144 1.66 .0962 -.00636 .07769

FISH_POP| -.02358 .03345 -.71 .4808 -.08914 .04197

FISH_SAV| .01320*** .00162 8.15 .0000 .01003 .01637

AQUATIC| .10656*** .03405 3.13 .0018 .03983 .17330

COST| -.02013*** .00218 -9.25 .0000 -.02439 -.01586

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -1091.87599

Restricted log likelihood -1441.96170

Chi squared [ 21 d.f.] 700.17142

Significance level .00000

McFadden Pseudo R-squared .2427843

Estimation based on N = 1321, K = 21

Inf.Cr.AIC = 2225.8 AIC/N = 1.685

Model estimated: Nov 14, 2012, 11:54:37

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -1365.7790 .2005******

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

NEITHER| -1.63728*** .33244 -4.92 .0000 -2.28886 -.98570

COM_FISH| .08871* .04991 1.78 .0755 -.00912 .18654

FISH_POP| -.02800 .07830 -.36 .7206 -.18147 .12546

FISH_SAV| .02596*** .00507 5.12 .0000 .01603 .03589

AQUATIC| .06163 .08284 .74 .4569 -.10074 .22400

|Nonrandom parameters in utility functions

COST| -.04209*** .00502 -8.39 .0000 -.05193 -.03226

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .29835 1.20020 .25 .8037 -2.05400 2.65069

NsCOM_FI| .14488 .12894 1.12 .2612 -.10785 .39761

NsFISH_P| .34534 .21084 1.64 .1014 -.06790 .75858

NsFISH_S| .03134 .02085 1.50 .1328 -.00952 .07221

NsAQUATI| .36594 .73244 .50 .6173 -1.06962 1.80150

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

COM_:NEI| -.09935 .17194 -.58 .5634 -.43634 .23765

FISH:NEI| -.02963 .33947 -.09 .9305 -.69498 .63573
7

FISH:COM| -.04728 .20356 -.23 .8163 -.44625 .35169

FIS0:NEI| .01206 .02925 .41 .6800 -.04526 .06939

FIS0:COM| -.04960** .02069 -2.40 .0165 -.09015 -.00904

FISH:FIS| .02064 .01917 1.08 .2817 -.01694 .05821

AQUA:NEI| .40595 .45000 .90 .3670 -.47604 1.28794

AQUA:COM| -.42749 .45389 -.94 .3463 -1.31711 .46212

AQUA:FIS| .30738 .40058 .77 .4429 -.47773 1.09249

AQU0:FIS| -.14265 .47791 -.30 .7653 -1.07934 .79405

|Standard deviations of parameter distributions

sdNEITHE| .29835 1.20020 .25 .8037 -2.05400 2.65069

sdCOM_FI| .17567 .12699 1.38 .1666 -.07323 .42457

sdFISH_P| .34982 .21724 1.61 .1073 -.07596 .77559

sdFISH_S| .06335** .02852 2.22 .0263 .00746 .11924

sdAQUATI| .77220 .60628 1.27 .2028 -.41608 1.96048

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

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

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

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

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

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

NEITHER| 1.00000 -.56553 -.08469 .19042 .52571

COM_FISH| -.56553 1.00000 -.06357 -.75333 -.75388

FISH_POP| -.08469 -.06357 1.00000 .41124 .42326

FISH_SAV| .19042 -.75333 .41124 1.00000 .57177

AQUATIC| .52571 -.75388 .42326 .57177 1.00000
8

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

Inland Region

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

|-> RESET $

|-> LOAD ; File="D:\Rob Johnston\Applied Research and Development\Wetland\SP

and BT 2009\Data Analysis\all_reg_model_data_w_wgts_071112.lpj"$

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

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

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -2281.59847

Estimation based on N = 2180, K = 6

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

Model estimated: Nov 14, 2012, 11:54:52

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.= 2180, skipped 0 obs

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

| Standard Prob. 95% Confidence
9

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

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

NEITHER| .17051 .18556 .92 .3581 -.19318 .53421

COM_FISH| .00698 .01645 .42 .6713 -.02526 .03922

FISH_POP| .00240 .02572 .09 .9256 -.04801 .05282

FISH_SAV| .01392*** .00204 6.82 .0000 .00992 .01792

AQUATIC| .07093*** .02590 2.74 .0062 .02017 .12169

COST| -.01935*** .00167 -11.57 .0000 -.02263 -.01608

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

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

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

Normal exit: 59 iterations. Status=0, F= 1820.480

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -1820.47987

Restricted log likelihood -2380.09089

Chi squared [ 21 d.f.] 1119.22205

Significance level .00000

McFadden Pseudo R-squared .2351217

Estimation based on N = 2180, K = 21

Inf.Cr.AIC = 3683.0 AIC/N = 1.689

Model estimated: Nov 14, 2012, 12:17:47

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -2281.5985 .2021******

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.53149*** .56418 -6.26 .0000 -4.63727 -2.42572

COM_FISH| .02335 .02773 .84 .3997 -.03100 .07770

FISH_POP| .01007 .04182 .24 .8097 -.07189 .09203

FISH_SAV| .01597*** .00393 4.06 .0000 .00826 .02368

AQUATIC| .04864 .04478 1.09 .2774 -.03913 .13641

|Nonrandom parameters in utility functions

COST| -.03280*** .00284 -11.55 .0000 -.03836 -.02723

|Diagonal values in Cholesky matrix, L.

NsNEITHE| 5.18984*** .89828 5.78 .0000 3.42925 6.95043

NsCOM_FI| .09955 .06150 1.62 .1055 -.02099 .22009

NsFISH_P| .03658 .13236 .28 .7823 -.22285 .29601

NsFISH_S| .01864** .00753 2.48 .0133 .00388 .03339

NsAQUATI| .06182 .12751 .48 .6278 -.18810 .31174

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

COM_:NEI| -.07666 .05058 -1.52 .1296 -.17579 .02247

FISH:NEI| .03644 .07329 .50 .6190 -.10721 .18009

FISH:COM| -.11299 .10825 -1.04 .2966 -.32516 .09917
10

FIS0:NEI| -.00369 .00839 -.44 .6599 -.02013 .01275

FIS0:COM| -.03288*** .00781 -4.21 .0000 -.04818 -.01757

FISH:FIS| -.02142** .00975 -2.20 .0281 -.04053 -.00231

AQUA:NEI| -.31014*** .08659 -3.58 .0003 -.47985 -.14043

AQUA:COM| .02081 .10435 .20 .8419 -.18371 .22534

AQUA:FIS| -.11070 .12644 -.88 .3813 -.35853 .13712

AQU0:FIS| -.04576 .13452 -.34 .7338 -.30942 .21791

|Standard deviations of parameter distributions

sdNEITHE| 5.18984*** .89828 5.78 .0000 3.42925 6.95043

sdCOM_FI| .12565** .06115 2.05 .0399 .00580 .24549

sdFISH_P| .12423 .08723 1.42 .1544 -.04674 .29520

sdFISH_S| .04359*** .00554 7.86 .0000 .03273 .05446

sdAQUATI| .33881*** .08224 4.12 .0000 .17762 .50000

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

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

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

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

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

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

NEITHER| 1.00000 -.61013 .29332 -.08464 -.91539

COM_FISH| -.61013 1.00000 -.89959 -.54587 .60718

FISH_POP| .29332 -.89959 1.00000 .51644 -.42059

FISH_SAV| -.08464 -.54587 .51644 1.00000 .13394

AQUATIC| -.91539 .60718 -.42059 .13394 1.00000
11

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

Pacific Region

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

|-> RESET $

|-> LOAD ; File="D:\Rob Johnston\Applied Research and Development\Wetland\SP

and BT 2009\Data Analysis\all_reg_model_data_w_wgts_071112.lpj"$

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

|-> REJECT; region = 3 $

|-> REJECT; region = 2 $

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

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

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -681.72007

Estimation based on N = 720, K = 6

Inf.Cr.AIC = 1375.4 AIC/N = 1.910

Model estimated: Nov 14, 2012, 12:17:49

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

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

| Standard Prob. 95% Confidence
12

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

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

NEITHER| .23678 .22627 1.05 .2954 -.20672 .68027

COM_FISH| .04028 .03004 1.34 .1799 -.01860 .09917

FISH_POP| .07344 .04629 1.59 .1126 -.01728 .16416

FISH_SAV| .01667*** .00162 10.32 .0000 .01351 .01984

AQUATIC| .17065*** .04766 3.58 .0003 .07723 .26407

COST| -.01727*** .00307 -5.62 .0000 -.02329 -.01124

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -537.50676

Restricted log likelihood -779.04062

Chi squared [ 21 d.f.] 483.06772

Significance level .00000

McFadden Pseudo R-squared .3100401

Estimation based on N = 720, K = 21

Inf.Cr.AIC = 1117.0 AIC/N = 1.551

Model estimated: Nov 14, 2012, 12:23:37

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -681.7201 .2115******

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

NEITHER| .18797 .49719 .38 .7054 -.78651 1.16244

COM_FISH| .11072 .08606 1.29 .1982 -.05794 .27939

FISH_POP| .16039 .12739 1.26 .2080 -.08929 .41007

FISH_SAV| .03645*** .00683 5.34 .0000 .02307 .04982

AQUATIC| .31963** .13338 2.40 .0166 .05820 .58106

|Nonrandom parameters in utility functions

COST| -.02177*** .00534 -4.07 .0000 -.03224 -.01129

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .01702 6.00963 .00 .9977 -11.76163 11.79568

NsCOM_FI| .31874 .23203 1.37 .1695 -.13602 .77351

NsFISH_P| .01653 .38841 .04 .9660 -.74473 .77780

NsFISH_S| .02457 .02761 .89 .3736 -.02955 .07868

NsAQUATI| .00153 1.38260 .00 .9991 -2.70831 2.71138

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

COM_:NEI| -.04717 1.03635 -.05 .9637 -2.07838 1.98403

FISH:NEI| -.03537 1.65838 -.02 .9830 -3.28574 3.21500
13

FISH:COM| -.23252 .33193 -.70 .4836 -.88309 .41804

FIS0:NEI| -.00836 .05617 -.15 .8817 -.11845 .10174

FIS0:COM| -.04305*** .01621 -2.66 .0079 -.07483 -.01127

FISH:FIS| .01778 .03789 .47 .6388 -.05648 .09205

AQUA:NEI| -.04427 2.44070 -.02 .9855 -4.82796 4.73941

AQUA:COM| -.28650 .35391 -.81 .4182 -.98014 .40714

AQUA:FIS| -.02449 1.04592 -.02 .9813 -2.07446 2.02548

AQU0:FIS| -.05496 .70314 -.08 .9377 -1.43309 1.32317

|Standard deviations of parameter distributions

sdNEITHE| .01702 6.00963 .00 .9977 -11.76163 11.79568

sdCOM_FI| .32221* .18467 1.74 .0810 -.03974 .68417

sdFISH_P| .23578 .25612 .92 .3573 -.26620 .73776

sdFISH_S| .05332*** .01546 3.45 .0006 .02303 .08361

sdAQUATI| .29608 .26966 1.10 .2722 -.23244 .82461

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

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

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

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

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

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

NEITHER| 1.00000 -.14641 -.15003 -.15671 -.14952

COM_FISH| -.14641 1.00000 -.95360 -.77577 -.93532

FISH_POP| -.15003 -.95360 1.00000 .84317 .97091

FISH_SAV| -.15671 -.77577 .84317 1.00000 .69160

AQUATIC| -.14952 -.93532 .97091 .69160 1.00000
14

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

National Survey Version

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

|-> RESET $

|-> LOAD ; File="D:\Rob Johnston\Applied Research and Development\Wetland\SP

and BT 2009\Data Analysis\all_reg_model_data_w_wgts_071112.lpj"$

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 $

|-> REJECT; w_var = -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

; Wts = w_var

$

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

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

Start values obtained using MNL model

Dependent variable Choice

Weighting variable W_VAR

Log likelihood function -815.04713

Estimation based on N = 789, K = 6

Inf.Cr.AIC = 1642.1 AIC/N = 2.081

Model estimated: Nov 14, 2012, 12:23:39

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.= 789, skipped 0 obs

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

| Standard Prob. 95% Confidence
15

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

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

NEITHER| -.01031 .21894 -.05 .9625 -.43941 .41880

COM_FISH| .05326* .02778 1.92 .0552 -.00118 .10770

FISH_POP| .01302 .04298 .30 .7619 -.07121 .09726

FISH_SAV| .01347*** .00195 6.92 .0000 .00966 .01729

AQUATIC| .01242 .04334 .29 .7745 -.07253 .09736

COST| -.01770*** .00277 -6.38 .0000 -.02314 -.01226

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

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

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

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

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

Random Parameters Logit Model

Dependent variable CHOICE

Weighting variable W_VAR

Log likelihood function -667.47691

Restricted log likelihood -864.47725

Chi squared [ 21 d.f.] 394.00067

Significance level .00000

McFadden Pseudo R-squared .2278838

Estimation based on N = 789, K = 21

Inf.Cr.AIC = 1377.0 AIC/N = 1.745

Model estimated: Nov 14, 2012, 12:28:22

Constants only must be computed directly

Use NLOGIT ;...;RHS=ONE$

At start values -815.0471 .1811******

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| -1.49361*** .48000 -3.11 .0019 -2.43439 -.55284

COM_FISH| .07350 .06015 1.22 .2218 -.04440 .19140

FISH_POP| .19524* .11069 1.76 .0778 -.02171 .41219

FISH_SAV| .02116*** .00665 3.18 .0015 .00811 .03420

AQUATIC| -.10694 .17007 -.63 .5295 -.44028 .22640

|Nonrandom parameters in utility functions

COST| -.03294*** .00578 -5.70 .0000 -.04427 -.02162

|Diagonal values in Cholesky matrix, L.

NsNEITHE| .08915 1.29659 .07 .9452 -2.45212 2.63041

NsCOM_FI| .02810 .46660 .06 .9520 -.88642 .94262

NsFISH_P| .53312 .48289 1.10 .2696 -.41332 1.47956

NsFISH_S| .05068 .03377 1.50 .1333 -.01550 .11686

NsAQUATI| .80904** .39228 2.06 .0392 .04019 1.57789

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

COM_:NEI| .30365 .18788 1.62 .1061 -.06459 .67189

FISH:NEI| .13757 .55148 .25 .8030 -.94331 1.21845
16

FISH:COM| -.32723 .64754 -.51 .6133 -1.59638 .94191

FIS0:NEI| .05259** .02679 1.96 .0497 .00008 .10509

FIS0:COM| .02042 .03441 .59 .5529 -.04702 .08785

FISH:FIS| -.03767 .02453 -1.54 .1247 -.08575 .01041

AQUA:NEI| .10157 .47154 .22 .8295 -.82264 1.02577

AQUA:COM| .24749 .71827 .34 .7304 -1.16029 1.65528

AQUA:FIS| -.29335 .60552 -.48 .6281 -1.48015 .89345

AQU0:FIS| -.78390 .78590 -1.00 .3185 -2.32423 .75643

|Standard deviations of parameter distributions

sdNEITHE| .08915 1.29659 .07 .9452 -2.45212 2.63041

sdCOM_FI| .30494* .17425 1.75 .0801 -.03658 .64647

sdFISH_P| .64049** .30347 2.11 .0348 .04570 1.23528

sdFISH_S| .08467*** .02698 3.14 .0017 .03179 .13756

sdAQUATI| 1.19443* .72258 1.65 .0983 -.22180 2.61066

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

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

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

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

Cor.Mat.| NEITHER COM_FISH FISH_POP FISH_SAV AQUATIC

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

NEITHER| 1.00000 .99575 .21479 .62103 .08503

COM_FISH| .99575 1.00000 .16680 .64061 .10376

FISH_POP| .21479 .16680 1.00000 -.36007 -.29203

FISH_SAV| .62103 .64061 -.36007 1.00000 -.18082

AQUATIC| .08503 .10376 -.29203 -.18082 1.00000