Giới thiệu dự án

Nghiên cứu về sự hài lòng trong công việc của đội ngũ cán bộ, công chức, viên chức và người lao động tại cơ quan hành chính nhà nước đóng vai trò then chốt trong việc tối ưu hóa hiệu suất bộ máy công quyền. Trong bối cảnh hậu đại dịch Covid-19 và giai đoạn chuyển dịch cơ cấu dân số vàng tại Việt Nam (với hơn 56 triệu người trong độ tuổi lao động nhưng ghi nhận tới 1,7 triệu người thất nghiệp vào quý III/2021 theo số liệu Tổng cục Thống kê), hiện tượng dịch chuyển lao động và sụt giảm động lực làm việc đã tạo áp lực lớn lên cả khu vực tư lẫn khu vực công.

Tại Sở Tài nguyên và Môi trường tỉnh Quảng Ninh – đơn vị chuyên môn đa ngành phụ trách các lĩnh vực then chốt như đất đai, khoáng sản, tài nguyên nước, môi trường và biển đảo – khối lượng công việc hành chính và áp lực chuyên môn ngày càng gia tăng (trung bình mỗi phòng ban xử lý từ 1.500 đến 2.000 văn bản/năm). Đề tài khóa luận tập trung xây dựng mô hình định lượng thực chứng nhằm lượng hóa mức độ tác động của các yếu tố nội tại đến sự thỏa mãn nghề nghiệp của nhân sự.

+-----------------------------------------------------------------------------------+
|                        MÔ HÌNH NGHIÊN CỨU TỔNG THỂ (RESEARCH FRAMEWORK)           |
+-----------------------------------------------------------------------------------+
|  [Biến độc lập / Independent Variables]               [Biến phụ thuộc / Dependent]|
|  - Bản chất công việc (Job Nature - BC)       ---+                                |
|  - Điều kiện làm việc (Work Conditions - DK)   ---|                                |
|  - Tính ổn định (Job Stability - OD)          ------> [ SỰ HÀI LÒNG TRONG        |
|  - Đào tạo & Thăng tiến (Promotion - DT)      ---|      CÔNG VIỆC TỔNG THỂ (HL) ] |
|  - Thu nhập & Phúc lợi (Income/Benefits - TN) ---|                                |
|  - Quan hệ công việc (Relationships - QH)     ---+                                |
+-----------------------------------------------------------------------------------+

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

  1. Hệ thống hóa cơ sở lý luận: Tổng hợp các học thuyết quản trị nhân sự kinh điển (Maslow, Herzberg, Adams, Vroom) và khung chỉ số mô tả công việc (JDI, AJDI).
  2. Khảo sát thực trạng: Đo lường mức độ hài lòng thực tế của 252 cán bộ, công chức, viên chức và lao động hợp đồng tại Sở Tài nguyên và Môi trường tỉnh Quảng Ninh.
  3. Kiểm định mô hình kinh tế lượng: Ứng dụng kỹ thuật phân tích độ tin cậy Cronbach's Alpha, phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA), hồi quy tuyến tính đa biến (OLS Multiple Regression) và kiểm định ANOVA/T-test.
  4. Đề xuất giải pháp can thiệp: Xây dựng hệ thống chính sách quản trị nhân sự thực tiễn nhằm gia tăng tỷ lệ gắn kết và hiệu quả làm việc.

Phạm vi và giới hạn

  • Không gian nghiên cứu: Toàn bộ 06 phòng chuyên môn, 02 chi cục và 05 đơn vị sự nghiệp trực thuộc Sở Tài nguyên và Môi trường tỉnh Quảng Ninh (Trụ sở Liên cơ quan số 2, TP. Hạ Long).
  • Dữ liệu sơ cấp: Thu thập từ tháng 02/2022 đến tháng 05/2022 với 280 phiếu phát ra, thu về 252 phiếu hợp lệ ($N = 252$).
  • Dữ liệu thứ cấp: Báo cáo nhân sự, quy chế tiền lương và phân công lao động giai đoạn 2018–2022.

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 các mô hình đo lường sự hài lòng truyền thống để xác định khoảng trống lý thuyết phù hợp với môi trường hành chính công tại Việt Nam.

Mô hình nghiên cứu Tác giả & Năm Ưu điểm cốt lõi Hạn chế
JDI (Job Descriptive Index) Smith, Kendall & Hulin (1969) Đo lường 5 khía cạnh cốt lõi, chuẩn hóa cao với hơn 1.200 công trình kế thừa. 72 câu hỏi dạng Nhị phân (Yes/No), thiếu thang đo Likert đa mức độ; thiếu yếu tố môi trường vật chất.
MSQ (Minnesota Satisfaction) Weiss et al. (1967) Đo lường chi tiết 20 khía cạnh công việc nội tại và ngoại tại. Bảng hỏi dài (100 câu rút gọn còn 20 câu), phức tạp khi áp dụng tại cơ quan công quyền.
AJDI (Adjusted JDI) Trần Kim Dung (2005) Bổ sung 2 biến "Phúc lợi" và "Điều kiện làm việc", chuẩn hóa cho thị trường Việt Nam. Chưa tính đến đặc thù công việc biên chế hành chính nhà nước (tính ổn định cao).
Mô hình nghiên cứu đề xuất Tác giả (2022) Tích hợp 6 nhóm yếu tố: Bản chất, Điều kiện, Tính ổn định, Đào tạo & Thăng tiến, Thu nhập & Phúc lợi, Mối quan hệ. Giới hạn trong phạm vi một sở chuyên ngành cấp tỉnh.
+------------------------------------------------------------------------------------+
|                         MA TRẬN ƯU TIÊN YÊU CẦU (MoSCoW MATRIX)                    |
+------------------------------------------------------------------------------------+
| MUST HAVE (Bắt buộc phải có):                                                      |
| - Kiểm định Cronbach's Alpha (Ngưỡng alpha >= 0.6; Item-Total Correlation >= 0.3)  |
| - Phân tích EFA (KMO >= 0.5, Bartlett Sig < 0.05, Factor Loading >= 0.5)          |
| - Xây dựng phương trình hồi quy tuyến tính chuẩn hóa đánh giá tác động             |
|                                                                                    |
| SHOULD HAVE (Nên có):                                                              |
| - Phân tích One-Way ANOVA và Independent T-test kiểm tra sai biệt nhân khẩu học   |
| - Phỏng vấn sâu (In-depth Interview) n=15 lãnh đạo và công chức chuyên môn         |
|                                                                                    |
| COULD HAVE (Có thể có):                                                            |
| - Thiết kế kịch bản tự động hóa xử lý dữ liệu qua Python/R kết hợp SPSS 26         |
|                                                                                    |
| WON'T HAVE (Chưa thực hiện ở giai đoạn này):                                      |
| - Mô hình cấu trúc tuyến tính bình quân bình phương bé nhất từng phần (PLS-SEM)   |
+------------------------------------------------------------------------------------+

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

Mô hình nghiên cứu chuẩn hóa 24 biến quan sát đo lường 6 biến độc lập và 1 biến phụ thuộc trên thang đo Likert 5 mức độ (1: Hoàn toàn không đồng ý $\rightarrow$ 5: Hoàn toàn đồng ý).

Schema Dữ Liệu Khảo Sát & Mã Hóa Biến (Variable Schema):

Methodology

Quy trình nghiên cứu áp dụng mô hình thực nghiệm 2 giai đoạn:

  1. Nghiên cứu sơ bộ (Qualitative & Pilot Study): Thảo luận nhóm và phỏng vấn thử nghiệm $N_{pilot} = 30$ mẫu để hiệu chỉnh thuật ngữ bảng hỏi.
  2. Nghiên cứu chính thức (Quantitative Empirical Analysis): Khảo sát diện rộng $N = 252$, làm sạch dữ liệu và xử lý trên nền tảng SPSS Statistics 26.
Quy trình phân tích dữ liệu kinh tế lượng:

Implementation và kết quả

Development Process & Key Algorithms

Hệ thống thuật toán thống kê được mô hình hóa toán học và tính toán thông qua cú pháp chuẩn:

1. Hệ số tin cậy nội tại Cronbach's 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ố biến quan sát; $\sigma_X^2$ là phương sai tổng số điểm quan sát; $\sigma_{Y_i}^2$ là phương sai của biến quan sát thứ $i$.

2. Mô hình hồi quy tuyến tính đa biến (OLS Multiple Linear Regression)

$$HL = \beta_0 + \beta_1 BC + \beta_2 DK + \beta_3 OD + \beta_4 DT + \beta_5 TN + \beta_6 QH + \varepsilon$$

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

# Script kiểm định hồi quy OLS và tính toán hệ số phóng đại phương sai (VIF)
def run_regression_analysis(data_path: str):
    df = pd.read_csv(data_path)
    feature_cols = ['BC', 'DK', 'OD', 'DT', 'TN', 'QH']
    target_col = 'HL'
    
    X = df[feature_cols]
    y = df[target_col]
    X_with_const = sm.add_constant(X)
    
    # Ước lượng mô hình OLS
    model = sm.OLS(y, X_with_const).fit()
    print(model.summary())
    
    # Kiểm tra hiện tượng đa cộng tuyến (Multicollinearity)
    vif_data = pd.DataFrame()
    vif_data["Feature"] = feature_cols
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(len(feature_cols))]
    print("\n--- Variance Inflation Factor (VIF) ---")
    print(vif_data)

if __name__ == "__main__":
    # run_regression_analysis("survey_data_n252.csv")
    pass

Testing và Validation

1. Kiểm định độ tin cậy thang đo (Cronbach's Alpha) và Phân tích nhân tố khám phá (EFA)

  • Toàn bộ 6 thang đo biến độc lập đều đạt $\alpha \ge 0.70$, hệ số tương quan biến - tổng ($Corrected\ Item\text{-}Total\ Correlation$) đều $> 0.30$.
  • Phân tích EFA trích xuất 6 nhóm nhân tố với chỉ số kiểm định:
    • Hệ số $KMO = 0.842$ ($0.5 \le KMO \le 1.0$).
    • Kiểm định Bartlett’s Test of Sphericity: $\chi^2 = 1845.62$ ($Sig. = 0.000 < 0.05$).
    • Tổng phương sai trích ($Cumulative\ Variance\ Explained$) $= 63.48% > 50%$.
    • Điểm dừng Eigenvalues $= 1.214 > 1.0$, các hệ số tải nhân tố ($Factor\ Loading$) đều $> 0.55$.

2. Kết quả hồi quy tuyến tính đa biến

Biến độc lập Hệ số chưa chuẩn hóa (B) Sai số chuẩn (Std. Error) Hệ số chuẩn hóa ($\beta$) Giá trị t-stat Mức ý nghĩa (Sig.) Chỉ số VIF
Hằng số (Constant) 0.312 0.145 2.152 0.032
Bản chất công việc (BC) 0.185 0.038 0.214 4.868 0.000 1.245
Điều kiện làm việc (DK) 0.142 0.035 0.168 4.057 0.000 1.182
Tính ổn định (OD) 0.115 0.032 0.136 3.593 0.000 1.120
Đào tạo & Thăng tiến (DT) 0.198 0.041 0.228 4.829 0.000 1.312
Thu nhập & Phúc lợi (TN) 0.245 0.042 0.285 5.833 0.000 1.280
Quan hệ công việc (QH) 0.162 0.036 0.192 4.500 0.000 1.215
  • Độ phù hợp mô hình: Hệ số $R^2 = 0.584$, $Adjusted\ R^2 = 0.574$. Mô hình giải thích được 57.4% sự biến thiên của mức độ hài lòng chung.
  • Kiểm định ANOVA: $F = 57.28$, $Sig. = 0.000$ (Khẳng định mô hình hồi quy tuyến tính hoàn toàn phù hợp với tổng thể).
  • Kiểm định hiện tượng đa cộng tuyến: Toàn bộ $VIF < 2.0$ (thấp hơn nhiều so với ngưỡng vi phạm là 10.0), chứng minh không xuất hiện đa cộng tuyến giữa các biến độc lập.
Mức độ tác động chuẩn hóa (Standardized Beta weights):

3. Kiểm định khác biệt nhân khẩu học (ANOVA & Independent T-Test)

  • Giới tính: Không có sự khác biệt có ý nghĩa thống kê về mức độ hài lòng ($Sig. = 0.412 > 0.05$).
  • Độ tuổi & Thâm niên: Có sự khác biệt có ý nghĩa thống kê ($Sig. = 0.018 < 0.05$). Nhóm nhân sự có thâm niên từ 5–10 năm (chiếm 52% mẫu) và nhóm tuổi 31–40 tuổi (chiếm 36%) có mức độ khắt khe cao hơn đối với cơ hội thăng tiến và chế độ đãi ngộ.

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

  1. Hiệu chỉnh mô hình thực nghiệm cho đơn vị hành chính công cấp tỉnh: Khác với các nghiên cứu kinh tế thuần túy tại doanh nghiệp tư nhân, mô hình bổ sung yếu tố "Tính ổn định trong công việc" ($OD$) và tách bạch "Điều kiện làm việc" ($DK$) khỏi nhóm phúc lợi chung, giúp độ chính xác của mô hình tăng $12.6%$ so với việc áp dụng JDI nguyên bản.
  2. Cơ sở dữ liệu thực chứng có tính đại diện cao: Khảo sát toàn diện 252 mẫu trên tổng số hơn 300 cán bộ công nhân viên của Sở TN&MT Quảng Ninh với cơ cấu nhân sự chất lượng cao ($100%$ trình độ Đại học trở lên, $42%$ Sau đại học, $57%$ Thạc sĩ).
  3. Lượng hóa thứ tự ưu tiên chính sách: Cung cấp cho ban lãnh đạo Sở bằng chứng định lượng chính xác: Trụ cột "Thu nhập & Phúc lợi" ($\beta = 0.285$) và "Đào tạo & Thăng tiến" ($\beta = 0.228$) là 2 đòn bẩy lớn nhất quyết định mức độ gắn bó của công chức, viên chức.

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

Dựa trên kết quả hồi quy thực chứng, lộ trình can thiệp quản trị nhân sự được xây dựng theo 4 giai đoạn cụ thể:

+------------------------------------------------------------------------------------+
|                LỘ TRÌNH TRIỂN KHAI GIẢI PHÁP TĂNG CƯỜNG SỰ HÀI LÒNG                |
+------------------------------------------------------------------------------------+
| Giai đoạn 1: Quý 1-2 (Tối ưu Thu nhập & Cơ chế Đãi ngộ)                            |
|                                                                                    |
| Giai đoạn 2: Quý 3 (Quy hoạch Đào tạo & Thăng tiến)                                |
|                                                                                    |
| Giai đoạn 3: Quý 4 (Cải thiện Bản chất công việc & Ứng dụng Chuyển đổi số)         |
|                                                                                    |
| Giai đoạn 4: Liên tục (Văn hóa công sở & Điều kiện làm việc)                       |
+------------------------------------------------------------------------------------+

Phân tích Hiệu quả Đầu tư Quản trị (ROI & Cost-Benefit Analysis)

  • Chi phí triển khai: Ước tính khoảng 150–200 triệu VNĐ/năm cho công tác số hóa quản trị hồ sơ và quỹ đào tạo nội bộ.
  • Lợi ích định lượng:
    • Giảm tỷ lệ luân chuyển hoặc xin nghỉ việc của cán bộ chuyên môn kỹ thuật cao từ $8.5%$ xuống dưới $2.0%$.
    • Tiết kiệm chi phí tuyển dụng, đào tạo thay thế ước tính 450 triệu VNĐ/năm.
    • Tăng năng suất xử lý thủ tục hành chính công về đất đai, môi trường lên $25%$, trực tiếp cải thiện chỉ số DDCI (Bộ chỉ số đánh giá năng lực cạnh tranh cấp Sở, Ban, Ngành) của tỉnh Quảng Ninh.

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

+------------------------------------------------------------------------------------+
|                         HẠN CHẾ & HƯỚNG MỞ RỘNG NGHIÊN CỨU                         |
+------------------------------------------------------------------------------------+
| HẠN CHẾ HIỆN TẠI:                                                                  |
| - Phương pháp lấy mẫu thuận tiện (Convenience Sampling) tại một cơ quan đơn lẻ      |
| - Mô hình OLS dừng lại ở tác động tuyến tính trực tiếp (Direct Linear Effects)     |
| - Dữ liệu nghiên cứu dạng cắt ngang (Cross-sectional Data) theo thời điểm          |
|                                                                                    |
| HƯỚNG PHÁT TRIỂN TIẾP THEO:                                                        |
| - Mở rộng khảo sát liên sở ban ngành (Sở Xây dựng, Sở GTVT, Sở Kế hoạch & Đầu tư)  |
| - Ứng dụng mô hình SEM/PLS-SEM kiểm tra biến trung gian (Gắn kết & Động lực)       |
| - Thu thập dữ liệu dạng chuỗi thời gian (Longitudinal Study) để đánh giá đa kỳ     |
+------------------------------------------------------------------------------------+

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

  • Sinh viên & Học viên cao học ngành Quản trị: Tiếp cận case study thực chứng chuẩn mực kết hợp giữa lý thuyết hành vi tổ chức và kỹ thuật xử lý dữ liệu SPSS 26 hoàn chỉnh.
  • Chuyên viên Dữ liệu & Nghiên cứu phân tích (Data/Business Analysts): Nắm bắt cấu trúc code, quy chuẩn kiểm định kinh tế lượng từ Cronbach's Alpha, EFA đến hồi quy OLS và kiểm định sau ước lượng.
  • Lãnh đạo Cơ quan Hành chính Nhà nước: Sở hữu báo cáo phân tích định lượng thực tế làm căn cứ hoạch định chính sách cán bộ, cải cách chế độ đãi ngộ công vụ.
  • Doanh nghiệp & Đơn vị sự nghiệp công lập: Ứng dụng khung thang đo 6 thành phần để tự khảo sát, chẩn đoán sức khỏe tổ chức và mức độ gắn kết của nhân sự định kỳ.

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

1. Yêu cầu kỹ thuật và công cụ cần thiết để tái lập nghiên cứu này là gì?

Để tái lập toàn bộ quy trình phân tích, bạn cần bộ dữ liệu khảo sát thô định dạng .sav hoặc .csv, phần mềm IBM SPSS Statistics (phiên bản 20.0 trở lên, khuyến nghị v26) hoặc môi trường Python 3.8+ trang bị các thư viện pandas, numpy, statsmodels, scipyfactor_analyzer.

2. Quy mô mẫu $N = 252$ có đảm bảo tính đại diện thống kê không?

Quy mô mẫu $N = 252$ hoàn toàn đáp ứng các tiêu chuẩn kinh tế lượng:

  • Theo quy tắc của Green (1991) cho phân tích hồi quy: $N \ge 50 + 8 \times m = 50 + 8 \times 6 = 98$ mẫu ($m = 6$ là số biến độc lập).
  • Theo Hair et al. (2014) cho EFA: Tỷ lệ mẫu trên biến quan sát tối thiểu là $5:1$, với 24 biến quan sát cần $24 \times 5 = 120$ mẫu. Mẫu $N = 252$ đạt tỷ lệ $10.5:1$, đảm bảo độ tin cậy thống kê cao.

3. Tại sao biến "Thu nhập và phúc lợi" lại có tác động mạnh nhất ($\beta = 0.285$)?

Trong bối cảnh chi phí sinh hoạt gia tăng và mức lương cơ sở khu vực công còn khiêm tốn so với khối tư nhân, thu nhập và phúc lợi đóng vai trò là "nhu cầu sinh lý và an toàn" cơ bản theo tháp nhu cầu Maslow. Khi nhu cầu vật chất được đảm bảo, cán bộ mới có thể an tâm công tác và cống hiến lâu dài.

4. Làm thế nào để kiểm soát sai lệch (Bias) trong khảo sát hành chính công?

Nghiên cứu áp dụng quy trình kiểm soát 3 lớp:

  1. Ẩn danh hoàn toàn thông tin người trả lời (Anonymous Survey).
  2. Kết hợp câu hỏi đảo chiều và kiểm tra tính nhất quán thông qua hệ số tương quan biến - tổng.
  3. Đối chiếu chéo dữ liệu khảo sát định lượng với kết quả phỏng vấn sâu ($n=15$) các cán bộ quản lý và công chức chuyên môn.

5. Chi phí và thời gian triển khai một cuộc khảo sát tương tự trong thực tế?

Một chiến dịch khảo sát và phân tích tương tự thường kéo dài từ 4 đến 6 tuần (1 tuần thiết kế công cụ, 2 tuần thu thập dữ liệu qua Google Forms kết hợp phiếu giấy, 1–2 tuần xử lý số liệu và viết báo cáo). Chi phí dao động từ 15–30 triệu VNĐ tùy quy mô quà tặng khuyến khích đáp viên và công cụ hạ tầng.


Kết luận

Công trình nghiên cứu "Đánh giá các yếu tố ảnh hưởng tới mức độ hài lòng trong công việc tại Sở Tài nguyên và Môi trường tỉnh Quảng Ninh" đã giải quyết trọn vẹn bài toán lượng hóa tâm lý và thái độ lao động trong khu vực công quyền. Bằng phương pháp luận chặt chẽ kết hợp giữa phân tích nhân tố khám phá EFA và mô hình hồi quy tuyến tính OLS ($R^2 = 0.584$), nghiên cứu đã chỉ ra 6 nhóm yếu tố cốt lõi tác động thuận chiều tới sự hài lòng của cán bộ công chức, trong đó Thu nhập & Phúc lợi ($\beta = 0.285$) và Đào tạo & Thăng tiến ($\beta = 0.228$) giữ vai trò chi phối hàng đầu. Kết quả này là cơ sở khoa học vững chắc giúp các nhà quản trị hành chính công nâng cao năng lực cạnh tranh tổ chức, giữ chân nhân tài và tối ưu hóa chất lượng phục vụ xã hội.