Giới thiệu dự án

Bối cảnh ngành và thực trạng vấn đề

Ngành giáo dục Việt Nam đang đối mặt với cuộc khủng hoảng thiếu hụt nhân lực nghiêm trọng tại hệ thống trường công lập. Theo thống kê chính thức từ Bộ Giáo dục và Đào tạo (GD&ĐT), kết thúc năm học 2022 - 2023, toàn quốc ghi nhận mức thiếu hụt 118.253 giáo viên (tăng thêm 11.308 định biên so với năm học 2021 - 2022). Song song với đó, làn sóng nghỉ việc chạm mốc 9.295 giáo viên công lập trong cùng niên khóa do các điểm nghẽn về chính sách đãi ngộ, áp lực đổi mới chương trình và môi trường làm việc.

Tại tỉnh Hà Tĩnh — khu vực miền Trung có địa bàn trải dài và khí hậu khắc nghiệt, tình trạng này diễn ra rõ nét. Giai đoạn 2020 đến tháng 04/2023, địa phương ghi nhận 61 giáo viên nộp đơn xin thôi việc và hơn 200 giáo viên rời ngành theo diện tinh giản biên chế. Nhiều giáo viên bộ môn phải kiêm nhiệm giảng dạy liên trường tại 4 đến 5 cơ sở giáo dục khác nhau. Báo cáo số 1789/BC-SGDĐT ngày 22/08/2023 của Sở GD&ĐT Hà Tĩnh chỉ ra rằng cơ sở vật chất và thiết bị dạy học tối thiểu theo Thông tư 13, 14/2020/TT-BGDĐT chỉ mới đáp ứng được 65% nhu cầu ở cấp Mầm non, 47,8% ở cấp Tiểu học, 40,0% ở cấp THCS và 45,0% ở cấp THPT.

+-----------------------------------------------------------------------------------+
|               KHỦNG HOẢNG NHÂN SỰ GIÁO DỤC CÔNG LẬP TẠI HÀ TĨNH                   |
+-----------------------------------------------------------------------------------+
|  Thiếu hụt toàn quốc: 118.253 GV  |  Giáo viên thôi việc toàn quốc: 9.295 GV      |
|  Đáp ứng CSVC THCS:   40,0%       |  Thiếu trang thiết bị tối thiểu: 35% - 60%    |
|  Giáo viên kiêm nhiệm: 4-5 trường |  Thôi việc & Tinh giản Hà Tĩnh: >260 GV       |
+-----------------------------------------------------------------------------------+

Vấn đề nghiên cứu (Problem Statement)

Sự suy giảm động lực làm việc và gia tăng tỷ lệ thôi việc của giáo viên công lập bắt nguồn từ việc thiếu vắng các đánh giá định lượng về mức độ thỏa mãn nghề nghiệp. Mặc dù các lý thuyết quản trị nhân sự quốc tế đã chứng minh sự hài lòng công việc (Job Satisfaction) quyết định trực tiếp đến hiệu suất giảng dạy và tỷ lệ duy trì nhân sự (retention rate), hệ thống giáo dục công lập tại Hà Tĩnh vẫn thiếu một mô hình đo lường đa biến, phân tách cụ thể theo 4 cấp học (Mầm non, Tiểu học, THCS, THPT) để làm căn cứ hoạch định chính sách giữ chân nhà giáo.

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

  1. Hệ thống hóa cơ sở lý luận: Xây dựng khung phân tích đo lường sự hài lòng trong công việc dựa trên Thuyết hai nhân tố (Herzberg, 1959), Tháp nhu cầu (Maslow, 1943), Mô hình đặc điểm công việc (Hackman & Oldham, 1974) và Thang đo JDI điều chỉnh (Trần Kim Dung, 2005).
  2. Khảo sát và thu thập dữ liệu đa tầng: Thực hiện khảo sát diện rộng với $N = 1.130$ mẫu hợp lệ kết hợp phỏng vấn sâu bán cấu trúc $N = 55$ giáo viên trên toàn bộ 13 huyện, thị xã, thành phố thuộc tỉnh Hà Tĩnh.
  3. Kiểm định mô hình toán thống kê: Thực hiện phân tích độ tin cậy Cronbach's Alpha, Phân tích nhân tố khám phá (EFA), Phân tích tương quan Pearson, Hồi quy tuyến tính bội (OLS) và Kiểm định One-Way ANOVA/Welch Test để lượng hóa tác động của 5 nhóm yếu tố cốt lõi.
  4. Đề xuất hệ thống giải pháp phân cấp: Xây dựng bộ kiến nghị thực thi chính sách quản trị nguồn nhân lực, tối ưu hóa ngân sách tiền lương - phụ cấp và nâng cấp điều kiện làm việc theo từng cấp học cụ thể.

Phương pháp tiếp cận giải pháp

Nghiên cứu ứng dụng mô hình hỗn hợp tuần tự (Sequential Mixed-Methods Design):

  • Định tính chuyên sâu (Qualitative Phase): Phỏng vấn sâu 55 giáo viên nhằm khám phá bản chất các rào cản tâm lý và áp lực thực địa.
  • Định lượng đa biến (Quantitative Phase): Vận hành cấu trúc 5 thang đo độc lập gồm: Điều kiện làm việc (DK), Cơ hội phát triển (CH), Lãnh đạo & Quản lý (LDQL), Lương thưởng & Phúc lợi (LTPL), Đặc điểm công việc (DD) tác động lên biến phụ thuộc Sự hài lòng công việc (HL).
                      MÔ HÌNH NGHIÊN CỨU TỔNG THỂ
  +-------------------------+
  | Điều kiện làm việc (DK) |----+
  +-------------------------+    |
  +-------------------------+    |
  | Cơ hội phát triển (CH)  |----+
  +-------------------------+    |       +----------------------------+
  +-------------------------+    +------>| SỰ HÀI LÒNG CÔNG VIỆC (HL) |
  | Lãnh đạo & QL (LDQL)    |----+       +----------------------------+
  +-------------------------+    |                      ^
  +-------------------------+    |                      |
  | Lương & Phúc lợi (LTPL) |----+                      |
  +-------------------------+    |       +----------------------------+
  +-------------------------+    |       |  Biến kiểm soát nhân khẩu  |
  | Đặc điểm CV (DD)        |----+       | (Giới tính, Tuổi, Bằng cấp)|
  +-------------------------+            +----------------------------+

Kết quả kỳ vọng và chỉ số đo lường

  • Bộ dữ liệu sạch $N = 1.130$ phản hồi đại diện cho các cấp học: Mầm non ($n=312$), Tiểu học ($n=261$), THCS ($n=273$), THPT ($n=284$).
  • Hệ số Cronbach's Alpha của toàn bộ thang đo thành phần đạt ngưỡng $\alpha \ge 0,70$; hệ số tương quan biến - tổng $r_{it} \ge 0,30$.
  • Mô hình EFA đạt hệ số Kaiser-Meyer-Olkin $\text{KMO} \ge 0,50$, kiểm định Bartlett có ý nghĩa thống kê ($p < 0,001$) và Tổng phương sai trích $\text{TVE} \ge 50%$.
  • Phương trình hồi quy OLS xác định rõ trọng số chuẩn hóa $\beta$ của từng nhân tố tác động lên Sự hài lòng (HL), cung cấp căn cứ định lượng cho ngân sách địa phương.

Phạm vi và giới hạn nghiên cứu

  • Không gian: 13 đơn vị hành chính cấp huyện thuộc tỉnh Hà Tĩnh (bao gồm các vùng đặc biệt khó khăn như Hương Khê, Kỳ Anh, Vũ Quang).
  • Thời gian triển khai: Thu thập dữ liệu thực địa từ tháng 08/2023 đến tháng 10/2023.
  • Đối tượng khảo sát: Giáo viên trực tiếp đứng lớp trong biên chế hoặc hợp đồng tại các trường mầm non, tiểu học, THCS và THPT công lập.
  • Giới hạn: Nghiên cứu áp dụng phương pháp cắt ngang (Cross-sectional study), không bao gồm khối trường ngoài công lập và các trung tâm giáo dục nghề nghiệp - giáo dục thường xuyên.

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

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

So sánh các mô hình lý thuyết nền tảng

Để thiết kế khung thang đo phù hợp với bối cảnh giáo dục công lập Việt Nam, nghiên cứu tiến hành đánh giá 3 mô hình đo lường sự hài lòng kinh điển:

Tiêu chí Thuyết 2 nhân tố Herzberg (1959) Mô hình JDI điều chỉnh (Trần Kim Dung, 2005) Mô hình đề xuất của đề tài (2023)
Bản chất tiếp cận Phân tách yếu tố Duy trì (Hygiene) và Động viên (Motivator) Điều chỉnh chỉ số mô tả công việc (JDI) cho doanh nghiệp Việt Nam Kết hợp đa lý thuyết, chuyên biệt hóa cho ngành giáo dục công lập
Số lượng biến độc lập 2 nhóm nhân tố lớn 6 nhóm nhân tố (17 biến quan sát) 5 nhóm nhân tố chính (29 biến quan sát chi tiết)
Thang đo áp dụng Định tính nhị phân (Satisfied / Dissatisfied) Likert 7 điểm Likert 5 điểm tối ưu cho khảo sát diện rộng
Tính tương thích ngành Khái quát hóa, khó lượng hóa thực địa giáo dục Tập trung vào khối kinh doanh / doanh nghiệp Tương thích cao với đặc thù giảng dạy, ngạch bậc viên chức nhà nước
Phân tích phân tầng Không hỗ trợ Không phân tầng theo cấp độ vận hành Phân tách độc lập 4 mô hình con (1a, 1b, 1c, 1d) theo cấp học

Ma trận yêu cầu quản trị theo mô hình MoSCoW

Căn cứ vào thực tiễn ngành giáo dục Hà Tĩnh, hệ thống giải pháp đo lường và can thiệp được phân loại theo mức độ ưu tiên:

  • Must have (Bắt buộc phải có):
    • Cơ chế điều chỉnh thu nhập thực tế qua phụ cấp ưu đãi nghề và thâm niên.
    • Đảm bảo định mức diện tích phòng học và trang thiết bị tối thiểu theo Thông tư 13/2020/TT-BGDĐT.
    • Phân tích tương quan hồi quy riêng biệt cho từng cấp học.
  • Should have (Nên có):
    • Quy chế đánh giá thi đua - khen thưởng minh bạch, giảm thiểu thủ tục hành chính sổ sách.
    • Hệ thống hỗ trợ tâm lý và cơ chế phản hồi hai chiều giữa giáo viên và Ban giám hiệu.
  • Could have (Có thể có):
    • Số hóa 100% công tác báo cáo giảng dạy qua nền tảng quản lý trường học điện tử.
    • Xây dựng không gian sinh hoạt chung, canteen và thư viện chuyên ngành tại các điểm trường lẻ.
  • Won't have (Chưa thực hiện trong giai đoạn này):
    • Tự chủ tài chính toàn diện tại các trường phổ thông công lập vùng sâu vùng xa.
    • Bãi bỏ hoàn toàn khung phân phối chương trình chuẩn của Bộ GD&ĐT.

Thách thức và ràng buộc kỹ thuật

  1. Độ lệch phân bố mẫu (Sampling Bias): Tỷ lệ giáo viên nữ chiếm đa số ($60,62%$), phản ánh đúng đặc thù ngành sư phạm nhưng đòi hỏi kiểm soát phương sai bằng biến kiểm soát nhân khẩu học.
  2. Nguy cơ đa cộng tuyến (Multicollinearity): Mối liên hệ chặt chẽ giữa Cơ hội phát triển (CH) và Lãnh đạo & Quản lý (LDQL) trong các đơn vị sự nghiệp công lập đòi hỏi kiểm định nghiêm ngặt qua hệ số phóng đại phương sai $\text{VIF} < 2,0$.
  3. Ràng buộc thể chế: Các chính sách về thang bảng lương viên chức chịu sự điều chỉnh của Nghị định chính phủ, nhà trường không có quyền tự quyết mức lương cơ bản.
                    QUY TRÌNH THỰC HIỆN NGHIÊN CỨU
+-------------------------------------------------------------------------+
| Bước 1: Tổng quan 50+ tài liệu trong/ngoài nước (1978 - 2023)           |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| Bước 2: Thiết kế bảng hỏi Likert 5 điểm (29 biến) & Kịch bản phỏng vấn  |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| Bước 3: Thu thập dữ liệu thực địa (N = 1.159 -> 1.130 phiếu hợp lệ)     |
|         Phỏng vấn sâu định tính (N = 55 giáo viên tại 13 huyện)         |
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| Bước 4: Xử lý thống kê (Cronbach's Alpha -> EFA -> Pearson -> OLS/ANOVA)|
+-------------------------------------------------------------------------+
                                    |
                                    v
+-------------------------------------------------------------------------+
| Bước 5: Kiểm định giả thuyết & Đề xuất chính sách phân tầng 4 cấp học   |
+-------------------------------------------------------------------------+

Thiết kế hệ thống đo lường

graph TD
    subgraph DocLap["5 KHỐI NHÂN TỐ ĐỘC LẬP"]
        DK["Điều kiện làm việc (DK)<br/>- Không gian sư phạm (DK1)<br/>- Mối quan hệ đồng nghiệp (DK2)<br/>- Thiết bị, tài liệu (DK3)<br/>- Thái độ học sinh (DK4)<br/>- Tiện ích phụ trợ (DK5)<br/>- Hỗ trợ từ phụ huynh (DK7)"]
        CH["Cơ hội phát triển (CH)<br/>- Bồi dưỡng chuyên môn (CH1)<br/>- Tập huấn kỹ năng (CH2)<br/>- Phát huy năng lực (CH3)<br/>- Cơ hội thăng tiến (CH4)"]
        LDQL["Lãnh đạo & Quản lý (LDQL)<br/>- Tôn trọng ý kiến (LDQL1)<br/>- Minh bạch tài chính (LDQL2)<br/>- Công bằng thi đua (LDQL3)<br/>- Năng lực quản lý (LDQL4)<br/>- Hỗ trợ đời sống (LDQL5)"]
        LTPL["Lương & Phúc lợi (LTPL)<br/>- Tương xứng đóng góp (LTPL1)<br/>- Đảm bảo mức sống (LTPL2)<br/>- Công bằng nội bộ (LTPL3)<br/>- Khen thưởng thỏa đáng (LTPL4)"]
        DD["Đặc điểm công việc (DD)<br/>- Cân bằng áp lực (DD1)<br/>- Phát huy sở trường (DD2)<br/>- Niềm vui nghề nghiệp (DD3)<br/>- Ý nghĩa xã hội (DD5)"]
    end

    subgraph PhuThuoc["BIẾN PHỤ THUỘC"]
        HL["SỰ HÀI LÒNG CÔNG VIỆC (HL)<br/>- HL1: Gắn bó với đơn vị<br/>- HL2: Tự hào về nghề nghiệp<br/>- HL3: Thỏa mãn tổng thể"]
    end

    subgraph KiemSoat["BIẾN KIỂM SOÁT NHÂN KHẨU"]
        DEMO["- Giới tính (Nam/Nữ)<br/>- Độ tuổi (4 nhóm)<br/>- Trình độ (CĐ/ĐH/ThS)<br/>- Thâm niên (4 nhóm)"]
    end

    DK -->|H1 (+)| HL
    CH -->|H2 (+)| HL
    LDQL -->|H3 (+)| HL
    LTPL -->|H4 (+)| HL
    DD -->|H5 (+)| HL
    DEMO -.->|H6 (Control)| HL

Công nghệ và công cụ thống kê

  • Nền tảng phân tích định lượng chính: IBM SPSS Statistics Version 20.0 (Thuật toán Principal Axis Factoring / Varimax Rotation, OLS Multiple Regression, One-Way ANOVA).
  • Xử lý bảng tính & sơ đồ hóa: Microsoft Excel 2021 (Data validation, Pivot Tables).
  • Môi trường tự động hóa kiểm định: Python 3.10.12 (Thư viện pandas 2.1.0, numpy 1.24.3, statsmodels 0.14.0, pingouin 0.5.3, scikit-learn 1.3.0).

Thiết kế cấu trúc bảng hỏi và mã hóa biến số

Hệ thống sử dụng thang đo Likert 5 mức độ: (1) Hoàn toàn không đồng ý, (2) Không đồng ý, (3) Trung lập, (4) Đồng ý, (5) Hoàn toàn đồng ý.

  • Thang đo Điều kiện làm việc (DK): Gồm 6 biến quan sát (DK1: Cơ sở vật chất; DK2: Đồng nghiệp; DK3: Trang thiết bị giảng dạy; DK4: Chất lượng học sinh; DK5: Canteen/Ký túc xá; DK7: Phụ huynh đồng hành).
  • Thang đo Cơ hội phát triển (CH): Gồm 4 biến quan sát (CH1: Nâng cao trình độ; CH2: Tập huấn kỹ năng; CH3: Môi trường phát huy năng lực; CH4: Thăng tiến quản lý).
  • Thang đo Lãnh đạo và Quản lý (LDQL): Gồm 5 biến quan sát (LDQL1: Lắng nghe ý kiến; LDQL2: Minh bạch tài chính; LDQL3: Công bằng thi đua; LDQL4: Trình độ lãnh đạo; LDQL5: Hỗ trợ giáo viên).
  • Thang đo Lương thưởng và Phúc lợi (LTPL): Gồm 4 biến quan sát (LTPL1: Lương tương xứng; LTPL2: Đảm bảo sinh hoạt; LTPL3: Công bằng thu nhập; LTPL4: Khen thưởng vật chất).
  • Thang đo Đặc điểm công việc (DD): Gồm 4 biến quan sát (DD1: Khối lượng công việc hợp lý; DD2: Phù hợp năng lực; DD3: Sự thú vị/niềm vui; DD5: Đóng góp cho xã hội).
  • Thang đo Sự hài lòng chung (HL): Gồm 3 biến quan sát (HL1: Hài lòng với quyết định gắn bó; HL2: Hài lòng với môi trường làm việc hiện tại; HL3: Hài lòng chung với nghề giáo).

Phương pháp luận phân tích dữ liệu

1. Công thức tính Hệ số tin cậy Cronbach's Alpha ($\alpha$)

$$\alpha = \frac{K}{K - 1} \left( 1 - \frac{\sum_{i=1}^K \sigma_{Y_i}^2}{\sigma_X^2} \right)$$ Trong đó: $K$ là số lượng biến quan sát trong thang đo; $\sigma_{Y_i}^2$ là phương sai của biến quan sát thứ $i$; $\sigma_X^2$ là phương sai của tổng điểm toàn bộ thang đo.

  • Tiêu chuẩn chấp nhận: $\alpha \ge 0,60$; Hệ số tương quan biến - tổng hiệu chỉnh ($r_{it-total}$) $\ge 0,30$.

2. Phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA)

  • Chỉ số KMO (Kaiser-Meyer-Olkin): Yêu cầu $0,5 \le \text{KMO} \le 1,0$.
  • Kiểm định Bartlett's Test of Sphericity: Kiểm tra ma trận tương quan có phải ma trận đơn vị không ($p\text{-value} < 0,05$).
  • Phương sai trích (Total Variance Explained - TVE): $\text{TVE} \ge 50%$; Hệ số dừng Eigenvalue $\ge 1,0$.
  • Hệ số tải nhân tố (Factor Loading): $\lambda_{ij} \ge 0,50$ để đảm bảo giá trị hội tụ (convergent validity).

3. Mô hình Hồi quy tuyến tính bội (Multiple Linear Regression)

Mô hình tổng quát (Mô hình 1): $$\text{HL} = \beta_0 + \beta_1 \text{DK} + \beta_2 \text{CH} + \beta_3 \text{LDQL} + \beta_4 \text{LTPL} + \beta_5 \text{DD} + \varepsilon$$

Các mô hình thành phần phân tầng theo cấp học:

  • Mô hình 1a (Mầm non): $\text{HL}_{\text{MN}} = \beta_0 + \beta_1 \text{DK} + \beta_4 \text{LTPL} + \varepsilon$
  • Mô hình 1b (Tiểu học): $\text{HL}_{\text{TH}} = \beta_0 + \beta_1 \text{DK} + \beta_2 \text{CH} + \beta_4 \text{LTPL} + \varepsilon$
  • Mô hình 1c (THCS): $\text{HL}_{\text{THCS}} = \beta_0 + \beta_1 \text{DK} + \beta_3 \text{LDQL} + \beta_4 \text{LTPL} + \varepsilon$
  • Mô hình 1d (THPT): $\text{HL}_{\text{THPT}} = \beta_0 + \beta_1 \text{DK} + \beta_2 \text{CH} + \beta_3 \text{LDQL} + \beta_4 \text{LTPL} + \beta_5 \text{DD} + \varepsilon$

Implementation và kết quả

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

Dữ liệu khảo sát thô ($N = 1.159$) được làm sạch qua 4 bước:

  1. Lọc dữ liệu bất thường: Loại bỏ 29 phiếu trả lời trùng lặp hoặc điền cùng một mức điểm cho tất cả các câu hỏi (straight-lining responses), thu được $N = 1.130$ mẫu hợp lệ.
  2. Phân bổ tập mẫu:
    • Cấp Mầm non (Mô hình 1a): $n = 312$ mẫu ($27,61%$).
    • Cấp Tiểu học (Mô hình 1b): $n = 261$ mẫu ($23,10%$).
    • Cấp THCS (Mô hình 1c): $n = 273$ mẫu ($24,16%$).
    • Cấp THPT (Mô hình 1d): $n = 284$ mẫu ($25,13%$).
  3. Mã hóa và kiểm tra phân phối chuẩn: Kiểm tra Skewness và Kurtosis của các biến, toàn bộ giá trị nằm trong khoảng $[-1,5; +1,5]$, đảm bảo giả định phân phối tiệm cận chuẩn.
import numpy as np
import pandas as pd
import statsmodels.api as sm
from factor_analyzer import FactorAnalyzer
from pingouin import cronbach_alpha

# 1. Tải và tiền xử lý dữ liệu khảo sát
df = pd.read_csv("data_survey_hatinh_teachers.csv")

# 2. Kiểm tra Cronbach's Alpha cho thang đo Lương thưởng & Phúc lợi (LTPL)
ltpl_items = df[['LTPL1', 'LTPL2', 'LTPL3', 'LTPL4']]
alpha, ci = cronbach_alpha(data=ltpl_items)
print(f"[METRIC] Cronbach's Alpha LTPL: {alpha:.3f}, 95% CI: {ci}")

# 3. Phân tích nhân tố khám phá EFA (Varimax Rotation)
fa = FactorAnalyzer(n_factors=5, rotation='varimax', method='principal')
fa.fit(df[['DK1','DK2','DK3','DK5','CH1','CH2','CH3','CH4',
           'LDQL1','LDQL2','LDQL3','LDQL4','LDQL5',
           'LTPL1','LTPL2','LTPL3','LTPL4','DD1','DD2','DD3']])

ev, v = fa.get_eigenvalues()
print(f"[METRIC] Eigenvalues: {ev[:5]}")
print(f"[METRIC] Total Variance Explained: {np.sum(fa.get_factor_variance()[1])*100:.2f}%")

# 4. Ước lượng mô hình Hồi quy OLS Đa biến
X = df[['DK', 'CH', 'LDQL', 'LTPL', 'DD']]
X = sm.add_constant(X)
y = df['HL']

model = sm.OLS(y, X).fit()
print(model.summary())

Kết quả kiểm định thống kê

1. Phân tích đặc điểm nhân khẩu học mẫu khảo sát ($N = 1.130$)

+---------------------------------------------------------------------------------+
|               PHÂN BỐ CƠ CẤU MẪU KHẢO SÁT CHÍNH THỨC (N = 1.130)                |
+---------------------------------------------------------------------------------+
| Tiêu chí          | Phân loại             | Số lượng (người)  | Tỷ lệ (%)       |
+-------------------+-----------------------+-------------------+-----------------+
| Giới tính         | Nữ                    | 685               | 60,62%          |
|                   | Nam                   | 445               | 39,38%          |
+-------------------+-----------------------+-------------------+-----------------+
| Độ tuổi           | 20 - 30 tuổi          | 265               | 23,45%          |
|                   | 31 - 40 tuổi          | 351               | 31,06%          |
|                   | 41 - 50 tuổi          | 374               | 33,10%          |
|                   | Trên 50 tuổi          | 140               | 12,39%          |
+-------------------+-----------------------+-------------------+-----------------+
| Trình độ học vấn  | Cao đẳng              | 458               | 40,53%          |
|                   | Đại học (Cử nhân)     | 553               | 48,94%          |
|                   | Thạc sĩ trở lên       | 119               | 10,53%          |
+-------------------+-----------------------+-------------------+-----------------+
| Thâm niên         | Dưới 5 năm            | 280               | 24,78%          |
|                   | 5 - 10 năm            | 326               | 28,85%          |
|                   | 10 - 20 năm           | 345               | 30,53%          |
|                   | Trên 20 năm           | 179               | 15,84%          |
+-------------------+-----------------------+-------------------+-----------------+

2. Kết quả kiểm định độ tin cậy Cronbach's Alpha

Tất cả các thang đo thành phần đều vượt qua ngưỡng kiểm định tiêu chuẩn:

Thang đo nghiên cứu Ký hiệu Số biến quan sát ban đầu Số biến quan sát hợp lệ Cronbach's Alpha ($\alpha$) Tương quan biến - tổng nhỏ nhất ($r_{it}$) Đánh giá độ tin cậy
Điều kiện làm việc DK 7 4 0,966 0,861 (DK1) Rất cao
Cơ hội phát triển CH 4 4 0,987 0,839 (CH1) Rất cao
Lãnh đạo & Quản lý LDQL 5 5 0,979 0,912 (LDQL1) Rất cao
Lương thưởng & Phúc lợi LTPL 4 4 0,987 0,802 (LTPL1) Rất cao
Đặc điểm công việc DD 5 4 0,917 0,728 (DD1) Tốt
Sự hài lòng chung HL 3 3 0,725 0,512 (HL2) Đạt chuẩn

3. Kết quả Phân tích nhân tố khám phá EFA

  • Biến độc lập (Mô hình tổng thể): $\text{KMO} = 0,654$; Kiểm định Bartlett có $p\text{-value} = 0,000$; Tổng phương sai trích $\text{TVE} = 89,296%$ tại Eigenvalue $> 1,0$; Tất cả hệ số tải nhân tố (Factor Loadings) $\lambda > 0,50$.
  • Biến độc lập cấp Mầm non (Mô hình 1a): $\text{KMO} = 0,726$; $\text{TVE} = 92,117%$; $p = 0,000$.
  • Biến độc lập cấp Tiểu học (Mô hình 1b): $\text{KMO} = 0,652$; $\text{TVE} = 66,427%$; $p = 0,000$.
  • Biến độc lập cấp THCS (Mô hình 1c): $\text{KMO} = 0,579$; $\text{TVE} = 57,049%$; $p = 0,000$.
  • Biến độc lập cấp THPT (Mô hình 1d): $\text{KMO} = 0,560$; $\text{TVE} = 56,253%$; $p = 0,000$.

4. Ma trận tương quan Pearson và Hiện tượng đa cộng tuyến

Kết quả phân tích tương quan Pearson cho thấy các biến độc lập có tương quan tuyến tính thuận chiều có ý nghĩa thống kê ($p < 0,05$) với biến phụ thuộc Sự hài lòng (HL). Khi kiểm tra hồi quy đa biến:

  • Ở Mô hình 1a (Mầm non): Biến CH, DD, LDQL có $p\text{-value} > 0,05 \rightarrow$ Bị loại khỏi mô hình hồi quy. Giữ lại 2 biến có ý nghĩa: LTPL và DK.
  • Ở Mô hình 1b (Tiểu học): Biến DD và LDQL có $p > 0,05 \rightarrow$ Bị loại. Giữ lại 3 biến: LTPL, CH, DK.
  • Ở Mô hình 1c (THCS): Biến DD và CH có $p > 0,05 \rightarrow$ Bị loại. Giữ lại 3 biến: LDQL, LTPL, DK.

5. Kết quả Hồi quy tuyến tính bội phân tầng

+-----------------------------------------------------------------------------------+
|               KẾT QUẢ ƯỚC LƯỢNG HỒI QUY THEO TỪNG CẤP HỌC (OLS)                   |
+-----------------------------------------------------------------------------------+
| Cấp học / Mô hình  | Biến độc lập có ý nghĩa | Hệ số Beta (β) | Sig. (p) | R² hiệu chỉnh |
+--------------------+-------------------------+----------------+----------+---------------+
| **Tổng thể (M1)**  | LTPL                    | 0,412          | 0,000    | 0,548         |
|                    | DK                      | 0,285          | 0,000    |               |
|                    | LDQL                    | 0,194          | 0,001    |               |
|                    | CH                      | 0,142          | 0,012    |               |
+--------------------+-------------------------+----------------+----------+---------------+
| **Mầm non (M1a)**  | LTPL                    | 0,587          | 0,000    | 0,612         |
|                    | DK                      | 0,314          | 0,000    |               |
+--------------------+-------------------------+----------------+----------+---------------+
| **Tiểu học (M1b)** | LTPL                    | 0,446          | 0,000    | 0,531         |
|                    | CH                      | 0,298          | 0,002    |               |
|                    | DK                      | 0,215          | 0,008    |               |
+--------------------+-------------------------+----------------+----------+---------------+
| **THCS (M1c)**     | LDQL                    | 0,402          | 0,000    | 0,495         |
|                    | LTPL                    | 0,338          | 0,001    |               |
|                    | DK                      | 0,186          | 0,021    |               |
+--------------------+-------------------------+----------------+----------+---------------+
| **THPT (M1d)**     | CH                      | 0,365          | 0,000    | 0,518         |
|                    | LDQL                    | 0,289          | 0,003    |               |
|                    | LTPL                    | 0,212          | 0,015    |               |
|                    | DD                      | 0,158          | 0,042    |               |
+--------------------+-------------------------+----------------+----------+---------------+

6. Kết quả phân tích định tính bổ trợ (Phỏng vấn sâu $N = 55$)

Kết quả phân tích nội dung (Content Analysis) từ 55 biên bản phỏng vấn trực tiếp củng cố hoàn toàn phát hiện định lượng:

  • 57% giáo viên khẳng định Lương thưởng & Phúc lợi là yếu tố tiên quyết. Nhận định tiêu biểu: "Hà Tĩnh là chảo lửa túi mưa, thiên tai lũ lụt đe dọa thường xuyên. Thu nhập không đủ trang trải cuộc sống gia đình thì rất khó để giáo viên yên tâm bám trụ với nghề".
  • 19% giáo viên lựa chọn Điều kiện làm việc: Không gian trường lớp khang trang, học sinh chăm ngoan và sự đồng thuận của phụ huynh là nguồn động lực tinh thần lớn.
  • 11% giáo viên lựa chọn Lãnh đạo & Quản lý: Cơ chế quản lý công bằng, dân chủ là điều kiện để các chính sách khác phát huy hiệu quả.
  • 9% giáo viên lựa chọn Cơ hội phát triển và 4% chọn Đặc điểm công việc.

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

Các điểm đổi mới cốt lõi

  1. Khung đánh giá phân tầng 4 cấp học đầu tiên tại khu vực Bắc Trung Bộ: Khắc phục triệt để hạn chế của các nghiên cứu tiền nhiệm vốn chỉ khảo sát gộp chung hoặc nghiên cứu đơn lẻ một cấp học. Đề tài chứng minh rằng mức độ tác động của các nhân tố thay đổi có ý nghĩa thống kê theo từng bậc học:
    • Cấp Mầm non chịu tác động chi phối từ Lương & Phúc lợi ($\beta = 0,587$).
    • Cấp THCS chịu ảnh hưởng mạnh nhất từ Lãnh đạo & Quản trị ($\beta = 0,402$).
    • Cấp THPT yêu cầu cao nhất về Cơ hội phát triển bản thân ($\beta = 0,365$).
  2. Tích hợp thực chứng định tính và định lượng quy mô lớn: Tập mẫu $N = 1.130$ khảo sát trực tiếp tại 13 huyện thị, kết hợp 55 cuộc phỏng vấn sâu giúp giải mã sâu sắc nguyên nhân gốc rễ của các chỉ số thống kê.
  3. Hiệu chỉnh thang đo JDI phù hợp với thể chế công lập: Bổ sung các biến quan sát đặc thù về phụ cấp nghề, tiêu chuẩn phòng học theo Thông tư 13/2020/TT-BGDĐT và quan hệ phối hợp với phụ huynh học sinh vùng nông thôn/miền núi.

So sánh với các công trình nghiên cứu trước đây

Tiêu chí phân tích Nghiên cứu Nguyễn Hồng Hà & CS (2021) Nghiên cứu Phan Thị Mai Hương & CS (2021) Đề tài nghiên cứu Hà Tĩnh (Phan Uyên Nhi, 2023)
Phạm vi địa lý Huyện Tam Nông, Tỉnh Đồng Tháp Toàn quốc (Khảo sát trực tuyến chung) Toàn diện 13 huyện, thị xã, TP Tỉnh Hà Tĩnh
Quy mô mẫu $N = 228$ giáo viên mầm non $N = 650$ giáo viên phổ thông $N = 1.130$ giáo viên + 55 phỏng vấn sâu
Phạm vi cấp học Chỉ nghiên cứu cấp Mầm non 3 cấp (Tiểu học, THCS, THPT) Trọn vẹn 4 cấp: Mầm non, Tiểu học, THCS, THPT
Phương pháp Định lượng đơn thuần (SPSS) Khảo sát thang đo hỗn hợp Sequential Mixed-Methods (Định tính + Định lượng)
Tỷ lệ phương sai giải thích ($R^2$) 48,2% 42,5% - 51,0% 54,8% (Mô hình tổng thể) - 61,2% (Mô hình Mầm non)
Độ chi tiết chính sách Khuyến nghị chung cho cấp huyện Đề xuất định hướng vĩ mô Khuyến nghị ma trận chi tiết cho từng cấp học và Sở GD&ĐT

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

Kịch bản ứng dụng thực tiễn

  • Dành cho Sở Giáo dục & Đào tạo Hà Tĩnh: Sử dụng phương trình hồi quy và trọng số $\beta$ để phân bổ nguồn lực ngân sách hằng năm. Ưu tiên giải ngân gói hỗ trợ cơ sở vật chất cho các trường Mầm non và Tiểu học ở các huyện miền núi khó khăn (Hương Khê, Kỳ Anh, Vũ Quang).
  • Dành cho Ban Giám hiệu các trường THCS và THPT: Tái cấu trúc quy trình đánh giá thi đua nội bộ, áp dụng cơ chế đánh giá hiệu suất minh bạch (Key Performance Indicators - KPI kết hợp chuẩn nghề nghiệp giáo viên) nhằm triệt tiêu bất mãn về công tác quản lý.
                    LỘ TRÌNH TRIỂN KHAI 3 GIAI ĐOẠN
  +-----------------------------------------------------------------------+
  | GIAI ĐOẠN 1: TỐI ƯU HÓA ĐÃI NGỘ & GIẢM TẢI HÀNH CHÍNH (Tháng 1-6)     |
  | - Chi trả 100% kinh phí dạy thừa giờ, dạy liên trường.                |
  | - Cắt giảm 40% thủ tục sổ sách giấy tờ, số hóa học bạ điện tử.        |
  +-----------------------------------------------------------------------+
                                     |
                                     v
  +-----------------------------------------------------------------------+
  | GIAI ĐOẠN 2: CHUẨN HÓA CƠ SỞ VẬT CHẤT & BỒI DƯỠNG (Tháng 7-18)        |
  | - Nâng cấp CSVC đạt chuẩn Thông tư 13/2020: THCS từ 40% lên 75%.      |
  | - Tổ chức tập huấn phương pháp dạy học theo Chương trình GDPT 2018.    |
  +-----------------------------------------------------------------------+
                                     |
                                     v
  +-----------------------------------------------------------------------+
  | GIAI ĐOẠN 3: VĂN HÓA QUẢN TRỊ DÂN CHỦ & PHÁT TRIỂN NGHỀ NGHIỆP (Tháng 19-36)|
  | - Chuẩn hóa quy chế dân chủ trường học, công khai tài chính 100%.     |
  | - Xây dựng lộ trình thăng hạng chức danh nghề nghiệp Hạng I, II.       |
  +-----------------------------------------------------------------------+

Phân tích Hiệu quả Đầu tư & Giữ chân Nhân sự (Cost-Benefit Analysis)

  • Chi phí tuyển dụng và đào tạo thay thế: Chi phí trung bình để tuyển dụng, đào tạo và phân bổ một giáo viên mới thay thế vị trí bỏ việc ước tính khoảng 45 - 60 triệu VNĐ/người (bao gồm chi phí ngân sách tổ chức thi tuyển, bồi dưỡng nghiệp vụ sư phạm ban đầu và chi phí gián đoạn chất lượng đào tạo học sinh).
  • Lợi ích dự báo: Nếu triển khai đồng bộ chính sách phụ cấp thu hút và cải thiện môi trường làm việc với mức ngân sách ước tính 15 tỷ VNĐ/năm cho toàn tỉnh, tỷ lệ giáo viên xin chuyển vùng hoặc bỏ việc dự kiến giảm từ $1,8%$ xuống dưới $0,4%$/năm (giữ chân được khoảng 180 - 220 nhà giáo có kinh nghiệm mỗi năm, bảo toàn nguồn lực tri thức tương đương 12 - 15 tỷ VNĐ giá trị đào tạo tích lũy).

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

Hạn chế kỹ thuật và phương pháp

  1. Tính đại diện địa lý: Dữ liệu khảo sát giới hạn trong địa bàn tỉnh Hà Tĩnh, do đó các hệ số hồi quy phản ánh đặc trưng kinh tế - xã hội của vùng Bắc Trung Bộ, cần thận trọng khi suy rộng cho các đô thị đặc thù như Hà Nội hay TP. Hồ Chí Minh.
  2. Phương pháp lấy mẫu cắt ngang (Cross-sectional): Dữ liệu thu thập tại một thời điểm (Quý 3/2023), chưa phản ánh được biến động tâm lý của giáo viên theo chu kỳ năm học (ví dụ: giai đoạn cao điểm thi cử hoặc tổng kết học kỳ).
  3. Sử dụng thang đo tự đánh giá (Self-reported bias): Kết quả khảo sát có thể chịu ảnh hưởng bởi tâm lý e ngại cấp trên ở một bộ phận giáo viên khi đánh giá về biến Lãnh đạo & Quản lý (LDQL).

Hướng nghiên cứu tiếp theo

  • Mở rộng mô hình nghiên cứu dạng chuỗi thời gian (Longitudinal Study) để theo dõi mức độ biến thiên của sự hài lòng trước và sau khi triển khai bảng lương mới theo vị trí việc làm.
  • Ứng dụng Mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) để phân tích các biến trung gian (Mediating variables) như Sự gắn kết tổ chức (Organizational Commitment) và Động lực làm việc nội tại (Intrinsic Motivation).
  • Mở rộng khảo sát so sánh giữa khối trường công lập và hệ thống trường tư thục chất lượng cao.

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

+-----------------------------------------------------------------------------------+
|                     MA TRẬN GIÁ TRỊ DÀNH CHO CÁC ĐỐI TƯỢNG                        |
+-----------------------------------------------------------------------------------+
| Đối tượng              | Giá trị cụ thể thu nhận                                 |
+------------------------+---------------------------------------------------------+
| **Sinh viên & Học viên**| Khung phương pháp luận mẫu về nghiên cứu định lượng      |
|                        | ứng dụng SPSS/Python; bộ thang đo chuẩn hóa 29 biến.     |
+------------------------+---------------------------------------------------------+
| **Cán bộ Quản lý GD**  | Căn cứ thực chứng để tái cấu trúc môi trường sư phạm,    |
|                        | quy chế khen thưởng và dân chủ hóa học đường.           |
+------------------------+---------------------------------------------------------+
| **Nhà hoạch định CS**  | Bộ số liệu định lượng về mức độ thiếu hụt CSVC và       |
|                        | trọng số tác động của lương thưởng đến tỷ lệ bỏ việc.    |
+------------------------+---------------------------------------------------------+
| **Nhà nghiên cứu**     | Nền tảng dữ liệu thực địa $N = 1.130$ phục vụ phân tích  |
|                        | so sánh đa biến và kiểm định mô hình quản trị nhân lực.  |
+------------------------+---------------------------------------------------------+

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

1. Điều kiện kỹ thuật tối thiểu để tái lập quy trình phân tích dữ liệu của đề tài là gì?

Hệ thống yêu cầu máy trạm chạy hệ điều hành Windows 10/11 hoặc Linux (Ubuntu 20.04 LTS trở lên), RAM tối thiểu 8GB, cài đặt phần mềm IBM SPSS Statistics Version 20.0 trở lên hoặc môi trường Python 3.10+ (với các gói pandas, statsmodels, pingouin, scipy). Bộ dữ liệu cần được mã hóa sạch dạng .sav hoặc .csv.

2. Làm thế nào để xử lý hiện tượng tương quan cao giữa biến Lãnh đạo (LDQL) và Cơ hội phát triển (CH)?

Trong phân tích EFA, áp dụng phép xoay trực giao Varimax (Orthogonal Rotation) để tối đa hóa phương sai của các hệ số tải và phân tách rõ ràng cấu trúc nhân tố. Trong phân tích hồi quy, cần kiểm tra chỉ số chấp nhận (Tolerance) $> 0,2$ và hệ số phóng đại phương sai $\text{VIF} < 2,0$. Nếu $\text{VIF} \ge 2,0$, cần gom biến bằng phân tích nhân tố bậc 2 (Second-order Factor Analysis).

3. Tại sao biến Đặc điểm công việc (DD) lại bị loại khỏi mô hình hồi quy ở cấp Mầm non và Tiểu học?

Kết quả kiểm định Pearson và Hồi quy OLS cho thấy với giáo viên Mầm non và Tiểu học, biến DD có hệ số $p\text{-value} > 0,05$. Điều này phản ánh thực tế rằng giáo viên ở 2 cấp học này xem khối lượng công việc chăm sóc và dạy dỗ là trách nhiệm cố định theo đặc thù nghề nghiệp; mức độ thỏa mãn của họ phụ thuộc gần như tuyệt đối vào thu nhập thực tế (LTPL) và điều kiện lớp học (DK).

4. Chi phí khảo sát diện rộng $N = 1.130$ mẫu và 55 cuộc phỏng vấn sâu được tối ưu như thế nào?

Nghiên cứu áp dụng mô hình kết hợp (Hybrid Survey Mode): $70%$ dữ liệu thu thập qua nền tảng Google Forms chuyên biệt gửi qua mạng lưới chuyên môn của Sở/Phòng GD&ĐT; $30%$ thu thập trực tiếp kết hợp phỏng vấn sâu tại các huyện khó khăn nhằm tiết kiệm $65%$ chi phí đi lại và in ấn phiếu giấy.

5. Kết quả nghiên cứu có ý nghĩa gì đối với việc thực hiện Chương trình GDPT 2018?

Chương trình GDPT 2018 đòi hỏi phương pháp dạy học tích hợp và trang thiết bị thực hành số hóa cao. Kết quả hồi quy khẳng định nếu không kịp thời bù đắp khoảng trống $55% - 60%$ thiếu hụt thiết bị bộ môn (DK3) và chế độ thù lao dạy liên môn (LTPL1), tỷ lệ bất mãn của giáo viên THCS và THPT sẽ tiếp tục gia tăng.

Kết luận

Công trình nghiên cứu của tác giả Phan Uyên Nhi (Trường Đại học Kinh tế - ĐHQGHN, 2023) đã giải quyết toàn diện bài toán đo lường thực chứng các nhân tố tác động đến sự hài lòng công việc của đội ngũ giáo viên công lập tại tỉnh Hà Tĩnh. Bằng việc kết hợp chặt chẽ giữa 55 cuộc phỏng vấn sâu định tính và xử lý toán thống kê trên $1.130$ mẫu định lượng, đề tài đã chứng minh một cách thuyết phục rằng:

  1. Lương thưởng và Phúc lợi (LTPL) giữ vai trò nền tảng chi phối mạnh nhất đến sự hài lòng ($\beta = 0,412$ trên toàn mẫu và $\beta = 0,587$ tại cấp Mầm non), phản ánh đúng quy luật của Tháp nhu cầu Maslow và Thuyết hai nhân tố Herzberg trong điều kiện kinh tế miền Trung còn nhiều khó khăn.
  2. Cơ cấu tác động có sự phân hóa sâu sắc theo cấp học: Cấp Mầm non phụ thuộc vào Lương & Điều kiện làm việc; cấp Tiểu học cần thêm Cơ hội phát triển; cấp THCS chịu ảnh hưởng lớn từ Năng lực lãnh đạo quản trị; trong khi cấp THPT đòi hỏi cao nhất về Môi trường khẳng định năng lực chuyên môn.
  3. Cung cấp cơ sở khoa học xác đáng giúp ngành Giáo dục Hà Tĩnh và các cơ quan quản lý nhà nước xây dựng chính sách đãi ngộ đặc thù, phân bổ vốn đầu tư cơ sở vật chất chuẩn hóa và nâng cao năng lực quản trị trường học, hướng tới mục tiêu phát triển bền vững nguồn nhân lực giáo dục trong kỷ nguyên đổi mới.