Tác động của thảm họa thiên nhiên đối với an ninh lương thực giai đoạn 2010-2016 tại vùng nông thôn Việt Nam: Kết quả từ phương pháp ước lượng bán tham số GAM

Bài viết Tác động của thảm họa thiên nhiên đối với an ninh lương thực giai đoạn 2010-2016 tại vùng nông thôn Việt Nam: Kết quả từ phương pháp ước lượng bán tham số GAM nghiên cứu mối quan hệ giữa thu nhập và an ninh lương thực (đo bằng thực phẩm tiêu thụ) ở cấp độ hộ gia đình, dưới tác động của các hiện tượng thời tiết khắc nghiệt giai đoạn 2010-2016. Thảm họa thiên nhiên, bao gồm: hạn hán, lũ lụt và bão, được tính bằng thời gian (tháng) từ lúc xảy ra thảm họa đến thời điểm điều tra. Mô hình ước lượng bán tham số GAM được sử dụng trên bộ số liệu điều tra mức sống dân cư VHLSS. Mời các bạn cùng tham khảo!

TÁC ĐỘNG CỦA THẢM HỌA THIÊN NHIÊN ĐỐI VỚI AN NINH LƯƠNG THỰC GIAI ĐOẠN 2010-2016 TẠI VÙNG NÔNG THÔN VIỆT NAM: KẾT QUẢ TỪ PHƯƠNG PHÁP ƯỚC LƯỢNG BÁN THAM SỐ GAM Trịnh Thị Hường1*, Michel Simioni2 Bộ mơn Tốn, Đại học Thương mại INRA, UMR-1110 MOISA, place Pierre Viala - Bât 26, 34 060 Montpellier Cedex 2, France Tóm tắt Bài báo nghiên cứu mối quan hệ thu nhập an ninh lương thực (đo thực phẩm tiêu thụ) cấp độ hộ gia đình, tác động tượng thời tiết khắc nghiệt giai đoạn 2010-2016 Thảm họa thiên nhiên, bao gồm: hạn hán, lũ lụt bão, tính thời gian (tháng) từ lúc xảy thảm họa đến thời điểm điều tra Mơ hình ước lượng bán tham số GAM sử dụng số liệu điều tra mức sống dân cư VHLSS Kết nghiên cứu cho thấy tượng lũ lụt thể tác động tức thời đến lượng thực phẩm tiêu thụ, hạn hán có tác động dài Kết báo đóng góp thêm chứng khoa học ảnh hưởng biến đối khí hậu đến an sinh xã hội Từ khóa: an ninh lương thực, bất thường thời tiết, thu nhập, lượng calo tiêu thụ, hồi quy bán tham số, VHLSS Giới thiệu Biến đổi khí hậu tượng khắc nghiệt thời tiết có ảnh hưởng đến nhiều lĩnh vực, đặc biệt hoạt động sản suất nông nghiệp lương thực Khu vực Châu Á Thái Bình Dương nằm khu vực chịu ảnh hưởng nặng nề Việt Nam nước chịu tác động mạnh thời tiết (FAO 2016) Từ đó, tượng khắc nghiệt thời tiết tác động mạnh đến an ninh lương thực – nhiệm vụ quan trọng 17 mục tiêu phát triển bền vững Liên Hợp Quốc Chính phủ Việt Nam cam kết (General 2015) Tại Việt Nam, đánh giá ảnh hưởng tượng khắc nghiệt thời tiết đến lĩnh vực kinh tế xã hội thời tiết thu hút quan tâm nhiều nhà hoạch định sách nhà khoa học (Arouri, Nguyen, and Youssef 2015; Narloch 2016; Diallo, Marchand, and Espagne 2019) Các nghiên cứu đánh giá tác động tầm vi mô, tức yếu tố cấp độ hộ gia đình thu nhập hay sản xuất nơng nghiệp Trong đó, tác giả khằng định thiên tai, bão, lũ lụt hạn hán có ảnh hưởng tiêu cực đến hoạt động sản xuất nông nghiệp, thu nhập chi tiêu hộ gia đình, bất bình đẳng thu nhập Tuy nhiên, chưa có đánh giá ảnh hưởng tượng thời tiết khắc nghiệt đến an ninh lương thực Đồng thời, nghiên cứu dừng lại việc đánh giá có/khơng có ảnh hưởng tượng thời tiết khắc nghiệt, mà chưa đưa 20 tác động theo thời gian khả phục hồi hộ gia đình sau thảm họa Xét góc độ vi mơ, có nhiều tài liệu ước lượng mối quan hệ tiêu dùng thực phẩm (đo Calo, Kcal) thu nhập hộ gia đình, dựa tảng số liệu Mức sống dân cư phương pháp tính đề xuất A Deaton1 Đánh giá mối quan hệ lượng thực phẩm tiêu thụ thu nhập có nhiều ý nghĩa Thứ nhất, kết nghiên cứu cho thấy tác động tăng trưởng thu nhập đến tiêu thụ thực phẩm, từ đó, hỗ trợ việc ban hành sách việc giải tình trạng thiếu dinh dưỡng cải thiện đa dạng phần ăn Thứ hai, đánh giá mối quan sử dụng để dự báo nhu cầu lương thực quốc gia, cung cấp thông tin an ninh lương thực tương lai Kết đánh giá có ý nghĩa quan trọng đất nước nông nghiệp Việt Nam Đồng thời, kết cần cập nhật liên tục để đưa chứng khoa học Nghiên cứu chúng tơi đóng góp thêm chứng khoa học cho mối quan hệ thu nhập hộ gia đình thực phẩm tiêu thụ, tác động tượng thời tiết khắc nghiệt yếu tố kinh tế xã hội Thứ nhất, chúng tơi đề xuất sử dụng mơ hình hồi quy ước lượng bán tham số cộng tuyến tổng quát (A generalized additive model (GAM)) để đánh giá mối quan hệ thu nhập lượng thực phẩm tiêu thụ Thứ hai, đưa thêm yếu tố thời gian, tức khoảng thời gian từ lúc chịu ảnh hưởng thảm họa đến điều tra, để xem xét tác động khả phục hồi sau thảm họa Nghiên cứu đánh giá từ số liệu Mức sống dân cư năm 2010, 2012, 2014 2016 Phương pháp nghiên cứu 2.1 Dữ liệu nghiên cứu Bộ liệu sử dụng để phân tích Điều tra mức sống dân cư 2010, 2012, 2014 2016 thực Tổng cục thống kê Việt Nam2 Đây khảo sát tiến hành 02 năm lần phạm vi nước phương pháp vấn trực tiếp chủ hộ (9000 hộ gia đình) cán chủ chốt xã (1000 xã nông thôn), từ năm 1993 Kết điều tra mức sống dân cư phủ quan sử dụng để đánh giá mức sống dân cư phục vụ hoạch định sách lập kế hoạch phát triển kinh tế - xã hội Bộ số liệu mang tính đại diện cấp vùng, cấp nơng thơn-thành thị cấp tỉnh Nghiên cứu hạn chế vùng nông thôn kết hợp phiếu điều tra cấp hộ cấp xã Sử dụng mã số đơn vị hành cấp xã, hộ gia đình gắn với thơng tin cấp xã tương ứng năm điều tra Cụ thể: Từ phiếu điều tra cấp hộ, trích xuất thơng tin kinh tế, xã hội hộ gia đình thu nhập hàng năm, thơng tin chủ hộ gia đình (giới tính, tuổi, dân tộc, trình độ giáo dục), số thành viên hộ gia đình nơi sống (vùng sinh thái) A Deaton trao giải Noel lĩnh vực kinh tế năm 2015 đóng góp ơng nghiên cứu tiêu dùng, nghèo đói an sinh xã hội https://en.wikipedia.org/wiki/Angus_Deaton Báo cáo toàn văn điều tra có website Tổng cục thống kê Việt Nam Ví dụ, kết Điều tra Mức sống dân cư 2016 có https://www.gso.gov.vn/du-lieu-va-so-lieu-thong-ke/2019/04/ket-qua-khao-sat-mucsong-dan-cu-viet-nam-nam-2016-2/ 21 Thực phẩm tiêu thụ đo lượng calo tiêu thụ bình quân/ngày sử dụng phương pháp tính truyền thống (Trinh, Simioni, and Thomas-Agnan 2018) Từ phiếu điều tra cấp xã, trích xuất thơng tin tượng thiên tai bất thường từ câu hỏi điều tra “Thiên tai, hỏa hoạn hay bệnh dịch xảy năm gần nhất” quan tâm đến kiện lũ lụt, bão hạn hán Đầu tiên, chúng tơi quan tâm xã có xảy tượng thời tiết bất thường hay không nghiên cứu của(Arouri, Nguyen, and Youssef 2015) Tiếp theo, từ thông tin ghi lại tháng/năm xảy kiện, chúng tơi tính số tháng từ lúc xảy tượng bất thường đến tháng tiến hành phiếu điều tra Bên cạnh đó, số thông tin sở hạ tầng xã như: 1) chợ 2) số dự án kết cấu hạ tầng xây mới, nâng cấp hay cải tạo tổng hợp 2.2 Phương pháp nghiên cứu Để đánh giá tác động thu nhập bình quân đến lượng calo tiêu thụ bình quân tác động tượng thời tiết bất thường yếu tố kinh tế xã hội, sử dụng mơ hình hồi quy ước lượng bán tham số cộng tuyến tổng quát (A generalized additive model (GAM)) Phương pháp Wood, Simon N phát triển (Wood 2017) Xét trường hợp mơ hình hai biến độc lập x z; biến phụ thuộc y mẫu quan sát i = 1, 2, , n Mơ hình cộng tuyến phát biểu sau y = α + f (x ) + f (z ) + ϵ , (1) Trong α hệ số chặn, f , f hàm trơn (smooth function) ϵ ∼ N(0, σ ) sai số Các hàm trơn (1) ước lượng thông qua hàm đa thức tuyến tính nhánh Ví dụ với δ hệ số đa thức sở {b (x), b (x), … , b (x)) f (x) = b (x)δ (2) b (x), j = 1, … , k hàm xây dựng từ nút x ∗ , j = 1, , k − (knots) tập giá trị biến x, x − x∗ ⎧ ∗ x ∗ < x ≤ x ∗ ∗ x − x ⎪ b (x) = x ∗ − x (3) ∗ ∗ ⎨ ∗ x < x < x ∗ ⎪x − x ⎩ khác Các hàm trơn f , f (1) (hay hệ số hàm trơn (2)) ước lượng thông qua phương pháp bình phương tối thiếu có phạt (the penalized least squares) kết hợp với ước lượng Spline3 Phương pháp kĩ thuật ước lượng xây dựng thành gói lệnh mgvc phần mềm R, thơng qua hàm ước lượng gam4 Các đường công Spline bao gồm đoạn đa thức kết nối với trơn nhẵn điểm nút Đây kĩ thuật thường dùng để làm trơn đường cong hồi quy trơn mơ hình cộng tuyến mở rộng (Green and Silverman, 1994) https://cran.r-project.org/web/packages/mgcv/mgcv.pdf 22 Sử dụng mô hình GAM trên, chúng tơi để xuất ước lượng mơ hình log(PCCI) = α + s(log(HHIN C )) + βX + s(log(HHINC), month ) + ϵ (4) Trong đó, i số hộ gia đình, PCCIlà lượng tiêu thụ calo bình quân/người/ngày, HHINClà thu nhập bình quân hộ gia đình, X thơng tin nhân học kinh tế xã hội hộ gia đình, month số tháng từ lúc xảy tượng bất thường k (như lũ lụt, hạn hán, bão) đến thời điểm điều tra ϵ sai số Kí hiệu s(x) hàm trơn (smooth function) PCCIvà HHINCđược chuyển đổi qua lơ ga rít mơ hình truyền thống ước lượng mối quan hệ (Trinh, Simioni, and Thomas-Agnan 2018) Ước lượng s(log(HHIN C ), month ) để đánh giá tác động phi tuyến đồng thời thu nhập tác động thiên tai5 Khi đó, ta nói tác động log(HHIN C ), month ) phi tuyến tác động biến cịn lại lên log(PCCI ) tuyến tính lên log(PCCI Mơ hình (4) thuộc dạng mơ hình bán tham số Các tính tốn nghiên cứu, thống kê mơ tả ước lượng mơ hình (4) tiến hành phần mềm Rstudio, phiên 4.0.2 Gói lệnh sử dụng mgvc Kết thảo luận 3.1 Đặc điểm thảm họa thiên tai vùng nông thôn Việt Nam Trong năm gần từ ngày diều tra, bão lũ lụt tượng xảy xã nhiều so với hạn hán Hiện tượng bão lũ lụt giao động khoảng 30-45% xã điều tra Đáng ý gia tăng thảm họa hạn hán vào năm 2016 so với năm trước Các xã vùng nơng thôn chứng kiến thảm họa thiên tai từ khoảng đến tháng trước thời điểm điều tra xu hướng thảm họa qua năm không đồng Số lượng cơng trình dân sinh xây dựng sửa chữa có xu hướng tăng, tương ứng với chương trình nâng cao hạ tầng để phát triển nơng thơn Số lượng xã có chợ (chợ hàng ngày, chợ phiên, chợ đầu mối) chưa cao, khoảng 25% mẫu quan sát 3.2 Đặc điểm nhân học hộ gia đình nơi sống Số lượng hộ gia đình nơng thơn mẫu điều tra có xu hướng giảm qua năm xu hướng thị hóa Thu nhập bình qn hộ gia đình tăng qua năm Về đặc điểm chủ hộ, tuổi chủ hộ cao trình độ giáo dục mức thấp Điều giải thích đặc điểm đa hệ gia đình Việt Nam chủ hộ thường hệ thứ gia đình ba hệ Qua năm, cỡ hộ (số thành viên hộ gia đình) có xu hướng giảm Nghiên cứu tập trung vùng nơng thơn nên số hộ gia đình Mơ hình tuyến tính dạng sau tác giả ước lượng cung cấp kết theo yêu cầu người đọc log(PCCI) = α + α log(HHNC ) + α log (HHINC) + + μ log (HHNC ) ∗ month + ϵ 23 βX + γ log(HHNC ) ∗ month chiếm đa số khu vực miền núi phía Bắc duyên hải miền Trung Lượng calo đầu người cao, khoảng 3500 calo/người/ngày có xu hướng tăng nhẹ vào năm 2016 sau giảm vào năm 2014 2012 Lượng calo bình quân cao mức trung bình năm mức calo tình từ thực phẩm mua, chưa qua chế biến thải bỏ (NIN 2011) Bảng 1: Thông tin thảm họa thiên tai xã nông thôn Việt Nam, gia đoạn 2010-2016 Năm 2010 Năm 2012 Năm 2014 Năm 2016 Không 53.87% 59.73% 56.64% 67.02% Có 46.13% 40.27% 43.36% 32.98% Khơng 71.45% 79.18% 80.03% 59.43% Có 28.55% 20.82% 19.97% 40.57% Xã bị bão năm gần Không 61.83% 61.44% 49.87% 56.94% Có 38.17% 38.56% 50.13% 43.06% Thời gian từ lúc thiên tai xảy đến lúc điều tra (tháng) Bão 6.3 (9.3) 6.4 (11) Chỉ tiêu Xã bị lũ lụt năm gần Xã bị hạn hán năm gần 7.7 (10.5) 6.7 (10.4) Lũ lụt (6.8) 3.9 (8.7) (7.4) 5.6 (9) Hạn hán 5.3 (8.7) 5.5 (9.7) 7.6 (9.9) 6.5 (10.5) Có 24.57% 25.90% 24.39% 23.72% Không 75.43% 74.10% 75.61% 76.28% 6.6 (2.7) 6.6 (2.6) 7.1 (2.6) (2.5) Xã có chợ hàng ngày, chợ phiên chợ đầu mối Số lượng công trình dân sinh xây dựng/sửa chữa năm Chú thích: Biến liên tục thể giá trị trung bình độ lệch chuẩn ngoặc kép, biến rời rạc thể tỉ lệ phần trăm mức độ Nguồn: Tác giả tính tốn từ số liệu Điều tra mức sống dân cư 2010 2012, 2014, 2016 Bảng 2: Đặc điểm nhân học hộ gia đình nơi sống Chỉ tiêu Năm 2010 Năm 2012 Năm 2014 Năm 2016 Số quan sát 5188 4638 3293 3958 71060.4 80180.2 85597.4 92399.5 Thu nhập bình quân hộ gia đình (nghìn đồng/năm) (54147.5) (60971.8) (65212.8) ( 84654.4) Tuổi chủ hộ (năm) 46.9 (13.9) 48.4 14.1) 50 (14.3) 51.1 (14.3) 24 3602 3477.5 3527.5 3668.2 (966.4) (1044.3) (1103.4) (1306.7) Nam 81.30% 81.22% 80.44% 76.07% Nữ 18.70% 18.78% 19.56% 23.93% Số thành viên hộ Có thành viên gia đình Có thành viên 15.38% 16.06% 19.07% 23.88% 18.47% 17.44% 18.43% 18.39% Có thành viê 31.92% 31.31% 30.09% 28.73% Có thành viên 18.14% 19.02% 17.49% 16.40% Có từ thành viên trở lên 16.09% 16.17% 14.91% 12.61% Dân tộc khác 25.69% 26.05% 23.60% 26% Kinh 74.31% 73.95% 76.40% 74% 59.08% 58.62% 53.66% 58.54% Trung học 38.63% 39.46% 43.46% 38.83% Đại học cao 2.29% 1.92% 2.88% 2.63% Đồng sông Hồng 14.11% 10.03% 13.48% 6.19% Trung du miền núi Phía Bắc 23.79% 22.19% 22.17% 14.88% Duyên hải miền Trung 31.05% 30.53% 33.62% 43.66% Tây Nguyên 8.46% 7.12% 7.74% 12.38% Đông Nam Bộ 7.02% 10.16% 8.02% 5.79% Đồng sông Cửu Long 15.57% 19.99% 14.97% 17.10% Lượng calo bình (Kcal/người/ngày) Giới tính Dân tộc chủ hộ qn đầu người Trình độ giáo Khơng cấp tiểu dục chủ hộ học Vùng sinh thái Chú thích: Biến liên tục thể giá trị trung bình độ lệch chuẩn ngoặc kép, biến rời rạc thể tỉ lệ phần trăm mức độ Nguồn: Tác giả tính tốn từ số liệu Điều tra mức sống dân cư 2010 2012, 2014, 2016 25 Bảng 3: Hệ số hồi quy mơ hình GAM Biến độc lập Hệ số hồi quy Sai số chuẩn Hệ số chặn 8.195 *** 0.053 *** Tuổi chủ hộ (năm) -0.002 *** -0.008 Giới tính chủ hộ (mặc định: Nữ) Nam Số thành viên hộ gia đình Có thành viên Có thành viê Có thành viên Có từ thành viên trở lên Dân tộc chủ hộ Kinh Trình độ giáo dục chủ hộ (mặc Trung học định: Không cấp tiểu học) Đại học cao Năm (mặc định: 2010) Năm 2012 Năm 2014 Năm 2016 0.003 0.013 -0.001 -0.034 *** -0.086 *** 0.006 -0.006 -0.092 *** -0.046 *** -0.034 *** -0.01 -0.008 0.031 0.05 -0.01 *** 0.066 *** -0.002 0.003 -0.017 -0.008 -0.009 -0.007 Chú thích: *, ** *** biểu diễn mức ý nghĩa 10%, 5% 1% Nguồn: Tác giả tính tốn từ số liệu Điều tra mức sống dân cư 2010 2012, 2014, 2016 3.3 Tác động thảm họa thiên nhiên lượng calo bình qn Kết ước lượng mơ hình GAM (4) minh họa Hình (phần tách động phi tham số) Bảng (phần tác động tham số) So với yếu tố mặc định biến độc lập Bảng 3, tác động có ý nghĩa âm Tăng tuổi chủ hộ hay tăng trình độ học vấn chủ hộ có tác động giảm lượng calo tiêu thụ Đồng thời, qua năm, so với năm 2010, tác động năm giảm Xu hướng giảm tương ứng với xu hướng chung xã hội Việt Nam đa dạng thực phẩm: Như tăng thực phẩm giàu chất đạm giảm thực phẩm giàu tinh bột (Trinh, Simioni, and Thomas-Agnan 2018; Nguyen and Hoang 2018) Hình biểu thị hàm trơn chuẩn hóa6 yếu tố phi tuyến đến lượng calo tiêu thụ giá trị dương hàm trơn chuẩn hóa thể tác động chiều giá trị âm thể tác động ngược chiều Biểu đồ a) thể tác động phi tuyến mối liên hệ Hàm trơn chuẩn hóa dựa vào ước lượng thực tế thuật toán Từ hàm trơn chuẩn hóa, tịnh tiến để xem tác động giá trị thực 26 thu nhập – lượng calo nghiên cứu trước (Trinh, Simioni, and Thomas-Agnan 2018) Kết có ý nghĩa kinh tế phát triển, đặc biệt khu vực nông thôn Tác động đồng thời bão thu nhập đến lượng calo khơng có ý nghĩa thống kê nên chúng tơi khơng phân tích thêm Giá trị hàm trơn chuẩn hóa lũ lụt (hình b)) có giá trị dương giá trị tương ứng hạn hán (hình c)) có tác động âm Chi tiết hơn, xét tác động khoảng thời gian từ lúc xảy kiện đến lúc điều tra, với mức thu nhập, số tháng tăng giá trị hàm trơn tăng, tức hộ gia đình có tăng lượng calo Tuy nhiên, số tháng gần 0, giá trị hàm trơn có xu hướng nhỏ 0, thể tác động tức thời lũ lụt đến an ninh lương thực Ngược lại, hạn hán có tác động dài hạn hơn, địa phương xảy hạn hán thời điểm điều tra hàm trơn có xu hướng 0, tức chưa có tác động tức thời, nhiên, thời điểm xảy tăng dần so với thời điểm điều tra, hàm trơn có xu hướng âm Hình 1: Ước lượng hàm mối liên hệ logarit lượng calo tiêu thụ logarit thu nhập tác động thời gian xảy thiên tai Kết ước lượng a), b), c) có mức ý nghĩa tương ứng 1%, 5% 1% Ước lượng d) khơng có ý nghĩa thống kê Nguồn: Tác giả tính toán từ số liệu Điều tra mức sống dân cư 2010 2012, 2014, 2016 Kết luận Nghiên cứu tập trung vào đánh giá tác động thu nhập hộ gia đình tượng bất thường thời tiết, yếu tố kinh tế xã hội đến lượng calo tiêu thụ khu vực nông thôn Việt Nam, giai đoạn 2010-2016 Các điều kiện kinh tế, giáo dục hạ tầng Việt Nam tăng trưởng qua năm có tác động tích cực đến lượng tiêu thụ calo Việt Nam chịu nhiều tác động bất thường thời tiết, đặc biệt khu vực nơng thơn Bão thu nhập hộ gia đình khơng có tác động có ý nghĩa thống kê đến lượng calo tiêu thụ Hiện tượng lũ lụt thể tác động tức thời đến lượng calo tiêu thụ hộ gia đình có xu hướng phục hồi thời gian thu nhập tăng lên Ngược lại, hạn hán có tác động dài hạn lượng calo tiêu thụ Tác động hạn hán đáng lo ngại, đặc biệt năm 2016 Việt Nam bị tác động El Nino (FAO 2016), tượng có tác động đến nơng nghiệp an ninh lương thực Kết báo đưa thêm 27 chứng cho sách an ninh lương thực tăng trưởng kinh tế bối cảnh biến đổi khí hậu thay đổi mạnh mẽ Chủ đề nên quan tâm nghiên cứu ứng dụng công cụ thống kê nâng cao Tài liệu tham khảo Arouri, Mohamed, Cuong Nguyen, and Adel Ben Youssef 2015 “Natural Disasters, Household Welfare, and Resilience: Evidence from Rural Vietnam.” World Development 70: 59–77 Diallo, Yoro, Sébastien Marchand, and Etienne Espagne 2019 “Impacts of Extreme Climate Events on Technical Efficiency in Vietnamese Agriculture,” no March FAO 2016 “‘El Nino’ Event in Vietnam: Agriculture Food Security and Livelihood Needs Assessment in Response to Drought and Sal Water Instrusion,” 75 http://www.fao.org/documents/card/en/c/7f0faa29-6ec6-47b8-a77cb7ff35c745da/?utm_source=facebook&utm_medium=social+media&utm_campai gn=fao+facebook General, UNG Assembly 2015 “Sustainable Development Goals.” SDGs Transform Our World 2030 Narloch, Ulf 2016 “The Varying Income Effects of Weather Variation: Initial Insights from Rural Vietnam.” The Varying Income Effects of Weather Variation: Initial Insights from Rural Vietnam, no July https://doi.org/10.1596/1813-9450-7764 Nguyen, Tuan T., and Minh V Hoang 2018 “Non-Communicable Diseases, Food and Nutrition in Vietnam from 1975 to 2015: The Burden and National Response.” Asia Pacific Journal of Clinical Nutrition 27 (1): 19–28 NIN, Viện dinh dưỡng quốc gia 2011 Tình Hình Dinh Dưỡng Việt Nam Năm 20092010 Nhà xuất Y học Trinh, Thi Huong, Michel Simioni, and Christine Thomas-Agnan 2018 “Assessing the Nonlinearity of the Calorie-Income Relationship: An Estimation Strategy – With New Insights on Nutritional Transition in Vietnam.” World Development 110 Wood, Simon N 2017 Generalized Additive Models An Introduction with R CRC press 28 ... Nguồn: Tác giả tính tốn từ số liệu Điều tra mức sống dân cư 2010 2012, 2014, 2016 3.3 Tác động thảm họa thiên nhiên lượng calo bình quân Kết ước lượng mơ hình GAM (4) minh họa Hình (phần tách động. .. thụ Tác động hạn hán đáng lo ngại, đặc biệt năm 2016 Việt Nam bị tác động El Nino (FAO 2016), tượng có tác động đến nông nghiệp an ninh lương thực Kết báo đưa thêm 27 chứng cho sách an ninh lương. .. f (1) (hay hệ số hàm trơn (2)) ước lượng thơng qua phương pháp bình phương tối thiếu có phạt (the penalized least squares) kết hợp với ước lượng Spline3 Phương pháp kĩ thuật ước lượng xây dựng