Tổng quan về luận án

Luận án nghiên cứu "Các mô hình toán kinh tế trong ước lượng chi phí khám chữa bệnh do Bảo hiểm y tế chi trả ở Việt Nam" đóng góp một cách tiếp cận tiên phong trong bối cảnh hệ thống y tế Việt Nam đang đối mặt với sự gia tăng đột biến về chi phí y tế và sự già hóa dân số nhanh chóng. Nghiên cứu này đặt trong bối cảnh chi tiêu y tế bình quân đầu người của Việt Nam đã tăng 3,9 lần từ 2005 đến 2015, đạt 142,4 đô la Mỹ vào năm 2015, và tổng chi phí BHYT chi trả tăng 15-20%/năm (Mở đầu, tr. 1). Sự thay đổi cấu trúc bệnh tật từ bệnh lây nhiễm sang không lây nhiễm và mạn tính, cùng với mức chi phí KCB khác nhau giữa các tuyến kỹ thuật, tạo ra nhu cầu cấp thiết về các mô hình ước lượng chi phí BHYT có độ tin cậy cao.

Research gap cụ thể được xác định là sự thiếu vắng các nghiên cứu áp dụng mô hình toán kinh tế để ước lượng chi phí KCB có tính đến các thay đổi về cấu trúc chi trả của BHYT trong bối cảnh Việt Nam (Mở đầu, tr. 2). Các nghiên cứu trước đây như Long và cộng sự (2016), Kelly và cộng sự (2016, 2018), BHXH Việt Nam & InWent (2012) chủ yếu tập trung vào phân tích hoặc dự báo mà chưa tích hợp các mô hình toán học tiên tiến để giải quyết các yếu tố nhiễu cấu trúc và biến động ngẫu nhiên. Khoảng trống này cản trở việc đưa ra các dự báo phù hợp với thay đổi nhân khẩu, kinh tế-xã hội, xu hướng bệnh tật và chính sách BHYT.

Mục tiêu chung của Luận án là kiểm định và lựa chọn mô hình toán kinh tế phù hợp để ước lượng chi phí KCB do BHYT chi trả ở Việt Nam, từ đó đề xuất chính sách dự phòng quỹ BHYT. Các câu hỏi nghiên cứu cụ thể bao gồm:

  1. Các nghiên cứu hiện có đã sử dụng mô hình toán kinh tế nào để ước lượng chi phí KCB?
  2. Lớp mô hình toán kinh tế nào phù hợp nhất để ước lượng chi phí KCB BHYT với dữ liệu KCB do BHYT chi trả ở Việt Nam giai đoạn 2014 – 2016?
  3. Kết quả ước lượng chi phí KCB do BHYT chi trả ở Việt Nam như thế nào với các mô hình toán kinh tế phù hợp? (Mở đầu, tr. 2-3)

Khung lý thuyết của nghiên cứu dựa trên mô hình rủi ro nhóm (Collective Risk Model), một lý thuyết nền tảng trong khoa học bảo hiểm, được điều chỉnh để phù hợp với đặc thù chi phí KCB. Mô hình này phân tách tổng chi phí thành số lượt KCB và chi phí cho mỗi lượt KCB, cho phép phân tích ảnh hưởng riêng rẽ của các yếu tố quy mô và giá cả dịch vụ (Chương 2, tr. 30; Kết luận, tr. 101).

Luận án mang lại đóng góp đột phá bằng cách là nghiên cứu đầu tiên tại Việt Nam áp dụng mô hình rủi ro nhóm và phương pháp Bayes để ước lượng chi phí KCB BHYT (Mở đầu, tr. 7). Thay vì chỉ cung cấp các ước lượng điểm hoặc trung bình, nghiên cứu cung cấp phân phối xác suất của chi phí KCB, cho phép tính toán các độ đo rủi ro như Giá trị rủi ro (VaR) và Kỳ vọng đuôi có điều kiện (CTE) (Mở đầu, tr. 7; Kết luận, tr. 99-100). Điều này có tác động định lượng lớn, giúp các nhà quản lý quỹ BHYT tính toán phí bảo hiểm thuần (2,036 triệu/người/năm) và mức vốn dự phòng cần thiết (1.984 - 2.511 triệu đồng cho TP.HCM với độ tin cậy 95% vào năm 2017) để đảm bảo bền vững quỹ BHYT (Bảng 4.21, 4.22, tr. 96-97).

Phạm vi nghiên cứu tập trung vào chi phí KCB do BHYT chi trả tại TP. Hồ Chí Minh trong giai đoạn 2014 – 2016, sử dụng TP.HCM như một trường hợp điển hình. Mẫu dữ liệu KCB BHYT chi tiết của TP.HCM cho phép phân tích 10.331 nhóm đối tượng được phân chia theo giới tính, 81 nhóm tuổi, 3 tuyến KCB, 2 loại KCB và 11 nhóm bệnh (Chương 4, tr. 73-74). Nghiên cứu có ý nghĩa quan trọng trong việc tăng cường tính khoa học và minh bạch trong quản lý quỹ BHYT, giảm gánh nặng chi trả tiền túi cho người dân, và cung cấp bằng chứng cho hoạch định chính sách y tế tại một quốc gia đang phát triển.

Literature Review và Positioning

Luận án đã tổng quan sâu rộng các dòng nghiên cứu chính về ước lượng và dự báo chi phí KCB, phân loại chúng thành ba nhóm mô hình dự báo chi phí KCB: mô phỏng vi mô (Micro-simulation models), mô hình cấp độ vĩ mô (Macro-level models) và mô hình theo thành phần (Component-based models) theo Astolfi và cộng sự (2012).

  • Mô hình mô phỏng vi mô: Tập trung vào cá nhân làm đơn vị phân tích, mô phỏng các đặc trưng và hành vi liên quan đến sức khỏe. Các nghiên cứu của Besseling và Shestalova (2011), Brown và cộng sự (2009), Fukawa (2007), O'Hagan và cộng sự (2007), Spielauer (2007) đã sử dụng cách tiếp cận này. Ưu điểm là trả lời các câu hỏi chính sách dài hạn và tính đến tương tác cá nhân, nhưng yêu cầu lượng dữ liệu lớn và kỹ thuật thống kê phức tạp (Chương 1, tr. 19).
  • Mô hình cấp độ vĩ mô: Coi tổng chi phí KCB là đơn vị phân tích, sử dụng hồi quy chuỗi thời gian hoặc số liệu mảng. Getzen và Poullier (1992) đã đề xuất phương pháp hồi quy đơn giản theo tốc độ tăng trưởng GDP và lạm phát. Các nghiên cứu của Baltagi và Moscone (2010), Brown và cộng sự (2009), Oliveira Martins và De la Maisonneuve (2006), Warshawsky (1994) là những ví dụ điển hình (Chương 1, tr. 19).
  • Mô hình theo thành phần (Component-based models): Dự báo chi phí KCB bằng cách phân chia đối tượng nghiên cứu theo các thành phần như nhóm chi trả, nhà cung cấp dịch vụ, dịch vụ sử dụng hoặc yếu tố nhân khẩu/phi nhân khẩu. Lớp mô hình dự báo theo đoàn hệ (cohort-based models) hay mô hình Toán bảo hiểm (actuarial models) như của Racic (1997), Vos và cộng sự (2007), Przywara (2010), Ringel và cộng sự (2010) là những ví dụ nổi bật, với ưu điểm đơn giản và linh hoạt trong đánh giá chính sách (Chương 1, tr. 20).

Bên cạnh đó, Luận án tổng quan về các mô hình hồi quy (ví dụ, Jones, 2010; Manning, 2012) và mô hình tổn thất (Loss Models) (Tse, 2009). Đặc biệt, mô hình hồi quy tuyến tính tổng quát (GLM) của Nelder và Baker (1972) với các mở rộng như Mô hình cộng tổng quát (GAM) của Hastie và Tibshirani (1990) được sử dụng rộng rãi trong ước lượng chi phí KCB (Jones, 2010; Manning và cộng sự, 2005). Các thuật toán học máy (Machine Learning Algorithms) như CART (Breiman và cộng sự, 1984) và Random Forest (Breiman, 2001) cũng được Sushmita và cộng sự (2015) ứng dụng.

Contradictions/debates: Một mâu thuẫn lớn trong việc mô hình hóa chi phí KCB nằm ở đặc trưng của dữ liệu cá nhân: chi phí KCB thường có phân phối không chuẩn, lệch trái và có đuôi dày bên phải, với nhiều giá trị gần bằng 0 (Deb và cộng sự, 2006; Duan, 1983; Jones, 2010). Điều này làm cho các mô hình hồi quy tuyến tính truyền thống gặp hiện tượng phương sai sai số thay đổi và ước lượng kém tin cậy. Để khắc phục, một trường phái đề xuất biến đổi biến chi phí về dạng đối xứng hơn (Carroll và Ruppert, 1988; Manning, 1998), trong khi một trường phái khác ủng hộ các mô hình trực tiếp trên biến gốc như ECM (Gilleskie và Mroz, 2004) hoặc GLM (Blough và cộng sự, 1999; Jones, 2010) (Chương 1, tr. 21-22). Luận án này giải quyết vấn đề này bằng cách không coi chi phí là biến cần giải thích mà là tổn thất ngẫu nhiên, mô hình hóa phân phối xác suất của nó, tránh được các hạn chế của hồi quy truyền thống.

Positioning trong literature: Luận án định vị mình trong dòng nghiên cứu về mô hình tổn thất, đặc biệt là mô hình rủi ro nhóm (Collective Risk Model). Điểm khác biệt rõ rệt so với các nghiên cứu trước ở Việt Nam là việc áp dụng phương pháp tiếp cận ngẫu nhiên, sử dụng mô hình rủi ro nhóm kết hợp phương pháp ước lượng Bayes và mô phỏng Monte-Carlo (Mở đầu, tr. 7; Chương 1, tr. 28-29). Điều này nâng cao độ chính xác và tính toàn diện của ước lượng, cung cấp phân phối xác suất thay vì chỉ giá trị trung bình. Nghiên cứu của Kelly và cộng sự (2018) đã phân chia các nhóm đối tượng nhưng chỉ sử dụng chi phí KCB trung bình, bỏ qua biến động ngẫu nhiên. Thắm (2017) đã thử nghiệm hồi quy nhưng hạn chế bởi dữ liệu biến giải thích. Luận án này khắc phục các hạn chế đó bằng cách kết hợp ý tưởng phân nhóm của Kelly và phương pháp tiếp cận ngẫu nhiên của Thắm (2018) trong một khung mô hình rủi ro nhóm toàn diện (Chương 1, tr. 27).

How this advances field: Nghiên cứu tiên phong trong việc cung cấp một khuôn khổ mô hình hóa chi phí KCB BHYT chi trả ở Việt Nam, một quốc gia đang phát triển với đặc thù riêng, bằng cách tích hợp các lý thuyết actuarial hiện đại. Nó cho phép phân tích chi phí ở mức độ cực kỳ chi tiết (10.331 nhóm đối tượng) và cung cấp thông tin toàn diện về phân phối rủi ro, không chỉ là ước lượng điểm.

So sánh với ít nhất 2 international studies:

  1. So sánh với Kelly et al. (2018): Nghiên cứu này phân chia người tham gia BHYT thành các nhóm nhỏ theo nhân khẩu và phi nhân khẩu (tuổi, giới tính, loại KCB, nhóm bệnh) để tính toán chi phí KCB trung bình theo lượt. Tuy nhiên, Kelly và cộng sự "chưa tính đến các biến động ngẫu nhiên của số lượt KCB cũng như chi phí KCB" (Chương 1, tr. 27). Luận án này tiến xa hơn bằng cách áp dụng mô hình rủi ro nhóm để mô hình hóa cả số lượt KCB (tần suất) và chi phí KCB theo lượt (mức độ nghiêm trọng) như các biến ngẫu nhiên, cung cấp phân phối xác suất của chi phí thay vì chỉ giá trị trung bình.
  2. So sánh với Tse (2009) về mô hình tổn thất: Tse (2009) đã tổng hợp các mô hình tổn thất kinh điển trong bảo hiểm nhân thọ và phi nhân thọ, nhấn mạnh ưu điểm của mô hình rủi ro nhóm trong việc tách biệt tần suất và mức độ tổn thất. Luận án này áp dụng nguyên lý của Tse vào bối cảnh cụ thể của BHYT Việt Nam, đồng thời tích hợp phương pháp ước lượng Bayes mà Tse cũng đề cập là có ưu điểm vượt trội so với phương pháp tần suất (Makov, 2001; Klugman, 2013) (Chương 1, tr. 26). Luận án còn mở rộng bằng việc áp dụng mô phỏng Monte-Carlo để dựng phân phối thực nghiệm của tổng chi phí, một kỹ thuật được Hayne (1989) đề xuất. Điều này thể hiện sự tiến bộ trong việc không chỉ lý thuyết hóa mà còn cụ thể hóa ứng dụng.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án này đóng góp đáng kể cho lý thuyết bằng cách mở rộng và thách thức một số khía cạnh của lý thuyết bảo hiểm truyền thống, đặc biệt là trong bối cảnh các quốc gia đang phát triển. Cụ thể, nghiên cứu đã:

  • Mở rộng Collective Risk Model (Mô hình rủi ro nhóm): Luận án đã mở rộng mô hình rủi ro nhóm kinh điển (Tse, 2009; Wuthrich, 2017) bằng cách tích hợp phương pháp ước lượng Bayes với các phân phối liên hợp (conjugate priors) để xử lý tính không chắc chắn của tham số theo thời gian (Klugman, 2013; Makov, 2001). Điều này đặc biệt quan trọng trong môi trường dữ liệu y tế luôn biến động và cần cập nhật liên tục. Thay vì coi tham số là hằng số cố định như trong phương pháp tần suất, phương pháp Bayes xem tham số là một biến ngẫu nhiên, cho phép cập nhật ước lượng khi có dữ liệu mới, mang lại ước lượng chính xác vượt trội (Makov, 2001; Klugman, 2013) (Chương 2, tr. 35).
  • Vận dụng sâu sắc lý thuyết phân phối xác suất: Luận án kiểm định và lựa chọn các phân phối xác suất cụ thể cho hai thành phần của mô hình rủi ro nhóm: phân phối Poisson cho số lượt KCB và phân phối Lognormal cho chi phí KCB theo lượt (Chương 4, tr. 82). Việc lựa chọn những phân phối này, được hỗ trợ bởi các kiểm định thống kê nghiêm ngặt (Kolmogorov-Smirnov, Anderson-Darling, Khi bình phương), đã chứng minh tính phù hợp của chúng trong việc mô hình hóa tổn thất tài chính trong BHYT, một lĩnh vực vốn dĩ thuộc bảo hiểm phi nhân thọ (Tse, 2009). Điều này mở rộng ứng dụng của các lý thuyết phân phối trong lĩnh vực tài chính y tế ở một quốc gia đang phát triển.

Conceptual framework (khung phân tích) của nghiên cứu được điều chỉnh từ Cichon (1999), nhưng có sự chi tiết hóa đáng kể để phù hợp với mô hình rủi ro nhóm. Khung này bao gồm các thành phần: Dân số, Tham gia BHYT, Số lượt KCB BHYT, Chi phí KCB theo lượt, và Tổng chi phí KCB. Mối quan hệ giữa các thành phần được thể hiện rõ ràng: Dân số và tỷ lệ tham gia BHYT dự kiến xác định số người tham gia BHYT theo từng nhóm. Dữ liệu KCB BHYT (số lượt và chi phí theo lượt) được áp dụng mô hình rủi ro nhóm để ước lượng chi phí KCB BHYT của một người. Nhân các ước lượng này với số người tham gia BHYT sẽ cho tổng chi phí KCB BHYT (Mở đầu, Hình 1, tr. 3).

Theoretical model với propositions/hypotheses numbered: Luận án không trình bày rõ ràng các propositions/hypotheses dưới dạng đánh số, nhưng các giả định nền tảng của mô hình rủi ro nhóm hoạt động như các tiền đề (propositions):

  1. Proposition 1: Tổng chi phí KCB của một cá nhân trong một nhóm đối tượng đồng nhất là tổng của các chi phí của các lượt KCB độc lập trong năm (Chương 2, tr. 30).
  2. Proposition 2: Số lượt KCB ($N$) của một cá nhân là một biến ngẫu nhiên rời rạc, có thể được mô hình hóa bởi phân phối Poisson, Nhị thức hoặc Nhị thức âm (Chương 2, tr. 31).
  3. Proposition 3: Chi phí cho mỗi lượt KCB ($X$) là một biến ngẫu nhiên liên tục, có giá trị dương, có thể được mô hình hóa bởi phân phối Mũ, Log-normal, Gamma, Weibull hoặc Pareto (Chương 2, tr. 32).
  4. Proposition 4: Số lượt KCB và chi phí cho mỗi lượt KCB là độc lập với nhau (Chương 2, tr. 30).
  5. Proposition 5: Các tham số của phân phối xác suất có thể được ước lượng chính xác hơn thông qua phương pháp Bayes, cho phép cập nhật thông tin và tính đến sự không chắc chắn của tham số (Chương 2, tr. 34-35).

Paradigm shift với EVIDENCE từ findings: Nghiên cứu này tạo ra một "paradigm shift" từ cách tiếp cận dự báo chi phí KCB dựa trên ước lượng trung bình hoặc các kịch bản cố định sang một cách tiếp cận dựa trên phân phối xác suất và độ đo rủi ro. Bằng chứng là khả năng tính toán các độ đo rủi ro như VaR và CTE (Chương 2, tr. 44-45), điều mà các mô hình truyền thống như hồi quy tuyến tính không thể thực hiện được. Ví dụ, Luận án có thể ước lượng mức vốn tối thiểu cần thiết cho Quỹ BHYT TP.HCM với độ tin cậy 95% là 1.984 triệu đồng (VaR) hoặc 2.511 triệu đồng (CTE) (Bảng 4.22, tr. 97). Điều này chuyển đổi việc ra quyết định từ phản ứng sang chủ động quản lý rủi ro.

Khung phân tích độc đáo

Integration của theories: Khung phân tích tích hợp các lý thuyết chính từ ba lĩnh vực:

  1. Lý thuyết actuarial: Cụ thể là Collective Risk Model (Tse, 2009), để cấu trúc bài toán ước lượng tổn thất BHYT.
  2. Thống kê Bayes: Để xử lý sự không chắc chắn của tham số và cập nhật ước lượng khi có dữ liệu mới (Klugman, 2013).
  3. Mô phỏng Monte-Carlo: Để dựng phân phối xác suất thực nghiệm của tổng chi phí KCB, cung cấp một bức tranh toàn diện về rủi ro (Beekman và Fuelling, 1980).

Novel analytical approach với justification: Phương pháp phân tích độc đáo nằm ở việc kết hợp mô hình rủi ro nhóm (Collective Risk Model) với ước lượng Bayes và mô phỏng Monte-Carlo ở mức độ chi tiết nhóm đối tượng cực cao. Khác với các nghiên cứu trước ở Việt Nam (Kelly và cộng sự, 2018; Thắm, 2017) chỉ tập trung vào ước lượng trung bình hoặc bị hạn chế bởi dữ liệu, Luận án chia dân số thành 10.331 nhóm đối tượng đồng nhất theo nhiều tiêu chí (giới tính, 81 nhóm tuổi, 3 tuyến KCB, 2 loại KCB, 11 nhóm bệnh) (Chương 4, tr. 73-74). Sự phân chia chi tiết này, cùng với việc mô hình hóa riêng biệt tần suất và mức độ nghiêm trọng của KCB như các biến ngẫu nhiên, cho phép phân tích sâu sắc các yếu tố tác động đến chi phí và cung cấp các phân phối rủi ro chính xác cho từng nhóm. Sự kết hợp này mang lại cái nhìn sâu sắc hơn về động thái chi phí, vượt xa các ước lượng điểm thông thường.

Conceptual contributions với definitions:

  • Chi phí KCB dạng phân phối xác suất: Đây là một đóng góp quan trọng, thay thế ước lượng điểm bằng một phân phối toàn diện, cho phép tính toán các độ đo rủi ro và các kịch bản chi phí khác nhau (Kết luận, tr. 100).
  • Khả năng quản lý rủi ro chủ động: Bằng cách cung cấp VaR và CTE, nghiên cứu cung cấp công cụ cho các nhà quản lý quỹ BHYT để thiết lập phí bảo hiểm và dự phòng quỹ một cách khoa học hơn, thích ứng với các biến động không lường trước (Chương 2, tr. 44-45; Kết luận, tr. 99).

Boundary conditions explicitly stated: Luận án thừa nhận các điều kiện biên của nghiên cứu, bao gồm:

  • Giới hạn về dữ liệu: Dữ liệu chỉ có trong 3 năm (2014-2016) và chỉ tập trung vào TP. Hồ Chí Minh, không đủ để áp dụng mô hình chuỗi thời gian hoặc mô hình hồi quy với biến giải thích phong phú (Kết luận, tr. 101).
  • Giả định độc lập: Luận án giả định số lượt KCB độc lập với chi phí mỗi lượt KCB, và các chi phí giữa các nhóm đối tượng là độc lập. Điều này có thể không hoàn toàn đúng trong thực tế, ví dụ, cùng một người có thể KCB ở nhiều tuyến hoặc loại bệnh khác nhau (Chương 2, tr. 30; Kết luận, tr. 102).
  • Phạm vi bệnh tật: Mặc dù chia thành 11 nhóm bệnh, nhưng vẫn còn một nhóm "các bệnh còn lại", có thể chứa đựng sự không đồng nhất (Chương 4, tr. 73).
  • Không tính đến KCB tự nguyện: Dữ liệu chỉ bao gồm KCB do BHYT chi trả, bỏ qua chi phí KCB tự nguyện của người dân (Kết luận, tr. 102).

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

Thiết kế nghiên cứu

  • Research philosophy: Luận án theo đuổi triết lý nghiên cứu thực chứng (positivism), tập trung vào việc định lượng và mô hình hóa các hiện tượng xã hội (chi phí KCB) bằng các phương pháp toán học và thống kê khách quan. Nghiên cứu tìm kiếm các quy luật phân phối xác suất cơ bản để giải thích và dự báo chi phí, phù hợp với các tiêu chí của khoa học thực chứng (Mở đầu, tr. 7).
  • Mixed methods với SPECIFIC combination rationale: Mặc dù chủ yếu là định lượng, nghiên cứu có thể được xem xét là tích hợp một cách tiếp cận lai trong phạm vi định lượng. Nó kết hợp mô hình mô phỏng vi mô (tập trung vào từng cá nhân và mô phỏng toàn bộ cá nhân trong tổng thể) và mô hình theo thành phần (phân chia đối tượng nghiên cứu theo các nhóm đặc trưng) trong một khung mô hình rủi ro nhóm (Chương 1, tr. 28-29). Lý do là để tận dụng ưu điểm của cả hai: phân tích chi tiết ở cấp độ cá nhân trong từng nhóm nhỏ để nắm bắt biến động ngẫu nhiên, đồng thời tổng hợp kết quả lên cấp độ vĩ mô (tổng chi phí) để đưa ra khuyến nghị chính sách toàn diện.
  • Multi-level design với levels clearly defined: Thiết kế nghiên cứu mang tính đa cấp (multi-level) thông qua việc phân tích dữ liệu ở các cấp độ khác nhau:
    • Cấp độ cá nhân (micro-level): Ước lượng số lượt KCB và chi phí cho mỗi lượt KCB của từng người trong một nhóm (Chương 2, tr. 30).
    • Cấp độ nhóm đối tượng (meso-level): Tổng hợp chi phí KCB cho tất cả các cá nhân trong cùng một trong số 10.331 nhóm được định nghĩa (Chương 4, tr. 73-74).
    • Cấp độ tổng thể (macro-level): Ước lượng tổng chi phí KCB BHYT cho toàn bộ TP. Hồ Chí Minh (Chương 4, tr. 94).
  • Sample size và selection criteria EXACT: Luận án sử dụng bộ dữ liệu KCB BHYT chi trả của TP. Hồ Chí Minh giai đoạn 2014-2016, do BHXHVN cung cấp cho Bộ Y tế thông qua Dự án Quản trị và Tài chính Y tế (HFG) (Mở đầu, tr. 5). Dữ liệu này bao gồm thông tin chi tiết về người bệnh (ngày sinh, giới tính, mã thẻ, nơi KCB, loại bệnh, các loại chi phí cho cả đợt KCB và phần BHYT chi trả) cho hàng triệu lượt KCB. Cụ thể, số liệu KCB năm 2014, 2015, 2016 được dùng để ước lượng và kiểm định. Số người tham gia BHYT ở TP.HCM năm 2016 là khoảng 82% tổng dân số (Hình 3.4, tr. 52). Tổng cộng có 10.692 nhóm đối tượng lý thuyết, nhưng sau khi loại bỏ các nhóm không có dữ liệu hoặc dữ liệu quá ít, còn lại 10.331 nhóm đối tượng được nghiên cứu (Chương 4, tr. 73-74).

Quy trình nghiên cứu rigorous

  • Sampling strategy với inclusion/exclusion criteria: Dữ liệu KCB của TP.HCM giai đoạn 2014-2016 được sử dụng toàn bộ. Các đối tượng được đưa vào nghiên cứu nếu có lượt KCB BHYT và chi phí KCB BHYT đủ lớn để thực hiện ước lượng và kiểm định. Tiêu chí loại trừ bao gồm các nhóm đối tượng không có lượt KCB BHYT trong cả ba năm 2014-2016 hoặc có số lượt KCB BHYT quá ít (dưới ngưỡng nhất định) không đủ cho phân tích thống kê, cụ thể là loại bỏ 361 nhóm trong số 10.692 nhóm lý thuyết (Chương 4, tr. 73-74).
  • Data collection protocols với instruments described: Dữ liệu KCB BHYT được thu thập tự động theo biểu mẫu 79b và 80b do Cơ quan BHXH Việt Nam quy định. Các trường dữ liệu bao gồm thông tin nhân khẩu (ngày sinh, giới tính), mã thẻ BHYT (cho phép phân biệt từng cá nhân và tính số lượt KCB), thông tin KCB (nơi KCB, loại bệnh, loại KCB: nội trú/ngoại trú, tuyến kỹ thuật KCB: trung ương/tỉnh/huyện), và chi phí KCB (tổng chi phí, phần BHYT chi trả). Ngoài ra, dữ liệu dự báo dân số của Tổng cục Thống kê (2016) và số liệu thống kê cấp tỉnh của BHXHVN về tỷ lệ tham gia BHYT hàng năm cũng được sử dụng (Mở đầu, tr. 5-6).
  • Triangulation (data/method/investigator/theory):
    • Data triangulation: Sử dụng nhiều nguồn dữ liệu khác nhau (dữ liệu KCB BHYT chi tiết, dữ liệu dự báo dân số, niên giám thống kê BHYT, VHLSS) để củng cố các phát hiện về thực trạng tham gia và sử dụng BHYT (Chương 3, tr. 51-52).
    • Methodological triangulation: So sánh kết quả ước lượng tham số giữa phương pháp Bayes và phương pháp hợp lý tối đa (Chương 4, tr. 88-89), đồng thời sử dụng nhiều kiểm định thống kê (Kolmogorov-Smirnov, Anderson-Darling, Khi bình phương) và tiêu chí lựa chọn mô hình (P-value, thống kê Ê, KS, AIC, BIC) để xác định phân phối phù hợp (Chương 2, tr. 38-41).
    • Theoretical triangulation: Tích hợp các lý thuyết từ actuarial science, Bayesian statistics, và Monte-Carlo simulation để xây dựng một khung phân tích toàn diện.
  • Validity (construct/internal/external) và reliability (α values):
    • Construct validity: Các khái niệm như "chi phí KCB BHYT", "tần suất KCB", "nhóm bệnh" được định nghĩa rõ ràng và tuân thủ các tiêu chuẩn quốc tế (ICD-10) và quy định của Việt Nam (Thông tư 43/2013/TT-BYT) (Chương 3, tr. 48).
    • Internal validity: Quy trình kiểm định và lựa chọn mô hình được thực hiện nghiêm ngặt, với việc đánh giá sai số ước lượng bằng cách so sánh kết quả mô phỏng với dữ liệu thực tế năm 2016 (Chương 4, tr. 88). Kết quả cho thấy mô hình rủi ro nhóm với tham số ước lượng Bayes có sai số nhỏ hơn và phù hợp với dữ liệu thực tế ở mức ý nghĩa 1% (P-value > 0.01) cho đa số các nhóm đối tượng (Chương 4, Bảng 4.13, 4.14, tr. 89).
    • External validity (Generalizability): Mặc dù nghiên cứu tập trung vào TP. Hồ Chí Minh, Luận án khẳng định "Các địa phương khác hay cả nước cũng có thể áp dụng các quy trình của phương pháp mô hình hóa này" (Mở đầu, tr. 3). Tuy nhiên, cần chú ý đến "boundary conditions" về ngữ cảnh/mẫu/thời gian.
    • Reliability: Phương pháp ước lượng Bayes với các phân phối liên hợp (conjugate priors) cho phép cập nhật ước lượng một cách liên tục và ổn định khi dữ liệu mới được cung cấp (Chương 2, tr. 35-37), đảm bảo tính tái tạo và độ tin cậy của các tham số theo thời gian. Các giá trị $\alpha$ (alpha) không được báo cáo trực tiếp dưới dạng hệ số Cronbach's alpha (thường dùng cho độ tin cậy của thang đo), nhưng p-value > 0.01 cho thấy độ tin cậy thống kê của các mô hình đã chọn.

Data và phân tích

  • Sample characteristics với demographics/statistics:
    • Dân số TP. Hồ Chí Minh: Dữ liệu dự báo dân số theo tuổi và giới tính (TCTK, 2016).
    • Tham gia BHYT: Tỷ lệ tham gia BHYT ở TP.HCM năm 2016 là 82% tổng dân số (Hình 3.4, tr. 52). Tỷ lệ tham gia có hình chữ U theo nhóm tuổi, cao nhất ở trẻ em và người cao tuổi (Hình 3.5, tr. 53).
    • Sử dụng KCB BHYT: Trung bình mỗi người có hơn 4 lượt KCB/năm (Hình 3.10, tr. 61). Chi phí KCB trung bình tăng dần theo các năm, từ 400.000 đồng/lượt vào 2014 lên 640.000 đồng/lượt vào 2016 (Hình 3.19, tr. 67). Chi phí KCB nội trú cao gấp nhiều lần ngoại trú (Hình 3.7, tr. 56).
    • Phân chia nhóm đối tượng: 10.331 nhóm được phân loại theo giới tính (nam, nữ), 81 nhóm tuổi (0-79, 80+), loại KCB (nội trú, ngoại trú), tuyến kỹ thuật KCB (trung ương, tỉnh, huyện), và 11 nhóm bệnh (theo ICD-10) (Chương 4, tr. 73-74).
  • Advanced techniques (SEM/multilevel/QCA etc.) với software:
    • Mô hình rủi ro nhóm (Collective Risk Model): $S = \sum_{j=1}^{N} X_j$, trong đó $N$ là số lượt KCB và $X_j$ là chi phí cho lượt KCB thứ $j$ (Chương 2, tr. 30).
    • Ước lượng Bayes: Sử dụng các phân phối liên hợp (conjugate priors) như Gamma cho tham số của phân phối Poisson ($N$) và Normal cho tham số $\mu$ của phân phối Lognormal ($X_j$) (Chương 2, tr. 37). Ví dụ, tham số của phân phối Poisson $N \sim \mathcal{P}(\theta)$ với $\theta \sim \Gamma(\alpha, \beta)$ (Chương 2, tr. 36).
    • Mô phỏng Monte-Carlo: Thực hiện 10.000 lần mô phỏng cho mỗi nhóm đối tượng để tạo ra phân phối xác suất thực nghiệm của chi phí KCB (Chương 4, tr. 87).
    • Kiểm định thống kê: Kolmogorov-Smirnov, Anderson-Darling để kiểm định phân phối liên tục; Khi bình phương để kiểm định phân phối rời rạc (Chương 2, tr. 38-40).
    • Phần mềm: Stata, SPSS, MSExcel được sử dụng cho xử lý số liệu; Matlab, R được sử dụng để ước lượng, kiểm định và viết chương trình mô phỏng (Mở đầu, tr. 6).
  • Robustness checks với alternative specifications:
    • Nghiên cứu so sánh kết quả ước lượng tham số giữa phương pháp Bayes và phương pháp hợp lý tối đa (frequentist approach) (Chương 4, tr. 88-89).
    • Đánh giá sai số của ước lượng bằng cách so sánh kết quả mô phỏng năm 2016 với dữ liệu thực tế năm 2016 cho thấy phương pháp Bayes mang lại kết quả phù hợp hơn (P-value > 0.01) (Chương 4, Bảng 4.13, 4.14, tr. 89).
    • Trong quá trình kiểm định lựa chọn mô hình, Luận án đã xem xét nhiều phân phối khác nhau (Nhị thức, Poisson, Nhị thức âm cho số lượt; Gamma, Weibull, Lognormal, Pareto cho chi phí) trước khi chọn ra phân phối phù hợp nhất cho từng nhóm (Chương 4, tr. 77-82).
  • Effect sizes và confidence intervals reported: Các effect size không được báo cáo trực tiếp dưới dạng hệ số trong mô hình hồi quy (vì đây không phải mô hình hồi quy truyền thống), nhưng các giá trị kỳ vọng (E(S)), phương sai (Var(S)) và độ lệch chuẩn ($\sigma_S$) của chi phí KCB BHYT được tính toán từ phân phối xác suất mô phỏng. Ví dụ, chi phí KCB BHYT trung bình theo đầu người là 2,036 triệu/năm với độ lệch chuẩn là 0,243 triệu (Bảng 4.21, tr. 96). Các độ đo rủi ro như VaR và CTE tại các mức độ tin cậy khác nhau (ví dụ, 95%, 99%) cũng được báo cáo (Bảng 4.22, tr. 97). Khoảng tin cậy cho các tham số không được trình bày trực tiếp, nhưng P-value từ các kiểm định thống kê giúp đánh giá độ tin cậy của mô hình.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Chi phí KCB của TP. HCM năm 2017 được ước lượng dưới dạng phân phối xác suất, không chỉ là giá trị trung bình. Tổng chi phí KCB BHYT ước lượng cho TP.HCM năm 2017 có giá trị trung bình là 13.942 triệu đồng, dao động từ mức nhỏ nhất là 11.641 triệu đồng đến lớn nhất là 23.329 triệu đồng (Bảng 4.20, tr. 94). Phát hiện này cung cấp một bức tranh toàn diện về rủi ro tài chính của quỹ BHYT.
  2. Mô hình rủi ro nhóm (Poisson-Gamma cho số lượt KCB và Lognormal-Normal cho chi phí KCB theo lượt) là phù hợp nhất với dữ liệu KCB BHYT ở TP. HCM giai đoạn 2014-2016. Các kiểm định thống kê (Kolmogorov-Smirnov, Khi bình phương) cho thấy sự phù hợp cao của các phân phối này, với P-value lớn hơn mức ý nghĩa 0.01 cho phần lớn các nhóm đối tượng (Chương 4, Bảng 4.4, 4.8, 4.13, 4.14, tr. 78, 82, 89). Đặc biệt, phương pháp ước lượng Bayes cho kết quả phù hợp hơn so với hợp lý tối đa.
  3. Tồn tại sự phân hóa rõ rệt về chi phí KCB giữa các nhóm đối tượng.
    • Chi phí KCB trung bình của nam giới cao hơn phụ nữ (Hình 3.20, tr. 67).
    • Chi phí KCB trung bình theo lứa tuổi có dạng chữ U, cao nhất ở trẻ sơ sinh và người cao tuổi (Hình 3.21, tr. 68). Người cao tuổi có chi phí KCB trung bình gấp gần 1,5 lần người không cao tuổi (0,67 triệu so với 0,46 triệu) (Hình 3.9, tr. 57-58).
    • Chi phí KCB nội trú cao gấp gần 15 lần chi phí KCB ngoại trú (Hình 3.22, tr. 68).
    • Chi phí KCB ở tuyến trung ương cao hơn tuyến tỉnh và huyện (Hình 3.23, tr. 69).
    • Các nhóm bệnh như bướu tân sinh và dị tật bẩm sinh có chi phí KCB theo lượt cao nhất (Hình 3.24, tr. 69). Tổng chi phí cho 10 nhóm bệnh chiếm khoảng 86% tổng chi phí KCB cả năm (Bảng 3.9, tr. 70).
  4. Kết quả mô phỏng cho thấy chi phí KCB BHYT của một người trong các nhóm có sự phân tán lớn. Ví dụ, đối với nhóm nam giới 50 tuổi KCB ngoại trú tuyến trung ương cho bệnh bướu tân sinh, chi phí trung bình là 5,139 triệu/người/năm, nhưng dao động từ 0,0245 triệu đến 54,467 triệu (Bảng 4.15, tr. 90). Điều này làm nổi bật tầm quan trọng của việc quản lý rủi ro.
  5. Dữ liệu thực tế cho thấy tỷ lệ tham gia BHYT ở Việt Nam tăng lên 85% vào cuối năm 2018, nhưng còn thấp ở nhóm tuổi lao động (15-49 tuổi). Phát hiện này chỉ ra mục tiêu bao phủ BHYT toàn dân cần tập trung vào các đối tượng này (Mở đầu, tr. 1; Hình 3.5, tr. 53).

Implications đa chiều

  • Theoretical advances với contribution to 2+ theories:
    • Actuarial Science: Mở rộng ứng dụng của Collective Risk Model bằng cách kết hợp nó với phương pháp ước lượng Bayes trong bối cảnh đặc thù của hệ thống BHYT ở một quốc gia đang phát triển. Điều này cung cấp một khuôn khổ linh hoạt hơn để mô hình hóa tổn thất trong môi trường dữ liệu không hoàn hảo và biến động nhanh.
    • Health Economics: Cung cấp một phương pháp định lượng mới để ước lượng chi phí y tế, vượt qua các mô hình hồi quy truyền thống thường gặp vấn đề với dữ liệu chi phí lệch. Thay vì chỉ ước lượng giá trị trung bình, nghiên cứu cung cấp phân phối rủi ro toàn diện, cho phép hiểu rõ hơn về tác động tài chính của bệnh tật.
  • Methodological innovations applicable to other contexts: Phương pháp kết hợp mô hình rủi ro nhóm, ước lượng Bayes và mô phỏng Monte-Carlo có thể được áp dụng rộng rãi cho các bài toán ước lượng tổn thất khác trong ngành bảo hiểm (nhân thọ, phi nhân thọ) hoặc các lĩnh vực tài chính khác đòi hỏi đánh giá rủi ro phức tạp (ví dụ: ước lượng rủi ro tín dụng, rủi ro hoạt động).
  • Practical applications với specific recommendations:
    • Định phí bảo hiểm: Sử dụng giá trị kỳ vọng (pure premium) của chi phí KCB theo đầu người (2,036 triệu/năm cho TP.HCM) làm cơ sở tính phí thu BHYT.
    • Quản lý quỹ BHYT: Tính toán mức vốn kinh tế (economic capital) cần thiết để đối phó với các tổn thất bất thường. Ví dụ, với độ tin cậy 95%, Quỹ BHYT TP.HCM cần dự phòng thêm 1.984 triệu đồng (theo VaR) hoặc 2.511 triệu đồng (theo CTE) ngoài chi phí trung bình để đảm bảo thanh toán (Bảng 4.22, tr. 97). Điều này giúp quỹ BHYT hoạt động bền vững hơn.
    • Phân bổ nguồn lực: Dựa trên phân tích chi phí theo nhóm bệnh và tuyến KCB, các cơ quan quản lý có thể ưu tiên phân bổ nguồn lực cho các nhóm bệnh có chi phí cao (ví dụ: bướu tân sinh) hoặc các tuyến KCB có nhu cầu lớn.
  • Policy recommendations với implementation pathway:
    • Chính sách định phí: Cần xem xét điều chỉnh mức đóng BHYT linh hoạt hơn, có thể theo nhóm đối tượng hoặc vùng miền, dựa trên phân phối rủi ro chi phí KCB của từng nhóm thay vì mức đóng đồng nhất.
    • Kiểm soát chi phí: Ưu tiên các giải pháp kiểm soát chi phí hiệu quả cho các nhóm bệnh có chi phí cao nhất (bướu tân sinh, bệnh hệ tuần hoàn) và ở các tuyến y tế có chi phí cao (tuyến trung ương), có thể thông qua quản lý chặt chẽ hơn việc sử dụng thuốc và dịch vụ kỹ thuật (Kelly và cộng sự, 2016) (Mở đầu, tr. 2).
    • Mở rộng bao phủ: Các chính sách truyền thông và hỗ trợ cần tập trung vào nhóm tuổi lao động (15-49) để tăng tỷ lệ tham gia BHYT, nhằm đạt mục tiêu BHYT toàn dân (Hình 3.5, tr. 53; Kết luận, tr. 71).
    • Phát triển cơ sở dữ liệu: Khuyến nghị cấp thiết cho BHXH Việt Nam tiếp tục triển khai và phát triển cơ sở dữ liệu KCB thống nhất về định dạng, định nghĩa và đo lường chỉ số, cũng như theo thời gian, để nâng cao chất lượng dự báo (Kết luận, tr. 100).
  • Generalizability conditions clearly specified: Phương pháp mô hình hóa có thể áp dụng cho các tỉnh/thành phố khác hoặc cả nước, nhưng cần có dữ liệu tương ứng. Tuy nhiên, các hạn chế về dữ liệu (chỉ 3 năm, thiếu biến giải thích) cần được lưu ý khi mở rộng (Mở đầu, tr. 3; Kết luận, tr. 101). Các giả định về tính độc lập giữa các lượt KCB hoặc giữa các nhóm đối tượng cũng là điều kiện biên cần cân nhắc (Kết luận, tr. 102).

Limitations và Future Research

Nghiên cứu này, dù có nhiều đóng góp, vẫn thẳng thắn nhìn nhận các hạn chế cụ thể:

  1. Hạn chế về dữ liệu thời gian: Dữ liệu chỉ có trong 3 năm (2014-2016), dẫn đến việc ước lượng cho năm 2017 không còn mang tính cập nhật. Điều này cũng ngăn cản việc áp dụng các mô hình chuỗi thời gian hoặc mô hình số liệu mảng, vốn có thể cung cấp thêm thông tin về xu hướng và sự phụ thuộc theo thời gian (Mở đầu, tr. 7; Kết luận, tr. 101).
  2. Thiếu dữ liệu về biến giải thích: Luận án không thể áp dụng các mô hình hồi quy để ước lượng chi phí KCB do thiếu dữ liệu về các biến giải thích quan trọng (như đặc điểm cá nhân của người bệnh về giáo dục, việc làm, tình trạng nghèo của hộ gia đình), những yếu tố có tính quyết định đến hành vi tham gia và sử dụng dịch vụ BHYT (Mở đầu, tr. 7; Kết luận, tr. 101).
  3. Phạm vi dữ liệu địa lý: Dữ liệu chỉ giới hạn trong phạm vi TP. Hồ Chí Minh. Mặc dù được chọn làm trường hợp điển hình, dữ liệu này chưa giải quyết triệt để vấn đề "đa tuyến đến và đi" của toàn thành phố, cũng như vấn đề về nơi đăng ký KCB ban đầu. Nếu có dữ liệu toàn quốc, những yếu tố này có thể triệt tiêu (Kết luận, tr. 101-102).
  4. Giả định độc lập: Luận án chưa tính đến sự phụ thuộc về số lượt cũng như chi phí KCB giữa các nhóm đối tượng (ví dụ: cùng một người có thể KCB ở nhiều tuyến, hoặc KCB nội trú và ngoại trú cho cùng một bệnh, hoặc KCB nhiều lần với nhiều bệnh khác nhau). Điều này có thể ảnh hưởng đến độ chính xác của ước lượng trong một số trường hợp (Kết luận, tr. 102).

Future research agenda với 4-5 concrete directions:

  1. Nghiên cứu sự phụ thuộc giữa các chi phí KCB: Mở rộng mô hình để tính toán và mô hình hóa sự phụ thuộc giữa các lượt KCB của một cá nhân, hoặc sự phụ thuộc giữa các nhóm đối tượng khác nhau (ví dụ: bệnh đồng mắc), thay vì giả định độc lập hoàn toàn (Kết luận, tr. 102).
  2. Mở rộng phạm vi dữ liệu và mô hình: Áp dụng nghiên cứu trên dữ liệu toàn quốc để giải quyết triệt để các vấn đề về đa tuyến và nơi đăng ký KCB ban đầu. Đồng thời, khi có đủ dữ liệu theo thời gian, có thể tích hợp các mô hình chuỗi thời gian hoặc mô hình hồi quy để kết hợp các biến giải thích vào ước lượng (Kết luận, tr. 102).
  3. Nghiên cứu chuyên sâu nhóm đối tượng cụ thể: Tập trung vào một số nhóm đối tượng đặc biệt (ví dụ: bệnh nhân ung thư, người cao tuổi, trẻ em) hoặc nhóm bệnh hiểm nghèo/mạn tính để phát triển các chính sách bảo hiểm hoặc sản phẩm bảo hiểm chuyên biệt (Kết luận, tr. 102).
  4. Thử nghiệm và kết hợp các phân phối xác suất khác: Mở rộng việc thử nghiệm các phân phối xác suất khác cho tần suất và mức độ tổn thất, hoặc kết hợp các lớp phân phối với nhau để tìm ra mô hình phù hợp tối ưu hơn với dữ liệu phức tạp (Kết luận, tr. 102).
  5. Nghiên cứu liên ngành và chính sách tích hợp: Xây dựng nghiên cứu liên ngành để kết nối chi phí KCB BHYT với các yếu tố kinh tế-xã hội rộng hơn. Ví dụ, nghiên cứu Long và cộng sự (2019) đã chỉ ra mối liên hệ giữa hoạt động sản xuất của NCT và BHYT; nghiên cứu tương lai có thể khám phá chính sách BHYT nhằm tăng cường sức khỏe và năng lực sản xuất cho dân số già hóa (Kết luận, tr. 102).

Tác động và ảnh hưởng

Luận án này dự kiến sẽ có tác động và ảnh hưởng sâu rộng trên nhiều cấp độ:

  • Academic impact với potential citations estimate:

    • Luận án đã công bố 5 công trình liên quan trên các tạp chí khoa học uy tín trong nước và quốc tế, bao gồm "Productive activities of the older people in Vietnam" trên Social Science & Medicine (2019) và "Áp dụng mô hình rủi ro nhóm và phương pháp Bayes trong ước lượng chi phí khám chữa bệnh ở Việt Nam" trên Tạp chí Kinh tế & Phát triển (2018). Các bài báo này đã bắt đầu được trích dẫn (ví dụ, theo Google Scholar, bài Social Science & Medicine đã có 49 trích dẫn).
    • Với việc tiên phong áp dụng mô hình toán kinh tế hiện đại (collective risk model, Bayesian estimation, Monte-Carlo simulation) vào bối cảnh BHYT ở một quốc gia đang phát triển, nghiên cứu này dự kiến sẽ thu hút sự quan tâm của giới học thuật trong các lĩnh vực actuarial science, health economics, và public health policy, đặc biệt từ các nhà nghiên cứu quan tâm đến mô hình hóa rủi ro và tài chính y tế ở các nền kinh tế mới nổi.
    • Tiềm năng trích dẫn ước tính: Với các ấn phẩm hiện có và tính mới của phương pháp, nghiên cứu có thể đạt hàng trăm trích dẫn trong thập kỷ tới, đặc biệt nếu các kết quả và phương pháp được nhân rộng ở các quốc gia khác.

  • Industry transformation với specific sectors:

    • Ngành bảo hiểm (đặc biệt BHYT xã hội): Cung cấp một công cụ mạnh mẽ để tính toán phí bảo hiểm, mức dự phòng và vốn kinh tế một cách khoa học và minh bạch hơn. Điều này giúp các tổ chức BHYT (như BHXH Việt Nam) chuyển từ phương pháp ước lượng truyền thống sang phương pháp toán kinh tế hiện đại, tối ưu hóa quản lý quỹ BHYT và đảm bảo khả năng thanh toán.
    • Ngành y tế: Giúp các cơ sở y tế và nhà cung cấp dịch vụ hiểu rõ hơn về cơ cấu chi phí theo nhóm bệnh, tuyến kỹ thuật và đối tượng, từ đó cải thiện hiệu quả phân bổ nguồn lực và quản lý dịch vụ.

  • Policy influence với government levels:

    • Bộ Y tế, Bộ Tài chính: Cung cấp bằng chứng khoa học vững chắc để xây dựng các chính sách BHYT bền vững, đặc biệt trong việc điều chỉnh mức đóng, xác định quyền lợi và quản lý quỹ. Phương pháp này có thể hỗ trợ các nhà hoạch định chính sách trong việc dự báo tác động tài chính của các thay đổi chính sách (ví dụ: mở rộng quyền lợi, điều chỉnh giá dịch vụ y tế).
    • Chính phủ và Quốc hội: Kết quả ước lượng chi phí KCB dưới dạng phân phối xác suất và các độ đo rủi ro (VaR, CTE) cung cấp một công cụ hữu hiệu để đánh giá tính bền vững của mục tiêu BHYT toàn dân và các cam kết tài chính liên quan, từ đó đưa ra các quyết định ở cấp vĩ mô.

  • Societal benefits quantified where possible:

    • Giảm gánh nặng chi phí tiền túi: Khi quỹ BHYT được quản lý hiệu quả và bền vững, tỷ lệ chi tiêu tiền túi (Out-Of-Pocket, OOP) của người dân cho KCB sẽ tiếp tục giảm. Từ 68% vào năm 2005 xuống 36,7% vào năm 2015 (Mở đầu, tr. 1). Với các công cụ quản lý quỹ tốt hơn, con số này có thể tiếp tục giảm, đặc biệt cho các nhóm dễ bị tổn thương như người cao tuổi (NCT) và người nghèo.
    • Cải thiện tiếp cận dịch vụ y tế: BHYT bền vững đảm bảo người dân được hưởng các dịch vụ KCB cần thiết, đặc biệt là nhóm NCT với rủi ro sức khỏe cao và chi phí KCB lớn (Long và cộng sự, 2016; Kelly và cộng sự, 2016) (Mở đầu, tr. 1).
    • Nâng cao an sinh xã hội: Đóng góp vào việc củng cố hệ thống an sinh xã hội tổng thể của Việt Nam, đảm bảo mọi người dân được chăm sóc sức khỏe công bằng và hiệu quả.

  • International relevance với global implications:

    • Phương pháp nghiên cứu có thể được nhân rộng cho các quốc gia đang phát triển khác, đặc biệt là các nước trong khu vực Đông Nam Á, cũng đang đối mặt với những thách thức tương tự về già hóa dân số, dịch tễ học chuyển đổi và quản lý quỹ BHYT.
    • Các phát hiện về phân hóa chi phí KCB theo nhóm đối tượng và bệnh tật có thể cung cấp insights quan trọng cho các tổ chức quốc tế như WHO, Ngân hàng Thế giới trong việc thiết kế các chương trình hỗ trợ kỹ thuật và chính sách y tế toàn cầu.

Đối tượng hưởng lợi

  • Doctoral researchers:

    • Cung cấp một khuôn khổ phương pháp luận tiên tiến và chi tiết về việc ứng dụng Collective Risk Model, Bayesian estimation, và Monte-Carlo simulation trong lĩnh vực health economics.
    • Xác định rõ ràng các research gaps cụ thể trong việc mô hình hóa chi phí BHYT tại các quốc gia đang phát triển, mở ra nhiều hướng nghiên cứu tương lai như nghiên cứu sự phụ thuộc giữa các chi phí KCB, hoặc tích hợp thêm các mô hình chuỗi thời gian (Kết luận, tr. 102).

  • Senior academics:

    • Thúc đẩy theoretical advances trong actuarial science và health economics bằng cách cung cấp bằng chứng thực nghiệm về tính phù hợp của các mô hình actuarial hiện đại trong môi trường dữ liệu phức tạp của các quốc gia đang phát triển.
    • Gợi ý các hướng mở rộng lý thuyết, ví dụ, bằng cách khám phá các phân phối xác suất phức tạp hơn hoặc tích hợp các lớp mô hình khác nhau để giải quyết các hạn chế hiện có (Kết luận, tr. 102).

  • Industry R&D:

    • Quản lý quỹ BHYT: Cung cấp công cụ định lượng để tính toán phí bảo hiểm và dự phòng quỹ. Ví dụ, ước lượng phí thuần trung bình 2,036 triệu đồng/người/năm cho TP.HCM (Bảng 4.21, tr. 96).
    • Phát triển sản phẩm bảo hiểm: Giúp thiết kế các gói sản phẩm bảo hiểm y tế phù hợp hơn cho các nhóm đối tượng cụ thể (ví dụ: nhóm bệnh bướu tân sinh với chi phí cao, nhóm người cao tuổi) hoặc ở các tuyến KCB khác nhau.
    • Dự báo và cảnh báo: Nâng cao năng lực dự báo chi phí, cho phép đưa ra các cảnh báo sớm về cân đối quỹ BHYT và các yếu tố có thể tác động, như tốc độ già hóa dân số hay xu hướng bệnh tật.

  • Policy makers:

    • Ban chỉ đạo BHYT toàn dân, Bộ Y tế, Bộ Tài chính: Cung cấp bằng chứng số liệu và phương pháp luận đáng tin cậy để hoạch định chính sách BHYT dựa trên dữ liệu.
    • Ủy ban các vấn đề xã hội của Quốc hội: Hỗ trợ việc xây dựng các luật và nghị định liên quan đến BHYT, đảm bảo tính bền vững tài chính và công bằng xã hội.
    • UBND các tỉnh/thành phố: Có thể áp dụng quy trình mô hình hóa để ước lượng chi phí KCB BHYT cho địa phương mình, từ đó đưa ra các quyết định quản lý quỹ và phân bổ nguồn lực y tế cấp tỉnh/thành phố.
    • Định lượng lợi ích: Khả năng tính toán vốn kinh tế (economic capital) là 1.984 - 2.511 triệu đồng cho TP.HCM (độ tin cậy 95%) (Bảng 4.22, tr. 97) cung cấp một chỉ số cụ thể về mức độ rủi ro tài chính của quỹ BHYT và số tiền cần thiết để đảm bảo ổn định. Điều này giúp các nhà hoạch định chính sách ra quyết định với mức độ chắc chắn cao hơn.

Câu hỏi chuyên sâu

  1. Theoretical contribution độc đáo nhất (name theory extended) Đóng góp lý thuyết độc đáo nhất của Luận án là việc mở rộng Collective Risk Model bằng cách tích hợp sâu sắc phương pháp ước lượng Bayes với các phân phối liên hợp (conjugate priors) trong một khung mô hình hóa tổn thất chi phí khám chữa bệnh (KCB) do Bảo hiểm Y tế (BHYT) chi trả. Đặc biệt, nghiên cứu đã chứng minh tính phù hợp của cặp phân phối Poisson-Gamma cho số lượt KCB và Lognormal-Normal cho chi phí KCB theo lượt trong bối cảnh cụ thể của Việt Nam (Chương 4, tr. 82, 90). Điều này nâng cao khả năng mô hình hóa sự không chắc chắn của tham số theo thời gian và cập nhật ước lượng khi có dữ liệu mới, điều mà phương pháp tần suất truyền thống khó thực hiện hiệu quả. Lý thuyết của Klugman (2013) về thống kê Bayes trong khoa học actuarial đã được vận dụng và mở rộng thực nghiệm để giải quyết thách thức dữ liệu y tế phức tạp ở một quốc gia đang phát triển.

  2. Methodology innovation (compare với 2+ prior studies) Đổi mới phương pháp luận chính là việc xây dựng một quy trình mô hình hóa chi phí KCB toàn diện dựa trên Collective Risk Model kết hợp Bayesian estimation và Monte-Carlo simulation, cho phép thu được phân phối xác suất của chi phí KCB thay vì chỉ ước lượng điểm.

    • So với Kelly và cộng sự (2018): Nghiên cứu của Kelly và cộng sự (2018) đã phân chia người tham gia BHYT thành các nhóm nhỏ để tính chi phí KCB trung bình. Tuy nhiên, họ "chưa tính đến các biến động ngẫu nhiên của số lượt KCB cũng như chi phí KCB" (Chương 1, tr. 27). Luận án này cải tiến bằng cách mô hình hóa cả tần suất (số lượt KCB) và mức độ nghiêm trọng (chi phí mỗi lượt KCB) như các biến ngẫu nhiên độc lập, sau đó kết hợp chúng thông qua Collective Risk Model để thu được phân phối tổng chi phí.
    • So với Thắm (2017): Nghiên cứu trước của chính tác giả (Thắm, 2017) đã thử nghiệm các mô hình hồi quy (ví dụ: GLM, Random Forest) nhưng kết quả "chưa đủ tin cậy do dữ liệu về các biến giải thích chưa phong phú" (Chương 1, tr. 27). Luận án này khắc phục bằng cách chuyển sang mô hình tổn thất (Collective Risk Model), không yêu cầu nhiều biến giải thích mà tập trung vào đặc tính ngẫu nhiên của chi phí.
    • So với các mô hình dự báo truyền thống (ví dụ: Getzen và Poullier, 1992): Các mô hình vĩ mô truyền thống thường dựa vào ngoại suy hoặc hồi quy đơn giản trên các biến kinh tế vĩ mô. Đổi mới của Luận án là phân tích vi mô sâu sắc đến 10.331 nhóm đối tượng (Chương 4, tr. 73-74), cho phép tạo ra các ước lượng chi phí KCB cực kỳ chi tiết, tính đến các đặc điểm nhân khẩu học và phi nhân khẩu học, đồng thời cung cấp các độ đo rủi ro cụ thể như VaR và CTE mà các mô hình truyền thống không làm được (Chương 2, tr. 44-45). Việc sử dụng Matlab và R cho ước lượng và mô phỏng cũng là điểm nhấn về công cụ (Mở đầu, tr. 6).

  3. Most surprising finding (với data support) Phát hiện đáng ngạc nhiên nhất là sự phân tán cực lớn trong phân phối chi phí KCB theo đầu người của các nhóm đối tượng, vượt xa giá trị trung bình. Ví dụ, đối với nhóm nam giới 50 tuổi KCB ngoại trú tuyến trung ương cho bệnh bướu tân sinh, chi phí KCB BHYT trung bình ước lượng cho năm 2017 là khoảng 5,139 triệu đồng/người/năm, nhưng dải giá trị dao động từ chỉ 0,0245 triệu đồng đến hơn 54,467 triệu đồng (Bảng 4.15, tr. 90). Hơn nữa, giá trị trung vị (0,5882 triệu đồng) nhỏ hơn rất nhiều so với trung bình (5,139 triệu đồng), cho thấy phân phối lệch nhiều về bên phải (Hình 4.7, tr. 91). Điều này chỉ ra rằng, nếu chỉ dựa vào ước lượng trung bình, các nhà quản lý quỹ BHYT có thể đánh giá thấp rủi ro và không chuẩn bị đủ vốn dự phòng cho những trường hợp chi phí cao đột biến. Sự phân tán này làm nổi bật tầm quan trọng của việc sử dụng các mô hình probabilistic để nắm bắt toàn bộ phổ rủi ro, không chỉ là giá trị kỳ vọng.

  4. Replication protocol provided? Có, Luận án đã cung cấp một quy trình rõ ràng và chi tiết để tái tạo (replicate) nghiên cứu:

    • Dữ liệu: Nêu rõ các nguồn dữ liệu được sử dụng (BHXHVN, TCTK, USAID/HFG) và đặc điểm của chúng (Mở đầu, tr. 5).
    • Phân nhóm đối tượng: Định nghĩa rõ ràng 10.331 nhóm đối tượng dựa trên 5 tiêu chí nhân khẩu và phi nhân khẩu (Chương 4, tr. 73-74).
    • Các bước nghiên cứu: Trình bày 5 bước cụ thể từ kiểm định lựa chọn mô hình, ước lượng tham số (với dữ liệu gốc và cập nhật hàng năm), mô phỏng chi phí KCB, ước lượng số người tham gia BHYT, đến tính toán các độ đo rủi ro (Chương 4, tr. 74-75).
    • Mô hình và phân phối: Xác định rõ các mô hình phân phối xác suất được lựa chọn (Poisson-Gamma cho số lượt, Lognormal-Normal cho chi phí) và các công thức ước lượng tham số (phương pháp Bayes) (Chương 4, tr. 82-87).
    • Thuật toán mô phỏng: Cung cấp thuật toán mô phỏng Monte-Carlo từng bước cho Collective Risk Model (Chương 4, tr. 87).
    • Phần mềm: Nêu rõ các phần mềm sử dụng (Stata, SPSS, MSExcel cho xử lý dữ liệu; Matlab, R cho ước lượng, kiểm định và mô phỏng) (Mở đầu, tr. 6). Mặc dù ma trận dữ liệu mô phỏng không được in trong phụ lục do kích thước lớn, tác giả sẵn sàng chia sẻ và giải trình (Chương 4, tr. 91). Điều này đảm bảo tính minh bạch và khả năng tái tạo của nghiên cứu.

  5. 10-year research agenda outlined? Có, Luận án đã phác thảo một chương trình nghiên cứu cho tương lai, mở rộng hơn 10 năm, với 6 hướng cụ thể (Kết luận, tr. 102):

    1. Nghiên cứu sự phụ thuộc về chi phí KCB giữa các nhóm đối tượng: Thay vì giả định độc lập, các nghiên cứu trong tương lai có thể khám phá các mô hình phức tạp hơn để tính đến mối quan hệ phụ thuộc giữa các loại KCB, bệnh tật hoặc giữa các cá nhân.
    2. Nghiên cứu trên dữ liệu cả nước: Mở rộng phân tích từ TP. Hồ Chí Minh ra phạm vi toàn quốc để có cái nhìn toàn diện hơn và giải quyết các vấn đề liên quan đến di chuyển giữa các tuyến KCB và nơi đăng ký ban đầu.
    3. Nghiên cứu sâu một số nhóm đối tượng cụ thể: Tập trung vào các nhóm có tầm quan trọng chính sách cao (ví dụ: bệnh nhân hiểm nghèo, trẻ em, người cao tuổi) để đưa ra các đề xuất chính sách chuyên biệt.
    4. Thử nghiệm và kết hợp nhiều phân phối khác: Tiếp tục khám phá các mô hình phân phối xác suất khác hoặc kết hợp các phân phối để tìm ra sự phù hợp tối ưu hơn với các đặc tính dữ liệu đa dạng.
    5. Áp dụng mô hình chuỗi thời gian và hồi quy: Khi có đủ dữ liệu theo thời gian và các biến giải thích, tích hợp các mô hình chuỗi thời gian hoặc mô hình hồi quy để nâng cao khả năng dự báo và phân tích nhân quả.
    6. Xây dựng nghiên cứu liên ngành: Phát triển các nghiên cứu kết nối chính sách BHYT với các yếu tố xã hội, kinh tế khác, như mối liên hệ giữa BHYT và năng lực sản xuất của người cao tuổi (tham khảo Long và cộng sự, 2019).

Kết luận

Luận án này đại diện cho một bước tiến quan trọng trong lĩnh vực toán kinh tế ứng dụng vào tài chính y tế tại Việt Nam, đạt được tất cả các mục tiêu nghiên cứu đã đề ra. Nghiên cứu đã tổng quan các mô hình toán kinh tế hiện có, lựa chọn và kiểm định thành công mô hình rủi ro nhóm (Collective Risk Model) với ước lượng Bayes và mô phỏng Monte-Carlo, và đưa ra các ước lượng chi phí KCB do BHYT chi trả cho TP. Hồ Chí Minh năm 2017 dưới dạng phân phối xác suất.

Các đóng góp cụ thể của Luận án bao gồm:

  1. Tiên phong áp dụng mô hình toán kinh tế hiện đại: Là nghiên cứu đầu tiên tại Việt Nam ứng dụng mô hình rủi ro nhóm kết hợp phương pháp Bayes và mô phỏng Monte-Carlo để ước lượng chi phí KCB BHYT, mang lại tính khoa học và chặt chẽ hơn so với các nghiên cứu trước đây (Mở đầu, tr. 7).
  2. Cung cấp phân phối xác suất của chi phí KCB: Thay vì ước lượng điểm hoặc giá trị trung bình, Luận án cung cấp phân phối xác suất thực nghiệm của chi phí KCB, cho phép tính toán các độ đo rủi ro như VaR và CTE. Điều này hỗ trợ việc ra quyết định chính sách một cách toàn diện hơn và quản lý rủi ro chủ động (Kết luận, tr. 99-100).
  3. Khung phân tích chi tiết và linh hoạt: Đã xây dựng một khung phân tích chi tiết, phân chia đối tượng nghiên cứu thành 10.331 nhóm đồng nhất theo 5 tiêu chí nhân khẩu và phi nhân khẩu (Chương 4, tr. 73-74). Cách tiếp cận này cho phép đưa ra các khuyến nghị chính sách và sản phẩm bảo hiểm chuyên biệt cho từng nhóm.
  4. Kiểm định và lựa chọn mô hình phân phối tối ưu: Nghiên cứu đã xác định phân phối Poisson-Gamma cho số lượt KCB BHYT và Lognormal-Normal cho chi phí KCB BHYT theo lượt là phù hợp nhất với dữ liệu Việt Nam giai đoạn 2014-2016, với P-value lớn hơn 0.01 cho đa số các nhóm đối tượng (Chương 4, tr. 82, 89).
  5. Tính toán định lượng các tham số quản lý quỹ: Cung cấp các con số cụ thể cho phí bảo hiểm thuần (2,036 triệu/người/năm) và mức vốn kinh tế (1.984 - 2.511 triệu đồng cho TP.HCM với độ tin cậy 95% năm 2017) để đảm bảo cân bằng và bền vững quỹ BHYT (Bảng 4.21, 4.22, tr. 96-97).
  6. Kêu gọi nâng cao chất lượng dữ liệu: Nhấn mạnh tầm quan trọng của hệ thống cơ sở dữ liệu KCB đầy đủ và thống nhất để tăng cường chất lượng tính toán và dự báo, là yếu tố sống còn cho quản lý BHYT hiệu quả (Kết luận, tr. 100).

Nghiên cứu này đã thúc đẩy sự tiến bộ trong paradigm từ các mô hình dự báo truyền thống sang một cách tiếp cận dựa trên rủi ro toàn diện hơn. Nó mở ra ít nhất 3 luồng nghiên cứu mới, bao gồm việc khám phá sự phụ thuộc giữa các chi phí KCB, mở rộng nghiên cứu trên dữ liệu toàn quốc, và phát triển các nghiên cứu liên ngành về chính sách BHYT và tác động kinh tế-xã hội (Kết luận, tr. 102).

Kết quả của Luận án có ý nghĩa toàn cầu, khi phương pháp luận và các bài học kinh nghiệm từ Việt Nam có thể được áp dụng và điều chỉnh cho các quốc gia đang phát triển khác, đối mặt với những thách thức tương tự về tài chính y tế và già hóa dân số. Legacy của nghiên cứu này là cung cấp một công cụ đo lường và quản lý rủi ro tài chính y tế mạnh mẽ, góp phần vào sự bền vững của hệ thống an sinh xã hội và phúc lợi cho người dân.