Nghiên cứu mức độ ảnh hưởng của các yếu tố lượng phân hữu cơ, lượng đạm bón và công lao động đến năng suất lúa
I MỞ ĐẦU Chuế Lưu xã nông huyện Hạ Hịa, tỉnh Phú Thọ, có truyền thống canh tác lúa nước từ lâu đời, thu nhập từ trồng lúa Năng suất lúa hộ xã năm gần xã tăng dần lên đầu tư nhiều yếu tố đầu vào Có nhiều yếu tố ảnh hưởng đến suất lúa, lượng phân bón, cơng lao động, chất lượng giống, tính chất đất canh tác, khí hậu, khoa học công nghệ áp dụng… Năng suất lúa tiêu quan trọng để đánh giá hiệu sản xuất lúa Do em tiến hành nghiên cứu mức độ ảnh hưởng yếu tố lượng phân hữu cơ, lượng đạm bón cơng lao động đến suất lúa để thấy ảnh hưởng yếu tố Từ đưa dự đoán, đề xuất nhằm tăng suất lúa tới mức cao Từ nâng cao giá trị sản xuất thu nhập cho người trồng lúa Đó lý em sử dụng mơ hình hồi quy tốn học để phân tích ảnh hưởng yếu tố đến suất lúa II PHƯƠNG PHÁP NGHIÊN CỨU Bước 1: Thu thập số liệu Tổng hợp số liệu điều tra từ hộ nông dân STT 10 11 12 Năng suất lúa 180 180 160 190 200 220 180 170 170 160 230 200 Lượng phân hữu 300 350 400 400 450 500 350 350 300 450 600 500 Lượng đạm 6 6.5 10 6.5 5.5 11 Công lao động 5 6 13 14 15 16 17 18 19 20 21 22 23 24 25 200 230 180 170 180 160 210 190 200 150 250 170 160 450 500 350 400 400 300 550 500 500 300 400 400 300 7.5 12 6 7.5 14 6.5 5 8 10 10 Bước 2: Thiết lập mối quan hệ lượng phân hữu bón, lượng đạm bón, cơng lao động với suất lúa qua mơ hình hồi quy sau, sử dụng mơ hình Cobb – Douglas: Yi = A X1a1 X2a X3a3 eui Lấy ln hai vế ta được: lnYi = lnA + a1 lnX1 + a2 lnX2 + a3 lnX3 + ui Hay : lnYi = a0 + a1 lnX1 + a2 lnX2 + a3 lnX3 + ui Trong đó: Yi: suất lúa (kg/ sào) X1: lượng phân hữu bón ( kg/sào) X2: lượng đạm bón (kg/ sào) X3: công lao động (người/ ngày/ sào) a0: hệ số tự do( = lnA) a1, a2 , a3 hệ số ảnh hưởng Xi đến đến Yi tương ứng ui:: sai số mơ hình III KẾT QUẢ NGHIÊN CỨU Qua điều tra 25 hộ ta thấy suất trung bình 187.6kg/ sào, lượng phân hữu bón trung bình 412kg/ sào, phân đạm bón trung bình 7.44kg/ sào, cơng lao động đầu tư trung bình 6.28 người/ngày/ sào Kết chạy mơ hình exel với độ tin cậy 95% ta được: stt 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 tong AP MP y 180 180 160 190 200 220 180 170 170 160 230 200 200 230 180 170 180 160 210 190 200 150 250 170 160 4690 x1 300 350 400 400 450 500 350 350 300 450 600 500 450 500 350 400 400 300 550 500 500 300 400 400 300 10300 0.45534 0.021683 x2 x3 6 6.5 10 6.5 5.5 11 7.5 12 6 7.5 14 6.5 186 25.21505 5.53273 5 6 8 5 8 10 10 157 29.87261 6.427888 SUMMARY OUTPUT Regression Statistics Multiple R 0.92083 R Square 0.847927 Adjusted R Square 0.826202 Standard Error 0.054965 Observations 25 ANOVA lny 5.192957 5.192957 5.075174 5.247024 5.298317 5.393628 5.192957 5.135798 5.135798 5.075174 5.438079 5.298317 5.298317 5.438079 5.192957 5.135798 5.192957 5.075174 5.347108 5.247024 5.298317 5.010635 5.521461 5.135798 5.075174 lnx1 5.703782 5.857933 5.991465 5.991465 6.109248 6.214608 5.857933 5.857933 5.703782 6.109248 6.39693 6.214608 6.109248 6.214608 5.857933 5.991465 5.991465 5.703782 6.309918 6.214608 6.214608 5.703782 5.991465 5.991465 5.703782 lnx2 1.791759 1.609438 1.791759 1.871802 2.079442 2.302585 1.871802 1.791759 1.704748 1.791759 2.397895 1.791759 2.014903 2.484907 1.791759 1.791759 2.079442 2.014903 2.197225 2.079442 1.94591 1.609438 2.639057 1.94591 1.871802 lnx3 1.609438 1.609438 1.386294 1.791759 1.791759 1.94591 1.386294 1.791759 1.609438 1.386294 2.197225 2.079442 2.079442 2.197225 1.609438 1.609438 1.791759 1.386294 2.079442 2.079442 2.302585 1.386294 2.302585 1.791759 1.609438 df 21 24 SS 0.353756 0.063445 0.417201 Coefficient s 4.122012 0.047619 0.219422 0.215177 Standard Error 0.405606 0.07904 0.06036 0.061995 Regression Residual Total Intercept lnx1 lnx2 lnx3 MS 0.117919 0.003021 F 39.03056 Significanc eF 9.06E-09 t Stat 10.1626 0.602468 3.635194 3.470849 P-value 1.46E-09 0.553314 0.001549 0.002284 Lower 95% 3.278508 -0.11675 0.093895 0.08625 Upper 95% 4.965516 0.211992 0.344948 0.344103 Lower 95.0% 3.278508 -0.11675 0.093895 0.08625 Sử dụng phương pháp OLS với độ tin cậy 95% ta kết sau: LnYi = 4.1220 + 0.0476 lnX1 + 0.2194 lnX2 + 0.2152 lnX3 Hay Yi = e4.1220 X10.0476 X20.2194 X30.2152 tkd(a1) = 0.6025 tkd(a2) = 3.6352 tkd(a3) = 3.4708 Hệ số tương quan R2 = 0.8479 Hệ số tương quan hiệu chỉnh bình phương: 0.8262 Bước 1: Kiểm định tham số ước lượng mơ hình Kiểm định a1 Giả thuyết H0 : a1 = Giả thuyết H1 : a1 ≠ Ở mức ý nghĩa α = 0.05, số bậc tự (n – k – 1) với n =25 tổng số mẫu, k = số biến độc lập mơ hình, tc với số bậc tự 21 có giá trị: tc(1 - α) = 2.080 Với: tkd(a1) = 0.6025 Upper 95.0% 4.965516 0.211992 0.344948 0.344103 Ta thấy |tkd(a1)| < tc(1 - α) , ta chấp nhận H 0, bác bỏ H1 có nghĩa hệ số hồi quy a1 khơng có ý nghĩa thống kê độ tin cậy 95% Hay mơ hình lượng bón phân hữu khơng ảnh hưởng đến suất lúa Kiểm định a2 Giả thuyết H0 : a2 = Giả thuyết H1 : a2 ≠ Ở mức ý nghĩa α = 0.05, tương tự ta có tc(1 - α) = 2.080 Với : tkd(a2) = 3.6352 Ta thấy |tkd(a2)| > tc(1 - α), ta chấp nhận H1 bác bỏ H0 Tức hệ số hồi quy a2 có ý nghĩa thống kê độ tin cậy 95% Hệ số có ý nghĩa tăng 1kg đạm bón suất tăng lên lượng MP2 = 5.53kg ( MP2 = a2*AP2) ( MP sản phẩm cận biên, AP sản phẩm bình quân) Kiểm định a3 Giả thuyết H0 : a3 = Giả thuyết H1 : a3 ≠ Ở mức ý nghĩa α = 0.05, tương tự ta có: tc(1 - α) = 2.080 Với : tkd(a3) = 3.4708 Ta thấy |tkd(a3)| > tc(1 - α), ta chấp nhận H1 bác bỏ H0 Tức hệ số hồi quy a3 có ý nghĩa thống kê độ tin cậy 95% Hệ số có ý nghĩa tăng cơng lao động suất tăng lên lượng MP3 = 6.43 kg ( MP3 = a3*AP3) Bước 4: Đánh giá độ chặt chẽ mơ hình Phân tích bảng ANOVA Regression = 0.353756 Residual = 0.063445 Total = 0.417201 Hệ số R2 = Regression / Total = 0.353756/0.417201 = 0.847927,mơ hình chặt chẽ Bước 5: Kiểm định mơ hình hồi quy (kiểm định R2) Giả thuyết H0 : R2= Giả thuyết H1 : R2 ≠ Dựa vào bảng kết phân tích ta có Fkd = 39.0306 Ta có Fc(1 – α) = 3.07, với số bậc tự 21(số bậc tự k n – k – 1) Ta thấy Fkd > Fc(1 – α), nên ta bác bỏ H0, chấp nhận H1, R2 ước lượng tin cậy hay có ý nghĩa thống kê độ tin cậy 95% Mơ hình giải thích 84.79% biến động suất lúa biến động yếu tố lượng phân đạm bón cơng lao động, 15.21% biến động suất yếu tố khác, hay tỷ lệ sai số mơ hình 15.21% IV KẾT LUẬN VÀ ĐỀ NGHỊ Qua phân tích mơ hình thấy suất lúa chịu ảnh hưởng lớn rõ rệt hai yếu tố lượng đạm bón cơng lao động Cịn lượng phân hữu khơng có ảnh hưởng Tuy nhiên lượng mẫu thống kê cịn chưa đủ lớn nên chưa thể khẳng định xác ảnh hưởng yếu tố mơ hình Vì ngồi yếu tố suất lúa chịu ảnh hưởng nhiều yếu tố rủi ro thời tiết, thiên tai, bệnh dịch hại, tính chất đất canh tác, chế chăm sóc, áp dụng khoa học kỹ thuật canh tác, khâu thu hoạch… Và tùy thuộc vào mục tiêu người nông dân trồng lúa tối đa hóa lợi nhuận hay tối đa hóa sản lượng, với giá yếu tố đầu vào đầu mà có lựa chọn cho phù hợp, tăng đầu tư đạm hay tăng công lao động Tuy nhiên nên đầu tư tới mức cụ thể mà tối ưu hóa lợi nhuận hay sản lượng mơ hình có dạng đồ thị sản phẩm cận biên có xu hướng giảm đầu tư tăng stt 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 tong AP MP y 180 180 160 190 200 220 180 170 170 160 230 200 200 230 180 170 180 160 210 190 200 150 250 170 160 4690 x1 300 350 400 400 450 500 350 350 300 450 600 500 450 500 350 400 400 300 550 500 500 300 400 400 300 10300 0.45534 0.021683 x2 x3 6 6.5 10 6.5 5.5 11 7.5 12 6 7.5 14 6.5 186 25.21505 5.53273 5 6 8 5 8 10 10 157 29.87261 6.427888 lny 5.192957 5.192957 5.075174 5.247024 5.298317 5.393628 5.192957 5.135798 5.135798 5.075174 5.438079 5.298317 5.298317 5.438079 5.192957 5.135798 5.192957 5.075174 5.347108 5.247024 5.298317 5.010635 5.521461 5.135798 5.075174 lnx1 5.703782 5.857933 5.991465 5.991465 6.109248 6.214608 5.857933 5.857933 5.703782 6.109248 6.39693 6.214608 6.109248 6.214608 5.857933 5.991465 5.991465 5.703782 6.309918 6.214608 6.214608 5.703782 5.991465 5.991465 5.703782 lnx2 1.791759 1.609438 1.791759 1.871802 2.079442 2.302585 1.871802 1.791759 1.704748 1.791759 2.397895 1.791759 2.014903 2.484907 1.791759 1.791759 2.079442 2.014903 2.197225 2.079442 1.94591 1.609438 2.639057 1.94591 1.871802 lnx3 1.609438 1.609438 1.386294 1.791759 1.791759 1.94591 1.386294 1.791759 1.609438 1.386294 2.197225 2.079442 2.079442 2.197225 1.609438 1.609438 1.791759 1.386294 2.079442 2.079442 2.302585 1.386294 2.302585 1.791759 1.609438 SUMMARY OUTPUT Regression Statistics Multiple R 0.92083 R Square 0.847927 Adjusted R Square 0.826202 Standard Error 0.054965 Observations 25 ANOVA df Regression Residual Total Intercept lnx1 lnx2 lnx3 21 24 SS 0.353756 0.063445 0.417201 Coefficient s 4.122012 0.047619 0.219422 0.215177 Standard Error 0.405606 0.07904 0.06036 0.061995 MS 0.117919 0.003021 F 39.03056 Significanc eF 9.06E-09 t Stat 10.1626 0.602468 3.635194 3.470849 P-value 1.46E-09 0.553314 0.001549 0.002284 Lower 95% 3.278508 -0.11675 0.093895 0.08625 Upper 95% 4.965516 0.211992 0.344948 0.344103 Lower 95.0% 3.278508 -0.11675 0.093895 0.08625 Upper 95.0% 4.965516 0.211992 0.344948 0.344103 RESIDUAL OUTPUT Observatio n 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 Predicted lny 5.133086 5.100421 5.09877 5.203579 5.254749 5.341898 5.109974 5.179658 5.113994 5.104378 5.42557 5.258545 5.30249 5.43598 5.140426 5.146785 5.24914 5.134033 5.352051 5.321668 5.340384 5.045065 5.481849 5.21984 5.150649 Residual s 0.059871 0.092536 -0.0236 0.043445 0.043569 0.05173 0.082983 -0.04386 0.021805 -0.0292 0.012509 0.039773 -0.00417 0.002099 0.052531 -0.01099 -0.05618 -0.05886 -0.00494 -0.07464 -0.04207 -0.03443 0.039611 -0.08404 -0.07547 10 lnx2 Line Fit Plot 5.5 4.5 lny lny lny lnx1 Line Fit Plot Predicted lny 5.6 5.8 6.2 6.4 6.6 lny lnx2 lnx3 Line Fit Plot 5.5 4.5 lny Predicted lny lny lnx1 5.5 4.5 lnx3 11 Predicted lny ... độ tin cậy 95% Mơ hình giải thích 84.79% biến động suất lúa biến động yếu tố lượng phân đạm bón cơng lao động, 15.21% biến động suất yếu tố khác, hay tỷ lệ sai số mơ hình 15.21% IV KẾT LUẬN VÀ... KẾT LUẬN VÀ ĐỀ NGHỊ Qua phân tích mơ hình thấy suất lúa chịu ảnh hưởng lớn rõ rệt hai yếu tố lượng đạm bón cơng lao động Cịn lượng phân hữu khơng có ảnh hưởng Tuy nhiên lượng mẫu thống kê chưa... đó: Yi: suất lúa (kg/ sào) X1: lượng phân hữu bón ( kg/sào) X2: lượng đạm bón (kg/ sào) X3: cơng lao động (người/ ngày/ sào) a0: hệ số tự do( = lnA) a1, a2 , a3 hệ số ảnh hưởng Xi đến đến Yi tương