4.2 .Kết quả thực nghiệm
4.2.1. Phân tích điểm nĩng dịch tả
4.2.1.1. Theo tháng
Dựa trên số liệu thống kê số ca bệnh tả theo tháng trong giai đoạn 2004 - 2010, cĩ thể thấy cao điểm của dịch tả xảy ra vào hai khoảng thời gian: tháng 3, 4, 5, 7 (tháng mưa ít) và tháng 10, 11 (tháng mưa nhiều). Tháng cĩ ít số ca bệnh nhất là tháng 1, 2, 8 và 9. Thế nhưng về mặt khơng gian, câu hỏi đặt ra là các ca bệnh thường xuất hiện ở khu vực nào? các ca bệnh phân bố tập trung thành cụm hay phân tán rải rác tồn vùng? Để trả lời cho hai câu hỏi này, nghiên cứu tiến hành phân tích điểm nĩng theo từng tháng với kết quả được thể hiện từ Hình 4.2 đến Hình 4.6. Theo đĩ, cĩ thể thấy các điểm nĩng về số ca bệnh tả thay đổi theo từng tháng, tuy nhiên thường tập trung quanh khu vực nội đơ bao gồm các quận Ba Đình, Hồn Kiếm, Hai Bà Trưng, Thanh Xuân, Đống Đa, Cầu Giấy. Đây là vùng tập trung dân cư đơng đúc, tiếp giáp với một số con sơng ơ nhiễm chảy qua địa bn.
Sóc Sơn
Mê Linh
Ba Vì
Son Tay ụng Anh
*#
Phc Thọ Đ an
Ba Vi Ph- ng
*#
Tp. Sơn Tây Tây Hoai Duc
*# Hồ Từ Liê m Long
Cu Ba ình Biê n
Hon
Hoi c GiÊyLang§ èng KiÕm Gia Lâm Thạ ch *# Đ aHai Bà
ThÊt Thanh Tr- ng Quốc Oai Xuân
Hoàng Ha Dong Mai *# Tp. Hà ụng Thanh Trì Ch- ơng Mü Thanh Oai Th- êng TÝn Hot Spot 2 Gi_Bin
Cold Spot - 99% Confidence Phĩ Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mü § øc Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Sóc Sơn
Mê Linh
Ba Vì
Son Tay Đ ơng Anh
*#
Phĩc Thä § an
Ba Vi Ph- ỵ ng
*#
Tp. Sơn Tây Tây Hoai Duc
*# Hồ Từ Liê m Long
CuBa ình Biê n
Hon
Hoi c GiấyLangĐ ống Kiếm Gia Lõm Th ch *#
Đ a Hai Bà ThÊt Thanh Tr- ng
Quèc Oai Xuâ n
Hồng Ha Dong Mai *# Tp. Hµ Đ ơng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 3 Gi_Bin
Cold Spot - 99% Confidence Phĩ Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mỹ Đ øc Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Mê Linh Mê Linh
Ba Vì Ba Vì
Son Tay Đ ông Anh Son Tay Đ ông Anh
*# *#
Phĩc Thä § an Phĩc Thä § an
Ba Vi Ph- ỵ ng Ba Vi Ph- ng
*# *#
Tp. Sơn Hoai Duc Tây Tp. Sn Tõy Tây Tây Hoai Duc
*# Hồ *# Hồ Từ Liê m Long Tõ Liª m Long
CầuBa Đ ình Biê n Cầu Ba Đ ình Biê n
Hồn Hồn
Hồi Đ ức GiÊyLang§ èng KiÕm Gia Lâm Hồi Đ ức GiÊyLang§ èng KiÕm Gia Lâm Thạ ch *# Đ a Hai Bà Th¹ ch *# Đ a Hai Bà
ThÊt Thanh Tr- ng ThÊt Thanh Tr- ng Quốc Oai Xuân Quốc Oai Xuõn
Hong Hoàng
Ha Dong Mai Ha Dong Mai
*# *# Tp. Hµ Tp. Hµ Đ ơng Đ ơng Thanh Thanh Trì Trì Ch- ơng Mỹ Ch- ơng Mỹ Thanh Thanh
Oai Th- ờng Oai Th- êng
TÝn TÝn
Hot Spot 4 Hot Spot 5
Gi_Bin Gi_Bin
Cold Spot - 99% Confidence Phĩ
Cold Spot - 99% Confidence Phĩ Cold Spot - 95% Confidence Xuyª n
Cold Spot - 95% Confidence Xuyê n
ứng Hòa ứng Hòa
Cold Spot - 90% Confidence Cold Spot - 90% Confidence Not Significant Mü § øc
Not Significant Mü § øc Hot Spot - 90% Confidence Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence Hot Spot - 99% Confidence
Hình 4.3. Phân tích điểm nĩng số ca bệnh tả tháng 4, 5
± ±
Sãc Sơn Mê Linh
Ba V×
Son Tay Đ ơng Anh
*#
Phĩc Thä Đ an
Ba Vi Ph- ng
*#
Tp. Sơn Hoai Duc Tây Tây
*# Hồ Từ Liê m Long
Cầu Ba ỡnh Biờ n
Hon
Hoi ức GiÊyLang§ èng KiÕm Gia Lâm Thạ ch *# § a Hai Bà Thất Thanh Tr- ng
Quốc Oai Xuâ n
Hoàng Ha Dong Mai *# Tp. H ụng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- ờng Tín Hot Spot 6 Gi_Bin Mờ Linh Ba Vỡ *#Son Tay Phc Thọ Đ an *#Ba Vi Ph- ng Tp. Sơn Tõy Hoi c Thạ ch ThÊt Qc Oai Ch- ¬ng Mü Hot Spot 7 Gi_Bin Sãc S¬n Đ ơng Anh
Hoai Duc Tây *# Hồ Từ Liê m Long
CầuBa ỡnh Biờ n
Hon
GiấyLangĐ ống Kiếm Gia Lâm
*# Đ a Hai Bà Thanh Tr- ng Xuâ n Ha Dong Hoàng Mai *# Tp. H ụng Thanh Trì Thanh Oai Th- ờng Tín
Cold Spot - 99% Confidence Ph
Cold Spot - 95% Confidence Xuyê n ứng Hòa
Cold Spot - 90% Confidence
Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Cold Spot - 99% Confidence
Phĩ
Cold Spot - 95% Confidence Xuyê n ứng Hòa
Cold Spot - 90% Confidence
Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
±
Mê Linh
Ba Vỡ
Son Tay Đ ông Anh
*#
Phc Thọ Đ an
*#Ba Vi Ph- ng
Tp. Sn Tõy Tây *#Hoai Duc Hồ
Từ Liê m Long
CuBa Đ ìnhHồn Biê n
Giấy
Hồi Đ øc LangĐ ống Kiếm Gia Lõm Thạ ch *# Đ a Hai Bà Thất Thanh Tr- ng
Quốc Oai Xuân
Hoàng Ha Dong Mai *# Tp. H ụng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- ờng Tín Hot Spot 8 Gi_Bin
Cold Spot - 99% Confidence Phĩ Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mỹ § øc Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Mê Linh
Ba Vì
Son Tay Đ ơng Anh
*#
Phĩc Thọ Đ an
Ba Vi Ph- ng
*#
Tp. Sn Tõy Tây Hoai Duc
*# Hồ Từ Liê m Long
CuBa ỡnhHon Biờ n
Giấy
Hoài Đ øc LangĐ ống Kiếm Gia Lâm Th¹ ch *# Đ a Hai Bà
Thất Thanh Tr- ng Quèc Oai Xuân
Hoàng Ha Dong Mai *# Tp. Hà Đ ụng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- ờng Tín Hot Spot 9 Gi_Bin
Cold Spot - 99% Confidence Ph Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mü § øc Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Hình 4.5. Phân tích điểm nĩng số ca bệnh tả tháng 8, 9
Sóc Sơn
Mê Linh
Ba Vì
Son Tay ụng Anh *#
Phc Thọ Đ an
*#Ba Vi Ph- ng
Tp. Sơn Hoai Duc Tây Tây
*# Hồ Từ Liê m Long
CầuBa ỡnh Hn Biờ n
Giy
Hoi ức Lang ng Kim Gia Lõm Thạ ch *# Đ aHai Bà Thất Thanh Tr- ng
Qc Oai Xn
Hồng Ha Dong Mai *# Tp. Hµ Đ ơng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 10 Gi_Bin
Cold Spot - 99% Confidence Phĩ
Cold Spot - 95% Confidence Xuyê n ứng Hòa
Cold Spot - 90% Confidence
Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
±
Sãc S¬n Mê Linh
Ba Vì
Son Tay Đ ông Anh
*#
Phc Thọ Đ an
Ba Vi Ph- ng
*#
Tp. Sơn Hoai Duc Tây Tây
*# Hồ Từ Liê m Long
Cu Ba ình Biê n
Hồn
Hồi Đ ức GiÊyLangĐ ống Kiếm Gia Lâm Th¹ ch *# Đ a Hai Bà Thất Thanh Tr- ng
Quốc Oai X n
Hồng Ha Dong Mai *# Tp. Hµ Đ ơng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 11 Gi_Bin
Cold Spot - 99% Confidence
Phĩ
Cold Spot - 95% Confidence Xuyê n ứng Hịa
Cold Spot - 90% Confidence
Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Hình 4.6. Phân tích điểm nĩng số ca bệnh tả tháng 10, 11
Theo số liệu thống kê số ca bệnh tả xuất hiện trong giai đoạn 2004 - 2010, cĩ thể thấy năm 2004 bắt đầu ghi nhận các ca mắc bệnh tả tại Tp. Hà Nội với số lượng chỉ 25 ca. Sau đĩ, dịch tả bùng phát liên tục từ năm 2007 đến 2010, đỉnh điểm là 2008. Thế nhưng về mặt khơng gian, tương tự như trường hợp theo tháng, câu hỏi đặt ra là các ca bệnh thường xuất hiện ở khu vực nào? các ca bệnh phân bố tập trung thành cụm hay phân tán rải rác tồn vùng? Để trả lời cho hai câu hỏi này, nghiên cứu tiến hành phân tích điểm nĩng theo từng năm với kết quả được thể hiện từ Hình 4.7 đến Hình 4.9. Theo đĩ, cĩ thể thấy các điểm nĩng về số ca bệnh tả thay đổi theo từng năm. Tuy nhiên, tương tự như tháng, các điểm nĩng thường tập trung quanh khu vực nội ơ bao gồm các quận Ba Đình, Hồn Kiếm, Hai Bà Trưng, Thanh Xuân, Đống Đa, Cầu Giấy. Đây là khu vực tập trung dân cư đơng đúc, tiếp giáp với sơng Hồng về phía Bắc.
Sãc S¬n
±
Mê Linh
Ba Vì
Son Tay Đ ông Anh
*#
Phĩc Thä § an
Ba Vi Ph- ỵ ng
*#
Tp. Sơn Tõy Tây Hoai Duc
*# Hồ Từ Liê m Long
CầuBa ỡnh Biờ n
Hon
Hoi c GiyLangĐ ống Kiếm Gia Lâm Thạ ch *# Đ aHai Bà
Thất Thanh Tr- ng Quốc Oai Xuõn
Hong Ha Dong Mai *# Tp. Hà Đ «ng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 2004 Gi_Bin
Cold Spot - 99% Confidence Ph Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mü § øc Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Sãc Sơn
Mờ Linh
Ba Vỡ
Son Tay ụng Anh
*#
Phc Thọ Đ an
*#Ba Vi Ph- ng
Tp. Sn Tõy Tây Hoai Duc
*# Hồ Từ Liê m Long
CầuBa ỡnh Biờ n
Hon
Hoi c GiyLangĐ ống Kiếm Gia Lâm Thạ ch *# Đ aHai Bà
Thất T hanh Tr- ng Qc Oai Xn
Hồng Ha Dong Mai *# Tp. Hµ Đ ơng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 2007 Gi_Bin
Cold Spot - 99% Confidence Phĩ Cold Spot - 95% Confidence Xuyê n
ứng Hòa Cold Spot - 90% Confidence
Not Significant Mü § øc Hot Spot - 90% Confidence
Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Mê Linh Mê Linh
Ba Vì Ba V×
Son Tay Đ ơng Anh Son Tay Đ ơng Anh
*# *#
Phĩc Thä § an Phĩc Thä § an
Ba Vi Ph- ỵ ng Ba Vi Ph- ỵ ng
*# *#
Tp. S¬n Hoai Duc Tây Tp. Sơn Hoai Duc T©y
T©y T©y
*# Hå *# Hồ
Từ Liê m Long Từ Liê m Long
CầuBa ình Biê n CầuBa Đ ình Biê n
Hoµn Hoµn
Hồi Đ ức GiÊyLangĐ ống Kiếm Gia Lâm Hoài Đ ức GiÊyLangĐ ống Kiếm Gia Lâm Th¹ ch *# Đ aHai Bà Th¹ ch *# Đ aHai Bà
Thất Thanh Tr- ng Thất Thanh Tr- ng Quốc Oai Xuâ n Quốc Oai Xuõ n
Hong Hoàng Ha Dong Ha Dong Mai Mai *# *# Tp. Hµ Tp. Hµ Đ ơng Đ ụng Thanh Thanh Trì Trì Ch- ơng Mỹ Ch- ơng Mỹ Thanh Thanh
Oai Th- êng Oai Th- êng
TÝn TÝn
Hot Spot 2008 Hot Spot 2009
Gi_Bin Gi_Bin
Cold Spot - 99% Confidence Phĩ Cold Spot - 99% Confidence Phĩ
Cold Spot - 95% Confidence Xuyª n Cold Spot - 95% Confidence Xuyê n
ứng Hòa ứng Hòa
Cold Spot - 90% Confidence Cold Spot - 90% Confidence
Not Significant Mü § øc Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence Hot Spot - 99% Confidence
Hình 4.8. Phân tích điểm nĩng số ca bnh t nm 2008, 2009
Sóc Sơn
Mê Linh
Ba Vì
Son Tay ụng Anh *#
Phc Thọ Đ an
Ba Vi Ph- ng *#
Tp. Sơn Hoai Duc Tây Tây
*# Hồ Từ Liê m Long
CầuBa ỡnhHn Biê n
Giy
Hoi c LangĐ ốngKiếm Gia Lâm Thạ ch *# aHai Bà
ThÊt Thanh Tr- ng Qc Oai Xn
Hồng Ha Dong Mai *# Tp. Hµ Đ ơng Thanh Trì Ch- ơng Mỹ Thanh Oai Th- êng TÝn Hot Spot 2010 Gi_Bin
Cold Spot - 99% Confidence Phĩ
Cold Spot - 95% Confidence Xuyê n ứng Hòa
Cold Spot - 90% Confidence
Not Significant Mü § øc
Hot Spot - 90% Confidence Hot Spot - 95% Confidence Hot Spot - 99% Confidence
Hình 4.9. Phân tích điểm nĩng số ca bệnh tả năm 2010
4.2.2.Xây dựng mơ hình hồi qui đa biến dự báo dịch tả trên địa bàn Tp. Hà Nội 4.2.2.1 Lựa chọn biến giải thích phát sinh dịch tả
Các kết quả phân tích điểm nĩng về ca bệnh tả theo tháng, năm, đều cho thấy các điểm nĩng thường tập trung tại những khu vực dân cư đơng đúc và nằm gần các con sơng. Từ nhận định trên kết hợp với các nghiên cứu đi trước về phân tích bệnh tả, nghiên cứu lựa chọn các biến giải thích phát sinh dịch tả , trên địa bàn Tp. Hà Nội như sau:
Theo tháng: các biến R, Sh, T, U, V lấy trung bình tháng; diện tích mặt nước (km2). Theo năm: các biến R, Sh lấy tổng theo năm; các biến T, U, V lấy trung bình năm; diện tích mặt nước (km2), dân số (nghìn người).
Do số ca mắc bệnh tả phân bố rất khơng đều theo tháng và theo năm nên nghiên cứu lựa chọn hàm hồi qui logarit để giải thích số ca bệnh tả (y) với dạng như sau:
Theo tháng: Logarit (y + 1) = α + β1 * R + β2 * Sh + β3 * T + β4 * U + β5* V + β6 * diện tích mặt nước + ε (sai số ngẫu nhiên)
Theo năm:Logarit (y + 1) = α + β1 * R + β2 * Sh + β3 * T + β4 * U + β5 * V + β6
* diện tích mặt nước + β7 * dân số + ε (sai số ngẫu nhiên) Trong đĩ: α là hệ số chặn, βi là hệ số hồi qui.
4.2.2.2. Mơ hình OLS
a, Theo tháng: Kết quả phân tích hồi qui OLS cho các tháng 3, 4, 5, 6, 7, 10, 11,
12 được thể hiện trong Bảng 4.2. Độ lệch chuẩn của phần dư (số ca bệnh thực tế - số ca bệnh mơ phỏng) được lần lượt tính tốn cho các tháng. Đối với các tháng 2, 8, 9, do số ca bệnh rất ít nên nghiên cứu khơng tìm ra được hàm hồi qui cĩ ý nghĩa thống kê để giải thích sự xuất hiện ca bệnh. Tháng 1 khơng cĩ ca bệnh nào nên khơng thiết lập hàm hồi quy. Bảng tổng hợp kết quả phân tích hồi qui OLS cho tháng được thể hiện ở bảng 4.2 và từ Hình 4.10 đến Hình 4.13. Từ các kết quả phân tích hồi qui OLS, cĩ thể rút ra một số nhận xét như sau: (1) Trong các tháng 3, 11, 12, yếu tố khí hậu giải thích được 55%, 38%, 32% số ca bệnh trên tồn khu vực; (2) Trong tháng 4 cao điểm về dịch bệnh, yếu tố mặt nước giải thích được 25% số ca bệnh trên tồn khu vực; và (3) Trong các tháng cịn lại bao gồm các tháng 5, 6, 7, 10, sự kết hợp của yếu
tố khí hậu và mặt nước giải thích được 72%, 41%, 57%, 55% số ca bệnh trên tồn khu vực.
Bảng 4.2. Tổng hợp kết quả phân tích hồi qui OLS theo tháng khu vực Hà Nội
Tháng Biến giải thích R2 *p_value
3 Hằng số, T,U,V 0.548761 p< 0,01 4 Hằng số, Mặt nước 0.250669 p< 0,01 5 Hằng số, Mặt nước, V 0.719093 p< 0,01 6 Hằng số, Mặt nước, R 0.414949 p< 0,01 7 Hằng số, Mặt nước, R, Sh,V 0.569390 p< 0,01 10 Hằng số, Mặt nước, Sh,T,V 0.549334 p< 0,01 11 Hằng số, R, Sh 0.380233 p< 0,01 12 Hằng số, Sh 0.324019 p< 0,01
*p_value: giá trị thống kê