Thiết kế thí nghiệm và xử lý kết quả bằng phần mềm thống kê IRRISTAT part 10

Từ kết cho thấy hồi quy lập ñược có ñộ tin cậy cao với phương trình tuyến tính bậc ñược xác ñịnh sau: Y = 4.555X + 4241.12 5.2 Hồi quy bôi tuyến tính Hồi quy bội tuyến tính ñược biểu diễn dạng phương trình tổng quát sau: Y = a1X1 + a2X2 + … + anXn + B Trong ñó: X1, X2, … Xn biến ñộc lập a1, a2, … hệ số biến ñộc lập B số Y biến phụ thuộc Ví dụ ta có tập hợp số liệu theo dõi suất phân tích hàm lượng N, P, K hạt từ 45 mảnh ruộng sau: PLOT Y14 G 1678 4265 2431 2431 4461 3110 4469 N% 0.9849 1.1714 1.0756 1.0435 1.2101 1.2084 1.1643 P% 0.0901 0.0926 0.0886 0.0807 0.0852 0.0845 0.0972 K% 0.3987 0.3814 0.4134 0.4027 0.3851 0.3746 0.3794 Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 104 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 4194 3379 7132 4359 3646 6917 6692 7028 7196 4220 7250 4915 4185 7463 7367 7860 7624 2830 3705 3280 2906 4041 3295 4147 4102 3509 5087 4353 3915 5122 4660 5150 5730 3744 5363 4628 4212 5063 1.2369 1.0247 1.4624 0.9954 1.0404 1.5309 1.5001 1.6112 1.6009 0.9392 1.6738 0.9525 0.8468 1.602 1.6678 1.5521 1.5707 0.8732 1.0743 1.0124 0.9288 1.0279 1.1412 1.1594 1.0733 0.8955 1.2949 1.0541 0.9689 1.4547 1.2749 1.3373 1.2832 0.9977 1.7496 1.0091 0.8938 1.6833 0.085 0.2901 0.2467 0.2872 0.2487 0.2868 0.2776 0.2852 0.3104 0.2558 0.1985 0.2712 0.2474 0.3213 0.1793 0.3058 0.2969 0.2582 0.2396 0.2665 0.2705 0.2492 0.243 0.2867 0.2624 0.2395 0.2548 0.2604 0.2469 0.3239 0.2454 0.3057 0.3059 0.2103 0.1617 0.2208 0.2008 0.2609 0.4373 0.3108 0.2545 0.3095 0.3137 0.3071 0.2506 0.3091 0.3168 0.3039 0.2474 0.3116 0.3051 0.3144 0.2295 0.3296 0.3078 0.3204 0.239 0.3016 0.3078 0.2925 0.2475 0.2924 0.2736 0.3815 0.2978 0.3597 0.3806 0.3195 0.2556 0.3158 0.3104 0.2943 0.2353 0.2993 0.2919 0.2983 Vào số liệu bảng IRRSTAT dạng sau: Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 105 Tạo thêm biến ñể ghi giá trị dư thừa giá trị phù hợp sau tính hồi qui ñể có ảnh sau: Tạo xong, save file lại ñể chạy hồi qui Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 106 Chọn Analysis Regression, viết lệnh cho Comand File xong Open ñể mở File số liệu ảnh dưới: Từ ảnh trên, ñưa biến GY14,N%, P%, K% vào hộp Analysis Variates Xong chọn Regression Model, chọn New chuyển biến từ Analysis Variates vào hộp Response Variate Independent Variates, chọn Add to the Model hộp Selection ñể có ảnh sau: Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 107 Chọn Option, Correclation, Compute and Test Residuals ảnh xong kích chuột vào OK ñể chạy mô hình Kết chạy mô hình: CORRELATION MATRIX AND STANDARD ERRORS FILE HQUYBOI 14/ 9/ 11:27 :PAGE SECTION 1 CORR AND RESID SES Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 108 Y14G N% P% K% 0.1584525E+04 0.8156641E+00 0.2676436E+00 0.4406525E+00 0.1723501E+00 0.7490746E-01 -0.4226677E+00 -0.3386678E+00 -0.5344221E+00 0.5054769E-01 Y14G N% P% K% REGRESSION OF Y14G (2 ) ON VARIATES FILE HQUYBOI 14/ 9/ 11:27 :PAGE TERM NO CONSTANT N% P% K% SOURCE REGRESSION RESIDUAL TOTAL COEFFICIENT -2284.148 4515.239 6554.451 38.24514 SS 0.8374954E+08 0.2672206E+08 0.1104716E+09 STD ERROR F-VALUE 1484.60 483.325 1922.42 2982.48 DF 41 44 2.367 87.274 11.625 0.000 MS 0.2791651E+08 651757.7 PRBF TOLERANCE 0.128 0.000 0.002 0.987 0.6571E-02 0.8852 0.7143 0.6517 F PRBF 42.833 0.0000 R-SQUARED - 75.8%, (ADJUSTED FOR D.F.- 74.0%) RECORD RECORD 42 HAS INFLUENTIAL X-VAL OBS= 45 HAS RESIDUAL= 2.6 SES OBS= 5363 5063 RES= RES= BOX PLOT OF STUDENTIZED RESIDUALS FROM LPLT= -2.570 NO.UPLT I -0 0.516 Kết cho ta phương trình tuyến tính với ba biến ñộc lập sau: Y = -2284,148 + 4515,239 N% + 6554,451 P% + 38,245 K% Với ñộ tin cậy phương trình hồi quy ñạt mức cao với xác suất nhận Ho nhỏ (0,000) Chứng tỏ phương trình có ý nghĩa cho việc suy ñoán từ yếu tố N, P, K suất Y Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT……… 109 ... 0.2556 0.3158 0. 3104 0.2943 0.2353 0.2993 0.2919 0.2983 Vào số liệu bảng IRRSTAT dạng sau: Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT …… 105 Tạo thêm biến... Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT …… 106 Chọn Analysis Regression, viết lệnh cho Comand File xong Open ñể mở File số liệu ảnh dưới: Từ ảnh trên, ñưa biến GY14,N%, P%, K% vào hộp... Variates vào hộp Response Variate Independent Variates, chọn Add to the Model hộp Selection ñể có ảnh sau: Trường ðại học Nông nghiệp Hà Nội – Giáo trình Thiết kế thí nghiệm xử lý kết IRRISTAT …… 107