tiểu luận kinh tế lượng nghiên cứu các nhân tố ảnh hưởng đến sản lượng lúa của 15 tỉnh thành giai đoạn 2007 2016

MỤC LỤC Lời mở đầu 1.1 Lí chọn đề tài 1.2 Mục tiêu nghiên cứu 2 Cơ sở lý thuyết 3 Phương pháp nghiên cứu 3.1 Mơ hình nghiên cứu 3.2 Biến số thước đo 3.3 Dữ liệu .4 Kết nghiên cứu thảo luận 4.1 4.1.1 Mô tả thống kê biến .5 4.1.2 Mô tả tương quan biến 4.2 Mô tả thống kê mô tả tương quan biến Kết ước lượng thảo luận Kết luận hàm ý sách 5.1 Kết luận 5.2 Hàm ý sách………………………………………………………………10 TÀILIỆUTHAMKHẢO………………………………………………………… 11 PHỤLỤC…………………………………………………………………………… 12 DANH MỤC BẢNG Bảng 3.2 Bảng miêu tả kỳ vọng Bảng 4.1 Mô tả thống kê biến Bảng 4.2 Mô tả tương quan biến Bảng 4.2a Bảng kết hồi quy MH0 Bảng 4.2b Bảng kết hồi quy MH1 Bảng 4.2c : Kết ước lượng Bảng 4.2d Bảng kết hồi quy MH2 1 Lời mở đầu 1.1 Lí chọn đề tài Như biết, kinh tế nước ta đa dạng phong phú, phụ thuộc vào điều kiện khí hậu ưu điểm vùng để phát triển Tuy nhiên, nhìn chung ngành nơng nghiệp ngày phát triển có giúp đỡ máy móc cơng nghiệp mà đặc biệt lúa nước Điều chứng Việt Nam nước xuất gạo đứng thứ giới Do vậy, nói, thị trường tiềm năng, vừa hội, vừa thách thức cho nhà sản xuất Làm để làm tăng suất chất lượng lúa, từ đưa giải pháp để khắc phục phát huy ưu điểm? Nhận thức tính cấp thiết việc trả lời câu hỏi trên, nhóm chúng em lựa chọn đề tài: “ Nghiên cứu nhân tố ảnh hưởng đến sản lượng lúa 15 tỉnh thành giai đoạn 2007-2016” để phân tích rõ yếu tớ có ảnh hưởng đến sản lượng Để từ đó, ta thấy tương quan yếu tố không chi việc sản xuất lúa nước nói riêng mà còn vấn đề lượng cung cầu khác kinh tế nói chung 1.2 Mục tiêu nghiên cứu Thông qua việc xem xét mức độ ảnh hưởng nhân tố đến sản lượng lúa 15 tinh thành gai đoạn 2007-2016, từ đưa đề xuất nhằm cải thiện sản lượng lúa đến mức có thể, để nâng cao giá trị sản lượng thu nhập cho người trồng lúa Bài báo cáo sử dụng mơ hình kinh tế lượng để thực mục tiêu nghiên cứu Các sớ liệu xử lí logarit hóa để xem xét biến động yếu tố ảnh hưởng đến sản lượng lúa Ở đây, chúng em sử dụng phần mềm STATA để phân tích biến động Do quy mô thời gian hạn chế nên đề tài nghiên cứu nhóm em chi tập trung đề cập đến yếu tớ ảnh hường đến lượng lúa như: Sản lượng, diện tích, suất trung bình, lượng mưa trung bình, sơ nắng trung bình độ ẩm khơng khí trung bình Do hạn chế mặt chun mơn, q trình thực tiểu luận khơng thể tránh khỏi sai sót, nhóm tác giả mong nhận đóng góp để tiểu luận hoàn thiện Cơ sở lý thuyết nhân tố ảnh hưởng đến sản lượng lúa địa phương Dựa theo lý thuyết mà nhóm tìm kiếm tổng hợp được, có nhân tố ảnh hưởng đến sản lượng lúa địa phương: 2.1 Diện tích Đất nơng nghiệp khu vực đất thích hợp cho sản xuất, canh tác nơng nghiệp Ngồi ra, diện tích đất nơng nghiệp còn phụ thuộc vào sách kinh tế q́c gia Những q́c gia có diện tích đất nơng nghiệp lớn có hội khả tăng sản lượng nông nghiệp tốt Biến lấy từ mơ hình Kavi Kumar đề tài nghiên cứu “Mức độ nhạy cảm nông nghiệp Ấn Độ đối với môi trường” năm 2009 2.2 Năng suất lúa trung bình Năng suất lúa trung bình tính tổng sản lượng lúa vụ chiêm vụ mùa vụ chia cho diện tích canh tác Năng suất lúa cao sản lượng lớn với diện tích canh tác 2.3 Lượng mưa Lượng mưa đại lượng thể mức độ mưa nhiều hay ít, đo độ sâu nước mưa thu bề mặt phẳng, đo đơn vị milimet Lượng mưa yếu tố quan trọng với suất lúa sản lượng lúa Lượng mưa khan trồng khơ héo, lượng mưa q nhiều trồng bị ngập úng, khơng phát triển được, theo “Đánh giá ảnh hưởng yếu tớ khí tượng thủy văn sản xuất nơng nghiệp đến suất lúa vùng đê bao lửng tinh An Giang” (nhóm tác giả Nguyễn Thị Mỹ Hạnh) năm 2012 “Nghiên cứu ảnh hưởng số yếu tố khí hậu tới việc bớ trí hệ thớng trớng huyện Thạch Thất, thành phố Hà Nội” năm 2014 Biến lấy từ mơ hình ba nhà khoa học thuộc đại học công nghệ Tshwane, Nam Phi 2.4 Số nắng Số nắng tổng thời lượng nắng tính tháng, năm tùy theo mục đích tổng hợp sớ liệu Sớ nắng có ý nghĩa quan trọng đối với lúa, đặc biệt vào kỳ trỗ bông, làm hạt…làm tăng giảm sản lượng lúa Các tháng có sớ nắng thấp năm vào đầu vụ lúa xuân nên diện tích gieo trồng lúa xuân thấp nhiều so với diện tích trồng lúa mùa, suất lúa xuân lại cao so với lúa mùa, thời kỳ trỗ bông, làm hạt lúa xuân vào tháng 4, tháng có sớ nắng cao thích hợp, đảm bảo đủ lượng tích ơn lúa so với thời kỳ ở lúa mùa tháng 10, tháng 11 có sớ nắng thấp Mặt khác, tùy theo số nắng mà thời vụ gieo trồng thay đổi theo năm (Theo luận văn thạc sĩ “Tác động biến đổi khí hậu đến ngành trồng lúa tinh Lào Cai đề xuất giải pháp ứng phó”, Dặng Thị Thanh Hoa, ĐHQG, 2013) 2.5 Độ ẩm Độ ẩm khơng khí ảnh hưởng lớn tới độ ẩm đất Để thỏa mãn sinh trưởng cây, đất cần độ ẩm sẵn có để trồng hút dễ dàng Mỗi lồi có kích thước rễ khác khả hút nước khác nhau, nhu cầu độ ẩm khác Biến lấy từ mô hình Kavi Kumar đề tài nghiên cứu “Mức độ nhạy cảm nông nghiệp Ấn Độ đối với môi trường” năm 2009 Phương pháp nghiên cứu 3.1 Mơ hình nghiên cứu Sau thực hồi quy liệu theo dạng: dạng hàm tuyến tính dạng hàm Cobb-Douglass Kết cho thấy hàm Cobb-Douglass có kết đáng tin cậy mặt thớng kê Hàm Cobb-Douglass có dạng: Yi = AXα1Xα2…Xkαkeui Trong đó: Yi biến phụ thuộc Xi ( i=1,2,…,k) biến độc lập αi ( i = 1,2,…,k) hệ số hồi quy A số, ui sai số Dựa số liệu sở lý thuyết nêu trên, nhóm tác giả đưa mơ hình hồi quy tổng thể: lnSL = + lnDT + lnNS + .lnLM + lnSGN + lnDM + ui (MH0) Trong đó: SL: Lượng lúa sản xuất trung bình 15 tinh thành (đơn vị: nghìn tấn) DT: Diện tích lượng lúa (đơn vị: nghìn ha) NS: Hiệu lao động trình sản xuất (đơn vị: tạ /ha) LM: Lượng mưa trung bình năm (đơn vị: mm) SGN: Sớ nắng trung bình khoảng thời gian năm (đơn vị: giờ) DM: Độ ẩm trung bình khoảng thời gian năm ( đơn vị: %) 3.2 Biến số thước đo Mô tả kỳ vọng: Theo lý thuyết nêu trên, diện tích trồng lúa tinh suất lúa TB có ảnh hưởng chiều đến sản lượng lúa TB tinh Từ nhóm tác giả đưa kỳ vọng thể bảng 3.2 sau đây: Bảng 3.2 Bảng miêu tả kỳ vọng STT Tên biến lnSL lnDT lnNS lnLM lnSGN Dấu kì vọng + + + + - lnDM - 3.3 Dữ liệu Dữ liệu biến số trích xuất từ bảng thớng kê trang web Tổng cục Thống kê mục Số liệu thống kê ( 01 Đơn vị hành chính, đất đai khí hậu, 06 Nông nghiệp, lâm nghiệp thủy sản) tổng hợp thành file liệu Kết nghiên cứu thảo luận 4.1 Mô tả thống kê mô tả tương quan biến 4.1.1 Mô tả thống kê biến Mô tả thống kê biến thể bảng 4.1 đây: Bảng 4.1 Mô tả thống kê biến Variable Obs Mean Std Dev Min lnSL 150 5.603873 0.919815 3.443618 lnDT 150 4.028302 0.878918 1.686399 lnNS 150 3.879045 0.184538 3.377588 lnLM 150 7.554225 0.293394 6.879561 lnSGN 150 7.538458 0.248198 6.813115 lnDA 150 4.397273 0.033083 4.31214 Nguồn: Nhóm tác giả tự tổng hợp từ liệu lấy Tổng cục thống kê phần mềm Stata 12 Max 7.106852 5.332235 4.131961 8.407601 7.985416 4.475061 với trợ giúp Dữ liệu bảng 4.1 cho biết: - Số quan sát mô hình 150, khơng có biến bị quan sát - Giá trị TB biến lnSL 5.6 suy sản lượng TB 15 tinh giai đoạn 2007-2016 e5.6 với giá trị xấp xi 270.4 nghìn với giá trị sản lượng nhỏ e3.44 xấp xi 31.19 nghìn giá trị sản lượng nhỏ e 7.1 xấp xi 1211.97 nghìn Điều cho thấy sản lượng lúa TB tinh ổn định qua năm tinh với có chênh lệch đáng kể có yếu tớ tác động diện tích, suất… 4.1.2 Mơ tả tương quan biến Mô tả tương quan biến thể bảng 4.2 đây: Bảng 4.2 Mô tả tương quan biến lnSL lnDT lnNS lnLM lnSGN lnDA lnSL 0.9801 0.3183 -0.0819 -0.3749 -0.0566 lnDT lnNS lnLM lnSGN lnDA 0.1249 -0.0973 -0.3404 -0.0646 0.068 -0.2611 0.028 -0.0834 0.5228 -0.1947 Dữ liệu bảng 4.2 cho biết:  Về mức độ tương quan biến độc lập đối với biến phụ thuộc:  r( lnSL, lnDT ) = 0.9801 > nên biến lnDT tương quan chiều với biến lnSL mức độ tương quan cao  r( lnSL, lnNS) = 0.3183 > nên biến lnNS tương quan chiều với biến lnSL mức độ tương quan trung bình  r( lnSL, lnLM) = -0,0819 < nên biến lnLM tương quan ngược chiều với biến lnSL mức độ tương quan thấp  r( lnSL, lnSGN) = -0,3749 < nên biến lnSGN tương quan ngược chiều với biến lnSL mức độ tương quan trung bình  r( lnSL, lnDA) = -0,0566 < nên biến lnDA tương quan ngược chiều với biến lnSL mức độ tương quan thấp  Về mức độ tương quan biến độc lập: nhìn chung khơng cao, chi có biến lnLM biến lnDA tương quan cao với r( lnLM, lnDA) = 0.5228 4.2 Kết ước lượng thảo luận Ta có mơ hình sau: ln(SL) = ln(DT) + ln(NS) + ln(LM) + ln(SGN) + ln(DA) + ei (mh0) Hồi quy mô hình (mh0) reg lnSL lnDT lnNS lnLM lnSGN lnDA Bảng 4.2a Bảng kết hồi quy MH0 lnSL lnDT lnNS lnLM lnSGN lnDA _cons Coef 1.00041 0.99579 -0.01023 0.0102 0.05146 -2.51475 Với mức ý nghĩa Std Err 0.0017 0.0078 0.0055 0.0062 0.0502 0.2241 t 587.44 127.29 -1.83 1.62 1.02 -11.22 P>t 0.000 0.000 0.069 0.107 0.308 0.000 [90% Conf Interval] 0.9970 1.00378 0.9803 1.01125 -0.0212 0.0008 -0.0022 0.02263 -0.0479 0.15086 -2.9578 -2.0716 : Biến lnDT có p-value < => Biến lnDT có ý nghĩa thớng kê Biến lnNS có p-value < => Biến lnNS có ý nghĩa thớng kê Biến lnLM có p-value < => Biến lnLM có ý nghĩa thớng kê Biến lnSGN có p-value > => Biến lnSGN khơng có ý nghĩa thớng kê, loại khỏi mơ hình Biến lnDA có p-value > => Biến lnDA khơng có ý nhĩa thớng kê, loại khỏi mơ hình Ta có mơ hình mới: ln(SL) = ln(DT) + ln(NS) + ln(LM) + ei (mh1) Hồi quy mô hình (mh1) reg lnSL lnDT lnNS lnLM Bảng 4.2b Bảng kết hồi quy MH1 lnSL lnDT Coef 0.9995 Std Err 0.0016 T 621.73 P>t 0.000 [90% Conf 0.9963 Interval] 1.0026 lnNS 0.9928 0.0076 129.99 0.000 0.9777 1.0079 lnLM -0.0081 0.0048 -1.68 0.094 -0.0175 0.0014 _cons -2.2124 0.0453 -48.8 0.000 -2.3020 -2.1228 Với mức ý nghĩa : Biến lnDT có p-value < => Biến lnDT có ý nghĩa thớng kê Biến lnNS có p-value < => Biến lnNS có ý nghĩa thớng kê Biến lnLM có p-value < => Biến lnLM có ý nghĩa thớng kê Ta có bảng kết ước lượng 4.2c sau: Bảng 4.2c : Kết ước lượng Variable LNDT LNNS LNLM _cons N r2 Kiểm định bỏ sót biến Kiểm định đa cộng tuyến Kiểm định phương sai sai số thay đổi Kiểm định phân phối nhiễu MH1 0.99946 0.00161 0.9928 0.00764 -0.00807 0.00479 -2.2124 0.04534 150 0.99966 Khơng bỏ sót biến MH2 (Robust) 0.99946 0.00033 0.9928 0.01129 -0.00807 0.00421 -2.2124 0.06119 150 0.99966 Khơng có đa cộng tuyến MH1 mắc khuyết tật phương sai sai số thay đổi Đã khắc phục MH2 Phân phối nhiễu chuẩn Như vậy, kết ước lượng mơ hình thể bảng 4.2b cho thấy mơ hình có ý nghĩa hệ số β1 = 0.9995 , β2 = 0.9928 , β3 = -0.0081 khác Tuy nhiên thực phép kiểm định thể cột MH1 bảng 4.2c , mơ hình bị mắc khuyết tật phương sai sai sớ thay đổi, chuyển sang mơ hình mơ hình hồi quy Robust khắc phục nhược điểm mơ hình Kết kiểm định MH2 thể cột MH2 bảng 4.2c Mơ hình Robust (mh2) reg lnSL lnDT lnNS lnLM, robust Bảng 4.2d Bảng kết hồi quy MH2 Robust 10 lnSL lnDT lnNS lnLM _cons Coef 0.9995 0.9928 -0.008 -2.212 Std Err 0.00033 0.01129 0.00421 0.06119 T 3053.3 87.96 -1.92 -36.16 P>t 0.000 0.000 0.057 0.000 [90% Conf Interval] 0.9988094 1.0001 0.9704946 1.01511 -0.0163935 0.00026 -2.33333 -2.09148 Kết hồi quy MH2 thể bảng 4.2d Ta có mơ hình: lnSL = 0.9995lnDT + 0.9928lnNS – 0.008lnLM + ui Nhận xét: β1 = -2.212 : β2 = 0.9995 : Khi diện tích trồng lúa tăng 1% sản lượng lúa tăng 0.9995% điều kiện yếu tố khác không đổi β3 = 0.9928 : Khi suất lúa tăng 1% sản lượng lúa tăng 0.9928% điều kiện yếu tố khác không đổi β4 = -0.008 : Khi lượng mưa tăng 1% sản lượng lúa giảm 0.008% điều kiện yếu tố khác không đổi Kết luận hàm ý sách 5.1 Kết luận Từ kiểm định ta rút sớ kết luận sau: Qua phân tích mơ hình thấy sản lượng lúa chịu ảnh hưởng rõ rệt từ nhân tố chủ yếu diện tích trồng lúa suất lúa, còn lượng mưa ảnh hưởng rõ ràng đến sản lượng lúa Ngoài ra, sản lượng lúa bị ảnh hưởng bởi bệnh dịch, chế độ chăm sóc, khoa học kỹ thuật canh tác thu hoạch Tuy nhiên nhân tố chi ảnh hưởng đến lượng nhỏ sản lượng, với lượng thơng tin tìm kiếm có hạn hạn chế trình độ thành viên nhóm nên chúng em không xét đến yếu tố này, làm giảm phức tạp mơ hình Mơ hình phù hợp với lý thuyết đưa ra, khơng bỏ sót biến, khơng có đa cộng tuyến bị mắc khuyết tật phương sai sai số thay đổi khắc phục mơ hình mơ hình Robust với phân phới chuẩn nhiễu 5.2 Hàm ý sách 11 Do biến suất tác động chiều đến sản lượng nên Chính phủ, Nhà làm sách tồn thể hộ nơng dân nước nên tập trung vào việc gia tăng sản lượng Cụ thể có biện pháp sau: - Gieo trồng với mật độ phù hợp - Chú ý chăm sóc vào thời kỳ lúa trổ - Thu hoạch thời điểm Ngồi ra, diện tích trồng lúa tác động chiều đến sản lượng lúa diện tích đất nơng nghiệp nước ta có hạn, với nước ta q trình cơng nghiệp hóa, đại hóa đất nước nên việc tác động vào diện tích trồng lúa cần phải có sách thích hợp lâu dài Tóm lại, thời gian tới nước ta cần phải xây dựng định hướng phát triển ngành lúa gạo Qua q trình thảo luận, nhóm chúng em xin phép tập hợp số giải pháp sau: - Phát triển sản xuất theo chiều sâu, nâng cao giá trị hạt lúa đơn vị diện tích đất sản xuất, tạo nhiều sản phẩm có giá trị gia tăng từ sản xuất lúa gạo, tiến đến xây dựng thương hiệu - Giảm diện tích gieo trồng lúa ở nơi sản xuất khơng có lợi thế, chuyển sang trồng khác có hiệu kinh tế - Tiến hành quy hoạch hoàn thiện cơng trình đê, hệ thớng giao thơng nội đồng nhằm đẩy mạnh ứng dụng giới hóa, quy hoạch đầu tư sở hạ tầng phát triển cụm công nghiệp trọng điểm để phục vụ tốt cho nhu cầu chế biến xuất lúa gạo - Xây dựng liên kết Cánh đồng lớn gắn kết sản xuất với thu mua, chế biến, tiêu thụ lúa gạo giúp nâng cao chất lượng, hạ giá thành sản xuất, nâng cao giá trị gạo - Lấy ứng dụng khoa học - công nghệ làm khâu đột phá, tăng cường việc ứng dụng tiến khoa học kỹ thuật vào sản xuất (áp dụng chương trình 3G3T, 1P5G,…) để nâng cao chất lượng giá trị sản phẩm nông nghiệp - Chọn lọc phát huy có hiệu chủ trương, sách phát triển nơng nghiệp ban hành - Xây dựng củng cố hệ thống hợp tác xã, tổ hợp tác vững mạnh nhằm xây dựng chuỗi liên kết chặt chẽ doanh nghiệp - nông dân Tăng cường tổ chức liên kết nông dân với nông dân, sản xuất theo chuỗi giá trị gia tăng 12 TÀI LIỆU THAM KHẢO GS.TS Nguyễn Quang Dong, PGS.TS Nguyễn Thị Minh, 2015, “ Giáo trình Kinh tế lượng”, NXB Đại học Kinh tế quốc dân PGS.TS Nguyễn Quang Dong, 1998, Bài tập Kinh tế lượng, NXB Khoa học kỹ thuật 13 PGS.TS Nguyễn Cao Văn, TS.Trần Thái Ninh, 2006, Giáo trình Lý thuyết xác suất thớng kê tốn, NXB Thớng kê https://www.gso.gov.vn/ (Tổng cục thớng kê) https://ongbien.vn/kien-thuc-nong-nghiep/cac-bien-phap-de-tang-nang-suat-lua-dongian-va-hieu-qua-14661dt.html 14 PHỤ LỤC Kiểm định MH1: Kiểm định bỏ sót biến: estat ovtest Ramsey RESET test using powers of the fitted values of lnSL Ho: model has no omitted variables F(3, 143) = 0.33 Prob > F = 0.8009 p-value > => Không bác bỏ Ho hay mô hình MH1 khơng bỏ sót biến Kiểm định đa cộng tuyển vif Variable lnDT lnNS lnLM Mean VIF VIF 1.03 1.02 1.02 1.02 1/VIF 0.97316 0.977911 0.984008 VIF < 10 => Mơ hình mh1 khơng có đa cộng tuyến Kiểm định phương sai sai số thay đổi estat hettest Breusch-Pagan / Cook-Weisberg test for heteroskedasticity Ho: Constant variance Variables: fitted values of LNSL chi2(1) = 10.15 Prob > chi2 = 0.0014 15 p-value < => Bác bỏ Ho hay mơ hình mh1 mắc khuyết tật phương sai sai sớ thay đổi Khắc phục: Dùng mơ hình Robust (mh2) reg lnSL lnDT lnNS lnLM, robust Robust lnSL lnDT lnNS lnLM _cons Coef 0.9995 0.9928 -0.008 -2.212 Std Err 0.00033 0.01129 0.00421 0.06119 t 3053.3 87.96 -1.92 -36.16 P>t 0 0.057 [90% Conf 0.9988094 0.9704946 -0.0163935 -2.33333 Interval] 1.0001 1.01511 0.00026 -2.09148 Kiểm định MH2: Kiểm định bỏ sót biến: estat ovtest Ramsey RESET test using powers of the fitted values of LNSL Ho: model has no omitted variables F(3, 143) = 0.33 Prob > F = 0.8009 Ta có: p-value = 0.8009 > => Khơng bác bỏ Ho hay mơ hình MH1 khơng bỏ sót biến Kiểm định đa cộng tuyến: Variable VIF 1/VIF LNDT LNNS LNLM Mean VIF 1.03 1.02 1.02 1.02 0.97316 0.977911 0.984008 16 VIF < 10 => MH1 khơng có đa cộng tuyến Kiểm định phân phối nhiễu: predict ei sktest ei Skewness/Kurtosis tests for Normality - joint Variable Obs Pr(Skewness) Pr(Kurtosis) adj chi2(2) Prob>chi2 ei 0.1967 2.49 150 p-value ei > 0.3771 0.2885 => Mơ hình MH2 có phân phối nhiễu chuẩn 17 ... “ Nghiên cứu nhân tố ảnh hưởng đến sản lượng lúa 15 tỉnh thành giai đoạn 2007- 2016 để phân tích rõ yếu tớ có ảnh hưởng đến sản lượng Để từ đó, ta thấy tương quan yếu tớ khơng chi việc sản. .. xuất lúa nước nói riêng mà còn vấn đề lượng cung cầu khác kinh tế nói chung 1.2 Mục tiêu nghiên cứu Thông qua việc xem xét mức độ ảnh hưởng nhân tố đến sản lượng lúa 15 tinh thành gai đoạn 2007- 2016, ... lượng lúa chịu ảnh hưởng rõ rệt từ nhân tớ chủ yếu diện tích trồng lúa suất lúa, còn lượng mưa ảnh hưởng rõ ràng đến sản lượng lúa Ngồi ra, sản lượng lúa bị ảnh hưởng bởi bệnh dịch, chế độ