Giới thiệu dự án

Thị trường giáo dục sớm (Early Childhood Education - ECE) tại Việt Nam trong những năm gần đây chứng kiến tốc độ tăng trưởng kép hàng năm (CAGR) đạt xấp xỉ 11.2% nhờ sự gia tăng của tầng lớp trung lưu và nhận thức ngày càng cao của phụ huynh về "giai đoạn vàng" 0–6 tuổi. Tuy nhiên, sự phát triển nhanh chóng này đi kèm với nhiều thách thức nghiêm trọng: các vụ việc vi phạm đạo đức nghề nghiệp tại cơ sở mầm non, sự thiếu minh bạch trong cam kết đào tạo và chất lượng dịch vụ không đồng đều khiến mức độ lo âu và yêu cầu kiểm soát chất lượng từ phụ huynh tăng vọt.

Đồ án khóa luận "Nghiên cứu các yếu tố ảnh hưởng đến sự hài lòng của phụ huynh đối với chất lượng dịch vụ giáo dục sớm tại Học viện Little Einsteins" được thực hiện nhằm giải quyết bài toán định lượng hóa sự hài lòng của khách hàng trong phân khúc giáo dục mầm non và phát triển tư duy.

+-------------------------------------------------------------------------------+
|                             BỐI CẢNH NGHIÊN CỨU                               |
|  • Thị trường ECE tăng trưởng mạnh (~11.2% CAGR)                              |
|  • Phụ huynh lo ngại rủi ro an toàn, bạo hành và chất lượng đào tạo           |
|  • Thiếu khung chuẩn định lượng đánh giá chất lượng dịch vụ giáo dục sớm      |
+---------------------------------------+---------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|                             GIẢI PHÁP ĐỀ XUẤT                                 |
|  • Mô hình tích hợp SERVPERF (Cronin & Taylor) + Cảm nhận giá cả (Fornell)   |
|  • Quy trình nghiên cứu thực nghiệm 2 giai đoạn (Sơ bộ N=100, Chính thức N=212)|
|  • Phân tích đa biến: Cronbach's Alpha, EFA, OLS Regression, ANOVA            |
+-------------------------------------------------------------------------------+

Problem Statement và Pain Points

Các trung tâm và học viện giáo dục sớm hiện nay gặp phải những vấn đề cốt lõi:

  • Thiếu mô hình đo lường chuyên biệt: Các đơn vị thường áp dụng công cụ đánh giá chung chung, không nắm bắt được đặc thù vô hình và đồng thời của dịch vụ giáo dục cho trẻ nhỏ (dưới 6 tuổi).
  • Độ vênh kỳ vọng - thực tế: Sự bất cân xứng thông tin giữa cam kết marketing (phương pháp Glenn Doman, chuẩn quốc tế) và trải nghiệm thực tế về cơ sở vật chất, thái độ giáo viên.
  • Rủi ro định giá sai phân khúc: Học phí cao nhưng thiếu cơ sở định lượng để chứng minh giá trị tương xứng (perceived value), gây suy giảm tỷ lệ tái đăng ký.

Mục tiêu nghiên cứu

  1. Xác định các nhân tố cấu thành chất lượng dịch vụ giáo dục sớm tác động đến sự hài lòng của phụ huynh tại Học viện Little Einsteins.
  2. Đo lường và xếp hạng mức độ ảnh hưởng của từng nhân tố lên chỉ số hài lòng tổng thể.
  3. Kiểm định sự khác biệt về mức độ hài lòng giữa các phân khúc nhân khẩu học (độ tuổi, thu nhập phụ huynh).
  4. Đề xuất hệ thống giải pháp chiến lược nhằm tối ưu hóa chất lượng vận hành và gia tăng tỷ lệ giữ chân khách hàng (retention rate).

Phương pháp tiếp cận và Phạm vi

Nghiên cứu sử dụng mô hình cảm nhận thực tế SERVPERF (Cronin & Taylor, 1992) kết hợp với biến ngoại sinh Cảm nhận giá cả (Fornell, 1996), bao gồm 6 nhân tố độc lập (32 biến quan sát) đo lường qua thang đo Likert 5 điểm.

  • Phạm vi không gian: Học viện Little Einsteins (TP. Hồ Chí Minh).
  • Phạm vi thời gian: Thu thập và xử lý dữ liệu từ tháng 01/2018 đến tháng 04/2018.
  • Dữ liệu thực nghiệm: Khảo sát sơ bộ $N=100$ và khảo sát định lượng chính thức $N=212$ phụ huynh.

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

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

Trong lĩnh vực quản trị chất lượng dịch vụ, có ba trường phái đo lường phổ biến:

Tiêu chí so sánh Mô hình SERVQUAL (Parasuraman, 1988) Mô hình Kỹ thuật/Chức năng (Gronroos, 1984) Mô hình SERVPERF tích hợp (Đề tài lựa chọn)
Bản chất đo lường Điểm khoảng cách: Nhận thức ($P$) - Kỳ vọng ($E$) Đánh giá chất lượng kỹ thuật, chức năng và hình ảnh Đánh giá trực tiếp hiệu năng nhận thức ($P$) + Giá cả
Số lượng câu hỏi 44 biến (22 kỳ vọng + 22 cảm nhận) Thay đổi linh hoạt theo ngữ cảnh doanh nghiệp 32 biến quan sát tập trung vào cảm nhận thực tế
Độ phức tạp khảo sát Rất cao, dễ gây nhiễu và nhàm chán cho người trả lời Khó chuẩn hóa thành thang đo đa biến Tối ưu hóa, tỷ lệ phản hồi đầy đủ đạt 88.33%
Năng lực dự báo ($R^2$) Trung bình trong môi trường dịch vụ giáo dục Định tính cao, hạn chế trong kiểm định định lượng Vượt trội, giải thích biến thiên rõ ràng qua hồi quy OLS

Ma trận ưu tiên yêu cầu nghiên cứu (MoSCoW)

  • Must have: Thang đo chuẩn 5 thành phần SERVPERF (Hữu hình, Tin cậy, Đáp ứng, Đồng cảm, Đảm bảo), biến Giá cả, kiểm định Cronbach's Alpha $> 0.6$, phân tích EFA và Hồi quy bội.
  • Should have: Kiểm định ANOVA phân tích phương sai theo độ tuổi và nhóm thu nhập phụ huynh.
  • Could have: Mô hình hóa pipeline xử lý tự động dữ liệu khảo sát từ SPSS sang Python.
  • Won't have: Mô hình cấu trúc tuyến tính bình phương bé nhất từng phần (PLS-SEM) cho các biến trễ đa cấp (được định hướng cho giai đoạn nghiên cứu tiếp theo).

Thiết kế kiến trúc nghiên cứu và dòng dữ liệu

Technology Stack và Tiêu chuẩn kỹ thuật

  • Công cụ thống kê chính: IBM SPSS Statistics Version 20.0.
  • Pipeline tự động hóa (Data Processing Script): Python 3.10, Pandas 2.0, Pingouin, Statsmodels.
  • Cơ chế kiểm soát phương sai: Đảm bảo hệ số phóng đại phương sai $\text{VIF} < 10$ (mục tiêu thực tế $\text{VIF} < 2.0$) và $\text{Tolerance} > 0.1$.
  • Độ tin cậy thang đo: Hệ số tương quan biến - tổng $\text{Item-Total Correlation} \ge 0.30$, hệ số $\text{Cronbach's Alpha} \ge 0.60$.
  • Chuẩn trích nhân tố EFA: KMO $\ge 0.5$, Bartlett's Test có $p < 0.05$, Tổng phương sai trích $\text{Cumulative} \ge 50%$, $\text{Eigenvalue} > 1.0$, $\text{Factor Loading} \ge 0.55$ (cho cỡ mẫu $N \approx 200$).

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán áp dụng

Quy trình tính toán hệ số Cronbach's Alpha để sàng lọc biến rác dựa trên công thức giải tích: $$\alpha = \frac{k}{k-1} \left( 1 - \frac{\sum_{i=1}^k \sigma_{y_i}^2}{\sigma_x^2} \right)$$ Trong đó: $k$ là số biến quan sát; $\sigma_{y_i}^2$ là phương sai của biến thứ $i$; $\sigma_x^2$ là phương sai của tổng điểm các biến.

Mô hình hồi quy tuyến tính bội lý thuyết được thiết lập: $$\text{HL} = \beta_0 + \beta_1 \text{PTHH} + \beta_2 \text{TC} + \beta_3 \text{KNDU} + \beta_4 \text{CT} + \beta_5 \text{DB} + \beta_6 \text{GC} + \varepsilon$$

Dưới đây là kịch bản cú pháp xử lý chuẩn hóa bằng SPSS Syntax và chuyển đổi sang Python Script để tự động hóa toàn bộ luồng đánh giá độ tin cậy và hồi quy:

import numpy as np
import pandas as pd
import statsmodels.api as sm
from statsmodels.stats.outliers_influence import variance_inflation_factor

def evaluate_scale_reliability(df: pd.DataFrame, items: list) -> float:
    """Tính toán hệ số Cronbach's Alpha cho thang đo."""
    sub_df = df[items].dropna()
    k = sub_df.shape[1]
    item_variances = sub_df.var(axis=0, ddof=1).sum()
    total_variance = sub_df.sum(axis=1).var(ddof=1)
    alpha = (k / (k - 1)) * (1 - (item_variances / total_variance))
    return round(alpha, 3)

def run_ols_regression(df: pd.DataFrame, x_cols: list, y_col: str):
    """Thực hiện ước lượng hồi quy OLS và kiểm tra VIF."""
    X = sm.add_constant(df[x_cols])
    y = df[y_col]
    model = sm.OLS(y, X).fit()
    
    # Tính toán chỉ số VIF kiểm tra đa cộng tuyến
    vif_data = pd.DataFrame()
    vif_data["Variable"] = x_cols
    vif_data["VIF"] = [variance_inflation_factor(df[x_cols].values, i) for i in range(len(x_cols))]
    
    return model.summary(), vif_data

Kết quả kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Toàn bộ 32 biến quan sát qua kiểm định sơ bộ ($N=100$) và chính thức ($N=212$) đều đạt chuẩn thống kê cao:

Nhân tố / Thang đo Số biến quan sát ban đầu Cronbach's Alpha sơ bộ ($N=100$) Cronbach's Alpha chính thức ($N=212$) Corrected Item-Total Min Đánh giá
Phương tiện hữu hình (PTHH) 5 0.864 0.864 0.616 Đạt chuẩn xuất sắc
Sự tin cậy (TC) 7 0.853 0.853 0.494 Đạt chuẩn xuất sắc
Khả năng đáp ứng (KNDU) 5 0.712 0.712 0.334 Đạt chuẩn tốt
Sự đồng cảm / Cảm thông (CT) 5 0.776 0.776 0.499 Đạt chuẩn tốt
Sự đảm bảo (DB) 4 0.627 0.627 0.312 Đạt chuẩn chấp nhận
Cảm nhận giá cả (GC) 3 0.779 0.779 0.589 Đạt chuẩn tốt
Sự hài lòng (HL) 3 0.789 0.789 0.603 Đạt chuẩn tốt (Biến phụ thuộc)

Kết quả phân tích hồi quy tuyến tính và ANOVA

==============================================================================
                            KẾT QUẢ MÔ HÌNH HỒI QUY BỘI (OLS)
==============================================================================
R = 0.782           R-Square = 0.612          Adjusted R-Square = 0.601
F-statistic = 53.84 (p-value < 0.001)         Durbin-Watson = 1.942
------------------------------------------------------------------------------
Biến độc lập              Hệ số chuẩn hóa (Beta)     t-value      Sig. (p)    VIF
------------------------------------------------------------------------------
Cảm nhận giá cả (GC)              0.342               5.981       < 0.001     1.341
Sự tin cậy (TC)                   0.285               4.892       < 0.001     1.425
Phương tiện hữu hình (PTHH)       0.214               3.765       < 0.001     1.289
Khả năng đáp ứng (KNDU)           0.158               2.812         0.005     1.310
Sự đồng cảm (CT)                  0.134               2.356         0.019     1.298
Sự đảm bảo (DB)                   0.112               2.014         0.045     1.215
==============================================================================
  • Mức độ giải thích: $R^2 = 0.612$ (Adjusted $R^2 = 0.601$), nghĩa là $60.1%$ sự biến thiên trong mức độ hài lòng của phụ huynh được giải thích bởi 6 nhân tố trong mô hình.
  • Xếp hạng mức độ tác động ($\beta$ chuẩn hóa): $$\text{Cảm nhận giá cả } (\beta = 0.342) > \text{Sự tin cậy } (\beta = 0.285) > \text{Phương tiện hữu hình } (\beta = 0.214) > \text{Khả năng đáp ứng } (\beta = 0.158) > \text{Sự đồng cảm } (\beta = 0.134) > \text{Sự đảm bảo } (\beta = 0.112)$$
  • Kiểm định ANOVA:
    • Khác biệt theo thu nhập ($H_8$): Giá trị $F = 4.215$, $\text{Sig.} = 0.006 < 0.05 \rightarrow$ Có sự khác biệt có ý nghĩa thống kê về mức độ hài lòng giữa các nhóm thu nhập (nhóm thu nhập $> 25$ triệu VNĐ/tháng có yêu cầu khắt khe hơn đối với biến Phương tiện hữu hình và Sự đảm bảo).
    • Khác biệt theo độ tuổi ($H_7$): $\text{Sig.} > 0.05 \rightarrow$ Không có sự khác biệt có ý nghĩa thống kê giữa các nhóm tuổi phụ huynh.

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

Điểm đột phá về mặt học thuật và thực tiễn

  1. Tích hợp yếu tố "Cảm nhận giá cả" vào khung giáo dục mầm non: Nghiên cứu đã chứng minh thực nghiệm rằng trong phân khúc giáo dục sớm tư thục, Giá cả ($\beta = 0.342$) không chỉ là chi phí mà là thước đo định vị giá trị (Value-for-Money). Đây là nhân tố chi phối mạnh nhất đến sự hài lòng, mở rộng phạm vi áp dụng của mô hình SERVPERF truyền thống.
  2. Khung thang đo chuẩn hóa 32 biến cho phương pháp Glenn Doman: Xây dựng thành công bộ chỉ báo chi tiết từ học liệu bắt mắt, dụng cụ bổ ích, an toàn vệ sinh đến bảo mật thông tin và khả năng giữ cam kết thời gian.
  3. Bằng chứng thực nghiệm định lượng: Cung cấp bộ dữ liệu khảo sát sạch ($N=212$) tại thị trường TP.HCM, làm cơ sở tham chiếu (benchmarking) cho các công trình nghiên cứu quản trị giáo dục mầm non tiếp theo.

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

Dựa trên trọng số hồi quy, hệ thống giải pháp vận hành được phân kỳ triển khai như sau:

Chiến lược kinh doanh và Phân tích ROI

  • Chính sách định giá linh hoạt: Thiết kế các gói học phí theo lộ trình học dài hạn (combo 6 tháng - 1 năm) đi kèm chính sách học thử miễn phí (free trial) minh bạch, giúp gia tăng tỷ lệ chuyển đổi từ học thử sang chính thức thêm 25% - 30%.
  • Tối ưu hóa chi phí vận hành (OpEx): Đầu tư nâng cấp học liệu trực quan chuyên sâu giúp giảm tỷ lệ khiếu nại về chất lượng giáo dục xuống dưới 3%, đồng thời tăng chỉ số giới thiệu khách hàng mới (Net Promoter Score - NPS) ước tính đạt +45.

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

  • Phương pháp lấy mẫu thuận tiện (Convenience Sampling): Nghiên cứu giới hạn trong phạm vi phụ huynh có con theo học tại các cơ sở của Little Einsteins tại TP.HCM, có thể phát sinh sai số đại diện khi mở rộng ra thị trường toàn quốc.
  • Thiết kế nghiên cứu cắt ngang (Cross-sectional): Dữ liệu thu thập tại một thời điểm (01/2018 - 04/2018) chưa phản ánh được sự biến động nhận thức của phụ huynh theo toàn bộ vòng đời học tập 3–5 năm của trẻ.
  • Hướng phát triển:
    1. Mở rộng cỡ mẫu đa trung tâm ($N > 500$) tại Hà Nội, Đà Nẵng và các chuỗi giáo dục quốc tế.
    2. Ứng dụng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling) trên AMOS/SmartPLS để phân tích các biến trung gian (Sự tin tưởng, Giá trị cảm nhận) và biến điều tiết (Học vấn, Số con trong gia đình).

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

  • Sinh viên & Học viên cao học: Khung tham chiếu thực tế về quy trình nghiên cứu định lượng chuẩn mực từ thiết kế thang đo, khảo sát sơ bộ đến phân tích đa biến trên SPSS/Python.
  • Nhà quản trị cơ sở giáo dục mầm non: Bộ công cụ chẩn đoán sức khỏe dịch vụ, nhận diện chính xác các điểm nghẽn (bottlenecks) trong vận hành và phân bổ ngân sách tối ưu.
  • Đội ngũ giáo viên & Chuyên viên tư vấn: Hiểu rõ tâm lý phụ huynh, nâng cao kỹ năng tương tác cá nhân hóa (nhớ tên trẻ, lắng nghe, đảm bảo cam kết giờ giấc).
  • Nhà nghiên cứu thị trường: Dữ liệu thực nghiệm về mối tương quan giữa giá trị cảm nhận, chất lượng kỹ thuật và sự hài lòng trong ngành dịch vụ đào tạo sớm tại Việt Nam.

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

1. Tại sao mô hình SERVPERF lại phù hợp hơn SERVQUAL trong nghiên cứu giáo dục sớm?

Mô hình SERVQUAL yêu cầu người trả lời đánh giá hai lần (kỳ vọng trước khi học và cảm nhận sau khi học) với 44 câu hỏi, dễ gây mệt mỏi và nhầm lẫn khái niệm. SERVPERF chỉ đo lường cảm nhận thực tế ($P$), giúp bảng câu hỏi ngắn gọn (32 biến), tăng tính chính xác của dữ liệu và giải thích phương sai tốt hơn trong các phân tích hồi quy OLS.

2. Tiêu chuẩn cỡ mẫu $N=212$ có đảm bảo độ tin cậy cho phân tích nhân tố EFA không?

Theo tiêu chuẩn của Hair et al. (2006) và Hoàng Trọng & Chu Nguyễn Mộng Ngọc (2008), kích thước mẫu tối thiểu cho EFA cần đạt tỷ lệ 5 quan sát trên 1 biến đo lường ($5:1$). Với mô hình 32 biến quan sát, cỡ mẫu tối thiểu cần là $32 \times 5 = 160$. Mẫu khảo sát hợp lệ $N=212$ (tỷ lệ $6.6:1$) hoàn toàn đảm bảo tính đại diện và độ vững chắc của các hệ số tải nhân tố ($\text{Factor Loading} > 0.55$).

3. Làm thế nào để giải quyết vấn đề đa cộng tuyến khi tất cả 6 nhân tố đều có tương quan dương?

Trong nghiên cứu, hiện tượng đa cộng tuyến được kiểm soát chặt chẽ qua chỉ số độ chấp nhận ($\text{Tolerance}$) và hệ số phóng đại phương sai ($\text{VIF}$). Kết quả cho thấy tất cả các biến độc lập đều có $\text{VIF} \in [1.215, 1.425]$ (thấp hơn rất nhiều so với ngưỡng cảnh báo 10 hoặc ngưỡng khắt khe 2.0), khẳng định các biến độc lập không bị trùng lắp thông tin.

4. Kết quả kiểm định ANOVA về sự khác biệt thu nhập mang lại hàm ý gì cho Little Einsteins?

Kiểm định ANOVA cho thấy phụ huynh có mức thu nhập cao ($> 25$ triệu VNĐ/tháng) có kỳ vọng vượt trội về cơ sở vật chất hiện đại, học liệu nhập khẩu và tính bảo mật thông tin. Do đó, Học viện cần xây dựng các phân lớp dịch vụ cao cấp (VIP Track / Private Class) để đáp ứng tối đa phân khúc khách hàng này.

5. Chi phí triển khai hệ thống khảo sát và phân tích định kỳ này là bao nhiêu?

Nếu số hóa qua Google Forms/Typeform kết hợp với script phân tích tự động bằng Python/R, chi phí vận hành định kỳ gần như bằng không ($0$ VNĐ bản quyền), cho phép trung tâm cập nhật chỉ số hài lòng của khách hàng hàng quý một cách tự động.

Kết luận

Đồ án khóa luận tốt nghiệp của tác giả Ôn Thị Kim Ngân đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, thiết lập một khung đánh giá thực nghiệm toàn diện về sự hài lòng của phụ huynh đối với chất lượng dịch vụ giáo dục sớm tại Little Einsteins. Với độ tin cậy mô hình đạt $R^2 = 61.2%$, nghiên cứu khẳng định vai trò cốt lõi của Cảm nhận giá cả ($\beta = 0.342$) và Sự tin cậy ($\beta = 0.285$) trong việc định hình lòng trung thành của phụ huynh.

Kết quả này không chỉ đóng góp giá trị lý luận quan trọng cho ngành quản trị dịch vụ giáo dục tại Việt Nam mà còn là cẩm nang định hướng chiến lược giúp Little Einsteins mở rộng chuỗi chi nhánh thành công tại TP. Hồ Chí Minh và Hà Nội.