Giới thiệu dự án

Trong kỷ nguyên kinh tế tri thức và cạnh tranh toàn cầu hóa, nguồn nhân lực đóng vai trò là tài sản chiến lược quyết định năng lực cạnh tranh cốt lõi của doanh nghiệp. Theo các khảo sát thị trường lao động tại Việt Nam trong lĩnh vực thương mại - dịch vụ kỹ thuật (đặc biệt là hệ thống phân phối và bảo dưỡng xe máy Honda), tỷ lệ biến động nhân sự bình quân dao động từ 22% đến 28%/năm, kéo theo chi phí tuyển dụng và tái đào tạo tăng cao gấp 1,5–2 lần chi phí tiền lương định kỳ.

Đề tài nghiên cứu "Các yếu tố ảnh hưởng đến kết quả làm việc của nhân viên tại Công ty TNHH Thương mại và Dịch vụ Khánh An" do tác giả Lâm Quang Khánh (Trường Đại học Ngân hàng TP. Hồ Chí Minh - HUB) thực hiện nhằm giải quyết trực tiếp bài toán tối ưu hóa năng suất lao động và giữ chân nhân sự chất lượng cao thông qua mô hình định lượng thực nghiệm.

Vấn đề thực tiễn và điểm nghẽn tổ chức (Problem Statement)

Thực trạng quản trị nhân sự tại Công ty TNHH Thương mại và Dịch vụ Khánh An bộc lộ các điểm nghẽn nghiêm trọng:

  • Phân công công việc thiếu minh bạch: Khối lượng công việc và bản mô tả công việc (Job Description) chưa chuẩn hóa; việc giao việc mang nặng tính chủ quan từ lãnh đạo phòng ban.
  • Bất bình đẳng trong thăng tiến: Cơ chế bổ nhiệm ưu tiên thâm niên thay vì năng lực thực tế, tạo tâm lý bất mãn cho nhân viên trẻ có năng lực.
  • Áp lực công việc mất cân bằng: Nhân viên kỹ thuật và kinh doanh đối mặt với chỉ tiêu (KPI) cao nhưng thiếu cơ chế kiểm soát sai sót hỗ trợ, dẫn đến căng thẳng (workplace stress) và suy giảm năng suất.
  • Chính sách đãi ngộ chưa tạo động lực: Lương và phúc lợi chưa phản ánh chính xác giá trị gia tăng mà nhân viên đóng góp.

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

  1. Xác định cụ thể 06 nhóm yếu tố cốt lõi tác động đến kết quả làm việc của nhân viên tại doanh nghiệp thương mại - dịch vụ.
  2. Đo lường mức độ ảnh hưởng (trọng số hồi quy $\beta$) của từng yếu tố lên biến phụ thuộc "Kết quả làm việc" thông qua dữ liệu định lượng.
  3. Xây dựng hệ thống giải pháp và hàm ý quản trị ứng dụng nhằm nâng cao hiệu suất làm việc tổng thể thêm tối thiểu 20% và giảm tỷ lệ luân chuyển lao động.

Phương pháp tiếp cận và kết quả kỳ vọng

  • Phương pháp tiếp cận: Kết hợp nghiên cứu định tính (phỏng vấn nhóm chuyên gia $n=10$) và nghiên cứu định lượng ($N=210$ bảng câu hỏi) thông qua phần mềm SPSS 22.0 với các kỹ thuật Cronbach's Alpha, EFA, OLS Regression và ANOVA/T-test.
  • Kết quả kỳ vọng: Mô hình đạt độ giải thích $R^2 \ge 50%$, hệ số phóng đại phương sai $\text{VIF} < 2.0$, cung cấp bằng chứng thực nghiệm chính xác để ban giám đốc tái cơ cấu chính sách nhân sự.
  • Phạm vi nghiên cứu: Nhân viên làm việc tại Công ty TNHH Thương mại và Dịch vụ Khánh An, thu thập dữ liệu từ tháng 09/2022 đến tháng 10/2022.

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

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

Để xây dựng cơ sở lý thuyết vững chắc, nghiên cứu đã tiến hành tổng quan các mô hình kinh điển: Tháp nhu cầu của Abraham Maslow (1943), Thuyết hai yếu tố của Frederick Herzberg (1959), Thuyết kỳ vọng của Victor Vroom (1964) và Chỉ số Mô tả Công việc JDI (Job Description Index) của Smith, Kendall & Hulin (1969).

Tiêu chí so sánh Mô hình Njambi (2014) Mô hình Đặng Thị Thùy Phương (2013) Mô hình Đề xuất Khánh An (2022)
Đối tượng nghiên cứu Nhân viên y tế (Amref Health Africa - Kenya) Nhân viên ngân hàng liên doanh (TP.HCM) Nhân viên đại lý thương mại & dịch vụ kỹ thuật
Cỡ mẫu ($N$) 412 quan sát 211 quan sát 210 quan sát
Yếu tố Áp lực CV Không tích hợp Không tích hợp Có tích hợp (theo Karasek & Jungwee Park)
Yếu tố Phúc lợi riêng Gộp chung trong thu nhập Tách biệt trong chế độ chính sách Tích hợp sâu theo bối cảnh SME Việt Nam
Phương pháp kiểm định OLS căn bản EFA + Hồi quy tuyến tính Cronbach's Alpha + EFA + OLS + ANOVA/T-test

Ma trận ưu tiên yêu cầu khảo sát (MoSCoW)

  • Must-have: Đo lường tác động của Lương & Phúc lợi (LP), Phong cách Lãnh đạo (LD), Quy trình Đào tạo (DT).
  • Should-have: Tích hợp biến Áp lực công việc (AL) và Đặc điểm công việc (CV).
  • Could-have: Phân tích sâu sai biệt theo phân khúc thu nhập và chức danh bằng ANOVA.
  • Won't-have (giai đoạn này): Mô hình hóa cấu trúc tuyến tính bậc cao PLS-SEM cho biến trung gian (Job Satisfaction).

Thiết kế hệ thống mô hình định lượng

Mô hình nghiên cứu thực nghiệm bao gồm 06 biến độc lập và 01 biến phụ thuộc với 24 biến quan sát chuẩn hóa theo thang đo Likert 5 mức độ ($1 = \text{Rất không đồng ý}$ đến $5 = \text{Rất đồng ý}$):

$$\text{KQ} = \beta_0 + \beta_1 \text{LP} + \beta_2 \text{MT} + \beta_3 \text{CV} + \beta_4 \text{AL} + \beta_5 \text{LD} + \beta_6 \text{DT} + \epsilon$$

+-------------------------------------------------------------------------+
|                       KHUNG MÔ HÌNH HỒI QUY TUYẾN TÍNH                   |
|                                                                         |
|  [H1] Lương và phúc lợi (LP) --------(+)--------\                       |
|  [H2] Môi trường làm việc (MT) ------(+)---------\                      |
|  [H3] Đặc điểm công việc (CV) -------(+)----------\ ---> [KẾT QUẢ       |
|  [H4] Áp lực công việc (AL) ---------(-)----------/     LÀM VIỆC (KQ)]  |
|  [H5] Lãnh đạo (LD) -----------------(+)---------/                      |
|  [H6] Quy trình đào tạo (DT) --------(+)--------/                       |
+-------------------------------------------------------------------------+

Công nghệ và công cụ xử lý

  • Nền tảng phân tích cốt lõi: IBM SPSS Statistics Version 22.0.
  • Ngôn ngữ mở rộng (Pipeline Validation): Python 3.10 với thư viện pandas (v2.0.3), statsmodels (v0.14.0), scikit-learn (v1.3.0).
  • Cơ chế kiểm soát dữ liệu: Thuật toán phát hiện dị biệt (Outlier Detection) dựa trên $Z\text{-score} > 3.0$ và khoảng cách Cook's Distance $< 1.0$.

Methodology (Phương pháp luận)

Quy trình nghiên cứu 2 giai đoạn đảm bảo tính giá trị nội tại (Internal Validity) và độ tin cậy ngoại tại (External Reliability):

Kế hoạch phân bổ và Quản trị rủi ro

Giai đoạn Nội dung công việc Tuần thực hiện Rủi ro tiềm ẩn Biện pháp giảm thiểu
Giai đoạn 1 Tổng quan tài liệu, phỏng vấn sâu $n=10$ Tuần 1 - 2 Thang đo không khớp thực tế Hiệu chỉnh bảng hỏi trực tiếp với HR Manager
Giai đoạn 2 Phát bảng hỏi ($N=210$), thu hồi dữ liệu Tuần 3 - 5 Tỷ lệ phản hồi thấp, sót dữ liệu Kết hợp Google Forms và bảng in giấy có đối soát
Giai đoạn 3 Xử lý Cronbach's Alpha, EFA, OLS trên SPSS Tuần 6 - 7 Vi phạm hiện tượng đa cộng tuyến Kiểm tra ma trận tương quan trước khi hồi quy
Giai đoạn 4 Kiểm định khác biệt ANOVA, viết kết luận Tuần 8 Sai số phân phối phi chuẩn Kiểm định Levene test và Kolmogorov-Smirnov

Implementation và kết quả

Development Process & Data Processing Pipeline

Dữ liệu khảo sát từ $N=210$ nhân viên được số hóa và chuẩn hóa theo quy trình xử lý thống kê nghiêm ngặt.

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

# 1. Pipeline tính hệ số Cronbach's Alpha chuẩn hóa
def calculate_cronbach_alpha(df_items: pd.DataFrame) -> float:
    k = df_items.shape[1]
    item_variances = df_items.var(axis=0, ddof=1).sum()
    total_variance = df_items.sum(axis=1).var(ddof=1)
    alpha = (k / (k - 1)) * (1 - (item_variances / total_variance))
    return float(alpha)

# 2. Xây dựng và kiểm định mô hình hồi quy đa biến OLS
def run_regression_pipeline(X: pd.DataFrame, y: pd.Series):
    X_const = sm.add_constant(X)
    model = sm.OLS(y, X_const).fit()
    
    # Tính toán Variance Inflation Factor (VIF)
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X.columns
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(X.shape[1])]
    
    return model.summary(), vif_data

Testing và Validation

1. Kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Tất cả 07 thang đo đều đạt chuẩn tin cậy với hệ số Cronbach's Alpha $> 0.70$ và hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $> 0.30$.

BẢNG TỔNG HỢP KIỂM ĐỊNH CRONBACH'S ALPHA
+-------------------------------+-------------------+----------------------------+
| Nhóm Thang đo                 | Số biến quan sát  | Hệ số Cronbach's Alpha     |
+-------------------------------+-------------------+----------------------------+
| Lương và phúc lợi (LP)        | 4                 | 0.832                      |
| Môi trường làm việc (MT)      | 4                 | 0.805                      |
| Đặc điểm công việc (CV)       | 3                 | 0.789                      |
| Áp lực công việc (AL)         | 3                 | 0.764                      |
| Lãnh đạo (LD)                 | 4                 | 0.868                      |
| Quy trình đào tạo (DT)        | 3                 | 0.812                      |
| Kết quả làm việc (KQ) - Phụ thuộc | 3             | 0.845                      |
+-------------------------------+-------------------+----------------------------+

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

  • Kiểm định KMO & Bartlett (Biến độc lập):
    • Hệ số $\text{KMO} = 0.842$ ($0.5 \le \text{KMO} \le 1.0$) $\rightarrow$ Mẫu hoàn toàn thích hợp cho phân tích nhân tố.
    • Kiểm định Bartlett: Chi-Square $\chi^2 = 1845.62$ ($p\text{-value} = 0.000 < 0.05$) $\rightarrow$ Các biến quan sát tương quan có ý nghĩa thống kê.
  • Giá trị trích xuất:
    • Tổng phương sai trích (Cumulative % Variance Explained) $= 63.85% > 50%$ tại điểm dừng Eigenvalue $= 1.21 > 1.0$.
    • Toàn bộ 21 biến quan sát độc lập trích thành đúng 6 nhóm nhân tố với trọng số tải nhân tố (Factor Loading) đều đạt chuẩn $\ge 0.58 > 0.50$.
  • Kiểm định EFA Biến phụ thuộc (KQ):
    • $\text{KMO} = 0.718$, $p\text{-value} = 0.000$, Phương sai trích đạt $68.42%$ tại Eigenvalue $= 2.05$.

Kết quả đạt được

Kết quả hồi quy tuyến tính bội (OLS)

MÔ HÌNH HỒI QUY: KQ = 0.115 + 0.312*LD + 0.285*LP + 0.214*CV + 0.186*DT + 0.142*MT - 0.118*AL + e
+--------------------+---------------------+------------------+------------+---------+-------+
| Biến độc lập       | Hệ số chưa chuẩn hóa | Hệ số Beta (β)   | Giá trị t  | Sig. (p)| VIF   |
+--------------------+---------------------+------------------+------------+---------+-------+
| (Hằng số / Const)  | 0.115               | --               | 0.542      | 0.588   | --    |
| Lãnh đạo (LD)      | 0.308               | 0.312            | 5.892      | 0.000   | 1.342 |
| Lương & Phúc lợi   | 0.274               | 0.285            | 5.418      | 0.000   | 1.285 |
| Đặc điểm CV (CV)   | 0.205               | 0.214            | 4.120      | 0.000   | 1.210 |
| Quy trình ĐT (DT)  | 0.179               | 0.186            | 3.654      | 0.000   | 1.198 |
| Môi trường LV (MT) | 0.138               | 0.142            | 2.845      | 0.005   | 1.156 |
| Áp lực CV (AL)     | -0.112              | -0.118           | -2.312     | 0.022   | 1.114 |
+--------------------+---------------------+------------------+------------+---------+-------+
| R = 0.790 | R-Square (R²) = 0.624 | Adjusted R² = 0.613 | F = 56.124 (Sig. = 0.000) | DW = 1.892 |
+---------------------------------------------------------------------------------------------+

Phân tích chẩn đoán mô hình:

  • Độ thích hợp mô hình: Hệ số $R^2 = 0.624$ (Hiệu chỉnh $= 0.613$), chứng minh rằng $61.3%$ sự biến thiên của kết quả làm việc được giải thích bởi 6 yếu tố trên.
  • Kiểm định F: $F = 56.124$ ($p = 0.000 < 0.001$) khẳng định mô hình hồi quy tuyến tính là hoàn toàn phù hợp.
  • Hiện tượng tự tương quan: Durbin-Watson $= 1.892$ (nằm trong khoảng $1.5 < DW < 2.5$), không có hiện tượng tự tương quan chuỗi bậc nhất.
  • Hiện tượng đa cộng tuyến: Toàn bộ hệ số $\text{VIF} \in [1.114; 1.342] < 2.0$, không vi phạm giả định đa cộng tuyến.
  • Kiểm định ANOVA: Phân tích phương sai một yếu tố chỉ ra có sự khác biệt có ý nghĩa thống kê về kết quả làm việc theo Vị trí công việc ($F = 4.125, p = 0.007 < 0.05$) và Mức thu nhập ($F = 3.842, p = 0.011 < 0.05$). Không có sự khác biệt theo Giới tính ($t = 0.452, p = 0.652$).

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

  1. Chuẩn hóa thang đo đặc thù ngành Dịch vụ Kỹ thuật ô tô - xe máy: Đã tích hợp thành công cấu trúc kiểm soát sai sót vào biến "Áp lực công việc (AL)" và "Đặc điểm công việc (CV)", giúp đo lường chính xác tâm lý người lao động làm việc theo quy chuẩn kỹ thuật nghiêm ngặt của Honda.
  2. Làm rõ tác động ngược chiều của Áp lực công việc: Khác với một số nghiên cứu trong ngành ngân hàng xem áp lực là đòn bẩy tích cực, nghiên cứu này chứng minh thực nghiệm tại Khánh An rằng áp lực quá tải tác động tiêu cực ($\beta = -0.118, p < 0.05$) đến kết quả hoàn thành nhiệm vụ.
  3. Phát hiện vai trò tối thượng của Năng lực Lãnh đạo: Yếu tố "Lãnh đạo" đạt trọng số ảnh hưởng cao nhất ($\beta = 0.312$), vượt qua cả "Lương và phúc lợi" ($\beta = 0.285$). Điều này chỉ ra rằng năng lực hỗ trợ, phản hồi và phong cách điều hành của quản lý trực tiếp là động lực then chốt quyết định hiệu suất.
  4. Bằng chứng số liệu đóng góp cho học thuật HRM: Bổ sung tập dữ liệu thực nghiệm $N=210$ tại thị trường lao động Việt Nam vào cơ sở dữ liệu nghiên cứu về sự thỏa mãn và hiệu quả lao động trong khối doanh nghiệp vừa và nhỏ (SMEs).

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

Dựa trên kết quả định lượng, ban lãnh đạo Công ty TNHH Thương mại và Dịch vụ Khánh An có thể triển khai lộ trình 4 giai đoạn tái cấu trúc nhân sự:

Phân tích Chi phí - Lợi ích (Cost-Benefit & ROI Analysis)

  • Chi phí triển khai dự kiến: 180.000.000 VNĐ (Bao gồm chi phí đào tạo quản lý, chuẩn hóa phần mềm đánh giá KPI và quỹ khen thưởng đột xuất).
  • Lợi ích kinh tế lượng hóa sau 12 tháng:
    • Giảm tỷ lệ nghỉ việc từ $25%$ xuống dưới $10%$ $\rightarrow$ Tiết kiệm 240.000.000 VNĐ chi phí tuyển dụng và đào tạo học việc.
    • Tăng năng suất lao động bình quân $18.5%$ thông qua việc giải tỏa áp lực và hoàn thiện kỹ năng $\rightarrow$ Đóng góp gia tăng doanh thu dịch vụ ước tính 450.000.000 VNĐ.
  • Tỷ suất sinh lời đầu tư (ROI):

$$\text{ROI} = \frac{(240.000.000 + 450.000.000) - 180.000.000}{180.000.000} \times 100% = 283.3%$$

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

Hạn chế nghiên cứu

  • Phương pháp lấy mẫu: Sử dụng kỹ thuật chọn mẫu thuận tiện phi xác suất ($N=210$), tập trung tại một doanh nghiệp thương mại dịch vụ cụ thể, giới hạn khả năng tổng quát hóa (Generalizability) cho toàn bộ ngành bán lẻ xe gắn máy toàn quốc.
  • Thời gian nghiên cứu: Dữ liệu cắt ngang (Cross-sectional data) thu thập trong 2 tháng (09/2022 - 10/2022) chưa theo dõi được biến động dài hạn (Longitudinal effect) của nhân sự qua các chu kỳ kinh doanh cao điểm (Tết/cuối năm).
  • Phương pháp phân tích: Mô hình OLS đo lường tác động trực tiếp, chưa phân tích các biến trung gian (Mediating variables) như "Sự gắn kết với tổ chức" hay biến điều tiết (Moderating variables) như "Văn hóa doanh nghiệp".

Định hướng phát triển

  • Ứng dụng mô hình cấu trúc tuyến tính bình phương bé nhất từng phần (PLS-SEM) để phân tích đa tầng các mối quan hệ nhân quả gián tiếp.
  • Mở rộng phạm vi nghiên cứu đa trung tâm (Multi-center study) trên 10+ đại lý Head Honda tại khu vực Đông Nam Bộ.

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

  • Sinh viên & Học viên cao học: Cung cấp tài liệu tham khảo chuẩn mực về cấu trúc khóa luận định lượng ngành Quản trị kinh doanh, minh họa chi tiết từng bước từ lập bảng hỏi đến chạy SPSS 22.0.
  • Chuyên viên Dữ liệu & Nhà phân tích Nhân sự (HR Data Analysts): Tham khảo quy trình tiền xử lý dữ liệu, bộ cú pháp kiểm định Cronbach's Alpha, EFA và các phép thử vi phạm mô hình kinh tế lượng.
  • Ban Giám đốc & Trưởng phòng Nhân sự: Sở hữu khung tham chiếu thực chứng để phân bổ ngân sách đào tạo và tái cơ cấu lương thưởng dựa trên mức độ ưu tiên: Lãnh đạo ($\beta=0.312$) $\rightarrow$ Lương thưởng ($\beta=0.285$) $\rightarrow$ Đặc điểm CV ($\beta=0.214$).
  • Nhà nghiên cứu học thuật: Tiếp cận cơ sở dữ liệu thực nghiệm đã kiểm định về mối quan hệ nghịch chiều giữa áp lực công việc và hiệu suất trong môi trường SME Việt Nam.

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

1. Cần chuẩn bị những điều kiện gì để triển khai nghiên cứu định lượng tương tự?

Cần xác định tập biến quan sát rõ ràng dựa trên các lý thuyết nền tảng (JDI, Herzberg), thiết kế bảng hỏi Likert 5 điểm, tối thiểu cỡ mẫu đạt tỷ lệ $5:1$ so với số biến quan sát (tối thiểu 120 mẫu cho 24 biến; khuyến nghị $\ge 200$ mẫu) và sử dụng phần mềm SPSS phiên bản 20.0 trở lên.

2. Vì sao biến "Áp lực công việc" lại có hệ số hồi quy mang dấu âm ($\beta = -0.118$)?

Trong môi trường kỹ thuật dịch vụ, áp lực quá mức từ chỉ tiêu thời gian sửa chữa và giao việc chủ quan gây ra căng thẳng tâm lý và mệt mỏi, làm tăng tỷ lệ lỗi kỹ thuật, dẫn đến suy giảm kết quả làm việc thực tế của nhân viên.

3. Làm thế nào để giải quyết vấn đề đa cộng tuyến khi hệ số tương quan giữa các biến độc lập quá cao?

Cần loại bỏ các biến quan sát có hệ số tương quan biến - tổng $< 0.30$, kiểm tra chỉ số $\text{VIF} < 2.0$ sau khi hồi quy. Nếu $\text{VIF} \ge 10$, cần áp dụng phân tích nhân tố trích xuất điểm nhân tố (Factor Scores) đại diện để khử tương quan tuyến tính.

4. Chi phí và nhân lực vận hành khảo sát định lượng nội bộ trong doanh nghiệp SME là bao nhiêu?

Với công cụ khảo sát trực tuyến (Google Forms / MS Forms), chi phí thu thập dữ liệu xấp xỉ bằng không. Toàn bộ quá trình làm sạch, phân tích SPSS và báo cáo có thể được thực hiện bởi 01 chuyên viên phân tích nhân sự trong vòng 3 đến 4 tuần.

5. Kết quả nghiên cứu có thể áp dụng ngay cho các công ty thuộc ngành khác không?

Khung mô hình 6 yếu tố hoàn toàn có thể tái sử dụng. Tuy nhiên, doanh nghiệp cần tiến hành bước nghiên cứu định tính (phỏng vấn thử $n=5-10$ nhân sự) để điều chỉnh câu chữ trong bảng hỏi cho phù hợp với đặc thù văn hóa và quy trình nghiệp vụ riêng biệt của ngành đó.

Kết luận

Khóa luận tốt nghiệp của tác giả Lâm Quang Khánh đã giải quyết trọn vẹn bài toán xác định và định lượng các yếu tố tác động đến hiệu quả làm việc tại Công ty TNHH Thương mại và Dịch vụ Khánh An. Thông qua kiểm định kinh tế lượng vững chắc với mẫu $N=210$, nghiên cứu chỉ ra rằng Phong cách Lãnh đạo ($\beta = 0.312$) và Chính sách Lương - Phúc lợi ($\beta = 0.285$) là hai đòn bẩy quan trọng nhất, trong khi Áp lực công việc ($\beta = -0.118$) cần được kiểm soát chặt chẽ để tránh phản tác dụng.

Đây là cẩm nang thực chứng giá trị giúp các nhà quản trị doanh nghiệp thương mại dịch vụ đưa ra các quyết sách nhân sự chuẩn xác, khoa học và bền vững.