BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH BÙI VĂN HOÀNG PHẢN ỨNG CỦA TƯƠNG QUAN CÁC CHỈ BÁO THỊ TRƯỜNG TÀI CHÍNH TRƯỚC CÁC CÚ SỐC DAO ĐỘNG: TRƯỜNG HỢP CỦA VIỆT NAM LUẬN VĂN THẠC SỸ KINH TẾ TP Hồ Chí Minh – Năm 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH BÙI VĂN HOÀNG PHẢN ỨNG CỦA TƯƠNG QUAN CÁC CHỈ BÁO THỊ TRƯỜNG TÀI CHÍNH TRƯỚC CÁC CÚ SỐC DAO ĐỘNG: TRƯỜNG HỢP CỦA VIỆT NAM LUẬN VĂN THẠC SỸ KINH TẾ TP Hồ Chí Minh – Năm 2016 BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC KINH TẾ TP HỒ CHÍ MINH BÙI VĂN HOÀNG PHẢN ỨNG CỦA TƯƠNG QUAN CÁC CHỈ BÁO THỊ TRƯỜNG TÀI CHÍNH TRƯỚC CÁC CÚ SỐC DAO ĐỘNG: TRƯỜNG HỢP CỦA VIỆT NAM Chuyên ngành: Tài Chính – Ngân Hàng Mã số: 60340201 LUẬN VĂN THẠC SỸ KINH TẾ NGƯỜI HƯỚNG DẪN KHOA HỌC: PGS TS NGUYỄN KHẮC QUỐC BẢO TP Hồ Chí Minh – Năm 2016 LỜI CAM ĐOAN Tôi xin cam đoan công trình nghiên cứu tôi, hướng dẫn người hướng dẫn nghiên cứu khoa học PGS.TS Nguyễn Khắc Quốc Bảo Các nội dung nghiên cứu kết đề tài trung thực chưa công bố công trình Những số liệu bảng biểu phục vụ cho việc phân tích, nhận xét, đánh giá tác giả thu thập từ nguồn khác có ghi phần tài liệu tham khảo Ngoài ra, luận văn sử dụng số nhận xét, đánh số liệu tác giả khác có thích nguồn gốc sau trích dẫn để dễ tra cứu, kiểm chứng Thành phố Hồ Chí Minh, ngày 12 tháng 09 năm 2016 Tác giả Bùi Văn Hoàng MỤC LỤC DANH MỤC BẢNG DANH MỤC HÌNH TÓM TẮT CHƯƠNG 1: MỞ ĐẦU 1.1 Giới thiệu 1.2 Mục tiêu câu hỏi nghiên cứu 1.3 Đối tượng phạm vi nghiên cứu 1.4 Phương pháp nghiên cứu 1.5 Đóng góp đề tài 1.6 Bố cục luận văn CHƯƠNG 2: TỔNG QUAN LÝ THUYẾT VÀ NGHIÊN CỨU TRƯỚC ĐÂY 2.1 Quan hệ thị trường cổ phiếu thị trường trái phiếu .7 2.2 Quan hệ thị trường ngoại hối thị trường trái phiếu 10 2.3 Quan hệ thị trường ngoại hối thị trường cổ phiếu .16 2.4 Quan hệ thị trường cổ phiếu, trái phiếu thị trường ngoại hối 21 CHƯƠNG 3: DỮ LIỆU VÀ PHƯƠNG PHÁP NGHIÊN CỨU 24 3.1 Phương pháp nghiên cứu 24 3.2 Lý thuyết đặc điểm kinh tế lượng sử dụng nghiên cứu 26 3.2.1 Thống kê mô tả kiểm định 26 3.2.1.1 Thống kê mô tả 26 3.2.1.2 Kiểm định tính dừng chuỗi liệu .29 3.2.2 Mô hình tự hồi quy vector - VAR 31 3.2.2.1 Mô hình ước lượng 31 3.2.2.2 Lựa chọn độ trễ tối đa 32 3.2.3 Mô hình dao động có điều kiện FIEGARCH 33 3.2.4 Mô hình tương quan có điều kiện động cDCC 38 3.2.5 Thuật toán xác định điểm chuyển đổi dao động PELT 40 3.3 Dữ liệu 44 CHƯƠNG 4: KẾT QUẢ NGUYÊN CỨU VÀ THẢO LUẬN 45 4.1 Thống kê mô tả .45 4.2 Kết ước lượng tiến trình VAR-FIEGARCH-cDCC 48 4.2.1 Mô hình tự hồi quy vector - VAR 48 4.2.2 Kết ước lượng tiến trình FIEGARCH-cDCC 50 4.3 Phân tích điểm chuyển đổi dao động 58 4.4 Ảnh hưởng cú sốc dao động tới tương quan động thị trường 64 CHƯƠNG 5: KẾT LUẬN 69 5.1 Kết nghiên cứu 69 5.2 Những hạn chế luận văn 71 5.3 Những gợi ý hướng nghiên cứu 72 TÀI LIỆU THAM KHẢO PHỤ LỤC DANH MỤC BẢNG Bảng 4.1 Kết thống kê mô tả liệu .45 Bảng 4.2 Kết xác định độ trễ tối ưu mô hình VAR 49 Bảng 4.3 Kiểm định tính ổn định mô hình 49 Bảng 4.4 Kết ước lượng mô hình VAR(1) 50 Bảng 4.5 Kết kiểm định chuẩn đoán phần dư VAR 51 Bảng 4.6 Kết ước lượng mô hình FIEGARCH(1,d,1) 52 Bảng 4.7 Kết ước lượng mô hình cDCC 53 Bảng 4.8 Thống kê mô tả chuỗi tương quan điều kiện động .58 Bảng 4.9 Mô tả thời điểm chuyển đổi dao động chuỗi phần dư VAR(1) 64 Bảng 4.10 Kết khảo sát tác động cú sốc dao động tới tương quan thị trường 68 DANH MỤC HÌNH Hình 4.1 Giá trị thực tế giá trị chênh lệch tỷ giá USD/VND, lãi suất số VN-Index giai đoạn từ tháng 4/2004 đến tháng 12/2015 46 Hình 4.2 Mật độ phân phối tỷ giá USD/VND, lãi suất số VN-Index giai đoạn từ tháng 4/2004 đến tháng 12/2015 48 Hình 4.3 Tương quan có điều kiện động cặp thị trường 54 Hình 4.4 Các điểm chuyển đổi dao động chuỗi phần dư VAR 61 Hình 4.5 Phản ứng tương quan có điều kiện động trước cú sốc thị trường 65 TÓM TẮT Bài viết phân tích biến động mối quan hệ thị trường cổ phiếu, trái phiếu thị trường ngoại hối Việt Nam giai đoạn từ tháng 04/2004 đến tháng 12/2015 Cụ thể hơn, viết trả lời cho câu hỏi có hay biến động đột ngột tương quan thị trường nhằm phản ứng lại với cú sốc thị trường? Trên sở đó, tác giả tiến hành kiểm tra để phát xem liệu biến động tức thời hay kéo dài? Bài viết sử dụng tiến trình ước lượng mô hình 𝑉𝐴𝑅(𝑝) − 𝐹𝐼𝐸𝐺𝐴𝑅𝐶𝐻 (1, 𝑑, 1) − 𝑐𝐷𝐶𝐶 thuật toán xác định điểm chuyển đổi dao động PELT kết hợp với mô hình hồi quy với biến giả Kết nghiên cứu khẳng định có tồn phụ thuộc lẫn thị trường cổ phiếu, trái phiếu ngoại hối Việt Nam Ngoài ra, mối tương quan không ổn định mà biến động theo thời gian Đặc biệt, cú sốc dao động tạo biến động đột ngột tương quan thị trường số thời điểm, điều chỉnh mang tính chất dài hạn Bài viết khuyến nghị nhà đầu tư nhà hoạch định sách nên quan tâm đến hiệu ứng lây lan dài hạn có hoạt động đầu tư xây dựng sách kinh tế trường hợp Việt Nam CHƯƠNG 1: MỞ ĐẦU 1.1 Giới thiệu Đối với nhiều nhà đầu tư hoạch định sách, việc thấu hiểu phụ thuộc giữ thị trường tài vô quan trọng Các thị trường bao gồm thị trường ngoại hối, thị trường cổ phiếu thị trường trái phiếu Đối với nhà đầu tư, cấu trúc tương quan thị trường sử dụng công cụ chủ chốt việc xây dựng điều chỉnh chiến lược đầu tư danh mục Còn người làm sách, việc phân tích kênh truyền dẫn thị trường lại quan trọng, nhằm đưa sách đắn hiệu quả, việc dự báo cách đầy đủ tác động mà sách họ gây cho kinh tế vĩ mô chúng thực thi Ngoài ra, việc phân tích trở nên quan trọng cấp thiết mà thực tế cho thấy rằng, mối tương quan không ổn định mà có xu hướng biến động theo thời gian Các nhà đầu tư có quan tâm lớn đến việc phát thay đổi lãi suất, từ giúp họ dự báo giá cổ phiếu Các nhà đầu tư quản lý vị danh mục cách hợp lý họ sử dụng thông tin thay đổi lãi suất báo đáng tin cậy cho chiều hướng giá thị trường Trong đó, nhà làm sách lại để ý tới tình thị trường cổ phiếu phản ứng lại với thay đổi lãi suất báo cho hoạt động thị trường cổ phiếu tương lai Các nhà phân tích nhà quản lý quỹ nghĩ nguyên tắc dựa chiến lược đầu tư để thu lợi suất vượt trội Trong hai thập kỷ qua, kinh tế toàn cầu trải qua giai đoạn khủng hoảng, đặc biệt kinh tế Có thể điểm qua số khủng hoảng gần đổ vỡ thị trường chứng khoán năm 1987 H0: rSP is I(1) ADF Statistics: -6.43235 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -2.56572 5% -1.94093 10% -1.61663 OLS Results y_1 dy_1 dy_2 RSS OBS Coefficient -0.775678 0.087288 0.045903 1.367692 137.000000 t-value -6.4324 0.83482 0.52935 Information Criteria (to be minimized) Akaike -1.724940 Shibata Schwarz -1.660999 Hannan-Quinn H0: rEX is I(1) -1.725872 -1.698956 ADF Statistics: -5.9691 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -2.56572 5% -1.94093 10% -1.61663 OLS Results y_1 dy_1 dy_2 RSS OBS Coefficient -0.877227 -0.103521 -0.107431 0.015328 137.000000 t-value -5.9691 -0.86044 -1.2508 Information Criteria (to be minimized) Akaike -6.216168 Shibata Schwarz -6.152227 Hannan-Quinn - -6.217100 -6.190184 TESTS : ======= ADF Test with lags Intercept and no time trend H0: rIR is I(1) ADF Statistics: -4.84467 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.4323 5% -2.86229 10% -2.56721 OLS Results y_1 Coefficient -0.497435 t-value -4.8447 dy_1 dy_2 Constant RSS OBS -0.187739 0.027069 -0.000094 0.006766 137.000000 -1.8084 0.31229 -0.15365 Information Criteria (to be minimized) Akaike -7.019090 Shibata Schwarz -6.933835 Hannan-Quinn H0: rSP is I(1) -7.020731 -6.984444 ADF Statistics: -6.47411 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.4323 5% -2.86229 10% -2.56721 OLS Results y_1 dy_1 dy_2 Constant RSS OBS Coefficient -0.787752 0.094885 0.050286 0.007009 1.361063 137.000000 t-value -6.4741 0.90262 0.57800 0.80482 Information Criteria (to be minimized) Akaike -1.715009 Shibata Schwarz -1.629754 Hannan-Quinn H0: rEX is I(1) -1.716650 -1.680363 ADF Statistics: -6.65976 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.4323 5% -2.86229 10% -2.56721 OLS Results y_1 dy_1 dy_2 Constant RSS OBS Coefficient -1.044426 0.007554 -0.051306 0.002598 0.014551 137.000000 t-value -6.6598 0.060513 -0.59248 2.6646 Information Criteria (to be minimized) Akaike -6.253387 Shibata Schwarz -6.168132 Hannan-Quinn - -6.255028 -6.218741 TESTS : ======= ADF Test with lags Intercept and time trend H0: rIR is I(1) ADF Statistics: -4.88688 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 dy_1 dy_2 Constant Trend RSS OBS Coefficient -0.508255 -0.180603 0.030854 -0.000096 -0.000011 0.006741 137.000000 t-value -4.8869 -1.7282 0.35462 -0.15669 -0.70663 Information Criteria (to be minimized) Akaike -7.008019 Shibata Schwarz -6.901450 Hannan-Quinn H0: rSP is I(1) -7.010560 -6.964712 ADF Statistics: -6.46594 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 dy_1 dy_2 Constant Trend RSS OBS Coefficient -0.790542 0.096441 0.051475 0.007034 -0.000088 1.359416 137.000000 t-value -6.4659 0.91390 0.58945 0.80506 -0.39999 Information Criteria (to be minimized) Akaike -1.701374 Shibata Schwarz -1.594805 Hannan-Quinn H0: rEX is I(1) -1.703915 -1.658067 ADF Statistics: -6.63905 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 dy_1 dy_2 Constant Trend RSS Coefficient -1.046689 0.009179 -0.050425 0.002604 0.000006 0.014544 t-value -6.6390 0.073170 -0.57973 2.6602 0.24416 OBS 137.000000 Information Criteria (to be minimized) Akaike -6.238992 Shibata Schwarz -6.132423 Hannan-Quinn - -6.241533 -6.195685 Phụ lục Kết chọn độ trễ tối ưu cho mô hình VAR(p) (Kết thu từ phần mềm EView ver 8.1) VAR Lag Order Selection Criteria Endogenous variables: REX RIR RSP Exogenous variables: C Date: 07/12/16 Time: 15:56 Sample: 2004M04 2015M12 Included observations: 130 Lag LogL LR FPE AIC SC HQ 10 964.7725 983.8987 991.7883 999.2909 1005.758 1016.715 1026.491 1029.234 1030.093 1035.070 1040.254 NA 37.07541 14.92952 13.85089 11.64115 19.21657* 16.69371 4.557759 1.387474 7.810594 7.895510 7.53e-11 6.44e-11* 6.55e-11 6.71e-11 6.98e-11 6.79e-11 6.73e-11 7.43e-11 8.46e-11 9.06e-11 9.68e-11 -14.79650 -14.95229* -14.93520 -14.91217 -14.87320 -14.90331 -14.91524 -14.81898 -14.69373 -14.63185 -14.57314 -14.73033* -14.68759 -14.47199 -14.25043 -14.01294 -13.84453 -13.65794 -13.36315 -13.03939 -12.77898 -12.52175 -14.76961 -14.84473* -14.74698 -14.64328 -14.52365 -14.47309 -14.40436 -14.22743 -14.02152 -13.87896 -13.73959 * indicates lag order selected by the criterion LR: sequential modified LR test statistic (each test at 5% level) FPE: Final prediction error AIC: Akaike information criterion SC: Schwarz information criterion HQ: Hannan-Quinn information criterion Phụ lục Kết tính ổn định mô hình VAR (Kết thu từ phần mềm Eviews ver 8.1) Roots of Characteristic Polynomial Endogenous variables: REX RIR RSP Exogenous variables: C Lag specification: 1 Date: 07/12/16 Time: 17:07 Root Modulus 0.380573 0.282888 -0.019347 0.380573 0.282888 0.019347 No root lies outside the unit circle VAR satisfies the stability condition Phụ lục Kết ước lượng mô hình VAR(1) (Kết thu từ package PCGive phần mềm OxMetrics ver 6.3) Ox Professional version 6.30 (Windows/U) (C) J.A Doornik, 1994-2011 PcGive 13.30 session started at 12:58:45 on 12-07-2016 SYS( 1) Estimating the system by OLS The dataset is: E:\HOANG DATA\Master\THESIS2\VIET BAO\Data\new12.in7 The estimation sample is: 2004(6) - 2015(12) URF equation for: rIR Coefficient rIR_1 0.382282 rSP_1 0.00233305 rEX_1 0.0537608 Constant U -0.000266689 sigma = 0.00721496 t-prob 0.0000 0.6962 0.3745 0.6747 Std.Error 1.107 0.08310 0.8409 0.008835 t-value 0.483 3.55 0.434 0.609 t-prob 0.6299 0.0005 0.6648 0.5433 t-value -0.846 -0.210 -0.385 2.77 t-prob 0.3989 0.8340 0.7007 0.0063 RSS = 1.364869782 URF equation for: rEX Coefficient rIR_1 -0.0967946 rSP_1 -0.00180254 rEX_1 -0.0334718 Constant U 0.00253163 sigma = 0.0103893 t-value 4.81 0.391 0.891 -0.421 RSS = 0.007027503754 URF equation for: rSP Coefficient rIR_1 0.534737 rSP_1 0.295304 rEX_1 0.365115 Constant U 0.00538327 sigma = 0.100549 Std.Error 0.07944 0.005963 0.06034 0.0006340 Std.Error 0.1144 0.008586 0.08689 0.0009129 RSS = 0.01457165693 log-likelihood 1056.39529 |Omega| 5.02724707e-011 R^2(LR) 0.241641 no of observations 139 -T/2log|Omega| 1648.09266 log|Y'Y/T| -23.4369651 R^2(LM) 0.0853012 no of parameters 12 F-test on regressors except unrestricted: F(9,323) = 4.33086 [0.0000] ** F-tests on retained regressors, F(3,133) = rIR_1 8.30200 [0.000]** rSP_1 4.19291 [0.007]** rEX_1 0.381592 [0.766] Constant U 3.05893 [0.031]* correlation of URF residuals (standard deviations on diagonal) rIR rSP rEX rIR 0.0072150 0.026478 0.12274 rSP 0.026478 0.10055 -0.13185 rEX 0.12274 -0.13185 0.010389 correlation between actual and fitted rIR rSP rEX 0.39657 0.29791 0.084537 Single-equation diagnostics using reduced-form residuals: rIR rIR rIR rIR rIR rIR rSP rSP rSP rSP rSP rSP rEX rEX rEX rEX rEX rEX Vector Vector Vector Vector Vector Vector : : : : : : : : : : : : : : : : : : Portmanteau(12): AR 1-7 test: ARCH 1-7 test: Normality test: Hetero test: Hetero-X test: Portmanteau(12): AR 1-7 test: ARCH 1-7 test: Normality test: Hetero test: Hetero-X test: Portmanteau(12): AR 1-7 test: ARCH 1-7 test: Normality test: Hetero test: Hetero-X test: Portmanteau(12): AR 1-7 test: Normality test: Hetero test: Hetero-X test: RESET23 test: Chi^2(11) F(7,128) F(7,125) Chi^2(2) F(6,132) F(9,129) Chi^2(11) F(7,128) F(7,125) Chi^2(2) F(6,132) F(9,129) Chi^2(11) F(7,128) F(7,125) Chi^2(2) F(6,132) F(9,129) Chi^2(99) F(63,335) Chi^2(6) F(36,560) F(54,636) F(18,359) = = = = = = = 26.484 [0.0055]** = 3.9704 [0.0006]** = 1.2872 [0.2619] = 210.14 [0.0000]** = 5.0034 [0.0001]** = 4.6081 [0.0000]** = 17.677 [0.0894] = 0.71270 [0.6613] = 3.1763 [0.0040]** = 16.676 [0.0002]** = 0.50654 [0.8025] = 0.37826 [0.9437] = 11.634 [0.3918] = 1.8834 [0.0774] = 0.049266 [0.9998] = 2037.3 [0.0000]** = 0.087163 [0.9975] = 0.091963 [0.9997] 132.81 1.3306 2078.5 1.4992 1.3907 1.7016 [0.0132]* [0.0594] [0.0000]** [0.0332]* [0.0377]* [0.0371]* Phụ lục Kết kiểm định chuẩn đoán phần dư ước lượng mô hình VAR(1) (Kết thu từ package PCGive phần mềm OxMetrics ver 6.3) Single-equation diagnostics using reduced-form residuals: rIR : Portmanteau(12): Chi^2(11) = 26.484 [0.0055]** rIR : AR 1-7 test: F(7,128) = 3.9704 [0.0006]** rIR : ARCH 1-7 test: F(7,125) = 1.2872 [0.2619] rIR : Normality test: Chi^2(2) = 210.14 [0.0000]** rIR : Hetero test: F(6,132) = 5.0034 [0.0001]** rIR : Hetero-X test: F(9,129) = 4.6081 [0.0000]** rSP : Portmanteau(12): Chi^2(11) = 17.677 [0.0894] rSP : AR 1-7 test: F(7,128) = 0.71270 [0.6613] rSP : ARCH 1-7 test: F(7,125) = 3.1763 [0.0040]** rSP : Normality test: Chi^2(2) = 16.676 [0.0002]** rSP : Hetero test: F(6,132) = 0.50654 [0.8025] rSP : Hetero-X test: F(9,129) = 0.37826 [0.9437] rEX : Portmanteau(12): Chi^2(11) = 11.634 [0.3918] rEX : AR 1-7 test: F(7,128) = 1.8834 [0.0774] rEX : ARCH 1-7 test: F(7,125) = 0.049266 [0.9998] rEX : Normality test: Chi^2(2) = 2037.3 [0.0000]** rEX : Hetero test: F(6,132) = 0.087163 [0.9975] rEX : Hetero-X test: F(9,129) = 0.091963 [0.9997] Vector Vector Vector Vector Vector Vector Portmanteau(12): AR 1-7 test: Normality test: Hetero test: Hetero-X test: RESET23 test: Chi^2(99) F(63,335) Chi^2(6) F(36,560) F(54,636) F(18,359) = = = = = = 132.81 1.3306 2078.5 1.4992 1.3907 1.7016 [0.0132]* [0.0594] [0.0000]** [0.0332]* [0.0377]* [0.0371]* Phụ lục Kết hồi quy mô hình FIEGARCH(1,d,1)-cDCC (Kết thu từ package G@ARCH phần mềm OxMetrics ver 6.3) Copyright for this package: S Laurent, 2007-2010 MG@RCH package version 1.3, object created on 4-04-2016 Copyright for this package: S Laurent, 2000-2010 G@RCH package version 6.1, object created on 4-04-2016 Copyright for this package: S Laurent, 2007-2010 Starting estimation process ******************** ** FIRST STEP ** ******************* -Estimating the univariate GARCH model for VrIR - ******************** ** SPECIFICATIONS ** ******************** Dependent variable : VrIR Mean Equation : ARMA (0, 0) model No regressor in the conditional mean Variance Equation : FIEGARCH (1, d, 1) model (Truncation order : 139) No regressor in the conditional variance Normal distribution Weak convergence (no improvement in line search) using numerical derivatives Log-likelihood = 430.896 Please wait : Computing the Std Errors Robust Standard Errors (Sandwich formula) Coefficient Std.Error t-value t-prob Cst(V) x 10^6 -1265554.600495 13870 -91.24 0.0000 d-Figarch 0.347425 0.037847 9.180 0.0000 ARCH(Phi1) 1.146765 0.085370 13.43 0.0000 GARCH(Beta1) -1.014973 0.017836 -56.91 0.0000 EGARCH(Theta1) 0.696591 1.7078 0.4079 0.6840 EGARCH(Theta2) 3.097866 0.68398 4.529 0.0000 No Observations Mean (Y) Skewness (Y) Log Likelihood : : : : 139 0.00000 0.10075 430.896 No Parameters Variance (Y) Kurtosis (Y) : : : 0.00005 14.22873 The sample mean of squared residuals was used to start recursion Estimated Parameters Vector : -1265554.600495; 0.347425; 1.146765;-1.014973; 0.696591; 3.097866 Elapsed Time : 0.31 seconds (or 0.00516667 minutes) -Estimating the univariate GARCH model for VrSP - ******************** ** SPECIFICATIONS ** ******************** Dependent variable : VrSP Mean Equation : ARMA (0, 0) model No regressor in the conditional mean Variance Equation : FIEGARCH (1, d, 1) model (Truncation order : 139) No regressor in the conditional variance Normal distribution Weak convergence (no improvement in line search) using numerical derivatives Log-likelihood = 150.044 Please wait : Computing the Std Errors Robust Standard Errors (Sandwich formula) Coefficient Std.Error t-value t-prob Cst(V) -5.188383 0.00012362-4.197e+004 0.0000 d-Figarch 1.917165 0.0024112 795.1 0.0000 ARCH(Phi1) -1.1677351.8911e-008-6.175e+007 0.0000 GARCH(Beta1) -0.731397 0.00014539 -5031 0.0000 EGARCH(Theta1) 0.131814 0.0035363 37.27 0.0000 EGARCH(Theta2) 0.359413 0.0020073 179.1 0.0000 No Observations Mean (Y) Skewness (Y) Log Likelihood : : : : 139 -0.00000 0.04877 150.044 No Parameters Variance (Y) Kurtosis (Y) : : : 0.00982 4.68673 The sample mean of squared residuals was used to start recursion Estimated Parameters Vector : -5.188383; 1.917165;-1.167735;-0.731397; 0.131814; 0.359413 Elapsed Time : 0.623 seconds (or 0.0103833 minutes) -Estimating the univariate GARCH model for VrEX - ******************** ** SPECIFICATIONS ** ******************** Dependent variable : VrEX Mean Equation : ARMA (0, 0) model No regressor in the conditional mean Variance Equation : FIEGARCH (1, d, 1) model (Truncation order : 139) No regressor in the conditional variance Normal distribution Weak convergence (no improvement in line search) using numerical derivatives Log-likelihood = 470.529 Please wait : Computing the Std Errors Robust Standard Errors (Sandwich formula) Coefficient Std.Error t-value Cst(V) x 10^4 0.059412 8.3744 0.007094 d-Figarch -0.205393 0.0045825 -44.82 ARCH(Phi1) 3.423879 0.20089 17.04 GARCH(Beta1) 0.992402 0.0076693 129.4 EGARCH(Theta1) -0.576960 0.067256 -8.579 EGARCH(Theta2) 1.023980 0.075626 13.54 No Observations Mean (Y) Skewness (Y) Log Likelihood : : : : 139 -0.00000 6.25436 470.529 No Parameters Variance (Y) Kurtosis (Y) : : : t-prob 0.9944 0.0000 0.0000 0.0000 0.0000 0.0000 0.00010 48.49811 The sample mean of squared residuals was used to start recursion Estimated Parameters Vector : 0.059412;-0.205393; 3.423879; 0.992402;-0.576960; 1.023980 Elapsed Time : 0.369 seconds (or 0.00615 minutes) ******************** ** SECOND STEP ** ******************* ************* ** SERIES ** ************* #1: VrIR #2: VrSP #3: VrEX ******************************* ** MG@RCH( 1) SPECIFICATIONS ** ******************************* Conditional Variance : Corrected Dynamic Correlation Model (Aielli) Multivariate Normal distribution Weak convergence (no improvement in line search) using numerical derivatives Log-likelihood = 1058.88 Please wait : Computing the Std Errors Robust Standard Errors (Sandwich formula) Coefficient Std.Error t-value t-prob rho_21 0.200525 0.092504 2.168 0.0319 rho_31 -0.032502 0.15638 -0.2078 0.8357 rho_32 -0.205260 0.081318 -2.524 0.0128 alpha 0.526357 0.18666 2.820 0.0055 beta 0.054426 0.056769 0.9587 0.3394 No Observations : 139 No Parameters : 23 No Series : Log Likelihood : 1058.884 Elapsed Time : 0.03 seconds (or 0.0005 minutes) Phụ lục Kết thống kê mô tả kiểm định số đặc tính kinh tế lượng chuỗi tương quan thu từ hồi quy mô hình FIEGARCH(1,d,1)-cDCC (Kết thu từ package G@ARCH phần mềm OxMetrics ver 6.3) TESTS : ======= ADF Test with lags Intercept and no time trend H0: CORR_VrIR_VrSP is I(1) ADF Statistics: -9.40889 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% 5% 10% -3.4323 -2.86229 -2.56721 OLS Results y_1 Constant RSS OBS Coefficient -0.789322 0.146807 3.512490 138.000000 t-value -9.4089 7.0569 Information Criteria (to be minimized) Akaike -0.803960 Shibata Schwarz -0.761536 Hannan-Quinn H0: CORR_VrIR_VrEX is I(1) -0.804372 -0.786720 ADF Statistics: -10.019 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.4323 5% -2.86229 10% -2.56721 OLS Results y_1 Constant RSS OBS Coefficient -0.849305 -0.031358 3.462051 138.000000 t-value -10.019 -2.2500 Information Criteria (to be minimized) Akaike -0.818424 Shibata Schwarz -0.776000 Hannan-Quinn H0: CORR_VrSP_VrEX is I(1) -0.818836 -0.801184 ADF Statistics: -10.2139 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.4323 5% -2.86229 10% -2.56721 OLS Results y_1 Constant RSS OBS Coefficient -0.882546 -0.144273 9.626463 138.000000 t-value -10.214 -5.3620 Information Criteria (to be minimized) Akaike 0.204231 Shibata Schwarz 0.246655 Hannan-Quinn - 0.203819 0.221471 TESTS : ======= ADF Test with lags Intercept and time trend H0: CORR_VrIR_VrSP is I(1) ADF Statistics: -9.47154 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 Constant Trend RSS OBS Coefficient -0.799774 0.148942 0.000365 3.483773 138.000000 t-value -9.4715 7.1289 1.0549 Information Criteria (to be minimized) Akaike -0.797543 Shibata Schwarz -0.733907 Hannan-Quinn H0: CORR_VrIR_VrEX is I(1) -0.798461 -0.771682 ADF Statistics: -9.98585 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 Constant Trend RSS OBS Coefficient -0.849684 -0.031408 -0.000071 3.460933 138.000000 t-value -9.9858 -2.2453 -0.20881 Information Criteria (to be minimized) Akaike -0.804120 Shibata Schwarz -0.740484 Hannan-Quinn H0: CORR_VrSP_VrEX is I(1) -0.805039 -0.778260 ADF Statistics: -10.1748 Asymptotic critical values, Davidson, R and MacKinnon, J (1993) 1% -3.96104 5% -3.41127 10% -3.12748 OLS Results y_1 Constant Trend RSS OBS Coefficient -0.882394 -0.144204 0.000088 9.624785 138.000000 t-value -10.175 -5.3394 0.15340 Information Criteria (to be minimized) Akaike 0.218683 Shibata 0.217764 Schwarz 0.282319 Hannan-Quinn 0.244543 - Ox Professional version 6.30 (Windows/U) (C) J.A Doornik, 1994-2011 Copyright for this package: S Laurent, 2000-2010 G@RCH package version 6.1, object created on 12-07-2016 Copyright for this package: S Laurent, 2000-2010 TESTS : ======= Database information -Sample: 2004(4) - 2015(12) (141 observations) Frequency: 12 Variables: Variable CORR_VrIR_VrSP CORR_VrIR_VrEX CORR_VrSP_VrEX Constant Trend #obs #miss type 139 double 139 double 139 double 141 double 141 double -0.75134 -0.76905 -0.84065 1 mean 0.18627 -0.036884 -0.16432 71 max 0.77679 0.68666 0.7467 141 std.dev 0.16261 0.15964 0.26497 40.702 Series #1/3: CORR_VrIR_VrSP Normality Test Skewness Excess Kurtosis Jarque-Bera - Statistic -1.1263 11.021 732.90 t-Test 5.4789 26.990 NaN P-Value 4.2791e-008 1.9578e-160 7.1339e-160 t-Test 1.2975 24.948 NaN P-Value 0.19446 2.2409e-137 1.2987e-131 t-Test 2.5895 4.0069 NaN P-Value 0.0096103 6.1521e-005 1.6123e-005 Series #2/3: CORR_VrIR_VrEX Normality Test Skewness Excess Kurtosis Jarque-Bera - Statistic -0.26672 10.188 602.75 Series #3/3: CORR_VrSP_VrEX Normality Test Skewness Excess Kurtosis Jarque-Bera - Statistic 0.53233 1.6362 22.071 Phụ lục 10 Kết xác định điểm chuyển đổi chuyễn phần dư thu từ mô hình VAR(1) (Kết thu từ package ‘changepoint’ phần mềm R-Studio) Change-point detection: on residual of USDVND return > require(ggplot2) > require(changepoint) > require(ggfortify) > cpEX cpIR cpSP ChangepointVIR plot(ChangepointVIR, cpt.width=1, xlab="Volatility shifts detection", ylab="VN-Index VAR residual") > cpts.ts(ChangepointVIR) [1] 2006.083 2008.083 2008.833 2009.833 2011.333 2011.583 2013.667 2014.083 2015.083 > ChangepointVIR Class 'cpt' : Changepoint Object ~~ : S4 class containing 12 slots with names date version data.set cpttype method test.stat pen.type pen.value minsegle n cpts ncpts.max param.est Created on : Fri Feb 19 03:38:41 2016 summary(.) : -Created Using changepoint version 2.2 Changepoint type : Change in variance Method of analysis : PELT Test Statistic : Normal Type of penalty : BIC with value, 9.868948 Minimum Segment Length : Maximum no of cpts : Inf Changepoint Locations : 21 45 54 66 84 87 112 117 129 -> ChangepointVEX plot(ChangepointVEX, cpt.width=1, xlab="Volatility shifts detection", ylab="USD/VND VA R residual") > cpts.ts(ChangepointVEX) [1] 2007.750 2009.250 2009.750 2010.083 2010.917 2011.083 > > ChangepointVEX Class 'cpt' : Changepoint Object ~~ : S4 class containing 12 slots with names date version data.set cpttype method test.stat pen.type pen.value minsegle n cpts ncpts.max param.est Created on : Fri Feb 19 03:38:41 2016 summary(.) : -Created Using changepoint version 2.2 Changepoint type : Change in variance Method of analysis : PELT Test Statistic : Normal Type of penalty : BIC with value, 9.868948 Minimum Segment Length : Maximum no of cpts : Inf Changepoint Locations : 41 59 65 69 79 81 -> ChangepointVSP plot(ChangepointVSP, cpt.width=1, xlab="Volatility shifts detection", ylab="USD/VND VA R residual") > cpts.ts(ChangepointVSP) [1] 2004.667 2004.833 2006.000 2009.333 > ChangepointVSP Class 'cpt' : Changepoint Object ~~ : S4 class containing 12 slots with names date version data.set cpttype method test.stat pen.type pen.value minsegle n cpts ncpts.max param.est Created on : Fri Feb 19 03:38:41 2016 summary(.) : -Created Using changepoint version 2.2 Changepoint type : Change in variance Method of analysis : PELT Test Statistic : Normal Type of penalty : BIC with value, 9.868948 Minimum Segment Length : Maximum no of cpts : Inf Changepoint Locations : 20 60 Phụ lục 11 Kết ước lượng hệ số biến giả D xác định ảnh hưởng cú sốc dao động đến tương quan động biến số Model 32: OLS, using observations 2004:07-2015:12 (T = 138) Dependent variable: CORR_VrIR_VrSP const CORR_VrIR_VrSP CORR_VrIR_VrSP_2 CORR_VrIR_VrSP_3 CORR_VrIR_VrSP_4 CORR_VrIR_VrSP_1 Coefficient 0.182551 −0.0387041 −0.209284 0.0719854 0.0242012 0.0780334 Std Error 0.0244261 0.0360352 0.057004 0.0477341 0.0404802 0.0861772 t-ratio 7.4736 −1.0741 −3.6714 1.5080 0.5979 0.9055 p-value