Mô hình dự báo số bệnh nhân lao bằng phương pháp arima tại thành phố hồ chí minh đến năm 2020

Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ Y TẾ ĐẠI HỌC Y DƯỢC THÀNH PHỐ HỒ CHÍ MINH - NGUYỄN THIỆN MINH MƠ HÌNH DỰ BÁO SỐ BỆNH NHÂN LAO BẰNG PHƯƠNG PHÁP ARIMA TẠI THÀNH PHỐ HỒ CHÍ MINH ĐẾN NĂM 2020 LUẬN VĂN THẠC SĨ Y TẾ CƠNG CỘNG TP Hồ Chí Minh - Năm 2018 Thơng tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� BỘ GIÁO DỤC VÀ ĐÀO TẠO BỘ Y TẾ ĐẠI HỌC Y DƯỢC TP HỒ CHÍ MINH - NGUYỄN THIỆN MINH MƠ HÌNH DỰ BÁO SỐ BỆNH NHÂN LAO BẰNG PHƯƠNG PHÁP ARIMA TẠI THÀNH PHỐ HỒ CHÍ MINH ĐẾN NĂM 2020 Ngành: Y Tế Công Cộng Mã số: 8720701 LUẬN VĂN THẠC SĨ Y TẾ CÔNG CỘNG Người hướng dẫn khoa học: PGS.TS ĐỖ VĂN DŨNG TP Hồ Chí Minh - Năm 2018 Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� i LỜI CAM ĐOAN Tôi xin cam đoan số liệu luận văn ghi nhận, nhập liệu phân tích cách trung thực Luận văn khơng có số liệu, văn bản, tài liệu Đại học Y Dược TP Hồ Chí Minh hay trường đại học khác chấp nhận để cấp văn đại học, sau đại học Luận văn khơng có số liệu, văn bản, tài liệu công bố trừ công khai thừa nhận Đề cương nghiên cứu chấp thuận mặt y đức nghiên cứu từ hội đồng duyệt đề cương Khoa Y tế công cộng ngày 28/06/2017 Ngày … tháng … năm 2018 Học viên NGUYỄN THIỆN MINH Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� ii MỤC LỤC MỤC LỤC ii DANH MỤC CÁC TỪ VIẾT TẮT iv DANH MỤC ĐỐI CHIẾU TỪ CHUYÊN MÔN TIẾNG VIỆT – TIẾNG NƯỚC NGOÀI vi DANH MỤC CÁC HÌNH, BIỂU ĐỒ vii DANH MỤC CÁC BẢNG viii ĐẶT VẤN ĐỀ Câu hỏi nghiên cứu: Mục tiêu tổng quát: Mục tiêu cụ thể: DÀN Ý NGHIÊN CỨU CHƯƠNG 1: TỔNG QUAN Y VĂN 1.1 Bệnh lao yếu tố nguy 1.2 Phương pháp phân tích chuỗi thời gian nghiên cứu y học 16 1.3 Ứng dụng phân tích chuỗi thời gian nghiên cứu lao phổi 26 CHƯƠNG 2: PHƯƠNG PHÁP NGHIÊN CỨU 29 2.1 Thiết kế nghiên cứu: 29 2.2 Dân số nghiên cứu: 29 2.3 Kỹ thuật chọn mẫu: 29 2.4 Tiêu chí chọn mẫu: 29 2.5 Thu thập kiện: 30 2.6 Liệt kê định nghĩa biến số: 32 2.7 Xử lý phân tích kiện: 35 2.8 Vấn đề y đức: 36 CHƯƠNG 3: KẾT QUẢ 37 3.1 Đặc tính mẫu nghiên cứu 37 3.2 Mối liên quan nhiệt độ lượng mưa với số bệnh nhân lao mắc TP.HCM 42 3.3 Chuyển đổi liệu số bệnh nhân lao 44 Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� iii 3.4 Các mơ hình dự báo 47 3.5 Đánh giá tính xác dự báo 53 3.6 Chọn mơ hình dự báo 56 3.7 Dự báo số bệnh nhân lao/100.000 dân từ năm 2018 đến năm 2020 59 CHƯƠNG 4: BÀN LUẬN 60 4.1 Bệnh lao thành phố Hồ Chí Minh 60 4.2 Mối liên quan nhiêt độ, lượng mưa trung bình số bệnh nhân lao mắc TP.HCM 62 4.3 Điểm gãy cấu trúc dự báo 72 4.4 Mơ hình dự báo 78 4.5 Kết dự báo 82 4.6 Những điểm mạnh hạn chế nghiên cứu 87 4.7 Những điểm nghiên cứu 87 KẾT LUẬN 89 ĐỀ XUẤT - KIẾN NGHỊ 90 TÀI LIỆU THAM KHẢO i PHỤ LỤC xiv Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� iv DANH MỤC CÁC TỪ VIẾT TẮT AIDS : Acquired immune deficiency syndrome AC (Hội chứng suy giảm miễn dịch mắc phải) : Autocorrelation (tự tương quan) ACF ANN : Autocorrelation function (hàm tự tương quan) : Artificial neural networks AR ARIMA : Autoregressive : Autoregressive integrated moving average (Trung bình động, đồng liên kết, tự hồi quy) AFB : Acid - Fast – Bacilli BCG BIOSIS EMBASE HIV : : : : IGRAs MA MAPE (virus gây hội chứng suy giảm miễn dịch người) : Interferon-Gamma Release Assays : Moving average : Mean absolute percentage error (Trung bình phần trăm sai số tuyệt đối) Bacille Calmette-Guerin Bioscience Information Service Excerpta Medica dataBASE Human immunodeficiency virus MEDLINE : Medical Literature Analysis and Retrieval System Online OLS : Ordinary least squares PAC : Partial autocorrelation (tự tương quan phần/ tự tương quan riêng phần) PACF : Partial autocorrelation function (Hàm tự tương quan phần/ tự tương quan riêng phần) QALY : Quality adjusted life year RMSE : Root-mean-square error (Căn bậc hai bình phương sai số) RSS SIGLE : Residual sum of squares : System for Information on Grey Literature in Europe Social Policy & Practice : Tuberculin skin test : Thành phố Hồ Chí Minh : World Health Organization TST TP.HCM WHO Thơng tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� v (Tổ chức y tế giới/ tổ chức sức khỏe giới) ZN Thông tin kết nghiên cứu : Ziehl-Neelsen Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� vi DANH MỤC ĐỐI CHIẾU TỪ CHUN MƠN TIẾNG VIỆT – TIẾNG NƯỚC NGỒI Biến nội sinh : Endogenous variable Biến ngoại sinh : Exogenous variable Biến ngoại sinh túy : Pure exogenous variables Biến giải thích : Explanatory variable Biến phản ứng : Response variable Chân trời dự báo : Forecast horizon Đồ thị tương quan : Correlogram Đồng kết hợp : Cointegration Hiệp phương sai : Covariance Hàm tự tương quan : Autocorrelation Hàm tự tương quan phần/ hàm tự tương quan riêng phần : Partial Autocorrelation Hồi quy giả mạo : Spurious regression Kinh tế lượng : Econometric Nghiệm đơn vị : Unit-root Phân tích chuỗi thời gian : Time series analysis Tính dừng/ tình ổn định : Stationary Ước lượng tuyến tính khơng lệch tốt Thơng tin kết nghiên cứu : Best unbiased linear estimator Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� vii DANH MỤC CÁC HÌNH, BIỂU ĐỒ Hình 1 Các bước thực ước lượng ARIMA theo phương pháp luận Box-Jenkins 21 Hình Các quốc gia có gánh nặng bệnh lao lao nhạy cảm, lao kháng đa thuốc, lao đồng nhiễm HIV phân loại WHO giai đoạn 2016-2020 60 Hình Số bệnh nhân lao phân loại theo tuổi giới tính Việt Nam, năm 2016 78 Biểu đồ Giới tính bệnh lao TP.HCM từ năm 2006 đến năm 2017 (n=144) 41 Biểu đồ Số bệnh nhân lao thể/ 100.000 dân số TP.HCM từ 2006 đến 2017 trước sau hiệu chỉnh theo ngày làm việc ước tính (n=144) 44 Biểu đồ 3 Số bệnh nhân lao giới nam/ 100.000 dân số nam TP.HCM từ 2006 đến 2017 trước sau hiệu chỉnh theo số ngày làm việc ước tính (n=144) 45 Biểu đồ Số bệnh nhân lao giới nữ/ 100.000 dân số nữ TP.HCM từ 2006 đến 2017 trước sau hiệu chỉnh theo số ngày làm việc (n=144) 46 Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� viii DANH MỤC CÁC BẢNG Bảng 1 Quy định ghi kết xét nghiệm AFB nhuộm ZN 15 Bảng Các dạng mơ hình lý thuyết đồ ACF đồ thị PACF 22 Bảng Nhiệt độ khơng khí trung bình lượng mưa trung bình theo tháng TP.HCM từ năm 2006 đến năm 2016 (n=132) 37 Bảng Số ngày làm việc ước tính tháng TP.HCM từ năm 2006 đến năm 2016 (n=132) 38 Bảng 3 Dân số chung, dân số nam, dân số nữ theo năm TP.HCM từ năm 2006 đến năm 2016 (n=132) 39 Bảng Số bệnh nhân lao thể theo giới tính theo năm TP.HCM từ năm 2006 đến năm 2017 (n=132) 40 Bảng Số bệnh nhân lao /100.000 dân theo giới tính TP.HCM từ năm 2006 đến năm 2016 (n=132) 41 Bảng Kết ước lượng hồi quy tuyến tính đơn biến nhiệt độ mơi trường trung bình theo tháng với số bệnh nhân lao mắc số bệnh nhân lao /100.000 dân TP.HCM từ năm 2006 đến năm 2016 (n=132) 42 Bảng Kết ước lượng hồi quy tuyến tính đơn biến lượng mưa trung bình theo tháng với số bệnh nhân lao mắc số bệnh nhân lao /100.000 dân TP.HCM từ năm 2006 đến năm 2016 (n=132) 43 Bảng Điểm gãy cấu trúc bệnh nhân lao/100.000 dân hiệu chỉnh theo ngày làm việc từ năm 2006 đến năm 2017 (n=144) 47 Bảng Tham số mơ hình hồi quy bệnh nhân lao/100.000 dân hiệu chỉnh theo ngày làm việc từ năm 2006 đến năm 2017 (n=132) 47 Bảng 10 Hệ số mơ hình Holt Winters khơng có yếu tố mùa 48 Bảng 11 Kiểm định Augmented Dickey-Fuller (ADF) với chuỗi thời gian bệnh nhân lao /100.000 dân TP.HCM 49 Bảng 12 Một số mô hình ARIMA bệnh nhân lao/100.000 dân TP.HCM có AIC thấp (n=120) 50 Thông tin kết nghiên cứu Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� reg d.lao time *tạo mơ hình hồi quy predict et, resid *trích phần dư mơ hình, vào biến “et” twoway line et time *Vẽ biểu đồ phần dư mơ hình theo thời gian **_Kiểm định tính dừng phần dư dfuller et, drift regress lags(0) * kiểm định ADF với độ trễ =0 corrgram et, lags(10) *tính AC, PAC, Q-test, P-value Q-test với 10 độ trễ ac et *vẽ biểu đồ ACF phần dư pac et * vẽ biểu đồ PACF phần dư *Phần 6: Tạo dự báo forecast create mơ_hình *tạo mơ hình tên “mơ_hình” regress lao thudam thuduong *thực hồi quy estimates store ước_lượng_1 *lưu thông tin ước lượng vào “ước_lượng_1” tsappend, add(10) *Mở rộng thời thời khoảng thêm 10 giai đoạn/kỳ forecast estimates ước_lượng_1, names(biến_mới) *Ước tính dự báo từ “ước_lượng_1” forecast solve, prefix(F_) begin(tq(2017m1)) periods(7) double static *Giải mơ hình để tình dự báo với biến có phần đầu tên biến thêm “F_”; tháng năm 2017, dự báo với thời đoạn twoway (tsline lao) (tsline F_lao) *vẽ biểu đồ cho biến gốc biến dự báo hiểu đồ Phụ lục Marco tạo câu lệnh cho Microsoft Excel 2013 dùng để chuyển đổi ngày âm lịch sang dương lịch ngược lại Public Lday, Lmonth As Byte, Lyear As Integer, isLeap, LunarInfo Sub LNI() LunarInfo = Array( _ &H3C4BD8, &H624AE0, &H4CA570, &H3854D5, &H5CD260, &H44D950, &H315554, &H5656A0, &H409AD0, &H2A55D2, &H504AE0, &H3AA5B6, &H60A4D0, &H48D250, &H33D255, &H58B540, &H42D6A0, &H2CADA2, &H5295B0, &H3F4977, _ &H644970, &H4CA4B0, &H36B4B5, &H5C6A50, &H466D40, &H2FAB54, &H562B60, &H409570, &H2C52F2, &H504970, &H3A6566, &H5ED4A0, &H48EA50, &H336A95, &H585AD0, &H442B60, &H2F86E3, &H5292E0, &H3DC8D7, &H62C950, _ &H4CD4A0, &H35D8A6, &H5AB550, &H4656A0, &H31A5B4, &H5625D0, &H4092D0, &H2AD2B2, &H50A950, &H38B557, &H5E6CA0, &H48B550, &H355355, &H584DA0, &H42A5B0, &H2F4573, &H5452B0, &H3CA9A8, &H60E950, &H4C6AA0, _ &H36AEA6, &H5AAB50, &H464B60, &H30AAE4, &H56A570, &H405260, &H28F263, &H4ED940, &H38DB47, &H5CD6A0, &H4896D0, &H344DD5, &H5A4AD0, &H42A4D0, &H2CD4B4, &H52B250, &H3CD558, &H60B540, &H4AB5A0, &H3755A6, _ &H5C95B0, &H4649B0, &H30A974, &H56A4B0, &H40AA50, &H29AA52, &H4E6D20, Thông tin kết nghiên cứu Tuân thủ Luật Sở hữu trí tuệ Quy định truy cập tài liệu điện tử Ghi rõ nguồn tài liệu trích dẫn Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� &H39AD47, &H5EAB60, &H489370, &H344AF5, &H50EA50, &H3D6A58, &H6256A0, &H4AAAD0, &H46C960, &H2ED954, &H54D4A0, &H3EDA50, &H5EA5D0, &H4A92B0, &H32AAB5, &H58A950, &H3C55D9, &H624BA0, &H4CA5B0, &H375176, &H307934, &H546AA0, &H3EAD50, &H2A5B52, &H48D260, &H32EA65, &H56D520, &H40DAA0, &H6249D0, &H4CA4D0, &H37D0B6, &H5AB250, &H54B5A0, &H3E55D0, &H2A55B2, &H5049B0, &H33B255, &H586D20, &H40AD60, &H2D4B63, &H4C54B0, &H3768A6, &H5ADA50, &H445AA0, &H4052E0, &H28D2E3, &H4EC950, &H38D557, &H5856A0, &H42A6D0, &H2C55D4, &H5252B0, &H34B6A6, &H5AAD50, &H4655A0, &H2EAB64, &H2AB173, &H4E6930, &H386B37, &H5E6AA0, &H42A570, &H2E52E4, &H50D160, &H3AE958, &H5A56D0, &H462AE0, &H30A9D4, &H54A2D0, &H4EB520, &H38D727, &H5EADA0, &H4A55B0, &H2EB2B4, &H54A950, &H3CB559, &H626B20, &H484570, &H326574, &H5852B0, &H406950, &H3BAAA7, &H5EA6D0, &H4A4AE0, &H35A2E5, &H52D550, &H3E5ABB, &H6256A0, &H4C96D0, &H30D4B5, &H56B290, &H40B550, &H2A6D52, &H609570, &H4A49B0, &H34A975, &H5A64B0, &H3E2B60, &H28AB61, &H4C9570, &H384AE6, &H54DA90, &H405B50, &H2C36D3, &H502AE0, &H4AC950, &H32D556, &H58B4A0, &H42B690, &H6425B0, &H4E92B0, &H36AAB6, &H5C6950, &H4055A0, &H2AAB73, &H522570, &H3A5377, &H346D56, &H585AA0, &H42AB50, &H2E56D4, &H4CD260, &H36EAA6, &H5AD550, &H465AA0, &H2AA9B3, &H50A4D0, &H3AD2B7, &H5EB250, 2199 */ End Sub Sub lunar(d, m, y) Dim DiffADate, Counter, I, Temp &H5A4970, &H3696D5, &H2A7552, &H42B4A0, &H5C5270, &H504B60, &H2D56A3, &H44B520, &H3AA577, &H525370, &H2FA6A4, &H5ED4A0, &H3CA9B8, &H54A570, &H48AD50, &H60D520, &H3ED150, &H362DB5, &H4CAD50, &H2A7953, &H5AA550, &H3949B6, &H504DA0, &H446A90, &H5CD160, &H3A93D7, &H2E5D94, &H4674A0, &H6052B0, &H544AE0, &H30ADA5, &H48B540, &H4464B0, &H2C74A3, &H5C92E0, _ &H4E56A0, &H38A7A7, &H2CBAA4, &H50AD50, &H466930, _ &H38A6E6, &H5EA4E0, &H5256D0, &H3C4AFB, &H2EDD25, _ &H5EA4B0, &H48AA50, &H3E49E8, &H60C970, &H54AAD0, _ &H46D950, &H325D55, &H62A930, &H4AB490, &H4052B0, _ &H332AD5, &H582B60, &H4ADA90, &H355AA6, &H28E952, _ &H5A45B0, &H44A2B0, &H385766, &H5C5370, &H505AA0, _ &H42D2A0, &H2DE2A4, &H5E4AB0, &H46A8D0, &H3B9567, _ &H2CBA94, &H526B50, &H46E4A0, &H2EED25, &H6092D0, _ &H5255B0, &H3E25FA, &H31B2A5, &H54AD50, &H4A6950, _ &H3CA570, &H2864D2, &H5695D0, &H404AD0, &H33D556)'' /* Years 2100- DiffADate = DateDiff("d", #1/31/1900#, CDate(d & "-" & m & "-" & y)) Counter = -1 Lyear = 1900 For I = Lyear To 2199 Temp = YearDays(I) Counter = Counter + Temp If Counter >= DiffADate Then Counter = Counter - Temp Exit For End If Lyear = Lyear + Thông tin kết nghiên cứu Tuân thủ Luật Sở hữu trí tuệ Quy định truy cập tài liệu điện tử Ghi rõ nguồn tài liệu trích dẫn Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� Next Leap = LeapMonth(Lyear) isLeap = "" Lmonth = For I = To 12 If Leap > And I = Leap + And isLeap = "" Then isLeap = "(N)" Lmonth = Lmonth - I = I - Temp = LeapDay(Lyear) Else Temp = MonthDays(Lyear, I) End If If isLeap = "(N)" And I Leap Then isLeap = "" Counter = Counter + Temp If Counter >= DiffADate Then Counter = Counter - Temp Exit For End If Lmonth = Lmonth + Next Lday = DiffADate - Counter End Sub Function LeapMonth(y) LNI If y >= 1900 Then LeapMonth = LunarInfo(y - 1900) And &HF Else LeapMonth = End Function ' Function LeapDay(y) LNI If LunarInfo(y - 1900) And &HF Then If LunarInfo(y - 1900) And &H10000 Then LeapDay = 30 Else LeapDay = 29 Else LeapDay = End If End Function ' Function MonthDays(y, m) LNI MonthMask = Array(32768, 16384, 8192, 4096, 2048, 1024, 512, 256, 128, 64, 32, 16) If LunarInfo(y - 1900) And MonthMask(m - 1) Then MonthDays = 30 Else MonthDays = 29 End Function ' Thông tin kết nghiên cứu Tuân thủ Luật Sở hữu trí tuệ Quy định truy cập tài liệu điện tử Ghi rõ nguồn tài liệu trích dẫn Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� Function YearDays(y) Dim I LNI MonthMask = Array(32768, 16384, 8192, 4096, 2048, 1024, 512, 256, 128, 64, 32, 16) YearDays = 348 For I = To 11 If LunarInfo(y - 1900) And MonthMask(I) Then YearDays = YearDays + Next YearDays = YearDays + LeapDay(y) End Function 'Ham so chuyen doi tu Duong lich sang Am lich dang ngay,thang,nam Public Function TransLu(d, m, y) Call lunar(m, d, y) TransLu = Lday & "/" & Lmonth & isLeap & "/" & Lyear End Function 'Ham so chuyen doi tu Duong lich sang Am lich dang tu mot o Public Function TransLu1(NT As Date) Call lunar(Day(NT), Month(NT), Year(NT)) TransLu1 = Lday & "-" & Lmonth & isLeap & "-" & CanchiV(Lyear - 0) End Function Public Function TransSolar(d, m, y) As Date'Ngay thang nam am lich sang duong lich Dim iSd As Date iSd = DateSerial(y, m, d) - 70 Do iSd = iSd + Loop Until TransLu(Day(iSd), Month(iSd), Year(iSd)) = d & "/" & m & "/" & y TransSolar = iSd End Function Thông tin kết nghiên cứu Tuân thủ Luật Sở hữu trí tuệ Quy định truy cập tài liệu điện tử Ghi rõ nguồn tài liệu trích dẫn Bản quyền tài liệu thuộc Thư viện Đại học Y Dược Tp.HCM.� Phụ lục: Kết AIC, BIC, RMSE mơ hình ARIMA mơ hình kết hợp Kết đánh giá mơ hình AIC, BIC, RMSE số bệnh nhân lao thể STT Mơ hình AIC BIC RMSE 10 11 12 13 14 15 16 17 18 19 20 21 SARIMA(0.1.0)(0.1.0.12) SARIMA(0.1.0)(0.1.1.12) SARIMA(0.1.1)(0.1.0.12) SARIMA(0.1.1)(0.1.1.12) SARIMA(0.1.2)(0.1.0.12) SARIMA(0.1.2)(0.1.1.12) SARIMA(1.1.0)(0.1.0.12) SARIMA(1.1.0)(0.1.1.12) SARIMA(1.1.1)(0.1.0.12) SARIMA(1.1.1)(0.1.1.12) SARIMA(1.1.2)(0.1.0.12) SARIMA(1.1.2)(0.1.1.12) SARIMA(2.1.0)(0.1.0.12) SARIMA(2.1.0)(0.1.1.12) SARIMA(2.1.1)(0.1.0.12) SARIMA(2.1.1)(0.1.1.12) SARIMA(2.1.2)(0.1.0.12) SARIMA(2.1.2)(0.1.1.12) SARIMA(3.1.0)(0.1.0.12) SARIMA(3.1.0)(0.1.1.12) SARIMA(3.1.1)(0.1.0.12) 493,00 429,40 398,50 341,50 371,90 315,70 438,60 376,20 383,40 323,10 373,80 316,30 409,70 345,50 379,50 318,80 375,70 320,60 394,10 335,00 378,80 498,10 434,50 403,60 349,20 382,10 325,90 446,20 386,40 393,60 333,30 386,60 326,50 420,00 355,80 392,30 331,50 391,00 338,50 406,90 350,30 394,20 169,89 173,92 72,21 74,10 114,22 62,54 166,45 90,62 89,60 61,16 112,32 62,34 138,93 76,18 111,71 62,66 117,59 62,99 138,13 71,54 114,37 Thông tin kết nghiên cứu Tuân thủ Luật Sở hữu trí tuệ Quy định truy cập tài liệu điện tử Ghi rõ nguồn tài liệu trích dẫn RMSE (SD) 116,95 84,29 65,34 59,91 97,62 51,05 126,06 69,41 71,72 52,22 94,10 51,25 114,58 56,02 95,23 50,53 101,83 51,34 115,08 54,88 99,03 p-value (ADF test) 0,040 0,007 0,016 0,397