Khi dịch Covid tăng nhanh hơn hàm mũ ... và dự báo
Update 25/7: Dự phóng đỉnh dịch TP HCM
Tóm tắt: dữ liệu về số ca nhiễm Corona ở TP HCM, từ ngày 25/06 tới 25/07 được vẽ lại là đưa ra nhận xét về xu hướng tiến triển của dịch bệnh. Xu hướng hiện tại vẫn và tăng theo hàm mũ (cấp số nhân). Tuy nhiên tín hiệu tốt là hệ số mũ có xu hướng giảm, mặc dù khá chậm, thể hiện cho việc dịch bệnh đang được ngăn chặn.
Chú thích: vẫn có khả năng là dịch bệnh đang lây lan nhanh hơn hàm mũ, như đề cập ở cuối bài viết.
UPDATE 1- ngày 23/7/2021 - Đã cập nhật kết quả thực tế để kiểm tra lại kết quả dự báo. (Hình 7) UPDATE 2 - ngày 25/7/2021 - Dự phóng đỉnh dịch với 2 mô hình.
Từ ngày 25/06/2021, số ca nhiễm Corona trong ngày tại TP. HCM tăng đột biến lên 724 ca. Từ đó đến nay, xu hướng tăng vẫn là chủ đạo, mặc dù các biện pháp cách li, phong toả được thực hiện có vẻ quyết liệt hơn, cụ thể là chỉ thị 16 được áp dụng từ này 09/07.
Trước tiên hãy cùng nhìn lại đồ thị về số ca nhiễm cho tới thời điểm này, như ở Hình 1 bên dưới.
Khi chỉ nhìn bằng mắt thường vào các cột dữ liệu thực tế, bạn cũng có thể thấy là hình dáng của hàm mũ - tức là dịch bệnh đang lây lan theo cấp số nhân. (Nếu vẽ đồ thị theo lang logarithm thì bạn sẽ thấy hình dáng của một đường thẳng dốc lên.) Tính theo từng tuần, số ca đã tăng từ: 2958 lên 4721, sau một tuần lên 10295, rồi tuần sau đó lên 23913.😇
Ý nghĩa của hàm mũ: nếu tốc độ lây lan là cố định, ví dụ: nếu mỗi người nhiễm đều lây cho 1.3 người khác, trong một ngày, thì tổng số ca bệnh mỗi ngày sẽ nhân thêm 2.3. Khi đó ta nói dịch bệnh tăng theo cấp số nhân, cũng tức là tăng theo hàm mũ.
Dùng toán để mô phỏng thực tế
Còn khi xem xét kĩ hơn, ta có thể rút ra một vài thông tin chi tiết hơn về dịch bệnh. Ở đây mình sử dụng mô hình gần như tối giản với 3 tham số
Dùng thuật toán bình phương tối thiểu thông dụng (Least square), ta sẽ tính được các tham số A, B, C.
Trong 3 tuần vừa qua, tham số B vào khoảng 0.11, tức là mỗi ca bệnh trung bình sẽ lây cho e^B-1 ≈ 0.12 người khác, trong một ngày. Từ đó suy ra một người nhiễm, trong 2 tuần, có thể lây cho 0.12*14 = 1.68 người khác.
Xem xét dữ liệu quá khứ
Bây giờ hãy nhìn lại quá khứ và xem dịch bệnh tiến triển theo thời gian như thế nào. Ở Hình 2, mình đã dùng dữ liệu của 2 tuần trước đó để "giả dự đoán" , tức là ngoại suy ra số ca nhiễm, từ ngày 10/7 tới 16/7. Ta thấy kết quả dự đoán khá phù hợp, tuy có thấp hơn thực tế một chút ở những ngày cuối.
Dự phóng cho 7 ngày tới
Bây giờ hãy dùng dữ liệu từ 2 tuần vừa qua để dự báo cho tuần tới. Ở hình 3, mình ngoại suy để tính toán số ca nhiễm dự phóng cho 7 ngày tới. Nếu tình hình lây nhiễm không nhiều thay đổi nhiều thì số ca nhiễm vào ngày 23/7 là khoảng 55000 người.
Một tín hiệu tích cực đó là hệ số mũ "B" giảm từ 0.11 về 0.10 . [xem hình 2 và hình 3]. Tuy rằng khá ít, nhưng cho thấy tín hiệu về tính hiệu quả của việc tiêm vaccine, cũng như của 1 tuần áp dụng chỉ thị 16.
Mong rằng công tác chống dịch hiệu quả hơn nữa để kết quả thực tế sẽ thấp hơn kết quả dự đoán ở đây.
Phần kết
Trước tiên , mình mong dịch bệnh mau chóng được ngăn chặn.
Thứ hai, mình mong thông qua bài này giúp mọi người có dược cái nhìn tổng quan về dịch bệnh hiện tại và trong tương lai gần, để không quá hốt hoảng khi nhìn thấy con số ca nhiễm tăng lên 50000 hay 100000.
Thứ ba, mình mong các công cụ thống kê được sử dụng tốt hơn để đối phó với dịch bệnh, ví dụ như phân tích tác động của biến thể Delta tại Việt Nam. Và từ đó để mọi người có thể tự điều chỉnh hành vi sinh hoạt, kế hoạch làm việc,... cho tương lai gần.
Happy Ending ....
_______________________________________________________________________
Chú thích thêm
Nếu chỉ dùng mô hình hàm mũ f(x) = A e^{Bx} thì số ca bệnh thực tế cao hơn khá nhiều so với dự báo, và dẫn tới kết luận là dịch bệnh tăng nhanh hơn cả hàm mũ/cấp số nhân.
Điều này vẫn có thể đúng do khả năng lây nhiễm mạnh của biến thể Delta, cũng như công tác chống dịch không hiệu quả, ví dụ như lây nhiễm chéo trong khu cách li.
Nhưng cá nhân mình tin tưởng mô hình có tham số C hơn, vì cho fitting tốt hơn, và có lẽ công tác chống dịch đang dần có hiệu quả.
Nguồn dữ liệu
Update 23/07/2021 - Review lại kết quả dự báo ở Hình 3
Đã được một tuần kể từ khi bài viết được đăng, mình cập nhật thêm dữ liệu thực tế ở đây. Có thể thấy rằng tổng số ca nhiễm đã giảm rõ rệt so với đường ngoại suy màu cam, thể hiện việc dịch bệnh đang dần được kiểm soát.
Chúc cho công tác chống dịch hiệu quả và khoa học hơn nữa.
Update 25/07/2021 - Dự phóng cho đỉnh dịch
Khoa học - Công nghệ
/khoa-hoc-cong-nghe
Bài viết nổi bật khác
- Hot nhất
- Mới nhất