(Luận văn) climate change and income diversification in the mekong river delta, a panel data analysis

to UNIVERSITY OF ECONOMICS HO CHI MINH CITY VIETNAM tn gh ERASMUS UNVERSITY ROTTERDAM INSTITUTE OF SOCIAL STUDIES THE NETHERLANDS p ie w n VIETNAM –THE NETHERLANDS PROGRAMME FOR M.A IN DEVELOPMENT ECONOMICS ad lo th yj an lu la ip uy CLIMATE CHANGE AND INCOME n va m ll fu DIVERSIFICATION IN THE MEKONG RIVER tz n oi DELTA: A PANEL DATA ANALYSIS z ht vb BY k jm NGUYEN THI TUYET NGA om l.c gm an Lu n va MASTER OF ARTS IN DEVELOPMENT ECONOMICS re y te th HO CHI MINH CITY, December 2016 to UNIVERSITY OF ECONOMICS HO CHIMINHCITY VIETNAM tn gh INSTITUTE OF SOCIAL STUDIES THE HAGUE THE NETHERLANDS p ie w n VIETNAM - NETHERLANDS PROGRAMME FOR M.A IN DEVELOPMENT ECONOMICS ad lo th yj la ip uy an lu CLIMATE CHANGE AND INCOME n va DIVERSIFICATION IN THE MEKONG RIVER m ll fu DELTA: A PANEL DATA ANALYSIS tz n oi z vb A thesis submitted in partial fulfilment of the requirements for the degree of ht MASTER OF ARTS IN DEVELOPMENT ECONOMICS k jm om NGUYEN THI TUYET NGA l.c gm By an Lu n va y te PHAM KHANH NAM re Academic Supervisor: th HO CHI MINH CITY, December 2016 ACKNOWLEDGEMENT t to ng I would first like to thank my thesis supervisorDr Pham Khanh Nam of the Vietnam – hi The Netherlands Programme (VNP) at Ho Chi Minh City University of Economics He ep consistently allowed this paper to be my own work, but steered me in the right direction whenever he thought I needed it w n ad lo I would like to express my gratitude to the VNP officers who were involved in mythesis processby updating thesis schedule and providing good conditions for my u yj th research process Without their passionate participation, the thesis process could not have yi been successfully conducted pl Finally, thanks are also due to my classmates for providing me with unfailing support n ua al and continuous encouragement throughout my years of study and through the process of n without them Thank you va researching and writing this thesis This accomplishment would not have been possible m ll fu t Ho Chi Minh City, December 2016 n oi Nguyen ThiTuyetNga z z j ht vb k m gm om l.c an Lu n va y te re Page i ABSTRACT t to ng The main objective of this study is to analyze the impact of climate change and other hi socio – economic determinants on the income diversification strategy in Mekong River ep Delta The data set is drawn from the VHLSS 2010,2012 and 2014, while climatic data including temperature and precipitation are extracted from the statistics website of The w n Ministry of Agriculture and Rural Development (MARD) Salinity data is collected from lo ad the Vietnam Institute of Meteorology, Hydrology and Climate Change.Findings of the u yj th study show that farmers have tendency to diversify their activities to reduce risk of crops failure when there is increasing temperature in dry season and precipitation in wet season yi pl However, it is recognized that those relationships are non-linear Diversification behavior n ua al is discovered not to be sensitive with the salinity intrusion and other climate variables Regarding socio-economic determinants, the household labor ratio and land area holding n va are found to be positively correlated with the income diversification, while educational m ll fu qualification has the negative effect A male household head would more likely to diversify their activities to disperse risk of climate change than female head From this n oi result, many policies are recommended in order to support farmers to access an effective t diversification strategy, helping them to response actively to climate change z consequences z j ht vb k m gm om l.c an Lu n va y te re Page ii TABLE OF CONTENT t to ng Chapter Page hi ep Acknowledgement i Abstract ii w Table of content iii n ad lo List of tables v List of figures vi u yj th Introduction yi 1.1 Research problem pl 1.2 Research objective n ua al 1.3 Research scope .3 1.4 Thesis structure va n Literature review m ll fu 2.1 Theoretical review 2.1.1 Climate change .4 n oi 2.1.2 Impact of climate change t 2.1.3 Adaptation of people to climate change .5 z 2.1.4 Income diversification z j ht vb 2.1.4.1 Definition and classification of income diversification 2.1.4.2 Motivations of income diversification m k 2.1.4.3 Income diversification measurements gm 2.2 Empirical review 10 om l.c 2.2.1 Impact of temperature and precipitation variation .10 2.2.2 Impact of high salinity intrusion to income diversification 12 an Lu 2.2.3 Impact of socio-economic characteristics on income diversification 14 Research methodology 16 3.2.3 Variable description .20 Page iii y te 3.2.2 Model specification 18 re 3.2.1 Income diversification index 17 n 3.2 Methodology 17 va 3.1 Analytical framework 16 (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis 3.3 Data sources 24 3.4 Salinity measurement 26 t to Result and discussion 28 ng 4.1 Overview of the Mekong River Delta 28 hi ep 4.1.1 Geographical position and natural conditions 28 4.1.2 Socio – economic conditions 28 4.1.3 Impact of climate change on the Mekong River Delta .29 w n 4.2 Salinity intrusion in the Mekong River Delta .31 lo ad 4.3 Descriptive statistics of variables 34 u yj th 4.3.1 Dependent variable 34 4.3.2 Independent variables 35 yi pl 4.4 Empirical results 41 n ua al 4.4.1 Findings of the Poisson model 41 4.4.2 Findings of the Tobit model .44 n va 4.4.3 Interpretation 45 fu Conclusion 50 m ll 5.1Conclusion .50 n oi 5.2Policy implications 51 t 5.3Research limitations and research directions 52 z z vb Reference 53 j ht Appendix .57 k m gm om l.c an Lu n va y te re Page iv (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis LIST OF TABLES t to ng Table 3.1.Variable description 21 hi ep Table 4.1.Descriptive statistics 34 Table 4.2.Results of the panel Poisson model and panel Tobit model 42 w n ad lo yi u yj th pl n ua al n va m ll fu t n oi z z j ht vb k m gm om l.c an Lu n va y te re Page v (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to tn gh p ie gender | 1667806 0507943 3.28 0.001 0672255 2663357 age | 0012656 0015894 0.80 0.426 -.0018495 0043807 land_ha | 0380917 0051055 7.46 0.000 0280851 0480983 scaled_salinity | -.0035491 0032377 -1.10 0.273 -.0098949 0027968 dry_temp | 5.436609 8.158859 0.67 0.505 -10.55446 21.42768 wet_temp | -16.22262 14.00837 -1.16 0.247 -43.67852 11.23328 dry_precipitation | 0005638 0017785 0.32 0.751 -.0029221 0040496 wet_precipitation | 0080996 0028302 2.86 0.004 0025525 0136467 sqr_dry_temp | -.1024101 1506563 -0.68 0.497 -.397691 1928709 sqr_wet_temp | 2908038 2512028 1.16 0.247 -.2015447 7831522 sqr_dryprecipitation | -2.87e-06 5.07e-06 -0.57 0.572 -.0000128 7.07e-06 sqr_wetprecipitation | -.0000175 6.12e-06 -2.85 0.004 -.0000295 -5.47e-06 _cons | 153.7981 256.2877 0.60 0.548 -348.5166 656.1127 -+ -/lnalpha | -16.47073 284.1286 -573.3526 540.4112 -+ -alpha | 7.03e-08 00002 9.9e-250 5.0e+234 -LR test of alpha=0: chibar2(01) = 0.00 Prob >= chibar2 = 1.000 w n ad lo th yj uy est sto re 0: 1: 2: 3: log log log log likelihood likelihood likelihood likelihood = = = = -809.66846 -807.99749 -807.99727 -807.99727 n va Iteration Iteration Iteration Iteration an lu la ip xtpoisson diversity_index1 `control' scaled_salinity dry_temp wet_temp dry_precipitation wet_precipitation sqr_dry_temp sqr_wet_temp sqr_dr > yprecipitation sqr_wetprecipitation , fe note: groups (7 obs) dropped because of only one obs per group m ll fu Conditional fixed-effects Poisson regression Group variable: id Number of obs Number of groups 1,064 355 Obs per group: = avg = max = 3.0 Log likelihood = -807.99727 tz n oi = = Wald chi2(16) Prob > chi2 = = 3.35 0.9996 z -diversity_index1 | Coef Std Err z P>|z| [95% Conf Interval] -+ -hh_size | 0064639 0308044 0.21 0.834 -.0539115 0668394 hh_labor_ratio | 0523888 2197282 0.24 0.812 -.3782706 4830481 migration | 0300864 0482235 0.62 0.533 -.0644299 1246028 education | 0334917 1367977 0.24 0.807 -.234627 3016103 gender | 0677677 3751553 0.18 0.857 -.6675233 8030586 age | 0014319 0360825 0.04 0.968 -.0692886 0721523 land_ha | 0149515 0136784 1.09 0.274 -.0118576 0417607 scaled_salinity | -.0038368 0202574 -0.19 0.850 -.0435405 035867 dry_temp | 5.342075 9.54153 0.56 0.576 -13.35898 24.04313 wet_temp | -3.154381 24.05328 -0.13 0.896 -50.29795 43.98919 dry_precipitation | 0002433 0026125 0.09 0.926 -.0048771 0053636 wet_precipitation | 0026639 0040137 0.66 0.507 -.0052027 0105306 sqr_dry_temp | -.0986839 1764474 -0.56 0.576 -.4445144 2471467 sqr_wet_temp | 0585224 4325856 0.14 0.892 -.7893298 9063745 sqr_dryprecipitation | -8.24e-07 7.47e-06 -0.11 0.912 -.0000155 0000138 sqr_wetprecipitation | -5.41e-06 8.52e-06 -0.63 0.526 -.0000221 0000113 ht vb k jm om l.c gm an Lu n va est sto fe re hausman fe re Page 38 th Coefficients -| (b) (B) (b-B) sqrt(diag(V_b-V_B)) | fe re Difference S.E -+ y te Note: the rank of the differenced variance matrix (11) does not equal the number of coefficients being tested (16); be sure this is what you expect, or there may be problems computing the test Examine the output of your estimators for anything unexpected and possibly consider scaling your variables so that the coefficients are on a similar scale (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to tn gh p ie hh_size | 0064639 0043159 002148 0275512 hh_labor_r~o | 0523888 1683662 -.1159774 2001915 migration | 0300864 0419248 -.0118383 031546 education | 0334917 -.0751008 1085924 135512 gender | 0677677 1667806 -.0990129 3717008 age | 0014319 0012656 0001663 0360475 land_ha | 0149515 0380917 -.0231402 0126898 scaled_sal~y | -.0038368 -.0035491 -.0002877 019997 dry_temp | 5.342075 5.436609 -.0945336 4.947101 wet_temp | -3.154381 -16.22262 13.06824 19.55316 dry_precip~n | 0002433 0005638 -.0003205 0019136 wet_precip~n | 0026639 0080996 -.0054356 0028459 sqr_dry_temp | -.0986839 -.1024101 0037262 0918496 sqr_wet_temp | 0585224 2908038 -.2322814 3521753 sqr_drypre~n | -8.24e-07 -2.87e-06 2.04e-06 5.48e-06 sqr_wetpre~n | -5.41e-06 -.0000175 0000121 5.93e-06 -b = consistent under Ho and Ha; obtained from xtpoisson B = inconsistent under Ha, efficient under Ho; obtained from xtpoisson w n ad lo th Ho: yj Test: difference in coefficients not systematic uy la ip chi2(11) = (b-B)'[(V_b-V_B)^(-1)](b-B) = 8.34 Prob>chi2 = 0.6822 an lu ** Choose the random effect model: n va xtpoisson diversity_index1 `control' scaled_salinity dry_temp wet_temp dry_precipitation wet_precipitation sqr_dry_temp sqr_wet_temp sqr_dr > yprecipitation sqr_wetprecipitation, re cluster() Iteration Iteration Iteration Iteration 0: 1: 2: 3: log log log log likelihood likelihood likelihood likelihood -1894.1417 -1665.3907 -1656.6608 -1655.4867 -1655.1319 -1655.0524 -1655.0379 -1655.0345 -1655.0337 -1655.0335 -1655.0335 Log likelihood Wald chi2(16) Prob > chi2 = -1655.0335 = avg = max = 3.0 = = 128.99 0.0000 re y te th Page 39 n va -diversity_index1 | Coef Std Err z P>|z| [95% Conf Interval] -+ -hh_size | 0043159 0137782 0.31 0.754 -.0226888 0313206 hh_labor_ratio | 1683662 0905751 1.86 0.063 -.0091578 3458902 migration | 0419248 0364741 1.15 0.250 -.0295631 1134126 education | -.0751008 0187117 -4.01 0.000 -.111775 -.0384265 gender | 1667806 0507943 3.28 0.001 0672255 2663357 age | 0012656 0015894 0.80 0.426 -.0018495 0043807 land_ha | 0380917 0051055 7.46 0.000 0280851 0480983 scaled_salinity | -.0035491 0032377 -1.10 0.273 -.0098949 0027968 dry_temp | 5.436609 8.158859 0.67 0.505 -10.55446 21.42768 wet_temp | -16.22262 14.00837 -1.16 0.247 -43.67852 11.23328 dry_precipitation | 0005638 0017785 0.32 0.751 -.0029221 0040496 wet_precipitation | 0080996 0028302 2.86 0.004 0025525 0136467 sqr_dry_temp | -.1024101 1506563 -0.68 0.497 -.397691 1928709 an Obs per group: 1,071 362 Lu Random effects u_i ~ Gamma = = om Number of obs Number of groups l.c Random-effects Poisson regression Group variable: id gm likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood k log log log log log log log log log log log jm 0: 1: 2: 3: 4: 5: 6: 7: 8: 9: 10: ht Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration vb = = = = = = = = = = = Fitting full model: z -1655.2474 -1655.0338 -1655.0335 -1655.0335 tz = = = = n oi m ll fu Fitting Poisson model: (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to tn gh p ie sqr_wet_temp | 2908038 2512028 1.16 0.247 -.2015447 7831522 sqr_dryprecipitation | -2.87e-06 5.07e-06 -0.57 0.572 -.0000128 7.07e-06 sqr_wetprecipitation | -.0000175 6.12e-06 -2.85 0.004 -.0000295 -5.47e-06 _cons | 153.7981 256.2877 0.60 0.548 -348.5166 656.1127 -+ -/lnalpha | -16.47073 284.1286 -573.3526 540.4112 -+ -alpha | 7.03e-08 00002 9.9e-250 5.0e+234 -LR test of alpha=0: chibar2(01) = 0.00 Prob >= chibar2 = 1.000 w ** Test the significance of the variable: n ad lo th yj test `control' scaled_salinity dry_temp wet_temp sqr_dry_temp sqr_wet_temp sqr_dryprecipitation sqr_wet > precipitation wet_precipitation [diversity_index1]hh_size = [diversity_index1]hh_labor_ratio = [diversity_index1]migration = [diversity_index1]education = [diversity_index1]gender = [diversity_index1]age = [diversity_index1]land_ha = [diversity_index1]scaled_salinity = [diversity_index1]dry_temp = [diversity_index1]wet_temp = [diversity_index1]dry_precipitation = [diversity_index1]wet_precipitation = [diversity_index1]sqr_dry_temp = [diversity_index1]sqr_wet_temp = [diversity_index1]sqr_dryprecipitation = [diversity_index1]sqr_wetprecipitation = an lu n va chi2( 16) = Prob > chi2 = 128.99 0.0000 jm 8.21 0.0165 k chi2( 2) = Prob > chi2 = ht [diversity_index1]wet_precipitation = [diversity_index1]sqr_wetprecipitation = vb ( 1) ( 2) z test wet_precipitation sqr_wetprecipitation tz n oi m ll fu ( 1) ( 2) ( 3) ( 4) ( 5) ( 6) ( 7) ( 8) ( 9) (10) (11) (12) (13) (14) (15) (16) la ip uy dry_precipitation gm xtpoisson diversity_index1 scaled_salinity c.dry_temp##c.dry_temp c.dry_precipitation##c.dry_precipitation c.wet_pre > cipitation##c.wet_precipitation `control', re cluster() c.wet_temp##c.wet_temp -1655.2403 -1655.0268 -1655.0265 -1655.0265 = = = = = = = = = -1894.1399 -1665.3846 -1656.6542 -1655.4798 -1655.1249 -1655.0454 -1655.031 -1655.0275 -1655.0267 re = = = = n likelihood likelihood likelihood likelihood va log log log log an 0: 1: 2: 3: Lu Iteration Iteration Iteration Iteration om Fitting Poisson model: l.c ** Calculate Marginal effect: 0: 1: 2: 3: 4: 5: 6: 7: 8: log log log log log log log log log likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood likelihood th Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration Iteration y te Fitting full model: Page 40 (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to Iteration 9: Iteration 10: Iteration 11: tn log likelihood = -1655.0265 log likelihood = -1655.0265 log likelihood = -1655.0265 gh p ie Random-effects Poisson regression Group variable: id Number of obs Number of groups Random effects u_i ~ Gamma Obs per group: w n Log likelihood = = 1,071 362 = avg = max = 3.0 = = 129.00 0.0000 Wald chi2(16) Prob > chi2 = -1655.0265 ad lo th diversity_index1 | Coef Std Err z P>|z| [95% Conf Interval] + -scaled_salinity | -.0035484 0032375 -1.10 0.273 -.0098938 002797 dry_temp | 5.442579 8.160219 0.67 0.505 -10.55116 21.43631 | c.dry_temp#c.dry_temp | -.1025196 1506811 -0.68 0.496 -.3978492 19281 | wet_temp | -16.27622 14.02314 -1.16 0.246 -43.76107 11.20862 | c.wet_temp#c.wet_temp | 2917645 2514663 1.16 0.246 -.2011005 7846294 | dry_precipitation | 0005633 0017784 0.32 0.751 -.0029224 0040489 | c.dry_precipitation#c.dry_precipitation | -2.86e-06 5.07e-06 -0.57 0.572 -.0000128 7.07e-06 | wet_precipitation | 0080959 0028312 2.86 0.004 0025468 0136449 | c.wet_precipitation#c.wet_precipitation | -.0000175 6.12e-06 -2.85 0.004 -.0000295 -5.45e-06 | hh_size | 0043187 0137781 0.31 0.754 -.0226859 0313233 hh_labor_ratio | 168364 0905751 1.86 0.063 -.0091599 3458879 migration | 0419269 0364738 1.15 0.250 -.0295603 1134142 education | -.0750974 0187117 -4.01 0.000 -.1117716 -.0384232 gender | 1667777 0507943 3.28 0.001 0672228 2663326 age | 0012656 0015894 0.80 0.426 -.0018495 0043807 land_ha | 0380908 0051054 7.46 0.000 0280844 0480973 _cons | 154.4647 256.5812 0.60 0.547 -348.4252 657.3546 + -/lnalpha | -17.66624 255.5196 -518.4755 483.143 + -alpha | 2.13e-08 5.43e-06 6.7e-226 6.7e+209 LR test of alpha=0: chibar2(01) = 5.4e-05 Prob >= chibar2 = 0.497 yj an lu la ip uy n va tz n oi m ll fu z ht vb ** Draw marginal effect graph of Precipitation in the wet season: jm k marginsplot Number of obs : wet_precip~n = 137.14 3._at : wet_precip~n = 167.14 4._at : wet_precip~n = 197.14 5._at : wet_precip~n = 227.14 6._at : wet_precip~n = 257.14 7._at : wet_precip~n = 287.14 8._at : wet_precip~n = 317.14 th 2._at y te 107.14 re = n : wet_precip~n va 1._at 1,071 an Expression : Linear prediction, predict() dy/dx w.r.t : wet_precipitation = Lu Average marginal effects Model VCE : OIM om margins, dydx(wet_precipitation) at( wet_precipitation =(107.14(30)345.29)) l.c gm Variables that uniquely identify margins: wet_precipitation Page 41 (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to tn gh p ie | Delta-method | dy/dx Std Err z P>|z| [95% Conf Interval] + -wet_precipitation | _at | | 004355 0015416 2.82 0.005 0013335 0073764 | 0033075 0011894 2.78 0.005 0009762 0056387 | 00226 0008503 2.66 0.008 0005934 0039265 | 0012125 0005489 2.21 0.027 0001366 0022884 | 000165 0003869 0.43 0.670 -.0005934 0009234 | -.0008825 0005178 -1.70 0.088 -.0018975 0001325 | -.00193 0008103 -2.38 0.017 -.0035182 -.0003418 | -.0029775 0011468 -2.60 0.009 -.0052252 -.0007298 - w n ad lo th // Tobit model: / is not a valid command name r(199); yj ip uy an lu la xttobit diversity_index2 scaled_salinity dry_temp wet_precipitation sqr_dry_temp sqr_wet_temp sqr_dryprecipitati > on sqr_wetprecipitation `control' , ul(1) re cluster() wet_temp dry_precipitation n va Obtaining starting values for full model: 0: 1: 2: 3: log log log log likelihood likelihood likelihood likelihood = = = = 164.51368 173.39111 173.56978 173.56995 = = = = = -338.06347 -270.15085 -269.09481 -269.08819 -269.08818 0: 1: 2: 3: 4: log log log log log likelihood likelihood likelihood likelihood likelihood tz Iteration Iteration Iteration Iteration Iteration n oi Fitting full model: m ll fu Iteration Iteration Iteration Iteration z = = Random effects u_i ~ Gaussian Obs per group: Integration method: mvaghermite Integration pts = Log likelihood Wald chi2(16) Prob > chi2 = = 1,071 362 ht 3.0 gm = -269.08818 k jm = avg = max = 12 47.59 0.0001 l.c Number of obs Number of groups vb Random-effects tobit regression Group variable: id om -diversity_index2 | Coef Std Err z P>|z| [95% Conf Interval] -+ -scaled_salinity | 0003955 0020019 0.20 0.843 -.0035282 0043193 dry_temp | -5.094096 2.920753 -1.74 0.081 -10.81867 6304745 wet_temp | -4.931524 6.094254 -0.81 0.418 -16.87604 7.012994 dry_precipitation | -.0003053 0006958 -0.44 0.661 -.0016691 0010585 wet_precipitation | -.0022184 001155 -1.92 0.055 -.0044822 0000454 sqr_dry_temp | 0948993 0539556 1.76 0.079 -.0108516 2006503 sqr_wet_temp | 0878975 1094411 0.80 0.422 -.1266032 3023982 sqr_dryprecipitation | 1.14e-06 1.98e-06 0.58 0.563 -2.73e-06 5.02e-06 sqr_wetprecipitation | 4.81e-06 2.46e-06 1.96 0.050 -6.47e-09 9.63e-06 hh_size | -.0089625 0072813 -1.23 0.218 -.0232336 0053087 hh_labor_ratio | -.047354 0497081 -0.95 0.341 -.1447801 0500721 migration | -.0157082 0147405 -1.07 0.287 -.0445989 0131826 education | 0508484 0129138 3.94 0.000 0255379 076159 gender | -.1413014 0345128 -4.09 0.000 -.2089452 -.0736576 age | -.0010515 001133 -0.93 0.353 -.003272 0011691 land_ha | -.003443 0036456 -0.94 0.345 -.0105882 0037023 _cons | 138.7191 107.0625 1.30 0.195 -71.11951 348.5577 -+ -/sigma_u | 2459029 0138871 17.71 0.000 2186847 273121 /sigma_e | 1954196 0060676 32.21 0.000 1835274 2073118 an Lu n va re y te th Page 42 (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to -+ -rho | 6129133 0312161 5505751 6724636 -0 left-censored observations 775 uncensored observations 296 right-censored observations tn gh p ie ** Test the significance of the variable: w n test scaled_salinity dry_temp wet_temp dry_precipitation wet_precipitation sqr_dry_temp sqr_wet_temp sqr_dryprecipitation sqr_wetprecipitat > ion `control' ad lo [diversity_index2]scaled_salinity = [diversity_index2]dry_temp = [diversity_index2]wet_temp = [diversity_index2]dry_precipitation = [diversity_index2]wet_precipitation = [diversity_index2]sqr_dry_temp = [diversity_index2]sqr_wet_temp = [diversity_index2]sqr_dryprecipitation = [diversity_index2]sqr_wetprecipitation = [diversity_index2]hh_size = [diversity_index2]hh_labor_ratio = [diversity_index2]migration = [diversity_index2]education = [diversity_index2]gender = [diversity_index2]age = [diversity_index2]land_ha = th yj an lu la ip uy n va ( 1) ( 2) ( 3) ( 4) ( 5) ( 6) ( 7) ( 8) ( 9) (10) (11) (12) (13) (14) (15) (16) 47.59 0.0001 [diversity_index2]dry_temp = [diversity_index2]sqr_dry_temp = chi2( 2) = Prob > chi2 = tz ( 1) ( 2) n oi test dry_temp sqr_dry_temp m ll fu chi2( 16) = Prob > chi2 = 7.05 0.0295 z 3.86 0.1448 k gm chi2( 2) = Prob > chi2 = jm [diversity_index2]wet_precipitation = [diversity_index2]sqr_wetprecipitation = ht ( 1) ( 2) vb test wet_precipitation sqr_wetprecipitation Lu c.wet_temp##c.wet_temp an xttobit diversity_index2 scaled_salinity c.dry_temp##c.dry_temp c.dry_precipitation##c.dry_precipitation c.wet_preci > pitation##c.wet_precipitation `control', re cluster() om l.c ** Calculate Marginal effect: va Obtaining starting values for full model: likelihood likelihood likelihood likelihood = = = = 164.51055 173.38895 173.56764 173.56781 log likelihood = log likelihood = 173.56781 173.56781 y te log log log log re 0: 1: 2: 3: n Iteration Iteration Iteration Iteration Iteration 0: Iteration 1: Random-effects tobit regression Group variable: id th Fitting full model: Number of obs Number of groups = = 1,071 362 Page 43 (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis (Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis(Luận.văn).climate.change.and.income.diversification.in.the.mekong.river.delta a.panel.data.analysis to Random effects u_i ~ Gaussian tn Obs per group: = avg = max = 3.0 gh p ie Integration pts = 12 Log likelihood Wald chi2(16) Prob > chi2 = = 41.25 0.0005 Integration method: mvaghermite = 173.56781 w diversity_index2 | Coef Std Err z P>|z| [95% Conf Interval] + -scaled_salinity | 0004745 0013884 0.34 0.733 -.0022467 0031957 dry_temp | -3.97197 2.255179 -1.76 0.078 -8.392039 4481002 | c.dry_temp#c.dry_temp | 0739636 0416519 1.78 0.076 -.0076726 1555999 | wet_temp | -3.129728 4.598307 -0.68 0.496 -12.14224 5.882788 | c.wet_temp#c.wet_temp | 0556771 0825489 0.67 0.500 -.1061158 21747 | dry_precipitation | -.0002633 0005265 -0.50 0.617 -.0012952 0007686 | c.dry_precipitation#c.dry_precipitation | 9.19e-07 1.49e-06 0.62 0.538 -2.01e-06 3.84e-06 | wet_precipitation | -.0014003 0008433 -1.66 0.097 -.003053 0002525 | c.wet_precipitation#c.wet_precipitation | 3.04e-06 1.80e-06 1.69 0.091 -4.87e-07 6.57e-06 | hh_size | -.0065029 0053121 -1.22 0.221 -.0169145 0039087 hh_labor_ratio | -.0536531 0355123 -1.51 0.131 -.1232558 0159497 migration | -.0161259 0116518 -1.38 0.166 -.038963 0067112 education | 0291604 0085039 3.43 0.001 0124931 0458276 gender | -.0845108 0228651 -3.70 0.000 -.1293257 -.0396959 age | -.0005967 0007589 -0.79 0.432 -.0020842 0008907 land_ha | -.0002609 0027889 -0.09 0.925 -.0057271 0052054 _cons | 98.29724 80.98127 1.21 0.225 -60.42313 257.0176 + -/sigma_u | 1595978 0085215 18.73 0.000 142896 1762996 /sigma_e | 1643946 0044296 37.11 0.000 1557128 1730763 + -rho | 485198 032434 4220835 5486863 - n ad lo th yj an lu la ip uy n va tz n oi m ll fu left-censored observations 1,071 uncensored observations right-censored observations z **Draw marginal effect graph of Temperature in the dry season: ht vb k jm marginsplot Number of obs = : dry_temp = 26.6 4._at : dry_temp = 26.7 5._at : dry_temp = 26.8 6._at : dry_temp = 26.9 7._at : dry_temp = 27 8._at : dry_temp = 27.1 th 3._at y te 26.5 re 26.4 = n = : dry_temp va : dry_temp 2._at an 1._at Lu Expression : E(diversity_index2|diversity_index2