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

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

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

******************************************************************************************************* 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% Confidence6

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

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

******************************************************************************************************* 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% Confidence9

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

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

******************************************************************************************************* 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% Confidence12

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

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

******************************************************************************************************* 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% Confidence15

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

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