Tổng quan nghiên cứu

Dự báo và phân tích chuỗi thời gian đóng vai trò đặc biệt quan trọng trong các quyết định kinh tế lượng, tài chính và khoa học kỹ thuật. Trong thực tế, dữ liệu thực nghiệm thường biến động phức tạp dưới tác động đan xen của xu thế dài hạn, tính mùa vụ và các xung lượng ngẫu nhiên. Luận văn thạc sĩ khoa học chuyên ngành Lý thuyết xác suất và thống kê toán học (mã số 60460106) của tác giả Dương Nhật Thăng, dưới sự hướng dẫn khoa học của Phó giáo sư Tiến sĩ Phan Viết Thư tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, giải quyết trực tiếp bài toán xây dựng nền tảng toán học vững chắc để nhận dạng, ước lượng và dự báo các quá trình ngẫu nhiên.

Mục tiêu cụ thể của công trình là hệ thống hóa lý thuyết quá trình dừng, phát triển các mô hình tự hồi quy tích hợp trung bình trượt tuyến tính và mở rộng sang mô hình mùa vụ, đồng thời hiện thực hóa thuật toán xử lý dữ liệu trên phần mềm thống kê R. Phạm vi nghiên cứu bao quát từ các mô hình toán học giải tích đến phân tích thực nghiệm trên chuỗi dữ liệu thực tế, bao gồm 28 quan sát theo quý của Công ty Điện - Nước Washington giai đoạn 1980–1986 và 144 quan sát theo tháng về lượng hành khách hàng không quốc tế giai đoạn 1949–1960.

Ý nghĩa học thuật và ứng dụng của đề tài thể hiện qua việc tối ưu hóa sai số dự báo trung bình bình phương (MSE) và thiết lập khoảng tin cậy 95% có tính đến cấu trúc tương quan chuỗi. Nghiên cứu chứng minh rằng khi tính tự tương quan thời gian được kiểm soát chặt chẽ, độ lệch chuẩn ước lượng được hiệu chỉnh chuẩn xác, loại bỏ hoàn toàn các sai số dự báo dài hạn có thể gây tổn thất kinh tế nghiêm trọng.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Khung lý thuyết của luận văn được xây dựng trên nền tảng giải tích ngẫu nhiên và lý thuyết xác suất hiện đại, kế thừa hệ thống lý thuyết chuỗi thời gian trong tài chính của tác giả Ngai Hai Chan kết hợp với phương pháp luận kinh điển của Box và Jenkins. Công trình tập trung vào các khái niệm cốt lõi:

  1. Quá trình ngẫu nhiên và Định lý tính nhất quán Kolmogorov: Xác lập sự tồn tại của quá trình ngẫu nhiên thông qua họ các hàm phân bố hữu hạn chiều, tạo tiền đề toán học cho các phân tích tiệm cận.
  2. Tính dừng và cấu trúc hiệp phương sai: Phân biệt chặt chẽ giữa tính dừng mạnh (phân bố xác suất đồng thời bất biến qua phép tịnh tiến thời gian) và tính dừng yếu (kỳ vọng không đổi, hàm tự hiệp phương sai chỉ phụ thuộc vào độ trễ thời gian k). Quá trình nhiễu trắng đóng vai trò là khối cơ sở với kỳ vọng bằng 0 và phương sai không đổi.
  3. Họ mô hình tuyến tính tổng quát: Hệ thống hóa mô hình tự hồi quy AR(p), mô hình trung bình trượt MA(q), mô hình kết hợp ARMA(p,q), mô hình tích hợp ARIMA(p,d,q) và mô hình mùa vụ đa chiều SARIMA(p,d,q)x(P,D,Q)s thông qua toán tử dịch chuyển lùi B.
  4. Tính nhân quả và khả nghịch: Chứng minh điều kiện cần và đủ để một mô hình AR(p) có tính nhân quả là mọi nghiệm của đa thức đặc trưng phải nằm hoàn toàn bên ngoài đường tròn đơn vị phức.

Phương pháp nghiên cứu

Nghiên cứu sử dụng nguồn dữ liệu thứ cấp chuẩn hóa quốc tế gồm 2 tập dữ liệu điển hình: tập dữ liệu doanh thu năng lượng với kích thước mẫu N = 28 quan sát theo quý và tập dữ liệu vận tải hàng không với kích thước mẫu N = 144 quan sát liên tục trong 12 năm. Phương pháp chọn mẫu là phương pháp thu thập toàn bộ theo dòng thời gian liên tục, đảm bảo tính toàn vẹn cấu trúc tự tương quan.

Lý do lựa chọn phương pháp phân tích: Phân tích chuỗi thời gian miền thời gian được chọn vì cho phép bóc tách trực tiếp 4 thành phần cấu tạo gồm xu thế Tt, chu kỳ C, mùa vụ St và nhiễu ngẫu nhiên Nt. Quy trình phân tích kết hợp phương pháp bình phương cực tiểu (LSE), ước lượng hợp lý cực đại (MLE) và hệ phương trình Yule - Walker để ước lượng vector tham số. Kỹ thuật tiền xử lý sử dụng biến đổi Logarit để ổn định phương sai, bộ lọc trung bình trượt đối xứng 13 điểm, bộ lọc đa thức Spencer 15 điểm với vector trọng số chuẩn hóa chia cho hằng số 320, và phương pháp làm trơn hàm mũ đơn với các tham số điều chỉnh 0,2; 0,3; 0,6 nhằm tối thiểu hóa sai số MSE.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô hình hóa và kiểm định thực nghiệm trên phần mềm R đã mang lại 4 phát hiện quan trọng:

  1. Ổn định hóa phương sai qua biến đổi dữ liệu: Chuỗi hành khách hàng không với 144 điểm dữ liệu tăng trưởng từ 112 nghìn lượt (tháng 1/1949) lên 622 nghìn lượt (tháng 7/1960), tương đương mức tăng hơn 455%, bộc lộ phương sai mở rộng theo thời gian. Sau khi áp dụng phép biến đổi Logarit, biểu đồ hộp cho thấy phương sai được kiểm soát hoàn toàn ổn định qua 12 năm, loại bỏ tính không đồng nhất.
  2. Bóc tách chu kỳ mùa vụ năng lượng: Phân tích chuỗi 28 quý của Công ty Điện - Nước Washington chỉ ra quy luật tuần hoàn rõ rệt. Doanh thu luôn đạt đỉnh ở quý 1 (do nhu cầu sưởi ấm mùa đông bằng khí gas và nhiệt điện) và chạm đáy ở quý 3 với mức sụt giảm trung bình từ 30% đến 40% so với quý 1.
  3. Vai trò của nghiệm đơn vị và sai phân: Đối với các chuỗi không dừng có xu thế tuyến tính, sai phân bậc 1 (1 - B) triệt tiêu hoàn toàn xu thế mà không làm biến dạng cấu trúc tương quan vi mô. Việc sử dụng sai phân mùa bậc 12 đưa chuỗi vận tải hàng không về trạng thái dừng tiệm cận.
  4. Hiệu chỉnh khoảng tin cậy: Nghiên cứu chỉ ra rằng nếu xem các quan sát là độc lập cùng phân phối (khi hệ số tự tương quan bằng 0), khoảng tin cậy 95% có độ rộng xấp xỉ 2 chia cho căn bậc hai của n. Tuy nhiên, khi hệ số tương quan tiến tới 1, khoảng tin cậy thực tế mở rộng theo tỷ lệ căn bậc hai của (1 + phi) chia cho (1 - phi), chứng minh việc bỏ qua tương quan chuỗi sẽ làm sai lệch nghiêm trọng đánh giá rủi ro.

Thảo luận kết quả

Cơ chế đằng sau sự biến động của chuỗi thời gian thực tế bắt nguồn từ các yếu tố tự nhiên và kinh tế - xã hội. Đỉnh tiêu thụ năng lượng mùa đông tại vùng Tây Bắc nước Mỹ phản ánh điều kiện khí hậu đặc trưng, trong khi đỉnh vận tải hàng không vào tháng 7 và tháng 8 hàng năm (với lượng khách cao hơn khoảng 50% so với các tháng mùa đông) phản ánh chu kỳ nghỉ lễ mùa hè toàn cầu.

Kết quả dữ liệu được trình bày trực quan hóa thông qua hệ thống biểu đồ chuỗi thời gian đa lớp, biểu đồ hộp phân bố theo từng năm và đồ thị hàm tự tương quan mẫu (ACF). Ma trận dữ liệu 11 hàng 12 cột sau khi lọc xu thế cho phép tính toán chính xác chỉ số mùa vụ cho từng tháng. So với các mô hình hồi quy tuyến tính cổ điển vốn giả định sai số độc lập, mô hình SARIMA vượt trội hơn nhờ khả năng mô hình hóa sự phụ thuộc nội tại giữa tháng hiện tại với các tháng liền kề và cùng kỳ năm trước. So với mạng nơ-ron hồi quy phi tuyến, mô hình chuỗi thời gian tuyến tính Box - Jenkins trong luận văn đảm bảo tính tường minh tuyệt đối về mặt toán học và khả năng suy diễn thống kê chuẩn tắc.

Đề xuất và khuyến nghị

  1. Chuẩn hóa quy trình làm sạch và kiểm định tính dừng: Các cơ quan phân tích kinh tế vĩ mô và doanh nghiệp tài chính cần thực hiện bước chuyển đổi Logarit và kiểm định nghiệm đơn vị định kỳ 100% trước khi đưa dữ liệu vào dự báo, hoàn thành áp dụng trong vòng 6 tháng nhằm giảm thiểu sai số ngoại lai dưới mức 5%.
  2. Triển khai mô hình mùa vụ SARIMA vào kế hoạch sản xuất: Khuyến nghị các tổng công ty năng lượng, vận tải và chuỗi cung ứng áp dụng cấu trúc SARIMA đa chiều để lập kế hoạch phân bổ nguồn lực theo từng quý và tháng, đặt mục tiêu giảm thiểu chi phí tồn kho và điều độ vận hành ít nhất 15% trong vòng 4 quý tài chính.
  3. Tối ưu hóa tham số làm trơn hàm mũ: Bộ phận phân tích dữ liệu tại các doanh nghiệp thương mại nên thiết lập quy trình tự động dò tìm hệ số trơn alpha tối ưu trong khoảng từ 0,2 đến 0,6 dựa trên tiêu chí cực tiểu hóa sai số bình phương trung bình (MSE), triển khai đồng bộ trong thời gian 12 tháng.
  4. Xây dựng thư viện phân tích tự động trên nền tảng mã nguồn mở R: Đề xuất các viện nghiên cứu và trường đại học chuẩn hóa các tệp lệnh phân tích phần dư, ước lượng hàm tự tương quan riêng (PACF) và phân tích sai phân thành các gói công cụ nội bộ trong vòng 1 năm nhằm nâng cao năng lực nghiên cứu ứng dụng.

Đối tượng nên tham khảo luận văn

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Toán ứng dụng, Thống kê: Nắm vững phương pháp chứng minh giải tích về tính nhân quả, khả nghịch, tính dừng tiệm cận và kỹ thuật giải hệ phương trình sai phân tuyến tính bậc cao.
  2. Chuyên viên phân tích tài chính và quản trị rủi ro: Ứng dụng quy trình mô hình hóa chuỗi thời gian tài chính để phân tích lợi suất cổ phiếu, dự báo biến động giá tài sản và xác lập khoảng tin cậy danh mục đầu tư trong chu kỳ từ 5 đến 10 năm.
  3. Kỹ sư dữ liệu và chuyên viên kế hoạch vận hành doanh nghiệp: Sử dụng quy trình bóc tách xu thế - thời vụ bằng phần mềm R để dự báo nhu cầu phụ tải điện, doanh thu bán hàng và lượng hành khách theo 12 tháng trong năm.
  4. Giảng viên và nhà nghiên cứu kinh tế lượng: Sử dụng làm tài liệu tham khảo giảng dạy chuyên sâu về phương pháp luận Box - Jenkins, kỹ thuật lọc Spencer 15 điểm và ứng dụng thực hành R trong kinh tế số.

Câu hỏi thường gặp

Tính dừng của chuỗi thời gian là gì và tại sao bắt buộc phải kiểm tra? Tính dừng đảm bảo kỳ vọng và phương sai của chuỗi không đổi theo thời gian, đồng thời hàm tự hiệp phương sai chỉ phụ thuộc vào độ trễ. Nếu chuỗi không dừng, các ước lượng thống kê sẽ mất hiệu lực, dẫn đến hiện tượng hồi quy giả với độ tin cậy danh định 95% nhưng thực tế sai số dự báo bùng nổ vô hạn.

Khi nào nên áp dụng mô hình SARIMA thay vì ARMA thông thường? Mô hình SARIMA bắt buộc áp dụng khi chuỗi dữ liệu tồn tại đồng thời hai tầng tương quan: tương quan giữa các chu kỳ liền kề và tương quan mùa vụ cố định lặp lại sau mỗi chu kỳ s thời điểm (như 4 quý hoặc 12 tháng), giúp loại bỏ triệt để tính phi dừng theo mùa.

Phép biến đổi Logarit có tác dụng then chốt nào trong xử lý dữ liệu? Biến đổi Logarit giúp ổn định hóa phương sai khi biên độ dao động của chuỗi gia tăng tỷ lệ thuận với giá trị trung bình theo thời gian. Minh chứng trên chuỗi 144 tháng vận tải hàng không, phép lấy logarit đã chuyển mô hình mùa vụ dạng nhân phức tạp về dạng cộng tuyến tính thuần nhất.

Làm cách nào để xác định cấp số bậc p của AR và q của MA? Bậc p và q được xác định thông qua việc đối sánh đồ thị hàm tự tương quan (ACF) và hàm tự tương quan riêng (PACF). Quá trình AR(p) có PACF triệt tiêu sau độ trễ p, trong khi quá trình MA(q) có ACF triệt tiêu hoàn toàn sau độ trễ q, kết hợp chỉ số thông tin để chọn bậc tối ưu.

Vì sao cần loại bỏ một số điểm dữ liệu ban đầu khi mô phỏng quá trình AR(1)? Quá trình AR(1) xuất phát từ một giá trị ban đầu xác định chỉ đạt trạng thái dừng khi thời gian t tiến tới vô cùng (dừng tiệm cận). Việc loại bỏ khoảng 50 đến 100 quan sát khởi đầu giúp triệt tiêu hoàn toàn ảnh hưởng của điều kiện biên ban đầu, đảm bảo mẫu mô phỏng đại diện chính xác cho phân bố dừng.

Kết luận

  • Luận văn hệ thống hóa toàn diện cơ sở toán học xác suất của các quá trình ngẫu nhiên dừng, mô hình AR, MA, ARMA, ARIMA và mở rộng sang cấu trúc mùa vụ phức hợp SARIMA.
  • Thiết lập hoàn chỉnh các kỹ thuật tiền xử lý gồm biến đổi Logarit, sai phân đơn, sai phân mùa vụ và các bộ lọc làm trơn trung bình trượt đối xứng từ 3 điểm đến 15 điểm.
  • Khẳng định tính ưu việt của phương pháp ước lượng hợp lý cực đại và hệ phương trình Yule - Walker trong việc xác định chính xác vector tham số mô hình trên không gian mẫu thực nghiệm.
  • Minh chứng năng lực ứng dụng thực tiễn vượt trội thông qua việc giải quyết trọn vẹn 2 bài toán dự báo doanh thu năng lượng và vận tải hàng không trên môi trường lập trình R.
  • Nghiên cứu mở ra lộ trình mở rộng kéo dài từ 6 đến 12 tháng tiếp theo để phát triển các mô hình chuỗi thời gian tài chính phi tuyến tính và ứng dụng mạng nơ-ron hồi quy trong điều kiện thị trường biến động mạnh. Quý độc giả và các nhà nghiên cứu quan tâm hãy tiếp cận toàn văn công trình để khai thác trọn vẹn thuật toán và hệ thống mã nguồn R chuyên sâu.