Luận án Tiến sĩ Tính dễ tổn thương và chiến lược thích ứng biến đổi khí hậu trong sinh kế của người dân tộc Khmer ở đồng bằng sông Cửu Long

ỦY BAN NHÂN DÂN TỈNH TRÀ VINH TRƯỜNG ĐẠI HỌC TRÀ VINH ISO 9001:2015 NGUYỄN THỊ THUÝ LOAN TÍNH DỄ TỔN THƯƠNG VÀ CHIẾN LƯỢC THÍCH ỨNG BIẾN ĐỔI KHÍ HẬU TRONG SINH KẾ CỦA NGƯỜI DÂN TỘC KHMER Ở ĐỒNG BẰNG SÔNG CỬU LONG LUẬN ÁN TIẾN SĨ NGÀNH QUẢN LÝ KINH TẾ TRÀ VINH, NĂM 2023 ỦY BAN NHÂN DÂN TỈNH TRÀ VINH TRƯỜNG ĐẠI HỌC TRÀ VINH NGUYỄN THỊ THUÝ LOAN TÍNH DỄ TỔN THƯƠNG VÀ CHIẾN LƯỢC THÍCH ỨNG BIẾN ĐỔI KHÍ HẬU TRONG SINH KẾ CỦA NGƯỜI DÂN TỘC KHMER Ở ĐỒNG BẰNG SÔNG CỬU LONG Ngành: Quản lý kinh tế Mã ngành: 9310110 LUẬN ÁN TIẾN SĨ NGÀNH QUẢN LÝ KINH TẾ Người hướng dẫn khoa học: PGS TS DIỆP THANH TÙNG GS TS VÕ XUÂN VINH TRÀ VINH, NĂM 2023 LỜI CAM ĐOAN  Tôi tên Nguyễn Thị Thuý Loan xin cam đoan số liệu, kết thảo luận, phân tích luận án cơng trình nghiên cứu tơi, đảm bảo tính trung thực khơng trùng với cơng trình Ngày 06 tháng năm 2023 Nghiên cứu sinh Nguyễn Thị Thuý Loan i LỜI CẢM ƠN  Luận án thực theo chương trình đào tạo ngành Quản lý Kinh tế, bậc tiến sĩ trường Đại học Trà Vinh Tác giả xin cảm ơn sâu sắc đến Qúy Thầy/Cô Ban Giám hiệu, Khoa Kinh tế, Luật, Trung tâm Đào tạo Logistics Thương mại điện tử Phòng Đào tạo Sau đại học trường Đại học Trà Vinh hướng dẫn, giúp đỡ hỗ trợ tận tình trình học tập tác giả Đặc biệt, tác giả xin chân thành cảm ơn sâu sắc đến Thầy PGS TS Diệp Thanh Tùng Thầy GS TS Võ Xuân Vinh người hướng dẫn khoa học, tận tình hướng dẫn đóng góp nhiều ý kiến quý báu suốt trình thực luận án Tác giả xin chân thành cảm ơn tất Qúy Thầy/Cô thành viên hội đồng bảo vệ đề cương, chuyên đề tổng quan chuyên đề tự chọn có nhiều ý kiến đóng góp qúy báu giới thiệu cho tác giả nhiều tài liệu tham khảo hữu ích Tác giả xin trân trọng cảm ơn cán Ban Dân tộc, Sở Nông nghiệp Phát triển nông thôn, cán UBND huyện, xã hộ gia đình Khmer tỉnh Trà Vinh, Sóc Trăng Kiên Giang, cung cấp số liệu tạo điều kiện thuận lợi cho tác giả thu thập số liệu nghiên cứu Ngoài ra, tác giả xin chân thành cảm ơn đến gia đình, đồng nghiệp bạn bè giúp đỡ ủng hộ tinh thần cho tác giả nhiều thời gian qua Do hạn chế kiến thức nên luận án khơng tránh sai sót Tác giả mong nhận ý kiến góp ý q báu q Thầy/Cơ để hồn chỉnh Trà Vinh, ngày 06 tháng năm 2023 ii MỤC LỤC LỜI CAM ĐOAN i LỜI CẢM ƠN ii DANH MỤC KÝ HIỆU, CHỮ VIẾT TẮT .vii DANH MỤC BẢNG BIỂU viii DANH MỤC HÌNH x TÓM TẮT x ABSTRACT xii CHƯƠNG TỔNG QUAN VẤN ĐỀ NGHIÊN CỨU 1.1 LÝ DO CHỌN ĐỀ TÀI 1.1.1 Bối cảnh thực tiễn 1.1.2 Bối cảnh lý thuyết 1.2 MỤC TIÊU NGHIÊN CỨU 1.2.1 Mục tiêu chung 1.2.2 Mục tiêu cụ thể 1.3 CÂU HỎI NGHIÊN CỨU 1.4 ĐỐI TƯỢNG NGHIÊN CỨU VÀ ĐỐI TƯỢNG KHẢO SÁT 1.4.1 Đối tượng nghiên cứu 1.4.2 Đối tượng khảo sát 1.5 PHẠM VI NGHIÊN CỨU 1.5.1 Phạm vi nội dung 1.5.2 Phạm vi không gian 1.5.3 Phạm vi thời gian 1.6 PHƯƠNG PHÁP NGHIÊN CỨU 1.6.1 Phương pháp thu thập liệu 1.6.2 Phương pháp phân tích liệu 1.7 CẤU TRÚC CỦA LUẬN ÁN 1.8 TÓM TẮT CHƯƠNG 10 CHƯƠNG CƠ SỞ LÝ THUYẾT VÀ KHUNG PHÂN TÍCH 11 2.1 CƠ SỞ LÝ THUYẾT 11 2.1.1 Các khái niệm liên quan 11 2.1.2 Khung phân tích sinh kế bền vững 17 iii 2.1.3 Khung phân tích tính dễ bị tổn thương sinh kế biến đổi khí hậu 20 2.2 TỔNG QUAN CÁC NGHIÊN CỨU 25 2.2.1 Đánh giá tác động biến đổi khí hậu đến sinh kế 25 2.2.2 Chiến lược thích ứng biến đổi khí hậu 29 2.2.3 Tính dễ bị tổn thương sinh kế 35 2.2.4 Đánh giá chung nghiên cứu khoảng trống nghiên cứu 45 2.3 KHUNG PHÂN TÍCH ĐỀ XUẤT 53 2.4 TÓM TẮT CHƯƠNG 55 CHƯƠNG THIẾT KẾ NGHIÊN CỨU 56 3.1 QUY TRÌNH NGHIÊN CỨU 56 3.2 PHƯƠNG PHÁP NGHIÊN CỨU 58 3.2.1 Phương pháp thu thập liệu 58 3.2.2 Phương pháp phân tích liệu 59 3.3 TÓM TẮT CHƯƠNG 84 CHƯƠNG KẾT QUẢ NGHIÊN CỨU VÀ THẢO LUẬN 85 4.1 THỰC TRẠNG VỀ BIẾN ĐỔI KHÍ HẬU Ở VIỆT NAM VÀ KHU VỰC ĐỒNG BẰNG SÔNG CỬU LONG TRONG THỜI GIAN QUA 85 4.1.1 Thực trạng biến đổi khí hậu Việt Nam thời gian qua 85 4.1.2 Thực trạng biến đổi khí hậu Đồng sông Cửu Long 87 4.2 PHÂN TÍCH THỰC TRẠNG TIẾP CẬN CÁC NGUỒN TÀI SẢN SINH KẾ CỦA CÁC HỘ GIA ĐÌNH NGƯỜI KHMER Ở ĐỒNG BẰNG SÔNG CỬU LONG TRƯỚC TÁC ĐỘNG CỦA BIẾN ĐỔI KHÍ HẬU 90 4.2.1 Vốn người 90 4.2.2 Vốn tự nhiên 93 4.2.3 Vốn tài 98 4.2.4 Vốn vật chất 102 4.2.5 Vốn xã hội 103 4.3 PHÂN TÍCH THỰC TRẠNG ÁP DỤNG CÁC CHIẾN LƯỢC THÍCH ỨNG BIẾN ĐỔI KHÍ HẬU CỦA CÁC HỘ GIA ĐÌNH NGƯỜI KHMER Ở ĐỒNG BẰNG SÔNG CỬU LONG 104 4.3.1 Chiến lược thích ứng ứng phó 108 4.3.2 Chiến lược thích ứng chủ động 109 iv 4.4 ĐO LƯỜNG TÍNH DỄ BỊ TỔN THƯƠNG VỀ SINH KẾ CỦA NGƯỜI KHMER Ở ĐỒNG BẰNG SÔNG CỬU LONG 111 4.4.1 Thống kê mô tả số thành phần phụ 112 4.4.2 Thống kê mô tả thành phần 117 4.4.3 Mức độ tổn thương sinh kế theo số LVI LVI-IPCC hộ gia đình người Khmer đồng sông Cửu Long 119 4.5 PHÂN TÍCH TÁC ĐỘNG CỦA CÁC CHIẾN LƯỢC THÍCH ỨNG BIẾN ĐỔI KHÍ HẬU ĐẾN TÍNH DỄ BỊ TỔN THƯƠNG VẾ SINH KẾ CỦA CÁC HỘ GIA ĐÌNH KHMER Ở ĐỒNG BẰNG SƠNG CỬU LONG 121 4.5.1 Kết mơ hình 121 4.5.2 Thảo luận kết nghiên cứu 123 4.6 TÓM TẮT CHƯƠNG 133 CHƯƠNG KẾT LUẬN VÀ HÀM Ý CHÍNH SÁCH 135 5.1 KẾT LUẬN 135 5.2 HÀM Ý CHÍNH SÁCH 137 5.2.1 Định hướng sách nhằm giảm thiểu tính dễ bị tổn thương sinh kế hộ gia đình người Khmer đồng sông Cửu Long 137 5.2.2 Một số hàm ý sách giảm thiểu tính dễ bị tổn thương sinh kế hộ gia đình Khmer đồng sơng Cửu Long 138 5.3 ĐÓNG GÓP CỦA LUẬN ÁN 144 5.3.1 Đóng góp mặt lý thuyết 144 5.3.2 Đóng góp mặt thực tiễn 145 5.3.3 Tính luận án 146 5.4 HẠN CHẾ VÀ HƯỚNG NGHIÊN CỨU 148 5.4.1 Hạn chế nghiên cứu 148 5.4.2 Hướng nghiên cứu 149 5.5 TÓM TẮT CHƯƠNG 149 DANH MỤC TÀI LIỆU THAM KHẢO 151 A VĂN BẢN QUY PHẠM PHÁP LUẬT 151 B TÀI LIỆU TIẾNG VIỆT 152 C TÀI LIỆU TIẾNG NƯỚC NGOÀI 155 D TÀI LIỆU ĐIỆN TỬ 170 v DANH MỤC CƠNG TRÌNH CỦA TÁC GIẢ 172 PHỤ LỤC PHỤ LỤC BẢN HỎI HỘ GIA ĐÌNH KHMER 2.1 BẢN HỎI HỘ GIA ĐÌNH KHMER SƠ BỘ 12 2.2 BẢN HỎI HỘ GIA ĐÌNH KHMER CHÍNH THỨC 14 PHỤ LỤC PHỎNG VẤN Ý KIẾN CHUYÊN GIA 12 2.1 DANH SÁCH PHỎNG VẤN CHUYÊN GIA 12 2.2 PHIẾU PHỎNG VẤN CHUYÊN GIA 14 2.2.1 Phiếu vấn chuyên gia lấy ý kiến điều chỉnh 14 2.2.2 Phiếu vấn chuyên gia lấy ý kiến hàm ý sách 22 2.3 TỔNG HỢP KẾT QUẢ PHỎNG VẤN Ý KIẾN CHUYÊN GIA 24 PHỤ LỤC KẾT QUẢ KIỂM ĐỊNH GIẢ THUYẾT 31 3.1 KIỂM ĐỊNH ANOVA 31 3.2 KIỂM ĐỊNH CHI-SQUARE 39 3.3 KIỂM ĐỊNH T-TEST 43 PHỤ LỤC PHÂN TÍCH THÀNH PHẦN CHÍNH PCA 49 4.1 KẾT QUẢ PCA CHO CÁC THÀNH PHẦN PHỤ 49 4.2 TÍNH PCA CỦA 03 NHÂN TỐ CHÍNH CỦA LVI-IPCC (KHẢ NĂNG THÍCH ỨNG, ĐỘ NHẠY CẢM VÀ MỨC ĐỘ PHƠI NHIỄM) 52 PHỤ LỤC MƠ HÌNH HỒI QUY 54 5.1 MƠ HÌNH HỒI QUY 54 5.3 KIỂM ĐỊNH SỰ PHÙ HỢP CỦA MƠ HÌNH 55 vi DANH MỤC KÝ HIỆU, CHỮ VIẾT TẮT Ac: Adaptive capacity – Khả thích ứng BĐKH: Biến đổi khí hậu DBTT: Dễ bị tổn thương DFID: Department for International Development - Cơ quan Phát triển Quốc tế Vương quốc Anh DTTS: Dân tộc thiểu số ĐBSCL: Đồng sông Cửu Long Ex: Exposure – Mức độ phơi nhiễm IPCC: Intergovernmental Panel on Climate Change - Ủy ban Liên Chính phủ Biến đổi khí hậu LVI: Livelihood Vulnerability Index - Chỉ số dễ bị tổn thương sinh kế LVI-IPCC: Livelihood Vulnerability Index - Intergovernmental Panel on Climate Change - Chỉ số dễ bị tổn thương sinh kế tính theo định nghĩa Ủy ban Liên Chính phủ Biến đổi khí hậu PCA: Principal Component Analysis – Phân tích thành phần Se: Sensitivity – Mức độ nhảy cảm SLF: Sustainable Livelihood Framework - Khung sinh kế bền vững vii DANH MỤC BẢNG BIỂU Bảng 2.1 Các hoạt động thích ứng BĐKH sinh kế hộ gia đình ĐBSCL 14 Bảng 2.2 Các hoạt động thích ứng theo nhóm chiến lược thích ứng chủ động ứng phó cấp độ hộ gia đình khu vực ĐBSCL 15 Bảng 2.3 Tổng hợp lược khảo tài liệu nghiên cứu nước 46 Bảng 3.1 Cỡ mẫu điều tra 59 Bảng 3.2 Tổng hợp số đại diện cho năm nguồn tài sản sinh kế 60 Bảng 3.3 Tổng hợp yếu tố phụ yếu tố để tính LVI 66 Bảng 3.4 Phân nhóm thành phần theo yếu tố ảnh hưởng IPCC 71 Bảng 3.5 Diễn giải ý nghĩa dấu kỳ vọng biến mơ hình 79 Bảng 4.1 Tổng thiệt hại ước tính thiên tai Việt Nam từ năm 2015 đến tháng 10 năm 2022 86 Bảng 4.2 Thực trạng số đại diện nguồn vốn người hộ gia đình Khmer ĐBSCL 90 Bảng 4.3 Kết phân tích sâu Anova số lao động trung bình tỉnh 91 Bảng 4.4 Kết phân tích sâu Anova trình độ học vấn trung bình tỉnh 91 Bảng 4.5 Thực trạng nguồn vốn đất đai hộ Khmer ĐBSCL 94 Bảng 4.6 Kết phân tích sâu Anova diện tích đất nơng nghiệp tỉnh 94 Bảng 4.7 Hiện trạng sẵn có nguồn vốn tự nhiên 95 Bảng 4.8 Thực trạng ảnh hưởng BĐKH hoạt động sản xuất người dân tộc Khmer ĐBSCL 98 Bảng 4.9 Thu nhập hộ gia đình Khmer ĐBSCL 99 Bảng 4.10 Kết phân tích sâu Anova tổng thu nhập hộ Khmer tỉnh 99 Bảng 4.11 Thực trạng vay vốn hộ gia đình Khmer ĐBSCL 101 Bảng 4.12 Thực trạng sở hạ tầng địa phương 102 Bảng 4.13 Thực trạng nhà hộ gia đình Khmer ĐBSCL 102 Bảng 4.14 Kết phân tích sâu Anova khác biệt nhóm chiến lược thích ứng 106 Bảng 4.15 Kết kiểm định Chi-Square khác biệt nhóm chiến lược thích ứng 107 viii 41 42 tab CLAP GIOI, col chi2 + -+ | Key | | -| | frequency | | column percentage | + -+ | GIOI CLAP | | Total -+ + -1 | 70 45 | 115 | 21.54 44.55 | 27.00 -+ + -2 | 33 | 41 | 10.15 7.92 | 9.62 -+ + -3 | 222 48 | 270 | 68.31 47.52 | 63.38 -+ + -Total | 325 101 | 426 | 100.00 100.00 | 100.00 Pearson chi2(2) = 20.7709 Pr = 0.000 tab CLAP NHA, col chi2 + -+ | Key | | -| | frequency | | column percentage | + -+ | NHA CLAP | | Total -+ + -1 | 68 47 | 115 | 30.91 22.82 | 27.00 -+ + -2 | 29 12 | 41 | 13.18 5.83 | 9.62 -+ + -3 | 123 147 | 270 | 55.91 71.36 | 63.38 -+ + -Total | 220 206 | 426 | 100.00 100.00 | 100.00 Pearson chi2(2) = 12.5704 Pr = 0.002 3.3 Kiểm định T-test sdtest LVIce, by (GIOI) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 325 3615722 01093 1970439 3400694 3830749 | 101 4893136 0191349 1923035 4513504 5272767 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 1.0499 Ho: ratio = degrees of freedom = 324, 100 Ha: ratio < Pr(F < f) = 0.6067 Ha: ratio != 2*Pr(F > f) = 0.7865 43 Ha: ratio > Pr(F > f) = 0.3933 ttest LVIce, by (GIOI) Two-sample t test with equal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 325 3615722 01093 1970439 3400694 3830749 | 101 4893136 0191349 1923035 4513504 5272767 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | -.1277414 0223212 -.1716153 -.0838674 -diff = mean(0) - mean(1) t = -5.7229 Ho: diff = degrees of freedom = 424 Ha: diff < Pr(T < t) = 0.0000 Ha: diff != Pr(|T| > |t|) = 0.0000 Ha: diff > Pr(T > t) = 1.0000 sdtest LVIce, by (NHA) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 220 4437625 0118905 1763644 4203281 4671969 | 206 3364266 0149967 2152436 3068589 3659942 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 0.6714 Ho: ratio = degrees of freedom = 219, 205 Ha: ratio < Pr(F < f) = 0.0019 Ha: ratio != 2*Pr(F < f) = 0.0038 Ha: ratio > Pr(F > f) = 0.9981 ttest LVIce, by (NHA) unequal Two-sample t test with unequal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 220 4437625 0118905 1763644 4203281 4671969 | 206 3364266 0149967 2152436 3068589 3659942 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | 1073359 0191386 0697103 1449616 -diff = mean(0) - mean(1) t = 5.6084 Ho: diff = Satterthwaite's degrees of freedom = 396.925 Ha: diff < Pr(T < t) = 1.0000 Ha: diff != Pr(|T| > |t|) = 0.0000 Ha: diff > Pr(T > t) = 0.0000 sdtest LVIce, by (HONGHEO) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 372 3771347 0104432 2014217 3565993 39767 | 54 4932876 0254054 1866904 4423309 5442443 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 1.1640 Ho: ratio = degrees of freedom = 371, 53 Ha: ratio < Pr(F < f) = 0.7468 Ha: ratio != 2*Pr(F > f) = 0.5065 44 Ha: ratio > Pr(F > f) = 0.2532 ttest LVIce, by (HONGHEO) Two-sample t test with equal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 372 3771347 0104432 2014217 3565993 39767 | 54 4932876 0254054 1866904 4423309 5442443 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | -.1161529 0290726 -.1732972 -.0590086 -diff = mean(0) - mean(1) t = -3.9953 Ho: diff = degrees of freedom = 424 Ha: diff < Pr(T < t) = 0.0000 Ha: diff != Pr(|T| > |t|) = 0.0001 Ha: diff > Pr(T > t) = 1.0000 sdtest LVIce, by (TINDUNG) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 162 3483102 0165217 2102869 315683 3809374 | 264 418581 0119529 1942117 3950454 4421165 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 1.1724 Ho: ratio = degrees of freedom = 161, 263 Ha: ratio < Pr(F < f) = 0.8730 Ha: ratio != 2*Pr(F > f) = 0.2541 Ha: ratio > Pr(F > f) = 0.1270 ttest LVIce, by (TINDUNG) Two-sample t test with equal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 162 3483102 0165217 2102869 315683 3809374 | 264 418581 0119529 1942117 3950454 4421165 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | -.0702708 0200073 -.1095967 -.0309448 -diff = mean(0) - mean(1) t = -3.5122 Ho: diff = degrees of freedom = 424 Ha: diff < Pr(T < t) = 0.0002 Ha: diff != Pr(|T| > |t|) = 0.0005 Ha: diff > Pr(T > t) = 0.9998 sdtest LVIce, by (TCTT) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 156 4785005 0128633 1606626 4530905 5039105 | 270 3417983 0126848 2084329 3168241 3667725 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 0.5942 Ho: ratio = degrees of freedom = 155, 269 Ha: ratio < Pr(F < f) = 0.0002 Ha: ratio != 2*Pr(F < f) = 0.0004 45 Ha: ratio > Pr(F > f) = 0.9998 ttest LVIce, by (TCTT) unequal Two-sample t test with unequal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 156 4785005 0128633 1606626 4530905 5039105 | 270 3417983 0126848 2084329 3168241 3667725 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | 1367022 0180657 101184 1722204 -diff = mean(0) - mean(1) t = 7.5669 Ho: diff = Satterthwaite's degrees of freedom = 390.34 Ha: diff < Pr(T < t) = 1.0000 Ha: diff != Pr(|T| > |t|) = 0.0000 Ha: diff > Pr(T > t) = 0.0000 sdtest LVIce, by (tcxh) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 163 4727385 0140034 1787834 4450858 5003913 | 263 3417309 0124204 2014254 3172744 3661875 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 0.7878 Ho: ratio = degrees of freedom = 162, 262 Ha: ratio < Pr(F < f) = 0.0488 Ha: ratio != 2*Pr(F < f) = 0.0976 Ha: ratio > Pr(F > f) = 0.9512 ttest LVIce, by (tcxh) Two-sample t test with equal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 163 4727385 0140034 1787834 4450858 5003913 | 263 3417309 0124204 2014254 3172744 3661875 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | 1310076 0192481 093174 1688413 -diff = mean(0) - mean(1) t = 6.8062 Ho: diff = degrees of freedom = 424 Ha: diff < Pr(T < t) = 1.0000 Ha: diff != Pr(|T| > |t|) = 0.0000 Ha: diff > Pr(T > t) = 0.0000 sdtest LVIce, by (chudong) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 156 5207166 0128377 1603429 4953572 5460761 | 270 3174068 0114146 1875611 2949335 3398801 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 0.7308 Ho: ratio = degrees of freedom = 155, 269 Ha: ratio < Pr(F < f) = 0.0158 Ha: ratio != 2*Pr(F < f) = 0.0316 46 Ha: ratio > Pr(F > f) = 0.9842 ttest LVIce, by (chudong) unequal Two-sample t test with unequal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 156 5207166 0128377 1603429 4953572 5460761 | 270 3174068 0114146 1875611 2949335 3398801 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | 2033098 0171785 1695287 2370909 -diff = mean(0) - mean(1) t = 11.8351 Ho: diff = Satterthwaite's degrees of freedom = 365.373 Ha: diff < Pr(T < t) = 1.0000 Ha: diff != Pr(|T| > |t|) = 0.0000 Ha: diff > Pr(T > t) = 0.0000 sdtest LVIce, by (ungpho) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 385 3859425 0106559 2090845 3649912 4068938 | 41 4474086 0192556 1232961 4084916 4863257 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 2.8757 Ho: ratio = degrees of freedom = 384, 40 Ha: ratio < Pr(F < f) = 0.9999 Ha: ratio != 2*Pr(F > f) = 0.0001 Ha: ratio > Pr(F > f) = 0.0001 ttest LVIce, by (ungpho) unequal Two-sample t test with unequal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 385 3859425 0106559 2090845 3649912 4068938 | 41 4474086 0192556 1232961 4084916 4863257 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | -.0614661 0220074 -.1053861 -.0175461 -diff = mean(0) - mean(1) t = -2.7930 Ho: diff = Satterthwaite's degrees of freedom = 67.5907 Ha: diff < Pr(T < t) = 0.0034 Ha: diff != Pr(|T| > |t|) = 0.0068 Ha: diff > Pr(T > t) = 0.9966 sdtest LVIce, by (HTDP) Variance ratio test -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 247 4124188 0121656 1911967 3884568 4363807 | 179 363487 016134 2158576 3316486 3953255 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -ratio = sd(0) / sd(1) f = 0.7846 Ho: ratio = degrees of freedom = 246, 178 Ha: ratio < Pr(F < f) = 0.0394 Ha: ratio != 2*Pr(F < f) = 0.0787 47 Ha: ratio > Pr(F > f) = 0.9606 ttest LVIce, by (HTDP) unequal Two-sample t test with unequal variances -Group | Obs Mean Std Err Std Dev [95% Conf Interval] -+ -0 | 247 4124188 0121656 1911967 3884568 4363807 | 179 363487 016134 2158576 3316486 3953255 -+ -combined | 426 3918583 0098414 2031235 3725145 4112021 -+ -diff | 0489317 0202066 0091921 0886714 -diff = mean(0) - mean(1) t = 2.4216 Ho: diff = Satterthwaite's degrees of freedom = 354.93 Ha: diff < Pr(T < t) = 0.9920 Ha: diff != Pr(|T| > |t|) = 0.0160 48 Ha: diff > Pr(T > t) = 0.0080 PHỤ LỤC PHÂN TÍCH THÀNH PHẦN CHÍNH PCA 4.1 Kết PCA cho thành phần phụ pca S11 S12 S13 S14 S15, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 6501644 Rotation: (unrotated = principal) Rho = 0.5024 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 326637 168826 0.5024 0.5024 Comp2 | 157811 0466974 0.2427 0.7451 Comp3 | 111114 0748182 0.1709 0.9160 Comp4 | 0362955 0179881 0.0558 0.9718 Comp5 | 0183074 0.0282 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S11 | 0.0201 | 01897 S12 | 0.3896 | 1317 S13 | 0.1125 | 0355 S14 | 0.5150 | 08176 S15 | 0.7550 | 05558 - pca S21 S22 S23 S24 S25, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 9896879 Rotation: (unrotated = principal) Rho = 0.3902 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 386145 189229 0.3902 0.3902 Comp2 | 196916 00727724 0.1990 0.5891 Comp3 | 189639 0452739 0.1916 0.7808 Comp4 | 144365 071742 0.1459 0.9266 Comp5 | 072623 0.0734 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S21 | 0.2176 | 1782 S22 | -0.6377 | 06878 S23 | 0.3600 | 07148 S24 | 0.0794 | 1951 S25 | 0.6404 | 09001 - pca S31 S32 S33 S34, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 5502126 Rotation: (unrotated = principal) Rho = 0.3710 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 204113 0498732 0.3710 0.3710 49 Comp2 | 15424 043819 0.2803 0.6513 Comp3 | 110421 0289826 0.2007 0.8520 Comp4 | 0814384 0.1480 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S31 | 0.8617 | 03812 S32 | 0.2709 | 1258 S33 | 0.1394 | 08131 S34 | 0.4057 | 1009 - pca S41 S42 S43 S44 S45 S46, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 8138914 Rotation: (unrotated = principal) Rho = 0.3543 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 288323 0393781 0.3543 0.3543 Comp2 | 248945 104563 0.3059 0.6601 Comp3 | 144381 0547727 0.1774 0.8375 Comp4 | 0896085 0628613 0.1101 0.9476 Comp5 | 0267473 0108602 0.0329 0.9805 Comp6 | 0158871 0.0195 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S41 | 0.0068 | 01605 S42 | 0.3359 | 2001 S43 | 0.7229 | 09495 S44 | 0.5930 | 08831 S45 | 0.0726 | 02814 S46 | 0.0877 | 098 - pca S51 S52, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 35657 Rotation: (unrotated = principal) Rho = 0.6890 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 245676 134782 0.6890 0.6890 Comp2 | 110894 0.3110 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S51 | 0.0209 | 1108 S52 | 0.9998 | 00004846 - pca S61 S62, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 3120377 50 Rotation: (unrotated = principal) Rho = 0.7105 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 221712 131387 0.7105 0.7105 Comp2 | 0903255 0.2895 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S61 | 0.1292 | 08882 S62 | 0.9916 | 001508 - pca S71 S72 S73, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 2078756 Rotation: (unrotated = principal) Rho = 0.5858 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 121783 0425344 0.5858 0.5858 Comp2 | 0792487 072405 0.3812 0.9671 Comp3 | 00684375 0.0329 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S71 | 0.1320 | 007152 S72 | 0.7039 | 03987 S73 | 0.6979 | 03907 - pca S81 S82 S83 S84, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 7177683 Rotation: (unrotated = principal) Rho = 0.5844 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 419491 249942 0.5844 0.5844 Comp2 | 169549 0615491 0.2362 0.8207 Comp3 | 108 0872725 0.1505 0.9711 Comp4 | 0207278 0.0289 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S81 | 0.0687 | 1293 S82 | 0.7003 | 01126 S83 | 0.7095 | 01079 S84 | -0.0385 | 147 - pca S91 S92 S93 S94 S95 S96, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 4295052 Rotation: (unrotated = principal) Rho = 0.5816 51 Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 249783 138562 0.5816 0.5816 Comp2 | 11122 0784568 0.2590 0.8405 Comp3 | 0327637 0046476 0.0763 0.9168 Comp4 | 0281161 0240125 0.0655 0.9823 Comp5 | 00410358 00058496 0.0096 0.9918 Comp6 | 00351862 0.0082 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S91 | -0.0370 | 03139 S92 | 0.9910 | 0006493 S93 | 0.1047 | 0454 S94 | 0.0099 | 003541 S95 | -0.0741 | 09452 S96 | -0.0075 | 004225 - pca S101 S102 S103 S104, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 6548799 Rotation: (unrotated = principal) Rho = 0.4959 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 324769 177089 0.4959 0.4959 Comp2 | 14768 0259314 0.2255 0.7214 Comp3 | 121748 0610654 0.1859 0.9073 Comp4 | 0606829 0.0927 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + S101 | 0.5527 | 09726 S102 | 0.4415 | 09538 S103 | 0.4518 | 0791 S104 | 0.5436 | 05838 4.2 Tính pca 03 nhân tố lvi-ipcc (khả thích ứng, độ nhạy cảm mức độ phơi nhiễm) pca M1c M2c M3c M4c, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 1.150437 Rotation: (unrotated = principal) Rho = 0.8038 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 92475 833688 0.8038 0.8038 Comp2 | 0910618 0100418 0.0792 0.8830 Comp3 | 08102 0274142 0.0704 0.9534 Comp4 | 0536057 0.0466 1.0000 -Principal components (eigenvectors) 52 Variable | Comp1 | Unexplained -+ + M1c | 0.1443 | 08241 M2c | 0.9786 | 003585 M3c | 0.1345 | 0558 M4c | 0.0587 | 08389 - pca M5c M6c M7c M8c, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 6664765 Rotation: (unrotated = principal) Rho = 0.3643 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 242828 041687 0.3643 0.3643 Comp2 | 201141 0248192 0.3018 0.6661 Comp3 | 176322 130136 0.2646 0.9307 Comp4 | 0461855 0.0693 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + M5c | 0.9175 | 03137 M6c | 0.0322 | 1763 M7c | 0.1123 | 0487 M8c | 0.3802 | 1673 - pca M9c M10c, components(1) covariance Principal components/covariance Number of obs = 426 Number of comp = Trace = 3384791 Rotation: (unrotated = principal) Rho = 0.7584 -Component | Eigenvalue Difference Proportion Cumulative -+ -Comp1 | 256705 174932 0.7584 0.7584 Comp2 | 0817737 0.2416 1.0000 -Principal components (eigenvectors) -Variable | Comp1 | Unexplained -+ + M9c | 0.9992 | 0001262 M10c | 0.0393 | 08165 53 PHỤ LỤC MƠ HÌNH HỒI QUY 5.1 Mơ hình hồi quy reg LVIce GIOI HVAN LD TYLE HONGHEO dientich NHA THUNN THUPNN TINDUNG TCTT TCXH HTDP M9 M10 chudong ungpho Source | SS df MS -+ -Model | 13.2355721 17 778563064 Residual | 4.29956807 408 010538157 -+ -Total | 17.5351402 425 041259153 Number of obs F(17, 408) Prob > F R-squared Adj R-squared Root MSE = = = = = = 426 73.88 0.0000 0.7548 0.7446 10266 -LVIce | Coef Std Err t P>|t| [95% Conf Interval] -+ -GIOI | 0703706 0124974 5.63 0.000 0458033 0949379 HVAN | -.0077211 0014858 -5.20 0.000 -.0106419 -.0048003 LD | 0176419 0046878 3.76 0.000 0084266 0268571 TYLE | 0160932 008962 1.80 0.073 -.0015244 0337107 HONGHEO | 0264328 0156868 1.69 0.093 -.0044043 0572699 dientich | -.0024424 0006863 -3.56 0.000 -.0037914 -.0010933 NHA | -.0252397 0108298 -2.33 0.020 -.0465289 -.0039505 THUNN | -.0005662 0001398 -4.05 0.000 -.000841 -.0002915 THUPNN | 0007598 0001285 5.91 0.000 0005073 0010123 TINDUNG | 0393671 0106663 3.69 0.000 0183994 0603348 TCTT | -.0693518 0109771 -6.32 0.000 -.0909306 -.047773 TCXH | -.037559 0072589 -5.17 0.000 -.0518286 -.0232894 HTDP | -.0191055 0110014 -1.74 0.083 -.040732 002521 M9 | 537839 0427332 12.59 0.000 4538344 6218436 M10 | 1071765 0184901 5.80 0.000 0708286 1435243 chudong | -.1228537 0136247 -9.02 0.000 -.149637 -.0960704 ungpho | -.0684946 0196661 -3.48 0.001 -.1071541 -.029835 _cons | 4896093 0299541 16.35 0.000 4307256 5484929 - reg LVIce GIOI HVAN LD TYLE HONGHEO dientich NHA THUNN THUPNN TINDUNG TCTT TCXH HTDP M9 M10 chudong ungpho, beta Source | SS df MS -+ -Model | 13.2355721 17 778563064 Residual | 4.29956807 408 010538157 -+ -Total | 17.5351402 425 041259153 Number of obs F(17, 408) Prob > F R-squared Adj R-squared Root MSE = = = = = = 426 73.88 0.0000 0.7548 0.7446 10266 -LVIce | Coef Std Err t P>|t| Beta -+ -GIOI | 0703706 0124974 5.63 0.000 1475143 HVAN | -.0077211 0014858 -5.20 0.000 -.1483867 LD | 0176419 0046878 3.76 0.000 1185824 TYLE | 0160932 008962 1.80 0.073 0547496 HONGHEO | 0264328 0156868 1.69 0.093 0433463 dientich | -.0024424 0006863 -3.56 0.000 -.1418601 NHA | -.0252397 0108298 -2.33 0.020 -.0621684 THUNN | -.0005662 0001398 -4.05 0.000 -.1656304 THUPNN | 0007598 0001285 5.91 0.000 1599836 TINDUNG | 0393671 0106663 3.69 0.000 0941962 TCTT | -.0693518 0109771 -6.32 0.000 -.1646806 TCXH | -.037559 0072589 -5.17 0.000 -.1432437 HTDP | -.0191055 0110014 -1.74 0.083 -.0464809 M9 | 537839 0427332 12.59 0.000 3250297 M10 | 1071765 0184901 5.80 0.000 1497642 chudong | -.1228537 0136247 -9.02 0.000 -.2917246 ungpho | -.0684946 0196661 -3.48 0.001 -.0995678 _cons | 4896093 0299541 16.35 0.000 54 5.3 Kiểm định phù hợp mơ hình sktest LVIce Skewness/Kurtosis tests for Normality joint -Variable | Obs Pr(Skewness) Pr(Kurtosis) adj chi2(2) Prob>chi2 -+ LVIce | 426 0.1541 0.0214 7.10 0.0288 predict residual, resid sktest residual Skewness/Kurtosis tests for Normality joint -Variable | Obs Pr(Skewness) Pr(Kurtosis) adj chi2(2) Prob>chi2 -+ residual | 426 0.5575 0.9348 0.35 0.8378 vif Variable | VIF 1/VIF -+ -THUNN | 2.78 0.359590 dientich | 2.64 0.378237 chudong | 1.74 0.574162 LD | 1.65 0.605298 TYLE | 1.55 0.646494 ungpho | 1.36 0.735349 HVAN | 1.36 0.737069 TCXH | 1.28 0.784125 THUPNN | 1.22 0.821428 HTDP | 1.19 0.838935 NHA | 1.18 0.844580 GIOI | 1.14 0.875651 TCTT | 1.13 0.884519 M10 | 1.11 0.900240 M9 | 1.11 0.901123 HONGHEO | 1.10 0.908169 TINDUNG | 1.08 0.922635 -+ -Mean VIF | 1.45 estat hettest Breusch-Pagan / Cook-Weisberg test for heteroskedasticity Ho: Constant variance Variables: fitted values of LVIce chi2(1) Prob > chi2 = = 0.45 0.5034 55