4.4.1 Mụ hình dự báo
Sau khi kiờ̉m định, phõn tích, và đánh giá kờ́t quả ta có thờ̉ dùng hàm hụ̀i qui vừa tìm được đờ̉ tiờ́n hành dự báo lượng cõ̀u nhọ̃p khõ̉u urờ cho các năm kờ́ tiờ́p khi cho trước giá trị của các biờ́n giải thích. Có thờ̉ tiờ́n hành dự báo điờ̉m và dự báo khoảng cho lượng cõ̀u nhọ̃p khõ̉u urờ trung bình với xác xuṍt tin cọ̃y (1-a)100%.
Giả sử véctơ các giá trị của biờ́n giải thích đã biờ́t là: X0 = (P0; LT0, S0, DV0, 1)
Đờ̉ dự báo điờ̉m lượng cõ̀u nhọ̃p khõ̉u trung bình của năm kờ́ tiờ́p ta chỉ viợ̀c thay giá trị X0 vào hàm hụ̀i qui trờn ta có được giá trị lnURE = f(x0) = ∧y 0; và do đó lượng cõ̀u nhọ̃p khõ̉u của năm đó được dự báo sẽ bằng e^( ∧
y0).
Khoảng dự báo cho giá trị lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
∧
y 0- t0,025(21 –5).se(Y∧ 0/X0) ≤ ln[E(Y/X0)] ≤ y∧0 + t0,025(21 –5).se(Y∧
0/X0)
Trong đó t0,025(21 –5) là giá trị phõn phụ́i t có 16 bọ̃c tự do với xác xuṍt 95%, se(Y∧ 0/X0) là đụ̣ lợ̀ch chuõ̉n của giá trị trung bình của biờ́n phụ thuụ̣c. Đờ̉ tính se(
Y∧ 0/X0), trước hờ́t ta tính phương sai
)X0 Trong đú:
- Cov(β∧ ) là ma trọ̃n hiợ̀p phương sai của vectơ hợ̀ sụ́ hàm hụ̀i qui mõ̃u β∧ ,
- σ2 là phương sai sai sụ́ ngõ̃u nhiờn của mụ hình hụ̀i qui tụ̉ng thờ̉ được ước
lượng bằng phương sai sai sụ́ ngõ̃u nhiờn của mụ hình hụ̀i qui mõ̃u.
- X là ma trọ̃n các sụ́ liợ̀u thụ́ng kờ với n quan sát của các biờ́n giải thích có dạng: X = − − − 1 .... ... ... ... 1 .... 1 .... )1 ( 2 1 2 )1 ( 22 12 1) 1 ( 21 11 n k n n k k x x x x x x x x x .
4.4.2 Dự báo giá trị của các biờ́n giải thích
Khi dự báo lượng cõ̀u nhọ̃p khõ̉u urờ cho các năm 2007, 2008 và 2009 ta có thờ̉ cho trước các giá trị các biờ́n giải thích tuỳ theo phõn tích tình hình thực tờ́.
Chẳng hạn đụ́i với lượng cung urờ trong nước trong năm 2007 cụng suṍt hai nhà máy Phú Mỹ và Đạm Hà Bắc khụng có thay đụ̉i lớn, nhà máy Khí điợ̀n đạm Cà Mau và Nhà máy Đạm Ninh Bình còn trong quá trình xõy dựng; vì vọ̃y tụ̉ng lượng cung urờ trong nước năm 2007 có thờ̉ đạt 900.000 tṍn (bằng 93% tụ̉ng cụng suṍt). Từ quí 2/2008 theo kờ́ hoạch nhà máy Khí - Điợ̀n - Đạm Cà Mau bắt đõ̀u sản xuṍt và cho cụng suṍt ban đõ̀u 200.000 tṍn urờ/năm (tụ̉ng cụng suṍt thiờ́t kờ́ nhà máy này 800000 tṍn urờ/năm); ta có thờ̉ lṍy mức cung urờ trong nước năm 2008 khoảng 1,1 triợ̀u tṍn; nờ́u nhà máy Khí - Điợ̀n - Đạm Cà Mau năm 2009 sản xuṍt được 500.000 tṍn urờ, thì sẽ đưa mức cung trong nước năm 2009 lờn 1,4 triợ̀u tṍn.
Với giá trị tụ̉ng sản lượng lương thực ta có thờ̉ tiờ́p cọ̃n theo cách của Aysen Tanyeri-Abur và Parr Rosson (2002), tiờ́n hành dự báo các giá trị của biờ́n này ở các
năm tiờ́p theo bằng cách tạo ra chuụ̃i sụ́ liợ̀u mới của nó dựa vào mụ hình hụ̀i qui giữa biờ́n giải thích này với các biờ́n trờ̃ của nó mụ̣t năm và hai năm. Cụ thờ̉ là:
LTt = 1,512 + 0,826 LTt-1 + 0,169 LTt-2; (4-38)
Từ đó ta có thờ̉ dự báo tụ̉ng sản lượng lương thực của VN năm 2007, 2008 và 2009 tương ứng là: 40,86; 41,88 và 42,93 triợ̀u tṍn.
Với các giá trị thực của giá urờ trong ba năm tới, ta thṍy giá urờ nhọ̃p khõ̉u tháng 1/2007 lại tăng lờn ở mức 247 USD/tṍn FOB, tăng 12 USD/tṍn so với tháng 12/2006 do giá urờ thờ́ giới tăng 10-20 USD/tṍn; dự báo trong tháng 2 và tháng 3/2007 giá urờ thờ́ giới võ̃n còn nhích lờn vì nguụ̀n cung từ Trung Đụng sẽ ít đi do mụ̣t sụ́ nhà máy ngừng hoạt đụ̣ng. Giá urờ thị trường trong nước vụ Đụng xuõn 2006- 2007 do vọ̃y cũng tăng nhẹ; ở miờ̀n Bắc tăng 100-200đụ̀ng/kg, phụ̉ biờ́n ở mức 4.550- 4.660 đụ̀ng/kg; ở ĐBSCL giá tăng ở mức 4.600-4.700 đụ̀ng/kg, riờng Kiờn Giang ở mức 4.800-4.850 đụ̀ng/kg. Năm 2007 tình hình cung urờ thờ́ giới võ̃n hạn chờ́ bởi nhiờ̀u nhà máy sản xuṍt urờ của Ấn Đụ̣ và Indonexia đóng cửa; Ấn Đụ̣ bãi bỏ trợ cṍp cho nụng nghiợ̀p, Trung Quụ́c tăng thuờ́ xuṍt khõ̉u urờ lờn 30%, giá dõ̀u lửa tăng trở lại trờn 60USD/thùng …. Mặt khác giá urờ thường biờ́n đụ̣ng thṍt thường, vì vọ̃y đờ̉ dự báo giá urờ ta có thờ̉ lṍy 2 mức: mức giá thṍp và mức giá cao. Dự báo giá thực urờ trong năm 2007, 2008 và 2009 với mức thṍp tương ứng là 4,5 triợ̀u đụ̀ng/tṍn; 4,6 triợ̀uđụ̀ng/tṍn; 4,7 triợ̀u đụ̀ng/tṍn; và với mức cao tương ứng là 4,8 triợ̀u đụ̀ng/tṍn; 4,9 triợ̀u đụ̀ng/tṍn; 5,0 triợ̀u đụ̀ng/tṍn.
Bảng 4-24: Dự báo giá thực của urờ, sản lượng lương thực và lượng cung urờ trong nước các năm 2007-2009
2007 2008 2009 Mức giá thṍp Mức giá cao Mức giá thṍp Mức giá cao Mức giá thṍp Mức giá cao P (triợ̀u VND/tṍn) 4,5 4,8 4,6 4,9 4,7 5,0 LT(triợ̀u tṍn) 40,86 40,86 41,88 41,88 42,93 42,93 S (tṍn) 900.000 900.000 1.100.000 1.100.000 1.400.000 1.400.000
4.4.3 Dự báo lượng cõ̀u nhọ̃p khõ̉u urờ
1. Năm 2007:
X0 = (P0; LT0, S0, DV0, 1) = (4,5; 40,86; 900 000; 1; 1) Suy ra:
lnUREt = 8.473- 0.5381n(Pt) + 2.410ln(LTt) - 0.253ln(St) + 0.822 (DVt) =13,96
Vọ̃y lượng nhọ̃p khõ̉u urờ cho năm 2007 với mức giá thṍp dự báo là: 1.105.724 tṍn với:
Ta tính được:
Var(Y∧ 0/X0) = X0TCov(β∧ )X0 = 0,0306
Suy ra: se(Y∧ 0/X0) = 0,175. Tra bảng phõn phụ́i t được t0,025(16) = 2,12. Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,96- 2,12* 0,175≤ ln[E(Y/X0)] ≤ 13,96+ 2,12* 0,175 763.873 ≤ E(Y/X0) ≤ 1.600.562 (tṍn) b. Với phương án giá urờ cao
X0 = (P0; LT0, S0, DV0, 1) = (4,8; 40,86; 900 000; 1; 1) Cov(β∧ )
=
LOG(P) LOG(LT) LOG(S) DV C
LOG(P) 0.060299 -0.07531 -0.01399 0.007632 0.367467 LOG(LT) -0.07531 0.275646 -0.00962 -0.0527 -0.71765 LOG(S) -0.01399 -0.00962 0.0140233 -0.00702 -0.114339
DV 0.007632 -0.0527 -0.007015 0.053144 0.208949 C 0.367467 -0.71765 -0.114339 0.208949 3.304577
Suy ra: lnUREt = 13,92
Vọ̃y lượng nhọ̃p khõ̉u urờ cho năm 2007 với mức giá cao là: 1.068.103 Var(Y∧ 0/X0) = 0,0302 ; se(Y∧ 0/X0) =
0,174
Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,556 ≤ ln[E(Y/X0)] ≤ 14,292 Hay
739.895 ≤ E(Y/X0) ≤ 1.541.901 (tṍn)
2. Năm 2008:
a. Với phương án giá urờ thṍp
X0 = (P0; LT0, S0, DV0, 1) = (4,6; 41,88; 1.100.000; 1; 1) Suy ra :
lnUREt = 8,473- 0.5381n(Pt) + 2,410ln(LTt) – 0,253ln(St) + 0,822 (DVt) = 13,96
Lượng nhọ̃p khõ̉u urờ cho năm 2008 với mức giá thṍp dự báo là: 1.102.234 tṍn . Tính Var(Y∧ 0/X0) = 0,0369; Suy ra se(Y∧ 0/X0) = 0,192
Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,548 ≤ ln[E(Y/X0)] ≤ 14,362 734.486 ≤ E(Y/X0) ≤ 1.654.107 (tṍn) b. Với phương án giá urờ cao
X0 = (P0; LT0, S0, DV0, 1) = (4,9; 41,88; 1.100.000; 1; 1)
Suy ra: lnUREt = 13,92. Vọ̃y lượng nhọ̃p khõ̉u urờ cho năm 2008 với mức giá cao dự báo là: 1.065.508 tṍn
Var(Y∧ 0/X0) = 0,0360; se(Y∧ 0/X0) = 0,190
Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,519 ≤ ln[E(Y/X0)] ≤ 14,324 713.429 ≤ E(Y/X0) ≤ 1.591.340 (tṍn)
3. Năm 2009:
a. Với phương án giá urờ thṍp
(P0; LT0, S0, DV0, 1) = (4,7; 42,93; 1.400.000; 1; 1)
Suy ra: lnUREt = 13,492. Vọ̃y lượng nhọ̃p khõ̉u urờ cho năm 2009 với mức giá thṍp dự báo là: 1.088.142 tṍn.
Tính được var(Y∧ 0/X0) = 0,0456 ; Suy ra se(Y∧ 0/X0) = 0,214
Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,4905 ≤ ln[E(Y/X0)] ≤ 14,395 692.821 ≤ E(Y/X0 ≤ 1.709.031 (tṍn) b. Với phương án giá urờ cao
X0 = (P0; LT0, S0, DV0, 1) = (5,0; 42,93; 1.400.000; 1; 1)
Suy ra: lnUREt = 13,909; Vọ̃y lượng nhọ̃p khõ̉u urờ cho năm 2009 với mức giá cao là: 1.052.622 tṍn
Var(Y∧ 0/X0) = 0,0443 ; se(Y∧ 0/X0) = 0,21
Khoảng dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình (ln) với xác xuṍt tin cọ̃y 95% là khoảng:
13,463 ≤ ln[E(Y/X0)] ≤ 14,3 674.702 ≤ E(Y/X0 ≤ 1.642.225 (tṍn)
Bảng 4-25: Dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình cho các năm 2007- 2009 với mức giá thṍp
Dự báo 2007 2008 2009
Lượng cõ̀u nhọ̃p
Dự báo 2007 2008 2009
(tṍn) (tṍn) (tṍn)
Đụ̣ lợ̀ch chuõ̉n (ln) 0,175 0,192 0,214
Khoảng tin cọ̃y 95% (763.873; 1.600.562) (tṍn)
(734.486; 1.654.107) (tṍn)
(692.821; 1.709.031) (tṍn)
Như vọ̃y, nờ́u tình hình thị trường urờ thờ́ giới khụng có biờ́n đụ̣ng lớn dự báo lượng nhọ̃p khõ̉u urờ của VN trong các năm 2007, 2008 và 2009 tương ứng khoảng 1,1; 1,1 và 1,09 triợ̀u tṍn. Trường hợp, thị trường urờ thờ́ giới có những biờ́n đụ̣ng làm giá urờ tăng mạnh thì dự báo lượng nhọ̃p khõ̉u urờ của VN trong các năm 2007, 2008 và 2009 tương ứng khoảng 1,06; 1,06 và 1,05 triợ̀u tṍn.
Bảng 4-26: Dự báo lượng cõ̀u nhọ̃p khõ̉u urờ trung bình cho các năm 2007- 2009 với mức giá cao
Dự báo 2007 2008 2009
Lượng cõ̀u nhọ̃p khõ̉u TB 1.068.103 (tṍn) 1.065.508 (tṍn) 1.052.622 (tṍn) Đụ̣ lợ̀ch chuõ̉n (ln) 0,174 0,190 0,210
Khoảng tin cọ̃y 95% (739.895; 1.541.901)
(tṍn) (713.429; 1.591.340) (tṍn) (674.702; 1.642.225) (tṍn) Mụ̣t điờ̀u lưu ý là khi dự báo càng xa thì đụ̣ lợ̀ch chuõ̉n càng lớn; đờ̉ khắc phục điờ̀u này chúng ta có thờ̉ dự báo cho từng năm mụ̣t. Tức là, muụ́n dự báo lượng cõ̀u nhọ̃p khõ̉u urờ cho năm kờ́ tiờ́p chính xác hơn, sau mụ̣t năm hoạt đụ̣ng ta phõn tích và nhọ̃p lại sụ́ liợ̀u với cỡ mõ̃u tăng thờm 1, sau đó tiờ́n hành hụ̀i qui đờ̉ xác định lại hàm cõ̀u nhọ̃p khõ̉u urờ và dự báo tương tự.