Giới thiệu dự án

Mô hình Y học gia đình (Family Medicine) và Chăm sóc sức khỏe ban đầu (Primary Health Care - PHC) đóng vai trò xương sống trong hệ thống y tế toàn cầu. Theo Tổ chức Y tế Thế giới (WHO), việc củng cố mạng lưới y tế cơ sở giúp giảm thiểu tới 30.000 ca tử vong trẻ em mỗi ngày tại hơn 100 quốc gia. Tại Việt Nam, Quyết định số 935/QĐ-BYT và Kế hoạch nhân rộng giai đoạn 2016–2020 của Bộ Y tế đã xác định việc phát triển mạng lưới Phòng khám Bác sĩ gia đình (PKBSGĐ) là giải pháp chiến lược nhằm giảm tải bệnh viện tuyến trên, khắc phục tình trạng phân bổ nguồn lực không đồng đều và nâng cao trải nghiệm điều trị cho người dân.

+-------------------------------------------------------------------------------+
|                             BỐI CẢNH NGHIÊN CỨU                              |
|                                                                               |
|  [Hệ thống y tế cơ sở] ---> [Thí điểm PKBSGĐ] ---> [Thách thức đánh giá CLDV] |
|        (Quá tải)            (Bệnh viện Quận 2)      (Thiếu công cụ chuẩn hóa) |
+-------------------------------------------------------------------------------+

Vấn đề thực tiễn và bài toán đặt ra (Problem Statement)

Mặc dù PKBSGĐ Bệnh viện Quận 2 là một trong những đơn vị thí điểm sớm nhất tại TP. Hồ Chí Minh, ban quản lý bệnh viện đối mặt với những thách thức cụ thể:

  • Thiếu hụt bộ công cụ chuẩn hóa: Chưa có thang đo lượng hóa chính xác các chiều kích chất lượng dịch vụ (CLDV) phù hợp với nhóm bệnh nhân cao tuổi và bệnh mãn tính tại tuyến y tế quận/huyện.
  • Độ lệch kỳ vọng của bệnh nhân: Khoảng cách giữa chất lượng kỹ thuật lâm sàng (Technical Quality) và chất lượng chức năng/trải nghiệm phục vụ (Functional Quality) chưa được phân tách rõ ràng.
  • Điểm nghẽn vận hành (Operational Bottlenecks): Thủ tục hành chính và vệ sinh cơ sở vật chất chưa có chỉ số định lượng để làm căn cứ cải tiến liên tục (CQI - Continuous Quality Improvement).

Mục tiêu nghiên cứu (Project Objectives)

  1. Kiểm định психометрический (Psychometric properties): Đánh giá độ tin cậy nội bộ (Cronbach’s Alpha) và giá trị cấu trúc (Exploratory Factor Analysis - EFA) của thang đo CLDV và sự hài lòng hiệu chỉnh.
  2. Đo lường định lượng CLDV: Xác định điểm trung bình từng thành phần CLDV và chỉ số hài lòng chung của người bệnh.
  3. Mô hình hóa tác động: Xác định mối quan hệ tương quan, hồi quy đa biến giữa các chiều kích CLDV, đặc điểm nhân khẩu học với mức độ hài lòng tổng thể.

Phạm vi và giới hạn (Scope & Limitations)

  • Phạm vi không gian: Phòng khám Bác sĩ gia đình, Bệnh viện Quận 2, TP. Hồ Chí Minh.
  • Phạm vi thời gian: Thu thập dữ liệu thực địa từ tháng 04/2017 đến tháng 07/2017.
  • Giới hạn kỹ thuật: Nghiên cứu cắt ngang mô tả (Cross-sectional study), chưa đánh giá mối quan hệ nhân quả dài hạn theo thời gian (Longitudinal analysis).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Nghiên cứu so sánh và kế thừa các mô hình kinh điển nhằm lựa chọn khung lý thuyết tối ưu:

Tiêu chí Mô hình SERVQUAL (Parasuraman, 1988) Mô hình Kang & James (2004) Mô hình nghiên cứu đề xuất (Nghiên cứu này)
Cơ chế đo lường Đánh giá khoảng cách: Kỳ vọng (Expectation) - Cảm nhận (Perception) Phân tách Chất lượng Kỹ thuật, Chức năng và Hình ảnh tổ chức Đánh giá CLDV cảm nhận trực tiếp (SERVPERF biến thể), tích hợp yếu tố chi phí
Thang đo áp dụng 7 điểm Likert 7 điểm Likert 5 điểm Likert (Tối ưu cho người cao tuổi, học vấn đa dạng)
Ưu điểm Đầy đủ 5 khía cạnh cốt lõi Bao quát cả năng lực chuyên môn và hình ảnh Tối ưu hóa thời gian khảo sát, giảm mệt mỏi nhận thức (Cognitive fatigue)
Hạn chế Trắc nghiệm dài dòng, dễ gây nhiễu sai lệch Khó lượng hóa chính xác chất lượng kỹ thuật y khoa Cần tái kiểm định cấu trúc nhân tố thực nghiệm

Phân loại yêu cầu theo mô hình MoSCoW

  • Must have: Thang đo kiểm soát được hệ số Cronbach's Alpha $\ge 0.70$, hệ số KMO $\ge 0.80$, phân tích EFA trích xuất phương sai tích lũy $\ge 50%$.
  • Should have: Chuyển đổi thang đo 7 mức độ sang 5 mức độ Likert, loại bỏ các mục hỏi gây nhiễu ngữ nghĩa kép (Double-barreled questions).
  • Could have: Tích hợp kiểm định tương quan giữa nguồn thông tin tiếp cận và mức độ trung thành của bệnh nhân.
  • Won't have (lần này): Tích hợp hồ sơ bệnh án điện tử (EMR) để đo lường tự động thời gian chờ khám.

Thiết kế hệ thống đo lường và cấu trúc thang đo

Mô hình nghiên cứu phân rã CLDV thành các thành phần chức năng tác động trực tiếp đến sự hài lòng:

Công nghệ và công cụ xử lý dữ liệu

  • EpiData Entry v3.1: Thiết lập biểu mẫu nhập liệu kiểm soát logic (Check codes), thực hiện quy trình nhập liệu kép (Double Data Entry) nhằm triệt tiêu sai số do con người.
  • Stata/SE v14.0: Xử lý thống kê mô tả, thống kê suy luận, phân tích nhân tố khám phá (EFA) với phép quay Promax và mô hình hồi quy tuyến tính Stepwise.
  • Python 3.10 (Data Pipeline Validation): Hỗ trợ kiểm toán ma trận tương quan và trực quan hóa phân bố điểm số.

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

Nghiên cứu ứng dụng công thức tính cỡ mẫu ước lượng giá trị trung bình với độ chính xác cao:

$$n = Z_{1-\alpha/2}^2 \cdot \frac{\sigma^2}{d^2}$$

  • $Z_{1-\alpha/2} = 1.96$ (với mức ý nghĩa thống kê $\alpha = 0.05$).
  • Độ lệch chuẩn kỳ vọng $\sigma = 1.27$ (tham chiếu từ nghiên cứu tiền khả thi của Trần Văn Khanh, 2015).
  • Sai số biên chấp nhận được $d = 0.12$.
  • Cỡ mẫu lý thuyết tối thiểu: $n = 1.96^2 \cdot \frac{1.27^2}{0.12^2} \approx 431$ đối tượng. Dự phòng 5% hao hụt $\rightarrow$ Cỡ mẫu đích $N = 452$. Thực tế thu thập thành công $N = 468$ phiếu khảo sát hợp lệ.

Implementation và kết quả

Quy trình xử lý và làm sạch dữ liệu

Dữ liệu được làm sạch thông qua thuật toán chuẩn hóa tuyến tính (Linear Stretching) để chuyển đổi thang Likert 5 điểm về thang phần trăm (0% – 100%):

$$\text{Mean}_{\text{normalized}} = \frac{\text{Mean} - \text{Min}}{\text{Max} - \text{Min}} \times 100% = \frac{\text{Mean} - 1}{4} \times 100%$$

* Stata script for psychometric evaluation and regression modeling
clear
import delimited "survey_data_pkbsgd.csv", clear

* 1. Kiểm tra tính tin cậy nội bộ Cronbach's Alpha cho toàn thang đo CLDV
alpha c1-c21 c23-c25, item

* 2. Kiểm định KMO và Bartlett Test
factortest c1-c21 c23-c25

* 3. Phân tích nhân tố khám phá (EFA) với phép quay Promax
factor c1-c21 c23-c25, pcf mineigen(1.0)
rotate, promax(3) blanks(0.4)

* 4. Hồi quy tuyến tính đa biến đánh giá tác động CLDV lên Hài lòng
regress satisfaction cldv_mean
regress satisfaction f_dapung f_thuantien f_bacsi f_dieutri f_vesinh

Kiểm định chất lượng thang đo (Testing & Validation)

  1. Độ tin cậy thang đo ban đầu: Thang đo CLDV ban đầu gồm 25 biến quan sát đạt Cronbach's $\alpha = 0.87$. Mục hỏi C22 ("Có đường dây nóng, hòm thư góp ý") có hệ số tương quan biến - tổng (Item-rest correlation) xấp xỉ $0.00$. Việc loại bỏ biến C22 giúp nâng Cronbach's $\alpha$ lên 0.88 đối với thang đo 24 biến.
  2. Kiểm định tính phù hợp dữ liệu EFA:
    • Hệ số Kaiser-Meyer-Olkin (KMO) đạt 0.862 ($> 0.80$, chứng minh mẫu nghiên cứu rất phù hợp cho phân tích nhân tố).
    • Kiểm định Bartlett đạt giá trị $\chi^2 = 5375.65$ với bậc tự do $df = 378$ ($p < 0.001$), bác bỏ hoàn toàn giả thuyết các biến không tương quan.
  3. Cấu trúc nhân tố khám phá (EFA): Trích xuất thành công 6 nhân tố có Eigenvalue $> 1.05$, giải thích 59.9% tổng phương sai trích:
    • Nhân tố 1 (Hài lòng): 5 biến (C19, C26, C27, C28, C29) - $\alpha = 0.83$.
    • Nhân tố 2 (Đáp ứng của NVYT): 4 biến (C1, C5, C13, C14) - $\alpha = 0.84$.
    • Nhân tố 3 (Sự thuận tiện): 4 biến (C8, C10, C12, C17) - $\alpha = 0.84$.
    • Nhân tố 4 (Năng lực bác sĩ): 3 biến (C3, C15, C16) - $\alpha = 0.84$.
    • Nhân tố 5 (Đảm bảo điều trị lâu dài): 6 biến (C2, C4, C20, C21, C23, C25) - $\alpha = 0.85$.
    • Nhân tố 6 (Vệ sinh môi trường): 2 biến (C9, C11) - $\alpha = 0.84$.
+-------------------------------------------------------------------------------+
|                      KẾT QUẢ PHÂN TÍCH NHÂN TỐ EFA                            |
|                                                                               |
|  KMO Measure: 0.862  |  Bartlett's Test: p < 0.001  |  Total Variance: 59.9%  |
|                                                                               |
|  [Nhân tố 1] Hài lòng chung (5 items) -------------------------> alpha = 0.83 |
|  [Nhân tố 2] Đáp ứng của NVYT (4 items) -----------------------> alpha = 0.84 |
|  [Nhân tố 3] Sự thuận tiện (4 items) --------------------------> alpha = 0.84 |
|  [Nhân tố 4] Năng lực bác sĩ (3 items) ------------------------> alpha = 0.84 |
|  [Nhân tố 5] Đảm bảo điều trị lâu dài (6 items) ---------------> alpha = 0.85 |
|  [Nhân tố 6] Vệ sinh môi trường (2 items) ---------------------> alpha = 0.84 |
+-------------------------------------------------------------------------------+

Kết quả nghiên cứu thực nghiệm

Khảo sát trên $N = 468$ bệnh nhân (trong đó 64.7% là nữ, 74.8% từ 50 tuổi trở lên, 97.0% sử dụng Bảo hiểm Y tế và 62.4% mắc bệnh mãn tính kéo dài trên 3 năm) thu được các chỉ số thực chứng sau:

Chiều kích đánh giá Điểm TB (Likert 1-5) Độ lệch chuẩn (SD) Điểm chuẩn hóa (%)
Năng lực của bác sĩ 4.73 $\pm 0.43$ 93.3%
Đảm bảo theo dõi/điều trị lâu dài 4.52 $\pm 0.40$ 88.1%
Đáp ứng của nhân viên y tế 4.27 $\pm 0.75$ 81.7%
Sự thuận tiện 4.20 $\pm 0.66$ 79.9%
Vệ sinh cơ sở vật chất 3.66 $\pm 0.88$ 66.4%
Chất lượng dịch vụ tổng thể (CLDV) 4.34 $\pm 0.42$ 83.6%
Mức độ hài lòng của bệnh nhân 4.40 $\pm 0.54$ 85.0%

Mô hình hồi quy tuyến tính

  • Mối tương quan Pearson giữa CLDV và Hài lòng đạt $r = 0.69$ ($p < 0.001$).
  • Hệ số xác định $R^2 = 0.468$: CLDV giải thích trực tiếp 46.8% sự biến thiên của mức độ hài lòng.
  • Phương trình hồi quy đơn biến: $\text{Hài lòng} = \beta_0 + 0.89 \times \text{CLDV}$ (Khoảng tin cậy 95%: $0.81 - 0.98, p < 0.001$).
  • Thứ tự tác động chuẩn hóa ($\beta$): Đáp ứng NVYT ($\beta = 0.60$) $>$ Sự thuận tiện ($\beta = 0.54$) $>$ Năng lực bác sĩ ($\beta = 0.41$) $=$ Theo dõi lâu dài ($\beta = 0.41$) $>$ Vệ sinh ($\beta = 0.39$).

Đổi mới và đóng góp

  1. Hiệu chỉnh thang đo thực chứng cho y tế cơ sở Việt Nam: Chuyển đổi thành công thang đo 7 điểm (vốn dễ gây quá tải nhận thức cho người bệnh cao tuổi tại Việt Nam) sang thang đo 5 điểm với các thuộc tính psychometric đạt độ tin cậy vượt trội ($\alpha = 0.88$).
  2. Phát hiện cấu trúc nhân tố mới: Tách biệt thành công khía cạnh "Vệ sinh" thành một nhân tố độc lập (trước đây thường bị gộp chung vào Phương tiện hữu hình), phản ánh chính xác tác động tiêu cực của điều kiện hạ tầng đang sửa chữa đến trải nghiệm bệnh nhân.
  3. Lượng hóa mức độ trung thành y tế tuyến cơ sở: 92.7% bệnh nhân cam kết tiếp tục khám lâu dài ($4.53 \pm 0.66$) và 90.4% sẵn sàng giới thiệu phòng khám cho người thân ($4.54 \pm 0.76$). Kênh truyền miệng (Word-of-mouth) chiếm 68.6% nguồn tiếp cận và là nhóm có điểm hài lòng cao nhất ($p = 0.012$).

Ứng dụng thực tế và triển khai

Kịch bản can thiệp nâng cao chất lượng dịch vụ

Dựa trên kết quả định lượng, một lộ trình hành động can thiệp được thiết lập nhằm nâng cao chỉ số vận hành cho các PKBSGĐ tuyến quận/huyện:

+------------------------------------------------------------------------------------+
|                         LỘ TRÌNH TRIỂN KHAI VÀ CAN THIỆP                           |
|                                                                                    |
|  [Giai đoạn 1: 0 - 3 tháng]     [Giai đoạn 2: 3 - 6 tháng]   [Giai đoạn 3: 6 - 12 tháng]   |
|  Cải tạo hạ tầng vệ sinh,     Tối ưu hóa thủ tục hành     Số hóa EMR, xây dựng     |
|  sắp xếp luồng tiếp nhận      chính, đào tạo giao tiếp    hệ thống feedback tức thì|
+------------------------------------------------------------------------------------+
  • Giai đoạn 1 (0–3 tháng - Quick Wins): Tập trung cải tạo hệ thống nhà vệ sinh bệnh viện (chỉ tiêu đang có điểm thấp nhất: $3.03 \pm 1.31$) và bố trí lại bảng hướng dẫn đường đi ($4.44 \pm 0.73$).
  • Giai đoạn 2 (3–6 tháng - Process Optimization): Tối ưu quy trình nhận thuốc và xét nghiệm Bảo hiểm Y tế (thủ tục nhanh gọn chỉ đạt $3.60 \pm 1.30$, với 26.1% phản hồi tiêu cực). Triển khai đào tạo kỹ năng lắng nghe và đồng cảm cho điều dưỡng/nhân viên tiếp đón.
  • Giai đoạn 3 (6–12 tháng - Digital Transformation): Triển khai ứng dụng quản lý lịch hẹn thông minh của BSGĐ, giảm thời gian chờ đợi trung bình xuống dưới 15 phút.

Hạn chế và hướng phát triển

  • Hạn chế chọn mẫu: Sử dụng phương pháp chọn mẫu thuận tiện tại một trung tâm duy nhất (Bệnh viện Quận 2), có thể hạn chế tính đại diện khi suy rộng cho các vùng nông thôn hoặc cơ sở y tế tư nhân.
  • Nhiễu bối cảnh môi trường: Thời điểm khảo sát trùng với giai đoạn bệnh viện thi công cải tạo, gây sụt giảm điểm số tại tiêu chí vệ sinh và không gian chờ.
  • Hướng nghiên cứu tương lai:
    1. Ứng dụng mô hình Structural Equation Modeling (SEM) để kiểm tra các tác động trung gian (Mediating effects) của niềm tin thương hiệu.
    2. Mở rộng khảo sát đa trung tâm tại toàn bộ 24 quận/huyện trên địa bàn TP. Hồ Chí Minh.

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

  • Sinh viên & Học viên Y tế công cộng: Nắm vững phương pháp luận kiểm định thang đo định lượng (EFA, Cronbach's Alpha) và phương pháp xây dựng đề cương nghiên cứu chuẩn mực.
  • Nhà quản lý bệnh viện & Hoạch định chính sách: Sở hữu bằng chứng thực nghiệm để phân bổ ngân sách tối ưu (ưu tiên năng lực giao tiếp nhân viên và thủ tục BHYT trước khi đầu tư trang thiết bị đắt đỏ).
  • Bác sĩ gia đình & Nhân viên y tế: Nhận thức rõ tầm quan trọng của việc giải thích đơn thuốc và lập kế hoạch theo dõi định kỳ đối với bệnh nhân mãn tính.
  • Người bệnh: Được hưởng lợi trực tiếp từ các cải tiến vận hành, rút ngắn thời gian làm thủ tục và cải thiện môi trường khám chữa bệnh.

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

1. Tại sao nghiên cứu lại chuyển đổi từ thang đo 7 điểm Likert sang 5 điểm Likert?

Qua khảo sát thử nghiệm trên 30 đối tượng, nhóm bệnh nhân thực tế phần lớn là người cao tuổi ($\ge 50$ tuổi chiếm 74.8%) và có trình độ học vấn dưới THCS chiếm tỷ lệ đáng kể (43.6%). Thang đo 7 mức độ gây quá tải nhận thức, dẫn đến hiện tượng trả lời ngẫu nhiên hoặc từ chối tham gia. Thang đo 5 điểm giúp cải thiện độ chính xác và tính khả thi thu thập dữ liệu.

2. Biến quan sát C22 (Đường dây nóng, hòm thư góp ý) tại sao bị loại bỏ?

Trong phân tích độ tin cậy, biến C22 có hệ số tương quan biến - tổng hiệu chỉnh (Item-rest correlation) bằng $0.00$. Điều này phản ánh bệnh nhân phòng khám BSGĐ hiếm khi sử dụng hoặc không quan tâm đến hòm thư góp ý truyền thống. Việc loại bỏ C22 giúp cải thiện hệ số Cronbach's Alpha toàn diện của thang đo từ 0.87 lên 0.88.

3. Thành phần nào của chất lượng dịch vụ có tác động mạnh nhất đến sự hài lòng?

Phân tích hồi quy đơn biến cho thấy nhân tố Đáp ứng của nhân viên y tế (NVYT) có trọng số tác động chuẩn hóa cao nhất ($\beta = 0.60$), tiếp theo là Sự thuận tiện ($\beta = 0.54$) và Năng lực của bác sĩ ($\beta = 0.41$).

4. Dữ liệu nghiên cứu được kiểm soát sai lệch thông tin như thế nào?

Dữ liệu được kiểm soát qua 3 bước: (1) Chỉ chọn bệnh nhân đã khám tại phòng khám ít nhất 1 lần trong 3 tháng gần nhất để giảm sai lệch hồi tưởng (Recall bias); (2) Nhập liệu độc lập 2 lần bằng phần mềm EpiData v3.1; (3) Kiểm tra logic và làm sạch ma trận dữ liệu bằng Stata v14.0.

5. Yếu tố nhân khẩu học nào ảnh hưởng đáng kể nhất đến sự hài lòng?

Tuổi tác và trình độ học vấn có mối liên hệ mật thiết: Nhóm bệnh nhân $\ge 60$ tuổi có điểm hài lòng trung bình ($4.59 \pm 0.41$) cao hơn có ý nghĩa thống kê ($p < 0.001$) so với nhóm $< 50$ tuổi ($4.11 \pm 0.60$). Ngược lại, bệnh nhân có học vấn từ THPT trở lên có kỳ vọng khắt khe hơn, dẫn đến điểm hài lòng thấp hơn ($4.15 \pm 0.63, p < 0.001$).

Kết luận

Nghiên cứu của tác giả Nguyễn Trường Viên (Đại học Y Dược TP. Hồ Chí Minh) đã chuẩn hóa thành công bộ công cụ đánh giá Chất lượng Dịch vụ và Sự Hài lòng của người bệnh tại Phòng khám Bác sĩ gia đình Bệnh viện Quận 2. Với cỡ mẫu thực nghiệm $N = 468$, công trình đã chứng minh mô hình 5 chiều kích CLDV (với Cronbach's $\alpha = 0.88, KMO = 0.862, R^2 = 0.468$) sở hữu độ tin cậy và giá trị cấu trúc vững chắc. Kết quả nghiên cứu cung cấp cơ sở khoa học quan trọng cho các nhà quản lý y tế trong việc nâng cao năng lực phục vụ của điều dưỡng, đơn giản hóa thủ tục Bảo hiểm Y tế và cải tạo cơ sở hạ tầng, hiện thực hóa mục tiêu chăm sóc sức khỏe toàn diện và lấy người bệnh làm trung tâm.