Luận văn thạc sĩ HUS chuỗi thời gian

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - DƢƠNG NHẬT THĂNG CHUỖI THỜI GIAN LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2015 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - DƢƠNG NHẬT THĂNG CHUỖI THỜI GIAN Chuyên ngành: Lý thuyết xác suất thống kê toán học Mã số: 60460106 LUẬN VĂN THẠC SĨ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: PGS.TS PHAN VIẾT THƢ Hà Nội – 2015 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỤC LỤC MỤC LỤC LỜI MỞ ĐẦU CHƯƠNG GIỚI THIỆU VỀ CHUỖI THỜI GIAN 1.1 MÔ TẢ SƠ LƢỢC 1.2 SƠ LƢỢC VỀ KĨ THUẬT 10 1.2.1 XU HƢỚNG 11 1.2.1.1 Phƣơng pháp bình phƣơng cực tiểu 11 1.2.1.2 Bộ lọc (hay phƣơng pháp làm trơn) 12 VỊNG TUẦN HỒN MÙA 17 1.2.2 1.3 CHUYỂN ĐỔI DỮ LIỆU 18 1.4 VÍ DỤ 19 CHƯƠNG 27 LÝ THUYẾT XÁC SUẤT CỦA QUÁ TRÌNH NGẪU NHIÊN 27 2.1 GIỚI THIỆU 27 2.2 QUÁ TRÌNH NGẪU NHIÊN 27 2.3 HÀM TỰ TƢƠNG QUAN MẪU 29 2.4 CÁC VÍ DỤ 30 CHƯƠNG 34 MƠ HÌNH TRUNG BÌNH TRƢỢT TỰ HỒI QUY 34 3.1 GIỚI THIỆU 34 3.2 MƠ HÌNH TRUNG BÌNH TRƢỢT 34 3.3 MƠ HÌNH TỰ HỒI QUY 37 Mối quan hệ tính nhân tính dừng 37 3.3.1 3.3.2 Tiệm cận tĩnh 39 3.3.3 Định lý nhân 40 3.3.4 Cấu trúc hiệp phƣơng sai mơ hình AR 40 3.4 MƠ HÌNH ARMA 41 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 3.5 MƠ HÌNH ARIMA 44 3.6 MƠ HÌNH ARIMA MÙA 46 CHƯƠNG 48 ƢỚC LƢỢNG TRONG MIỀN THỜI GIAN 48 4.1 GIỚI THIỆU 48 4.2 CÁC ƢỚNG LƢỢNG MOMENT 48 4.3 ƢỚC LƢỢNG TRONG MƠ HÌNH TỰ HỒI QUY AR(p) 49 4.4 ƢỚC LƢỢNG CHO MƠ HÌNH TRUNG BÌNH TRƢỢT 51 4.5 ƢỚC LƢỢNG CHO MƠ HÌNH ARMA 53 4.6 ƢỚC LƢỢNG HỢP LÝ CỰC ĐẠI 54 4.7 HỆ SỐ TỰ TƢƠNG QUAN RIÊNG (PACF) 58 4.8 CHỌN LỰA BẬC 62 4.9 PHÂN TÍCH PHẦN DƢ 66 4.10 XÂY DỰNG MƠ HÌNH 67 CHƯƠNG 68 CÁC VÍ DỤ SỬ DỤNG R 68 5.1 GIỚI THIỆU 68 5.2 VÍ DỤ 68 5.3 VÍ DỤ 73 CHƯƠNG 80 DỰ BÁO 80 6.1 GIỚI THIỆU 80 6.2 DỰ ĐOÁN ĐƠN GIẢN 81 6.3 TIỆM CẬN BOX - JENKINS 83 6.4 VÍ DỤ VỀ TÍN PHIẾU KHO BẠC 84 KẾT LUẬN 89 TÀI LIỆU THAM KHẢO 90 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com LỜI MỞ ĐẦU Chuỗi thời gian đƣợc sử dụng nhƣ công cụ hữu hiệu để phân tích kinh tế, xã hội nhƣ nghiên cứu khoa học Nghiên cứu dự báo chuỗi thời gian ln tốn gây đƣợc ý nhà toán học, kinh tế, xã hội học, Các quan sát thực tế thƣờng đƣợc thu thập dƣới dạng chuỗi số liệu Làm để từ chuỗi số liệu khô khan, ta tìm đƣợc mơ hình hay quy luật q trình có đủ sở xác để phản ánh đƣợc chân thực liệu có (kiểm tra) đồng thời lại dự đốn cho thời điểm tƣơng lai chƣa xẩy ra?! Mà việc dự đoán đƣợc tƣơng lai có lẽ ln mong đợi thƣờng trực xã hội lồi ngƣời Chính tầm quan trọng việc phân tích chuỗi thời gian nhƣ vậy, nhiều tác giả nghiên cứu đề xuất cơng cụ để phân tích thời chuỗi thời gian nhƣ sử dụng công cụ thống kê hồi qui, phân tích Furie, mơ hình ARIMA BoxJenkins, Sau có nhiều ngƣời sử dụng mạng Nơron để xử lý tính chất phi tuyến chuỗi số liệu, tìm thấy sách chun khảo vấn đề thí dụ nhƣ Mandic Chambers “ Recurrent neural network and prediction” in vào năm 2001 Một hƣớng khác sử dụng khái niệm mờ để đƣa thuật ngữ “ Chuỗi thời gian mờ” Phƣơng pháp sử dụng chuỗi thời gian mờ đƣợc đƣa từ năm 1994 đến đƣợc tiếp tục nghiên cứu để làm tăng độ xác dự báo Nghiên cứu chuỗi thời gian đòi hỏi kiến thức rộng lớn Thống kê, xác suất, kinh tế khoa học Các phép tính chủ yếu dựa cac số liệu rời rạc, liên tục địi hỏi phải có thuật tốn hay phần mềm tính tốn cho xác suất thống kê chun dụng Có nhiều phần mềm nhƣ nhƣ Eview, S –plus, R, v.v Trong khuôn khổ luận văn học viên trƣờng khoa học nhƣ Đại học Khoa học Tự nhiên – Đại học Quốc Gia Hà Nội, việc nghiên cứu LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com dạng khác chuỗi thời gian, tính chất nó, phƣơng pháp thiết lập mơ hình nghiên cứu có tính chất tảng cần thiết Việc phát triển kĩ thuật dự đốn khác nhau, mơ hình cho chuỗi liệu khơng tuyến tính kinh tế, tài chính, v.v nên dành cho cho nghiên cứu sâu trƣờng chuyên ngành Kinh tế Chính thế, Luận văn tập trung vào chƣơng đầu “Chuỗi thời gian – Các ứng dụng tài với R S-plus“ – Time series Applications to Finance with R and S – Plus (R) tác giả Ngai Hai Chan – The University of Hong Kong Bao gồm Chƣơng 1: Giới thiệu chuỗi thời gian Chƣơng 2: Giới thiệu Lý thuyết xác suất trình ngẫu nhiên Chƣơng 3: Giới thiệu mơ hình trung bình trƣợt tự hồi quy Chƣơng 4: Bàn Ƣớc lƣợng miền thời gian Chƣơng 5: Trình bày ví dụ có sử dụng phần mềm R cho tính tốn đƣợc trình bày chƣơng trƣớc Chƣơng 6: Trình bày Dự đốn Luận văn đƣợc hồn thành dƣới hƣớng dẫn tận tình PGS TS Phan Viết Thƣ, tác giả xin bày tỏ lòng biết ơn chân thành thầy! Tác giả xin chân thành cảm ơn thầy giáo trƣờng ĐH Khoa học Tự Nhiên – Đại học Quốc Gia Hà Nội tham gia giảng dạy giúp đỡ em suốt trình học tập nâng cao trình độ kiến thức trƣờng năm học Cao học Tuy nhiên, điều kiện thời gian khả có hạn nên luận văn khơng thể tránh khỏi thiếu sót Tác giả kính mong thầy bạn góp ý kiến để đề tài đƣợc hồn thiện LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG GIỚI THIỆU VỀ CHUỖI THỜI GIAN 1.1 MÔ TẢ SƠ LƢỢC Chuỗi thời gian đƣợc sử dụng nhƣ cơng cụ hữu hiệu để phân tích kinh tế, xã hội nhƣ nghiên cứu khoa học Nghiên cứu chuỗi thời gian đòi hỏi có số lƣợng lớn quan sát cho đại lƣợng thích hợp để nghiên cứu mối liên hệ đại lƣợng Các quan sát đƣợc tiến hành đặn qua thời kì chẳng hạn theo tháng, quý, năm thời điểm đặc biệt nhƣ thời kỳ xảy khủng hoảng kinh tế Dãy quan sát ta gọi chuỗi thời gian Mơ hình CAPM tiếng mơ hình dao động ngẫu nhiên ví dụ mơ hình tài có chứa cấu trúc chuỗi thời gian Khi nghĩ đến chuỗi thời gian, thƣờng nghĩ đến tập hợp giá trị {X t : t 1, , n} số dƣới t T thời gian t mà mốc X t đƣợc theo dõi CHUỖI RỜI RẠC Chuỗi thời gian rời rạc nhƣ tập số T tập rời rạc (thí dụ, chuỗi báo cáo doanh thu cƣớc phí điện thoại hàng tháng bƣu điện từ tháng năm 2010 đến tháng 12 năm 2014) CHUỖI LIÊN TỤC Chuỗi thời gian gọi liên tục T khoảng liên tục Ví dụ: biểu đồ nhịp tim bệnh nhân hay Biểu đồ theo dõi dƣ trấn dƣới lịng đất cảnh báo sóng thần, động đất tháng CHUỖI LẶP LẠI Dữ liệu đại diện cho phép đo lƣờng lắp lại số thơng qua ngày khác Ví dụ, giám sát doanh thu tuần dựa vào số lƣợng khách hàng siêu thị theo thời gian CHUỖI KẾT HỢP Thay số đo chiều, X t véc-tơ với thành phần cấu thành đại diện cho chuỗi thời gian đơn lẻ Ví dụ, thu LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com thập tổ hợp chứa p phần tử đƣợc viết dƣới dạng X t ( X1t , , X pt )' với X it , i 1,2, , p đại diện cho thu thập phần tử tổ hợp Trong trƣờng hợp này, không ý đến cấu trúc tƣơng quan chuỗi với phần tử mà phải ý đến cấu trúc tƣơng quan chéo phần tử khác CHUỖI PHI TUYẾN, KHÔNG DỪNG VÀ TÍNH KHƠNG ĐỒNG NHẤT Nhiều chuỗi thời gian bắt gặp thực tế phi tuyến tính Trong vài trƣờng hợp chuyển đổi liệu, nhƣng thƣờng phải thiết lập hệ thống phức tạp để tính tốn cho khía cạnh khơng đƣợc quy chuẩn Ví dụ, tính chất khơng đối xứng doanh số bán hàng nghiên cứu mô hình GARCH Mặc dù khía cạnh nêu quan trọng, nhƣng ta bàn luận chủ yếu chuỗi thời gian chuẩn Sau hiểu rõ kĩ thuật độ khó việc phân tích chuỗi thời gian ngắt qng vơ hƣớng giải số khía cạnh khơng đƣợc quy chuẩn Trong thống kê cổ điển, thƣờng giả sử giá trị X độc lập Trong chuỗi thời gian, giá trị X thƣờng có tƣơng quan mục tiêu việc phân tích chuỗi thời gian nhằm sử dụng cấu trúc tƣơng quan chuỗi cho việc xây dựng mơ hình trúc tốt Ví dụ dƣới miêu tả điều quan điểm khoảng tin cậy CI (confidence interval) Ví dụ 1.1 Cho X t at Rõ ràng : E ( X t ) at víi at ~ N (0,1) i.i.d vµ var X t cov( X t , X t k ) E[( X t E[(at , , nh­ vËy )( X t )] k at )(at k at k )] | k | 1, k trường hợp khác LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com n Cho X X t / n th× t var n n n n2 Xt t var X t t n2 n t cov( X t X j ) t j DÔ thÊy X var X n(1 ) (n 1) n n2 (1 2 ) n n [(1 ) ] n n Từ X ~ N ( , X ) Vì khoảng tin cậy 95% (CI) cho X Nếu 2 X [(1 n X )2 khoảng tinh cậy (CI) trở thành X 1/2 ] n khớp với trường hợp n phân bố độc lập phân phối (i.i.d) Sự khác biệt khoảng tin cậy vµ viết dạng L( ) [(1 )2 1/2 ] n Bảng 1 Giá trị khác khoảng tin cậy với n = 50 Ví dụ sử dụng khoảng tin cậy cho việc sai cấu trúc CI gây tốn nhiều thời gian Cấu trúc tương quan thời gian đưa theo mơ hình giúp suy luận tốt trường hợp LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 1.2 SƠ LƢỢC VỀ KĨ THUẬT Tất kỹ thuật phân tích chuỗi thời gian giựa giả định có mẫu hình tiềm ẩn số liệu nghiên cứu với yếu tố ngẫu nhiên ảnh hƣởng lên hệ thống xét Công việc phân tích chuỗi thời gian nghiên cứu kỹ thuật để tách mẫu hình sử dụng nhƣ sở để dự báo cho tƣơng lai Nhìn chung, chuỗi thời gian đƣợc phân tích thành nhiều thành phần vĩ mô vi mô Cấu thành vĩ mô thƣờng đƣợc miêu tả thông qua xu hƣớng , thời vụ chu kỳ , cấu trúc vi mơ cần kết hợp nhiều phƣơng pháp phức tạp để miêu tả Một cách tổng quát có yếu tố cần nghiên cứu:  Xu ( Tt ): Đó thay đổi biến quan trắc Y xét thời gian dài  Chu kỳ tƣợng: (C) :Là thời gian mà tƣợng lặp lại phối hợp với xu ( Tt ) chu kỳ nhiều năm  Biến đổi theo mùa ( St ): Xét đến biến đổi tuần hoàn chu kỳ  Dao động ngẫu nhiên (I ) :Xét đến dạo động ngẫu nhiên xung quanh xu thế, làm ảnh hƣởng đến chu kỳ biến đổi theo mùa quan sát Hình 1 Các đặc trƣng chuỗi thời gian 10 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Tiếp theo, ta đặt mơ hình MA(13) vào liệu Wt W mong đợi có nhiều hệ số khơng > ddacc.2 ddacc.2$model$theta Ta có dự đốn 13 thông số MA(13) -0.67638240 -0.01548842 -0.19295582 -0.21266992 0.25524876 -0.34868415 0.06292423 -0.02186688 0.02991082 0.21487710 0.03356127 -0.60105664 0.62806027 > addacc.3 ddacc.3$model$theta > tsdiag(ddacc.3) Kết -0.4715024 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 -0.6914782 0.2065221 Chỉ có ƣớc lƣợng q Giá trị AIC Do 2log L(  ) -0.4715024, 8547.47 -0.6914782, 12 13 0.2065221 với 2log L(  ) 2( p q 1) với p 0,q 13 834.47 nên xác có thơng số độc q Thống kê Portmanteau đƣợc Hình 5.7 dƣới 77 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (Hình 5.7 Giả thuyết thống kê Portmanteau cho ddacc.3) Trong trƣờng hợp này, giá trị p thống kê Portmanteau cho mơ hình cụ thể dƣờng nhƣ nghèo (poor) (Hình 5.7) Ta thử nghiệm xa với mơ hình bao gồm khác không >ddacc.4 ddacc.4$model$theta Kết -0.6496481 0.0000000 0.0000000 0.0000000 0.0000000 -0.3113243 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 -0.5495682 0.4827663 Bây q thay Giả thuyết thống kê Hình 5.8 hợp lý Mơ hình tổng thể là: (1 B)(1 B12 ) X t 28.83 (1 0.6497 B 0.3113B6 +0.5496B12 0.4828B13 )Zt , 78 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (Hình 5.8 – Giả thuyết thống kê Portmanteau cho ddacc.4) Ngồi ra, để dự đốn, ta dùng lệnh sau R 79 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG DỰ BÁO 6.1 GIỚI THIỆU Khi theo dõi chuỗi thời gian Y1, ,Yn , thƣờng ý đến việc dự báo giá trị tƣơng lai Yn h Cần phân biệt hai loại dự báo, dự báo sau dự báo trƣớc Dự báo sau dự báo đƣợc tiến hành mà giá trị “tƣơng lai” đƣợc biết Nói đại khái giả sử ta biết đƣợc liệu 10 năm, ta lập mô hình dựa số liệu năm đầu, từ mơ hình dự đốn giá trị năm thứ 6, 7, ,10 (dù liệu thực tế năm 6,7, ,10 có!) Loại dự báo dùng nhƣ phƣơng tiện để kiểm tra liệu sẵn có từ đó, đánh giá mơ hình dự đốn đủ tốt hay chƣa Trong đó, dự báo trƣớc đƣa dự báo giá trị tƣơng lai dựa số liệu hiên Bây giả sử có chuỗi liệu Y1, ,Yn sử dụng Y1, ,YT (T n) để ƣớc lƣợng mô hình (Ví dụ ARMA(1,1)) sử dụng mơ hình để ƣớc tính YT 1, ,Yn Do so sánh với YT 1, ,Yn thực tế nên kiểm nghiệm đƣợc tính hợp lý mơ hình Tiếp theo ta dùng mơ hình dự đoán Yn 1, ,Yn h (h 0) chƣa biết Bây ta giới thiệu ba đại lƣợng quan trọng Đầu tiên hàm dự đoán Y n (h) (cách viết khác Y (n, h), Y n Y n Với PY n n h h h hay Ynn h ) định nghĩa nhƣ sau: E (Yn h | Yn , Yn 1, ) PnYn Ynn h h Y1, ,Yn hiểu tổ hợp tuyến tính PY n n h (6.1) a0 a1Yn anY1 Thành phần thứ hai sai số dự báo tƣơng ứng đƣợc định nghĩa en (h) Yn h Ynn h , (6.2) 80 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Thành phần thứ ba phƣơng sai sai số đƣợc cho Pnn h E (Yn Ynn h )2 | Yn , h (6.3) Chú ý trung bình phƣơng sai có điều kiện dựa liệu khứ 6.2 DỰ ĐOÁN ĐƠN GIẢN Nếu chuỗi thời gian sở có xu hƣớng đơn giản nhƣ thảo luận chƣơng 1, Yn đƣợc dự đốn phép ngoại suy đơn Đây đƣợc coi phƣơng pháp hữu ích cần dự đoán khoảng thời gian dài, chuỗi thời gian khơng phù hợp với mơ hình phức tạp Xét Yn mn X n Giả sử mn hằng, chuỗi dừng Do Ynn h PY n n n m h h n, m (6.4)  n , đƣợc ƣớc lƣợng cách làm trơn mũ hàm Ở đại lƣợng cuối cùng, m số mũ Nhớ lại từ chƣơng 1, xu hƣớng đƣợc ƣớc lƣợng thời điểm n đƣợc diễn tả tổ hợp lồi theo dõi dự đốn trƣớc t m Yt (1 Yt (1 )Y t t )m Yt (1  t (t 1, , n) m Y )m 1 Lập lại phép quy hồi, giá trị ƣớc tính thời điểm n mn trọng số trung t bình liệu m t (1 ) j Yt j (1 )t 1Y1 j Khi xu hƣớng khơng hằng, hàm dự đốn viết PY n n h an bn h , 81 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Ở an bn mức hệ số góc hàm xu hƣớng thời điểm n Đầu tiên ta dự đoán mức an cách làm trơn số mũ giá trị giá trị dự đoán cuối nhƣ sau t m Yt (1 Yt (1 )Y t )(at bt ) (t 2, , n) Đối với hệ số góc, trƣớc hết cân nhắc đến: Ynn PY n n an 2bn Mặc khác, Yn sẵn có, tính Ynn 21 Pn 1Yn a n b n Do Ynn Ynn 21 dự đoán giá trị tƣơng lai Yn b n  an a n (a t 1 2b n Do viết b t a t ) b t nhƣ sau b t (a t a t ) (1 cân giá trị ta có dƣới dạng tổ hợp lồi )b t (t 2, , n) Bằng cách kết hợp phép truy hồi dự đốn mức hệ số góc với việc cố định điều kiện ban đầu a Y2 b Y2 Y1 a n b n tính đƣợc thu đƣợc dự đoán Ynn h Tƣơng tự cách làm trơn mũ, giá trị n v¯ ei2 cực tiểu hóa đƣợc chọn tổng bình phƣơng sai số dự đoán i Trong thực hành v¯ nằm khoảng từ 0.1 đến 0.3 Kĩ phƣơng pháp Holt Winters cho xu hƣớng thời vụ nàycó thể đƣợc tìm Kendall and Ord (1990) 82 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com 6.3 TIỆM CẬN BOX - JENKINS Phƣơng pháp đề cập đến việc xây dựng mơ hình ARIMA phù hơp sử dụng mơ hình vào dự báo Thảo luận sâu tìm Abraham and Ledolter (1983) Cho Yt ARMA( p, q) BYt tuân theo quan hệ nhân mơ hình Zt khí Yt i Zt i , Ta có nhận định i sau Điều thứ hai có E (Ys | Yn , ) Ys Ysn s n, s n E (Zs | Yn , ) Zs s n, s n đƣợc Yn i Zn i , 1, i E(Zs | Yn , ) E(Zs | Zn , ) Thêm vào đó, hàm dự báo thỏa mãn Ynn h E (Yn h | Yn , ) i E (Z n h i | Yn , ) i i Zn h i , i Kéo theo sai số dự đoán trở thành en (h) Yn h Ynn h i Zn h i i Đặc biệt, phƣơng sai sai số dự đoán đƣợc thể h Pnn h E (en (h)2 | Yn , ) 2 Z i i Do vậy, miễn sai số dự đoán (sự cải tiến) đƣợc phân phối chuẩnkhoảng tin cậy xác suất (1 sau Ynn h z /2 ) cho giá trị dự đốn đƣợc xây dựng nhƣ Pnn h , Khi z biểu thị cho phân vị thứ biến ngẫu nhiên quy chuẩn 83 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Ví dụ 6.1 Giả sử mơ hình SARIMA(1,0,0) x(0,1,1)12 đƣợc đặt phù hợp vào liệu Yn Yn 12 (Yn Yn 13 ) Zn Zn 12 (6.5) Từ chỗ E(Zn+1 | Yn , ) = ta có Ynn E(Yn | Yn ,Yn 1, ) Yn 11 (Yn Yn 12 ) Z n 11 Từ Ynn E (Yn | Yn , Yn , ) E (Yn Yn 10 10 (Yn Yn 11 ) Z n (Ynn Yn 11 ) Zn 10 | Yn , ) Z n 10 (6.6) Q trình tiếp tục lặp lại đệ quy Ví dụ, Yn có sẵn, theo (6.6) ta có Ynn 21 Yn 10 (Yn Yn 11 ) Yn 10 (Ynn Yn 11 ) Ynn (Yn Ynn Z n 1 Z n 10 | Z n 10 Ynn Yn Ynn )  Đối với mơ hình ARMA( p, q) tổng qt, phƣơng trình tƣơng tự dùng để tính giá trị dự đốn Nhắc lại rằng, giá trị xấp xỉ cho Z n xuất dự báo đƣợc tính đệ quy với n p 1, p 2, băng cách giải Z n phƣơng trình ARMA với giả sử Zn (n p) 6.4 VÍ DỤ VỀ TÍN PHIẾU KHO BẠC Xét ví dụ Lãi suất tín phiếu kho bạc nhà nƣớc kì hạn tháng từ năm 1950 đến tháng năm 1988, liệu đƣợc lƣu file Ustbill.dat gồm 462 liệu Dùng lệnh R ta có biểu đồ chuỗi thời gian, hàm tự tƣơng quan ACF chuỗi chuyển đổi Logarit đƣợc cho Hình 6.1 Lệnh R nhƣ sau 84 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com > ttbill ttbill ttbill ttbill plot.ts(ttbill,main="Chuỗi ttbill - hình a") > acf(ttbill,main="ACF ttbill") > plot.ts(log(ttbill),main="Chuỗi Log(ttbill) - hình c") > acf(log(ttbill),main="ACF chuỗi Log(ttbill) - hình d") ( Hình 6.1 – Chuỗi thời gian, ACF tttbill Log(ttbill) ) Từ biểu đồ thấy liệu khơng dừng trung trung bình phƣơng sai Vì ta lấy sai phân liệu chuyển đổi “Log(ttbill)” đặt tên “dlntbill” kí hiệu Wt Tính ACF, PACF biểu diễn ba biểu đồ chuỗi dlntbill hình 6.2 > dlntbill plot.ts(dlntbill,main='Chuỗi phƣơng sai dlntbill - hình a') > acf(dlntbill,main='ACF dlntbill - hình b') > acf(dlntbill,30,type='partial',main='PACF dlntbill hình c') 85 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (Hình 6.2 Chuỗi thời gian, ACF, PACF với dlntbill) Trong hình 6.2, ta thấy có tƣơng quan mạnh độ trễ 1, có thẻ 17 Do tƣơng quan trễ 17 không đáng kể đƣợc tạo nhiều yếu tố (ví dụ phƣơng sai khơng dừng, ) Vì thế, ta bắt đầu mơ hình AR(6), sau MA(6) cuối ARMA(6,6) cho dlntbill Xét chuẩn đốn thơng tin AIC ARMA(6,6) MA(6) thuận lợi AR(6) (Xem biểu đồ chuẩn đốn thặng dƣ mơ hình AR(6) Hình 6.3) Vì ta chọn mơ hình AR(6) phù hợp dự đốn tƣơng lại dựa vào mơ hình 86 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com (Hình 6.3 – Biểu đồ chuẩn đoán thặng dƣ AR(6)) Ta chọn mơ hình AR(6) cho liệu Wt ( B)Yt Yt Wt Wt , Wt (1 B) X t , X t ( B) 0.45B 0.2B2 Zt , log(ttbill ) 0.09B3 0.04B4 0.01B5 0.21B6 Cụ thể lệnh R > d3tsdiag(d3) > d33 tfore tfore[1] for(i in 2:6){tfore[i] tfore [1] 5.820030 5.896966 5.960716 5.957511 5.899653 5.864232 87 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Lệnh dự đốn giá trị từ mơ hình Y t t lại ta có W Y t W với h 1, ,6 t h  Xt h W X t h Y t h W Lấy lại X t từ Wt :  t h e X t h exp(Y t h W  Lấy lại ttbill X t h ) h t W h h W tính ngƣợc h  Xt h Vẽ biểu đồ mơ tả giá trị dự đốn liệu thực tế năm 1988 cho Hình 6.4 nhƣ sau > ttbill tfore ts.plot(ttbill[457:462],tfore,lty=c(1:2)) > leg.names legend(locator(1),leg.names,lty=c(1,2)) t Dự đoán Thực tế Tháng Tháng Tháng Tháng Tháng Tháng 5.820030 5.896966 5.960716 5.957511 5.899653 5.864232 5.81 5.66 5.70 5.91 6.26 6.46 (Hình 6.4 – Biểu đồ giá trị thực giá trị dự đốn) Trong Hình 6.4, mẫu dự đốn cho tháng đầu năm 1988 hợp lý, nhiên tín phiếu bắt đầu đột ngột tăng cao vào tháng 5/1988 mẫu khơng đủ dự đốn kịp Ngồi ra, xu hƣớng chung giá trị đƣợc dự đoán theo 88 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com giá trị nhận đƣợc trƣớc Đây khiếm khuyết việc sử dụng mơ hình ARIMA để dự đốn chuỗi tài bất ổn KẾT LUẬN Luận văn bàn vấn đề phân tích chuỗi thời gian dừng Các kĩ thuật, mã lệnh thực hành với R vài lý thuyết trình ngẫu nhiên dừng đƣợc cố gắng trình bày rõ Tuy nhiên đối tƣợng mà luận văn đề cập đến chuỗi thời gian dừng, cịn chuỗi thời gian khơng dừng khơng quy chuẩn tài chƣa đƣợc đề cập tới hạn hẹp thời gian nhƣ trình độ cá nhân ngƣời viết luận văn Với cố gắng thân cộng với giúp đỡ thầy cô Trƣờng Đại học Khoa học Tự nhiên – ĐHQG Hà Nội, đặc biệt PGS.TS Phan Viết Thƣ, em hoàn thành xong luận văn Em xin chân thành cảm ơn tiếp đọc góp ý GS.TSKH Đặng Hùng Thắng, TS Nguyễn Thịnh, TS Nguyễn Hồng Hải, TS Trần Mạnh Cƣờng, TS Tạ Ngọc Ánh giúp đỡ em hoàn thành 89 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com TÀI LIỆU THAM KHẢO Ngai Hai Chan (2010) , “Time series Application to finance with R and S-plus – 2nd ed” , Published by John Wiley & Sons, Inc., Hoboken, New Jersey Nguyễn Thị Kim Loan, “Mơ hình chuỗi thời gian mờ dự báo chuỗi thời gian”, Luận văn thạc sỹ công nghệ thông tin, Đại Học Thái Nguyên Nguyễn Văn Hữu – Nguyễn Hữu Dƣ, “Phân tích thống kê dự báo”, NXB Đại học Quốc gia Hà Nội 90 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự - Hạnh phúc ============ GIẤY XÁC NHẬN CHỈNH SỬA LUẬN VĂN THẠC SĨ học viên cao học Dƣơng Nhật Thăng Tên là: DƢƠNG NHẬT THĂNG, tác giả luận văn với đề tài: “Chuỗi thời gian” đƣợc bảo vệ trƣớc Hội đồng chấm luận văn thạc sĩ trƣờng Đại học Khoa học Tự nhiên ngày 24 tháng 12 năm 2015 Theo góp ý Hội đồng, tơi xin bổ sung chỉnh sửa nội dung sau: Sửa thuật ngữ: ”Bằng phân bố” (trang 26) thành ” phân bố xác suất đồng thời chúng không phụ thuộc vào thời gian.” Sửa thuật ngữ “Hội tụ phân bố” (trang 48) thành “Hội tụ theo phân bố” Sửa tên đề mục “4.7 Ước lượng hệ số tự tương quan riêng PACF” thành “4.7 Hệ số tự tương quan riêng PACF” Tôi xin trân trọng đề nghị Hội đồng xác nhận việc chỉnh sửa, cho phép đƣợc làm thủ tục xin cấp Thạc sĩ khoa học Hà Nội, ngày 26 tháng 12 năm 2015 CÁN BỘ HƢỚNG DẪN KHOA HỌC HỌC VIÊN CAO HỌC PGS.TS Phan Viết Thƣ Dƣơng Nhật Thăng CHỦ TỊCH HỘI ĐỒNG CHẤM LUẬN VĂN THẠC SĨ KHOA HỌC GS.TSKH Đặng Hùng Thắng 91 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ... chứa cấu trúc chuỗi thời gian Khi nghĩ đến chuỗi thời gian, thƣờng nghĩ đến tập hợp giá trị {X t : t 1, , n} số dƣới t T thời gian t mà mốc X t đƣợc theo dõi CHUỖI RỜI RẠC Chuỗi thời gian rời rạc... CHƯƠNG GIỚI THIỆU VỀ CHUỖI THỜI GIAN 1.1 MÔ TẢ SƠ LƢỢC Chuỗi thời gian đƣợc sử dụng nhƣ công cụ hữu hiệu để phân tích kinh tế, xã hội nhƣ nghiên cứu khoa học Nghiên cứu chuỗi thời gian địi hỏi có... đƣa thuật ngữ “ Chuỗi thời gian mờ” Phƣơng pháp sử dụng chuỗi thời gian mờ đƣợc đƣa từ năm 1994 đến đƣợc tiếp tục nghiên cứu để làm tăng độ xác dự báo Nghiên cứu chuỗi thời gian đòi hỏi kiến