(Tiểu luận FTU) dự báo chuỗi số liệu về lượng du khách quốc tế đến việt nam

TRƯỜNG ĐẠI HỌC NGOẠI THƯƠNG KHOA KINH TẾ QUỐC TẾ ******************** TIỂU LUẬN DỰ BÁO KINH TẾ ĐỀ TÀI: DỰ BÁO CHUỖI SỐ LIỆU VỀ LƯỢNG DU KHÁCH QUỐC TẾ ĐẾN VIỆT NAM Giảng viên giảng dạy: ThS.Nguyễn Thúy Quỳnh Lớp tín : KTE418(1-1920).2_LT Danh sách thành viên nhóm 15: Nguyễn Phùng Linh Giang 1614410045 Phan Thị Huyền 1614410085 Lê Phương Thảo 1614410168 Trần Phương Thảo 1614410171 Hà Nội, tháng 12 năm 2019 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com ĐIỂM ĐÁNH GIÁ CÁ NHÂN Nguyễn Phan Thị Lê Phương Trần Phương Phùng Huyền Thảo Thảo Linh Giang Linh Giang 10 10 10 10 Phan Huyền 10 10 10 10 Lê Thảo 10 10 10 10 Trần Thảo 10 10 10 10 Điểm trung bình 10 10 10 10 cá nhân LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MỤC LỤC DANH MỤC BẢNG BIỂU LỜI MỞ ĐẦU CHƯƠNG 1: TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU CHƯƠNG 2: KHẢO SÁT DỮ LIỆU Những đặc tính chuỗi số liệu 2 Kiểm định tính mùa vụ Kruscal Wallis 10 Kiểm định tính dừng chuỗi phân tích giản đồ tự tương quan ACF giản đồ tự tương quan riêng phần PACF 12 CHƯƠNG 3: LỰA CHỌN PHƯƠNG PHÁP DỰ BÁO 18 San mũ Winters 18 ARIMA 19 2.1 Các bước dự báo 19 2.2 Kiểm định phân phối chuẩn 21 2.3 Kiểm định tự tương quan 22 2.4 Kiểm định phương sai sai số thay đổi 22 3 Lựa chọn mô hình tốt 23 CHƯƠNG 4: KẾT QUẢ DỰ BÁO 24 KẾT LUẬN 26 DANH MỤC TÀI LIỆU THAM KHẢO 27 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com DANH MỤC BẢNG BIỂU Bảng 2.1: Biểu đồ thể lượng du khách quốc tế đến Việt Nam từ tháng 1/2012 đến tháng 10/2019 Bảng 2.2: Biểu đồ mơ tả tính mùa vụ chuỗi visitor 10 Bảng 2.3: Kiểm định tính mùa vụ Kruscal Wallis 11 Bảng 2.4: Scaling Factors sau hiệu chỉnh yếu tố mùa vụ cho chuỗi 12 Bảng 2.5: Kiểm định tính dừng chuỗi sau hiệu chỉnh mùa vụ Unit Root Test 13 Bảng 2.6: Kiểm định tính dừng sai phân bậc chuỗi sau hiệu chỉnh mùa vụ Unit Root Test 13 Bảng 2.7: Giản đồ tự tương quan tự tương quan riêng phần sai phân bậc chuỗi visitorsa 15 Bảng 2.8: Kiểm định tính dừng cho sai phân bậc chuỗi sau hiệu chỉnh mùa vụ Unit Root Test 16 Bảng 2.9: Giản đồ tự tương quan tự tương quan riêng phần sai phân bậc chuỗi visitorsa 17 Bảng 3.1: Bảng kết san mũ Winters 18 Bảng 3.2: Bảng kết ước lượng để chọn bậc p, q tối ưu 20 Bảng 3.3: Bảng kết hồi quy mơ hình MA(1) 20 Bảng 3.4: Két dự báo cho chuỗi visitorsa mơ hình MA(1) 21 Bảng 3.5: Kiểm định phân phối chuẩn cho mơ hình MA(1) 22 Bảng 3.6: Bảng so sánh sai số dự báo mơ hình san mũ Winters mơ hình MA(1) 23 Bảng 4.1: Bảng kết giá trị dự báo cho chuỗi visitor từ tháng 11/2019 đến tháng 12/2020 25 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com LỜI MỞ ĐẦU Hiện nay, du lịch coi ngành kinh tế tổng hợp phát triển Nhiều quốc gia lấy du lịch làm ngành kinh tế trọng điểm góp phần thúc đẩy ngành kinh tế khác phát triển, tăng thu ngoại tệ, tạo nhiều việc làm nguồn thu nhập cho người dân địa từ nâng cao mức sống họ Việt Nam nước có kinh tế đà phát triển mạnh mẽ đẩy mạnh phát triển ngành du lịch xem giải pháp hữu hiệu để thúc đẩy kinh tế, tạo nhiều hội giao lưu, hội nhập quốc tế Xét vị trí địa lý lẫn điều kiện tự nhiên, với vị trí nằm trung tâm khu vực Đông Nam Á, thiên nhiên khí hậu nhiệt đới ơn hịa, người thân thiện văn hóa ẩm thực phong phú, Việt Nam xem điểm đến đầy thú vị dành cho khách du lịch quốc tế Do việc phân tích thống kê lượng khách du lịch quốc tế đến Việt Nam cần thiết để đánh giá quy mô, hiệu hoạt động ngành có sở để dự báo lập kế hoạch cho phát triển ngành du lịch năm Nhận thấy tầm quan trọng vấn đề này, nhóm nghiên cứu xin phép chọn đề tài : “Dự báo chuỗi số liệu lượng du khách quốc tế đến Việt Nam” làm đề tài nghiên cứu Mục tiêu nghiên cứu khảo sát chuỗi liệu nhằm rút đặc điểm chuỗi liệu, lựa chọn phương pháp dự báo phù hợp, đánh giá kết dự báo từ lựa chọn mơ hình dự báo phù hợp tốt cho chuỗi số liệu Trong nghiên cứu này, biến lựa chọn để dự báo visitor- số lượng khách du lịch nước ngồi đến Việt Nam Số liệu nhóm thu thập số liệu theo tháng thời gian từ tháng năm 2012 đến tháng 10 năm 2019, liệu trích nguồn từ CSDL Thống kê Du lịch Tổng cục du lịch Việt Nam Khoảng thời gian dự báo từ tháng 11 năm 2019 đến tháng 12 năm 2020 Nhóm nghiên cứu xin gửi lời cảm ơn chân thành đến giáng viên, Thạc sỹ Nguyễn Thúy Quỳnh tận tình dạy dỗ, bảo kiến thức môn Dự báo kinh tế phần mềm Eview để nhóm hồn thành tiểu luận Tuy nhiên, hạn chế kiến thức kinh nghiệm thực tiễn giới hạn thời gian, tiểu luận nhóm khơng thể LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com tránh khỏi thiếu sót nội dung lẫn hình thức Nhóm nghiên cứu mong nhận bảo, góp ý để tiểu luận hoàn thiện LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG 1: TỔNG QUAN TÌNH HÌNH NGHIÊN CỨU Hiện có số nghiên cứu dự báo lượng khách du lịch quốc tế đến Việt Nam nhiều phương pháp khác Nhóm nghiên cứu xin trích dẫn số nghiên cứu tiêu biểu Tác giả Mai Phương với luận văn nghiên cứu Các phương pháp thống kê nghiên cứu lượng khách du lịch phương pháp thống kê nghiên cứu lượng khách du lịch tiêu biểu sau: Phương pháp số tương đối số tuyệt đối, Dãy số thời gian, Hồi quy tương quan Trong bật phương pháp dãy số thời gian, từ số liệu thực tế qua việc sử dụng dãy số thời gian, tác giả tiến hành tính tốn tiêu, phân tích, tìm quy luật, kết luận số lượng khách du lịch dự đoán cho tương lai Trong nghiên cứu này, tác giả số phương pháp dự đoán thống kê ngắn hạn thường dùng du lịch là:  Dự đốn dựa vào lượng tăng (giảm) tuyệt đối trung bình: Điều kiện áp dụng : Áp dụng trường hợp dãy số có lượng tăng (giảm) tuyệt đối liên hồn xấp Mơ hình có dạng : Trong đó:- h: Thời hạn dự đốn (tầm xa dự đốn) - : Mức độ dùng làm gốc ngoại suy - : mức độ cuối dãy số, mức độ bình qn vài thời kỳ cuối dãy số - ̅ :Là lượng tăng (giảm) tuyệt đối bình quân : Là giá trị dự đoán điểm thứ n+h Do yêu cầu tài liệu đơn giản nên phương pháp phù hợp với việc nghiên cứu lượng khách du lịch, nhiên hạn chế mức độ xác chưa cao phụ thuộc vào mức độ đầu cuối, bỏ qua biến động tượng thời gian nghiên cứu  Dự đoán dựa vào hàm xu biến động thời vụ: LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Mơ hình dự đốn : Trong đó: = : Mức độ dự đốn thời gian thứ I thuộc năm j i: Tháng, quý; i=1, j: Năm: j=1, : Là mức độ dự đoán năm j : Là số thời vụ thời gian thứ i Trong ngành du lịch nay, phương pháp sử dụng tương đối phổ biến phán ảnh xác biến động tổng số lượng khách quốc tế đén Việt Nam qua năm Tuy nhiên phương pháp đòi hỏi phải có số liệu tháng, quý Tác giả Phạm Quang Hưng có viết Một số phương pháp dự báo du lịch tập trung đề cập đến hai phương pháp dự báo đơn giản ứng dụng nhiều phương pháp ngoại suy khơng đổi phương pháp ngoại suy bình quân giản đơn, số nhu cầu lượng khách du lịch quốc tế dự kiến đón làm đối tượng xác định cần dự báo  Phương pháp ngoại suy không đổi: Điều kiện áp dụng: Thường áp dụng bối cảnh tình hình tăng trưởng du lịch tương đối ổn định, bị tác động yếu tố độc lập (tác động bên ngoài) Các bước tiến hành: -Bước 1: Thu thập liệu lượng khách du lịch quốc tế đón giai đoạn hành khứ -Bước 2: Tính tỷ lệ % tăng giảm lượng khách qua năm tỷ lệ chênh lệch -Bước 3: Sử dụng tỷ lệ % thay đổi so với giai đoạn hành làm sở dự báo cho giai đoạn kỳ năm tới Công thức dự báo: Phương trình a: Ft = At - Phương trình b: Ft = At – × (At – 1/At – 2) LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Trong đó: F = Giá trị cần dự báo; A = giá trị thực tế đạt được; t = Mốc thời gian cần dự báo  Phương pháp ngoại suy bình quân giản đơn Điều kiện áp dụng: Khi muốn dự báo chắn với số tương đối an toàn bao hàm cà yếu tố khập khiễng, thất thường tăng trường, yếu tố không mong đợi Các bước tiến hành: Bước 1: Thu thập liệu nhu cầu du lịch cần dự báo (Lượng khách du lịch quốc tế đón được) đạt giai đoạn hành q khứ Bước 2: Tính bình qn tỷ lệ % tăng giảm qua năm Bước 3: Sử dụng tỷ lệ % thay đổi so với giai đoạn hành làm sở dự báo cho số nhu cầu du lịch giai đoạn kỳ năm tới Cơng thức dự báo: Phương trình: Ft = (At-1 + At-2 + At-3)/n Trong đó: F = Giá trị cần dự báo; A = Giá trị thực tế đạt được; t= Mốc thời gian cần dự báo; n = Khoảng chuỗi thời gian khứ Từ tổng quan tình hình nghiên cứu liên quan đến đề tài nghiên cứu cho thấy rằng, công tác dự báo vô cần thiết quan trọng Trên thực tế có nhiều phương pháp phổ biến để thực dự báo, nhiên tùy thuộc vào thực tế tình hình thực tế, mục tiêu dự báo khả sử dụng nguồn số liệu mà cá nhân, tổ chức định lựa chọn phương pháp riêng, phù hợp với khả điều kiện cho phép Qua trình tìm hiểu tiếp cận nhiều nguồn số liệu, nhóm nghiên cứu tập trung lựa chọn phương pháp san mũ Winters ARIMA để làm phương pháp dự báo cho lượng khách quốc tế đến Việt Nam khoảng thời gian tháng 11 năm 2019 đến tháng 12 năm 2020 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG 2: KHẢO SÁT DỮ LIỆU 2.1 Những đặc tính chuỗi số liệu Số liệu mà nhóm tác giả sử dụng lượng khách du lịch quốc tế đến Việt Nam từ tháng năm 2012 đến tháng 10 năm 2019, gồm 94 quan sát Số liệu nhóm tổng hợp từ nguồn CSDL Thống kê Du lịch Tổng cục Du lịch Việt Nam, thuộc Bộ Văn hóa, Thể thao Du lịch Website http://thongke.tourism.vn/ Mở cửa sổ Series: Visitor vào View/ Descriptive Statistics & Tests/ Stats Table Ta bảng mô tả số thống kê quan trọng sau:  Số quan sát (Observations): 94  Giá trị trung bình (Mean): 881892.4  Số trung vị (Median): 771968.5  Giá trị lớn (Maximum: 1618337  Giá trị nhỏ (Minimum): 417429.0  Độ lệch chuẩn (Std Dev.): 323506.6 Trên cửa sổ Series: Visitor vào View/ Graph/ Line and symbol, ta có biểu đồ mô tả số liệu: LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 2.5: Kiểm định tính dừng chuỗi sau hiệu chỉnh mùa vụ Unit Root Test Null Hypothesis: VISITORSA has a unit root Exogenous: Constant Lag Length: (Automatic - based on SIC, maxlag=11) Augmented Dickey-Fuller test statistic Test critical values: t-Statistic Prob.* 2.002222 0.9999 1% level -3.503049 5% level -2.893230 10% level -2.583740 Nguồn: Eview Từ bảng kết quả, ta thấy: P-value = 0.9999 > α = 0.05 → Khơng có sở bác bỏ H0 → Chuỗi visitorsa không dừng Ta thử lấy sai phân bậc chuỗi visitorsa Trên cửa sổ Series: Visitorsa vào lại View/ Unit Root Tests, phần Test for unit root in chọn 1st difference, ta được: Bảng 2.6: Kiểm định tính dừng sai phân bậc chuỗi sau hiệu chỉnh mùa vụ Unit Root Test Null Hypothesis: D(VISITORSA) has a unit root Exogenous: Constant Lag Length: (Automatic - based on SIC, maxlag=11) t-Statistic Prob.* Augmented Dickey-Fuller test statistic -11.41728 0.0001 Test critical values: 1% level -3.503049 5% level -2.893230 10% level -2.583740 Nguồn: Eview 13 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Từ bảng kết quả, ta thấy: P-value = 0.0001 < α = 0.05 → Có thể bác bỏ H0 →Chuỗi Visitorsa dừng sai phân bậc Để xác định giá trị p, q cho mô hình ARIMA, ta phải dựa vào giản đồ tự tương quan ACF tự tương quan riêng phần PACF Trên cửa sổ Series: Visitorsa vào View/ Correlogram Trên cửa sổ Correlogram Specification, chọn Correlogram of 1st difference độ trễ phần Lags to include 36 Ta có giản đồ tự tương quan tương quan riêng phần sau: 14 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 2.7: Giản đồ tự tương quan tự tương quan riêng phần sai phân bậc chuỗi visitorsa Nguồn: Eview Ta thấy, khơng có độ trễ vượt khỏi đường biên Vậy, ta thử lấy sai phân bậc chuỗi sau: Trên cửa sổ Series: Visitorsa vào lại View/ Unit Root Tests, phần Test for unit root in chọn 2nd difference, ta được: 15 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 2.8: Kiểm định tính dừng cho sai phân bậc chuỗi sau hiệu chỉnh mùa vụ Unit Root Test Null Hypothesis: D(VISITORSA,2) has a unit root Exogenous: Constant Lag Length: (Automatic - based on SIC, maxlag=11) t-Statistic Prob.* Augmented Dickey-Fuller test statistic -9.947685 0.0000 Test critical values: 1% level -3.505595 5% level -2.894332 10% level -2.584325 Nguồn: Eview Từ bảng kết quả, ta thấy: P-value = 0.0000 < α = 0.05 → Có thể bác bỏ H0 →Chuỗi Visitorsa dừng sai phân bậc Trên cửa sổ Series: Visitorsa vào View/ Correlogram Trên cửa sổ Correlogram Specification, chọn Correlogram of 2nd difference độ trễ phần Lags to include 36 Ta có giản đồ tự tương quan tương quan riêng phần sau: 16 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 2.9: Giản đồ tự tương quan tự tương quan riêng phần sai phân bậc chuỗi visitorsa Nguồn: Eview Từ giản đồ ACF PACF trên, ta chọn độ trễ cho AR 1,2,3,6 MA 17 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG 3: LỰA CHỌN PHƯƠNG PHÁP DỰ BÁO Nhóm tác giả định lựa chọn phương pháp san mũ Winters (chuỗi có tính xu mùa vụ), ARIMA (chuỗi dừng sai phân bậc 2, có tính mùa vụ) để dự báo cho chuỗi 3.1 San mũ Winters Xét thấy chuỗi visitor có tính xu mùa vụ, nhóm chọn phương pháp san mũ Winters để tiến hành báo chuỗi Kích đúp chuột trái vào chuỗi visitor -> Proc -> Exponential Smoothing -> Simple Exponential Smoothing -> Smoothing Method: Holts-Winters – Multiplicative -> Smoothed series: visitorwm -> OK Ta bảng kết sau: Bảng 3.1: Bảng kết san mũ Winters Parameters: Alpha 0.7800 Beta 0.0000 Gamma 0.0000 Sum of Squared Residuals 2.67E+11 Root Mean Squared Error 53297.97 End of Period Levels: Mean 1778629 Trend 10309.56 Seasonals: 2018M11 1.020312 2018M12 1.029115 2019M01 1.148441 2019M02 1.186161 2019M03 1.038712 2019M04 1.073662 2019M05 0.921972 2019M06 0.850721 2019M07 0.928087 2019M08 1.009870 2019M09 0.894196 2019M10 0.898751 Tính MAPE MRSE mơ hình theo cơng thức: 18 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com MAPE: genr mape = @mean(@abs(visitor-visitorwm)/visitor)  MAPE= 0.047684= 4.7% MRSE: genr mrse = @sqrt (@mean(visitor-visitorwm)^2)  Các MRSE= 3756.634 giá trị dự báo tính Excel theo cơng thức: visitorf=(Mean+h*Trend)*Seasonals (với h số kỳ sau giá trị san cuối cùng) 3.2 ARIMA 3.2.1 Các bước dự báo Xác định p, q giản đồ tự tương quan ACF giản đồ tự tương quan PACF: Kích đúp chuột trái vào chuỗi sau hiệu chỉnh yếu tố mùa vụ visitorsa -> vào View -> Correlogram, hộp thoại Correlogram Specification xuất Ở khung Correlogram of, chọn 2nd difference, ô Lags to include để mặc định độ trễ 36 Ta thu kết sau: Ở giản đồ PACF, độ trễ có vạch vượt đường biên 1, 2, 3, -> p=1, 2, 3, Ở giản đồ ACF, độ trễ có vạch vượt ngồi đường biên -> q=1 Ta ước lượng mơ hình với AR(1, 2, 3, 6) MA(1) Chạy lệnh: ls d(visitorsa,2) c ar(1)ar(2)ar(3)ar(6)ma(1) kết quả: 19 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 3.2: Bảng kết ước lượng để chọn bậc p, q tối ưu Variable Coefficient Std Error t-Statistic Prob C 351.6397 222.4514 1.580748 0.1179 AR(1) -0.229459 0.116149 -1.975557 0.0517 AR(2) -0.070108 0.136358 -0.514147 0.6086 AR(3) -0.113110 0.131320 -0.861328 0.3916 AR(6) -0.034837 0.119609 -0.291252 0.7716 MA(1) -0.972348 0.021734 -44.73811 0.0000 Nguồn: Eview Loại AR(1), AR(2), AR(3), AR(6) hệ số khơng có ý nghĩa thống kê Ước lượng với MA(1) Chạy lệnh: ls d(visitorsa,2) c ma(1) kết quả: Bảng 3.3: Bảng kết hồi quy mô hình MA(1) Variable Coefficient Std Error t-Statistic Prob C 468.6368 219.5264 2.134763 0.0355 MA(1) -0.987412 0.015813 -62.44226 0.0000 R-squared 0.588190 Mean dependent var 248.7000 Adjusted R-squared 0.583614 S.D dependent var 84158.71 S.E of regression 54305.92 Akaike info criterion 24.66415 Sum squared resid 2.65E+11 Schwarz criterion 24.71898 Log likelihood -1132.551 Hannan-Quinn criter 24.68628 Durbin-Watson stat 2.440068 F-statistic 128.5471 Prob(F-statistic) 0.000000 Inverted MA Roots 99 Nguồn: Eview Kết hồi quy cho thấy hệ số ước lượng có ý nghĩa thống kê mức ý nghĩa 5% 20 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Từ bảng kết ước lượng chọn Forecast -> OK Ta kết dự báo cho chuỗi visitorsa mẫu sau: Bảng 3.4: Két dự báo cho chuỗi visitorsa mô hình MA(1) 5,000,000 Forecast: VISITORSAF Actual: VISITORSA Forecast sample: 2012M01 2019M10 Adjusted sample: 2012M03 2019M10 Included observations: 92 Root Mean Squared Error 347368.4 Mean Absolute Error 270856.5 Mean Abs Percent Error 26.89489 Theil Inequality Coefficient 0.157260 Bias Proportion 0.551339 Variance Proportion 0.376506 Covariance Proportion 0.072155 4,000,000 3,000,000 2,000,000 1,000,000 -1,000,000 2012 2013 2014 2015 2016 VISITORSAF 2017 2018 2019 ± S.E Nguồn: Eview  MAPE = 26.89%  MRSE = 347368.4 3.2.2 Kiểm định phân phối chuẩn Cặp giả thuyết: H0: Biến có phân phối chuẩn H1: Biến khơng có phân phối chuẩn Từ cửa số ước lượng chọn View -> Residual Diagnostics -> Histogram Normality Test Ta kết quả: 21 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 3.5: Kiểm định phân phối chuẩn cho mơ hình MA(1) 14 Series: Residuals Sample 2012M03 2019M10 Observations 92 12 10 Mean Median Maximum Minimum Std Dev Skewness Kurtosis -3795.027 -1411.331 208429.5 -153957.2 53871.74 0.310684 5.024853 Jarque-Bera Probability 17.19682 0.000184 -150000 -100000 -50000 50000 100000 150000 200000 Nguồn: Eview Probability = 0.000184 < 5% -> bác bỏ H0 -> biến khơng có phân phối chuẩn 3.2.3 Kiểm định tự tương quan Cặp giả thuyết: H0: Mơ hình khơng có tự tương quan H1: Mơ hình có tự tương quan Từ cửa số ước lượng chọn View -> Residual Diagnostics -> Serial Correlation LM Test Lags to include: -> OK Prob Chi-Square(2) = 0.1161 > 5% -> khơng có sở bác bỏ H0 -> Mơ hình khơng có tự tương quan 3.2.4 Kiểm định phương sai sai số thay đổi Cặp giả thuyết: H0: Mơ hình khơng có phương sai sai số thay đổi H1: Mơ hình có phương sai sai số thay đổi Từ cửa số ước lượng chọn View -> Residual Diagnostics -> Heteroskedasticity Test -> ARCH -> OK 22 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Prob Chi-Square(1) = 0.9122 > 5% -> bác bỏ H0 -> Mơ hình có phương sai sai số thay đổi 3.3 Lựa chọn mơ hình tốt Lựa chọn mơ hình dựa sai số dự báo hai mơ hình Ta có bảng so sánh sai số dự báo mơ hình san mũ Winters mơ hình MA(1) sau: Bảng 3.6: Bảng so sánh sai số dự báo mơ hình san mũ Winters mơ hình MA(1) Mơ hình Winters MAPE 4.7% MRSE 3756.634 Mơ hình MA(1) 26.89% 347368.4 Nhóm xét thấy sai số mơ hình Winters nhỏ nên lựa chọn mơ hình Winters để dự báo chuỗi 23 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com CHƯƠNG 4: KẾT QUẢ DỰ BÁO Sau dự báo mẫu chọn mơ hình san mũ Winters, nhóm tiến hành dự báo mẫu cho chuỗi Đổi range mẫu: kích đúp chuột vào range: end date: 2020M12 Kích đúp chuột trái vào chuỗi visitor -> Proc -> Exponential Smoothing -> Simple Exponential Smoothing -> Smoothing Method: Holts-Winters – Multiplicative -> Smoothed series: visitorwm -> OK Ta bảng kết bảng 3.1 Các giá trị dự báo cho chuỗi visitor tính Excel theo cơng thức: visitorf=(Mean+h*Trend)*Seasonals Trong h số kỳ sau giá trị cuối chuỗi visitor (2019M10) Cuối nhóm thu giá trị dự báo cho chuỗi visitor từ tháng 11/2019 đến tháng 12/2020 sau: 24 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com Bảng 4.1: Bảng kết giá trị dự báo cho chuỗi visitor từ tháng 11/2019 đến tháng 12/2020 Đơn vị: lượt người h Month Mean Trend Seasonals visitorf 2019M11 1778629 10309.56 1.020312 1825275 2019M12 1.029115 1851633 2020M01 1.148441 2078170 2020M02 1.186161 2158656 2020M03 1.038712 1901027 2020M04 1.073662 1976060 2020M05 0.921972 1706382 2020M06 0.850721 1583282 2020M07 0.928087 1736836 10 2020M08 1.00987 1900297 11 2020M09 0.894196 1691849 12 2020M10 0.898751 1709733 13 2020M11 1.020312 1951503 14 2020M12 1.029115 1978950 25 LUAN VAN CHAT LUONG download : add luanvanchat@agmail.com KẾT LUẬN Chuỗi số liệu lượng du khách quốc tế đến Việt Nam khoảng thời gian từ tháng 1/2012 đến tháng 10/2019 chuỗi có tính mùa vụ, xu tăng, thân chuỗi không dừng sai phân bậc có dừng nên nhóm chọn phương pháp san mũ Winters phương pháp ARIMA để dự báo cho chuỗi Khi so sánh hai mô hình mơ hình san mũ Winters cho sai số dự báo thấp nên mơ hình phù hợp với chuỗi số liệu Cụ thể sai số MAPE mơ hình Winter 4,7% (