3.1. Dự báo xu hƣớng diễn biến giá xăng dầu trong thời gian tới
3.1.2. Dự báo định lượng về giá dầu thô thế giới trong ngắn hạn
Dựa theo số liệu giá dầu thô thế giới theo ngày từ tháng 7 năm 2014 đến nay, tác giả sẽ sử dụng mơ hình tự hồi quy kết hợp trung bình trượt ARIMA để dự báo giá cả của dầu thơ thế giới trong vịng một tuần tiếp theo.
ARIMA là một mơ hình sử dụng phương pháp Box- Jenkins để phân tích dữ liệu chuỗi thời gian đã được áp dụng trong nhiều nghiên cứu khoa học và có khả năng dự báo tốt trong ngắn hạn. Một chuỗi thời gian trung bình trượt kết hợp tự hồi quy ARIMA(p,d,q) là một chuỗi thời gian dừng sau khi được sai phân d lần, với p biểu thị số các số hạng tự hồi quy (AR), q là số các số hạng trung bình trượt (MA). Phương pháp Box- Jenkins gồm 4 bước:
Bước 1: Nhận dạng
Tìm các giá trị thích hợp của p, d và q thông qua biểu đồ tương quan correlogram và biểu đồ tương quan riêng phần.
Bước 2: Ước lượng
Sau khi xác định được p và q thích hợp, ta ước lượng các hệ số của các số hạng tự hồi quy và trung bình trượt trong mơ hình.
Bước 3: Kiểm tra chẩn đoán
Một kiểm định chẩn đoán đơn giản thường dùng là xem các phần dư ước lượng từ mơ hình có tính ngẫu nhiên thuần túy hay khơng; nếu có, ta có thể chấp nhận sự phù hợp này của mơ hình; nếu khơng ta phải lặp lại từ đầu.
Bước 4: Dự báo
Sau khi tìm được một mơ hình phù hợp ta có thể thực hiện dự báo với phương trình ước lượng.
Sau khi thu thập giá dầu thơ thế giới (tính theo giá bán FOB, đơn vị tính là USD/thùng) từ ngày 01/07/2014 đến 02/05/2015, ta kiểm định xem chuỗi dữ liệu này có đạt tính dừng (giá trị trung bình, phương sai và hiệp phương sai tại các độ trễ khác nhau không đổi) hay không thông qua kiểm định Augmented Dickey-Fuller.
Kết quả kiểm tra tính dừng của chuỗi giá dầu thơ thế giới :
Hình 3.1: Kiểm định tính dừng của chuỗi giá dầu thế giới
Null Hypothesis: SERIES10 has a unit root Exogenous: Constant
Lag Length: 1 (Automatic based on SIC, MAXLAG=14)
t-Statistic Prob.* Augmented Dickey-Fuller test statistic -1.524830 0.5192 Test critical values: 1% level -3.462253
5% level -2.875468 10% level -2.574271
Nguồn: Tính tốn của tác giả
Kết quả cho p-value>0,05 nên ở mức ý nghĩa 5% kết luận chuỗi này có nghiệm đơn vị (unit root), do đó nó khơng có tính dừng.
Để biến đổi chuỗi không dừng về dạng chuỗi dừng, thường ta sử dụng sai phân bậc một hoặc bậc hai của chuỗi. Nếu chuỗi dừng ở sai phân bậc 1, ta gọi chuỗi là chuỗi tích hợp bậc 1 và nếu chuỗi dừng ở sai phân bậc 2 thì chuỗi là tích hợp bậc 2.
Ta tạo chuỗi DEL là sai phân bậc 1 của chuỗi giá dầu thô thế giới với cơng thức:
DELt= Yt+1-Yt Trong đó Y là giá bán dầu thơ thế giới.
Kết quả kiểm định tính dừng của chuỗi DEL trên Eviews:
Hình 3.2: Kiểm định tính dừng của chuỗi sai phân bậc 1 của giá dầu thế giới giá dầu thế giới
Null Hypothesis: DEL has a unit root Exogenous: Constant
Lag Length: 0 (Automatic based on SIC, MAXLAG=14)
t-Statistic Prob.* Augmented Dickey-Fuller test statistic -17.57996 0.0000 Test critical values: 1% level -3.462253
5% level -2.875468 10% level -2.574271 *MacKinnon (1996) one-sided p-values.
Nguồn: Tính tốn của tác giả
Kết quả p-value <0,05 nên ở mức ý nghĩa 5% chuỗi này khơng có nghiệm đơn vị, do đó nó có tính dừng. Sau khi dùng lệnh Correlogram và đốn định dạng tự hồi quy theo hàm ACF và PACF, tác giả nhận định chuỗi DEL có thể là chuỗi tự hồi quy dạng AR(1).Kết quả correlogram đến bậc 25 của chuỗi DEL được thể hiện qua hình sau:
Hình 3.3: Correlogram của chuỗi sai phân bậc một của giá dầu thế giới
Autocorrelation Partial Correlation AC PAC Q-Stat Prob **|. | **|. | 1 -0.207 -0.207 8.9368 0.003 .|* | .|. | 2 0.098 0.058 10.964 0.004 .|. | .|. | 3 -0.036 -0.005 11.232 0.011 .|. | .|. | 4 -0.006 -0.020 11.239 0.024 .|. | .|. | 5 -0.021 -0.025 11.335 0.045 .|* | .|* | 6 0.092 0.090 13.148 0.041 .|. | .|. | 7 -0.026 0.011 13.299 0.065 .|. | .|. | 8 -0.014 -0.034 13.340 0.101 .|. | .|. | 9 0.050 0.048 13.889 0.126 .|. | .|. | 10 0.039 0.068 14.215 0.163 .|. | .|. | 11 -0.008 0.005 14.228 0.221 .|. | .|. | 12 -0.000 -0.018 14.228 0.286 .|* | .|* | 13 0.085 0.094 15.815 0.259 .|. | .|. | 14 -0.061 -0.019 16.636 0.276 .|. | .|. | 15 0.038 -0.002 16.965 0.321 .|. | .|. | 16 0.059 0.073 17.759 0.338 .|. | .|. | 17 0.001 0.034 17.759 0.404 .|. | .|. | 18 0.023 0.018 17.877 0.464 .|. | *|. | 19 -0.048 -0.066 18.402 0.496 .|. | .|. | 20 0.038 0.031 18.743 0.539 .|. | .|. | 21 -0.006 0.020 18.751 0.601 .|. | .|. | 22 0.055 0.028 19.454 0.617 *|. | *|. | 23 -0.086 -0.083 21.183 0.570 .|. | .|. | 24 0.054 0.025 21.867 0.587 .|. | .|. | 25 -0.030 0.002 22.077 0.631
Kết quả hồi quy chuỗi DEL theo dạng AR(1):
Hình 3.4: Kết quả hồi quy chuỗi sai phân bậc một của giá dầu thế giới
Variable Coefficient Std. Error t-Statistic Prob. C -0.242433 0.094990 -2.552195 0.0114 AR(1) -0.206789 0.068646 -3.012411 0.0029
R-squared 0.042790 Mean dependent var
- 0.242049 Adjusted R-squared 0.038074 S.D. dependent var 1.673461 S.E. of regression 1.641294 Akaike info criterion 3.838555 Sum squared resid 546.8509 Schwarz criterion 3.870975 Log likelihood -391.4519 Hannan-Quinn criter. 3.851668 F-statistic 9.074621 Durbin-Watson stat 1.976449 Prob(F-statistic) 0.002921
Inverted AR Roots -.21
Nguồn:Tính tốn của tác giả
Kết quả hồi quy cho thấy hệ số chặn và hệ số hồi quy của số hạng tự hồi quy bậc 1 đều có ý nghĩa thống kê ở mức ý nghĩa 5% (p-value <0,05), đồng thời kiểm định F cũng cho biết mơ hình phù hợp (p-value = 0,002921<0,05). Bước tiếp theo, ta cần kiểm định tính dừng của chuỗi phần dư của hồi quy trên. Sau khi tạo phần dư E, sử dụng kiểm định Augmented Dickey-Fuller ta thu được kết quả kiểm định nghiệm đơn vị của E như sau:
Hình 3.5:Kết quả kiểm định tính dừng của phần dƣ
Null Hypothesis: RESID01 has a unit root Exogenous: Constant
Lag Length: 0 (Automatic based on SIC, MAXLAG=14)
t-Statistic Prob.* Augmented Dickey-Fuller test statistic -14.08249 0.0000 Test critical values: 1% level -3.462253
5% level -2.875468 10% level -2.574271 *MacKinnon (1996) one-sided p-values.
Nguồn: Tính tốn của tác giả
Do p-value = 0,000<0,005 nên ở mức ý nghĩa 5% công nhận chuỗi phần dư E khơng có nghiệm đơn vị. Ngồi ra, kết quả kiểm tra correlogram của E cũng cho
kết quả khơng có tự tương quan hay tự tương quan riêng phần nào có ý nghĩa thống kê. Vậy, chuỗi phần dư E là hồn tồn ngẫu nhiên.
Hình 3.6: Correlogram của chuỗi phần dƣ
Autocorrelation Partial Correlation AC PAC Q-Stat Prob .|. | .|. | 1 0.012 0.012 0.0300 0.862 .|. | .|. | 2 0.055 0.055 0.6572 0.720 .|. | .|. | 3 -0.019 -0.020 0.7317 0.866 .|. | .|. | 4 -0.019 -0.021 0.8053 0.938 .|. | .|. | 5 -0.004 -0.002 0.8096 0.976 .|* | .|* | 6 0.090 0.092 2.5444 0.863 .|. | .|. | 7 -0.012 -0.015 2.5752 0.921 .|. | .|. | 8 -0.010 -0.021 2.5963 0.957 .|. | .|. | 9 0.060 0.066 3.3868 0.947 .|. | .|. | 10 0.051 0.056 3.9639 0.949 .|. | .|. | 11 -0.000 -0.010 3.9639 0.971 .|. | .|. | 12 0.016 0.003 4.0231 0.983 .|* | .|* | 13 0.079 0.089 5.4060 0.965 .|. | .|. | 14 -0.040 -0.038 5.7579 0.972 .|. | .|. | 15 0.042 0.021 6.1469 0.977 .|. | .|* | 16 0.073 0.076 7.3628 0.966 .|. | .|. | 17 0.019 0.022 7.4475 0.977 .|. | .|. | 18 0.015 -0.000 7.4978 0.985 .|. | .|. | 19 -0.039 -0.059 7.8466 0.988 .|. | .|. | 20 0.030 0.048 8.0570 0.991 .|. | .|. | 21 0.014 0.016 8.0997 0.995 .|. | .|. | 22 0.040 0.007 8.4764 0.996 *|. | *|. | 23 -0.070 -0.079 9.6223 0.993 .|. | .|. | 24 0.034 0.042 9.8903 0.995 .|. | .|. | 25 -0.001 0.005 9.8905 0.997
Nguồn: Tính tốn của tác giả
Vậy kết quả hồi quy chuỗi ARIMA thu được sai phân của giá dầu thô thế giới là chuỗi tự hồi quy AR(1) theo phương trình:
DEL = -0,242433282423 + [AR(1)=-0,206789114188] Sau khi biến đổi phương trình thu được:
Vậy, dự báo giá dầu thế giới trong thời gian một tuần tới là:
Bảng 3.1: Dự báo giá dầu thô trong ngắn hạn
Đơn vị: USD/thùng
Ngày 1 2 3 4 5 6 7
Vậy, kết quả ước lượng theo ARIMA cho thấy xu hướng ngắn hạn của giá dầu là giảm nhẹ quanh mốc 55USD/thùng. Như đã nói ở trên, thị trường dầu thế giới đang trải qua nhiều biến động phức tạp do đó các kết quả ước lượng chỉ mang tính tham khảo. Theo EIA, hiện nay thị trường đang ở tình trạng dư cung – là điều kiện để giá dầu duy trì ở mức thấp, tuy nhiên tình hình sản xuất dầu ở Mỹ lại có nhiều diễn biến phức tạp: lượng dầu dự trữ giảm đột ngột vào tuần đầu tháng 5/2015 khiến giá cả tăng trong một thời gian ngắn trước khi suy giảm nhanh chóng khi năng suất lọc dầu giảm từ 93% xuống 91,2%. Do Mỹ hiện là nước khai thác dầu nhiều nhất trên thế giới và cũng là một trong những nước có nhu cầu tiêu thụ lớn nhất, mọi động thái của thị trường dầu ở Mỹ đều có thể tác động ngay tức thì tới diễn biến gía cả trên thị trường thế giới. Ngoài ra, báo cáo ngày 06/05/2015 của EIA cho biết sản lượng khai thác dầu của một số nước như Trung Quốc, Việt Nam và Malaysia đã tăng trong những tháng đầu năm 2015, và được dự báo là sẽ tiếp tục tăng đến khoảng 100.000 thùng dầu/ ngày trong năm nay.
Tóm lại, thị trường dầu thế giới đang trải qua những biến động hết sức phức tạp và ảnh hưởng của những biến động đó lên giá dầu là rất khó xác định chính xác. Hơn nữa, do nước ta điều hành giá xăng dầu theo phương pháp cịn mang nặng tính can thiệp hành chính, các quyết định điều chỉnh giá có một độ trễ nhất định so với biến động giá thị trường nên thực tế có thể những điều chỉnh giá xăng dầu ở Việt Nam có thể đi ngược lại với diễn biến giá xăng dầu thế giới tại cùng một thời điểm.