Giới thiệu dự án

Trong bối cảnh toàn cầu hóa và cách mạng công nghiệp, ngành dệt may Việt Nam đóng vai trò xương sống với kim ngạch xuất khẩu chiếm trên 12-15% tổng kim ngạch cả nước. Tuy nhiên, đặc thù ngành là thâm dụng lao động lớn với tỷ lệ biến động nhân sự (turnover rate) trung bình dao động từ 15% đến 25%/năm. Sự dịch chuyển lao động liên tục giữa các khu công nghiệp tạo ra áp lực tài chính nghiêm trọng cho doanh nghiệp trong việc tuyển dụng, đào tạo bù đắp và suy giảm năng suất sản xuất.

Công ty May Hòa Thọ - Đông Hà (chi nhánh Tổng Công ty Cổ phần Dệt May Hòa Thọ tại Quảng Trị) sở hữu quy mô hơn 2.060 công nhân viên, tạo doanh thu trên 230 tỷ VNĐ (năm 2019). Dù năng lực sản xuất đạt 5,2 triệu sản phẩm/năm, doanh nghiệp đối mặt với bài toán nan giải: nhóm lao động có thâm niên dưới 6 năm chiếm tới 75,8% (trong đó thâm niên dưới 3 năm chiếm 31,1%), phản ánh tỷ lệ giữ chân lao động chưa bền vững.

+-------------------------------------------------------------------------------+
|                      BÀI TOÁN QUẢN TRỊ NHÂN SỰ TẠI DOANH NGHIỆP               |
|                                                                               |
|  [Quy mô 2.060 nhân sự] ---> [Thâm niên < 6 năm: 75,8%] ---> [Turnover cao]   |
|                                                                               |
|  Thách thức: Chi phí tuyển dụng tăng + Thất thoát kỹ năng + Gián đoạn chuyền  |
+-------------------------------------------------------------------------------+

Mục tiêu nghiên cứu được xác định cụ thể:

  1. Hệ thống hóa cơ sở lý luận về động viên và lòng trung thành nhân viên dựa trên các học thuyết kinh điển: Tháp nhu cầu Maslow, Thuyết hai nhân tố Herzberg, Thuyết kỳ vọng Vroom và Thuyết công bằng Adams.
  2. Xây dựng và kiểm định thang đo định lượng gồm 7 nhân tố độc lập (Lương, Điều kiện làm việc, Đồng nghiệp, Cấp trên, Khen thưởng, Phúc lợi, Cơ hội đào tạo & thăng tiến) tác động đến biến phụ thuộc Lòng trung thành của cán bộ công nhân viên (CBCNV).
  3. Đánh giá sự khác biệt về mức độ trung thành theo các đặc trưng nhân khẩu học (giới tính, độ tuổi, trình độ học vấn, vị trí công tác, thâm niên).
  4. Đề xuất ma trận giải pháp quản trị nguồn nhân lực (HRM) có khả năng nâng cao tỷ lệ gắn kết thêm 20-25% và giảm thiểu 15% tỷ lệ nghỉ việc tự nguyện.

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

  • Không gian: Các phòng ban gián tiếp và 20 chuyền may trực tiếp tại Nhà máy May Hòa Thọ - Đông Hà, Khu công nghiệp Nam Đông Hà, tỉnh Quảng Trị.
  • Thời gian: Khảo sát sơ cấp thu thập từ tháng 01/2020 đến tháng 04/2020 kết hợp phân tích chuỗi dữ liệu thứ cấp giai đoạn 2017–2019.
  • Giới hạn: Dữ liệu cắt ngang (cross-sectional data) với dung lượng mẫu hiệu dụng $n = 190$.

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

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

Nhiều nghiên cứu trước đây tiếp cận lòng trung thành của người lao động dưới các góc độ khác nhau nhưng bộc lộ nhiều khoảng trống khi áp dụng vào môi trường sản xuất thâm dụng lao động phổ thông tại miền Trung.

Tiêu chí phân tích Mô hình Trần Kim Dung (2003, 2005) Mô hình Vũ Khắc Đạt (2008) Mô hình nghiên cứu Hòa Thọ (Đề xuất)
Đối tượng khảo sát Nhân viên văn phòng đa ngành tại TP.HCM Nhân viên Vietnam Airlines Miền Nam Công nhân may trực tiếp & văn phòng xí nghiệp
Cấu trúc biến số Gộp chung chính sách đãi ngộ và môi trường Gộp Lương + Phúc lợi; Gộp Đồng nghiệp + ĐK làm việc Tách biệt độc lập 7 nhân tố: Lương, ĐKLV, Đồng nghiệp, Cấp trên, Khen thưởng, Phúc lợi, Thăng tiến
Độ nhạy phân tầng Thấp với lao động chân tay Trung bình (ngành dịch vụ hàng không) Cao (phân bổ tỷ lệ 82,1% công nhân may : 17,9% văn phòng)
Hạn chế kỹ thuật Tồn tại nghịch lý tương quan âm giữa lương và nỗ lực Thang đo gộp làm mờ tác động đơn lẻ của chính sách phúc lợi Cần loại bỏ biến rác có hệ số tải chéo trước khi hồi quy

Phân tích yêu cầu theo khung MoSCoW:

  • Must-have (Bắt buộc): Đo lường tính chuẩn xác của hệ số Cronbach's Alpha ($\alpha \ge 0.70$), hệ số KMO ($\ge 0.50$), tổng phương sai trích ($TVE \ge 50%$) và mức ý nghĩa thống kê ($Sig. < 0.05$).
  • Should-have (Nên có): Kiểm định ANOVA và Independent-Samples T-Test để phát hiện sai biệt giữa các nhóm độ tuổi và thâm niên.
  • Could-have (Có thể có): Mở rộng đánh giá tương quan giữa thu nhập bình quân (tăng từ 5,8 triệu lên 7,3 triệu VNĐ/tháng) với chỉ số thỏa mãn cá nhân.
  • Won't-have (Chưa thực hiện): Mô hình hóa cấu trúc tuyến tính SEM (Structural Equation Modeling) đa biến trung gian trong khuôn khổ khóa luận này.

Thiết kế hệ thống nghiên cứu

Mô hình nghiên cứu được xây dựng dưới dạng hệ thống xử lý dữ liệu thống kê đa tầng, chuyển hóa phản hồi định tính qua thang đo Likert 5 điểm thành các vector trọng số định lượng.

Technology Stack và công cụ phân tích:

  • Phần mềm chính: IBM SPSS Statistics Version 20.0 (Core Engine).
  • Ngôn ngữ bổ trợ mô phỏng: Python 3.10 với các thư viện pandas, statsmodels, factor_analyzer, scikit-learn.
  • Thang đo: Likert 5 mức độ (1: Rất không đồng ý $\rightarrow$ 5: Rất đồng ý).
  • Phương trình hồi quy tổng quát: $$LOY = \beta_0 + \beta_1 LUONG + \beta_2 DKLV + \beta_3 DONGNGHIEP + \beta_4 CAPTREN + \beta_5 KHENTHUONG + \beta_6 PHUCLOI + \beta_7 COHOITHANGTIEN + \epsilon$$

Methodology

Nghiên cứu áp dụng quy trình 2 giai đoạn chuẩn mực:

  1. Giai đoạn định tính: Nghiên cứu tài liệu kết hợp phỏng vấn sâu chuyên gia quản trị nhân sự tại Tổng Công ty Dệt May Hòa Thọ để hiệu chỉnh danh mục 38 biến quan sát ban đầu.
  2. Giai đoạn định lượng: Chọn mẫu phân tầng tỷ lệ (Stratified Random Sampling) trên tổng thể $N = 2.060$ nhân viên. Cỡ mẫu tối thiểu theo Hair et al. (1998) là $38 \times 5 = 190$ mẫu. Thực tế phát 210 phiếu, thu về 190 phiếu hợp lệ (đạt tỷ lệ 90,48%).
+-----------------------------------------------------------------------------------+
|                            MA TRẬN QUẢN LÝ RỦI RO NGHIÊN CỨU                     |
+--------------------------+------------+-----------+-------------------------------+
| Rủi ro tiềm ẩn           | Mức độ     | Tác động  | Giải pháp xử lý               |
+--------------------------+------------+-----------+-------------------------------+
| Sai lệch phân bố mẫu     | Trung bình | Nghiêm trọng | Lấy mẫu phân tầng theo đúng tỷ|
|                          |            |           | lệ xưởng sản xuất (82,1%)     |
| Thiên lệch xã hội (SDB)  | Cao        | Vừa       | Ẩn danh hoàn toàn bảng hỏi    |
| Tải chéo giữa các nhân tố| Cao        | Nghiêm trọng | Áp dụng phép xoay Varimax và   |
|                          |            |           | loại biến có Factor Load < 0.5|
+--------------------------+------------+-----------+-------------------------------+

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu được thực thi tuần tự qua các sprint phân tích số liệu:

[Sprint 1: Thu thập & Mã hóa] 
[Sprint 2: Kiểm định Độ tin cậy Alpha Lần 1] 
[Sprint 3: Khám phá Nhân tố EFA] 
[Sprint 4: Hồi quy & Phân tích Đa biến]

Đoạn mã Python mô phỏng logic xử lý Cronbach's Alpha và trích xuất nhân tố tương đương SPSS:

import numpy as np
import pandas as pd
from factor_analyzer import FactorAnalyzer
from factor_analyzer.factor_analyzer import calculate_kmo, calculate_bartlett_sphericity

def cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số tin cậy Cronbach's Alpha chuẩn xác."""
    item_vars = df_items.var(axis=0, ddof=1)
    total_var = df_items.sum(axis=1).var(ddof=1)
    n_items = df_items.shape[1]
    return (n_items / (n_items - 1)) * (1 - (item_vars.sum() / total_var))

def evaluate_factor_suitability(df_matrix: pd.DataFrame):
    """Kiểm tra điều kiện KMO và Bartlett's Test of Sphericity."""
    chi_square_value, p_value = calculate_bartlett_sphericity(df_matrix)
    kmo_all, kmo_model = calculate_kmo(df_matrix)
    return {
        "KMO": round(kmo_model, 3),
        "Bartlett_ChiSq": round(chi_square_value, 2),
        "Bartlett_pValue": p_value
    }

Testing và validation

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

Toàn bộ các nhóm thang đo ban đầu đều đạt độ tin cậy cao ($\alpha > 0.75$), chứng minh tính nhất quán nội tại của bảng hỏi:

Nhóm nhân tố Số biến ban đầu Cronbach's Alpha Lần 1 Biến bị loại do Item-Total < 0.3 Cronbach's Alpha Sau lọc
Lương 5 0.946 Không 0.946
Điều kiện làm việc 4 0.820 dklv2 0.808
Đồng nghiệp 3 0.753 Không 0.753
Cấp trên 4 0.806 Không 0.806
Khen thưởng 3 0.910 Không 0.910
Phúc lợi 5 0.800 phucloi1 0.795
Cơ hội thăng tiến 6 0.826 cohoithangtien1 0.812

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

  • EFA Lần 1: Phát hiện các biến cohoithangtien1, dklv2, phucloi1 có hệ số tương quan biến-tổng hoặc hệ số tải (Factor Loading) $< 0.30$, tiến hành loại bỏ.
  • EFA Lần 2 (26 biến quan sát còn lại):
    • Hệ số Kaiser-Meyer-Olkin (KMO) $= 0.824 > 0.50$.
    • Kiểm định Bartlett's Test of Sphericity: $\chi^2 = 2451.02$, $df = 351$, $Sig. = 0.000 < 0.05$ (bác bỏ giả thuyết $H_0$ về ma trận tương quan đơn vị).
    • Trích xuất 7 nhân tố có Eigenvalue $> 1.0$.
    • Tổng phương sai trích (Total Variance Explained) $= 70,131% > 50%$ (mô hình giải thích được 70,131% biến thiên của tập dữ liệu).
+-------------------------------------------------------------------------------+
|                      KẾT QUẢ PHÂN TÍCH EFA LẦN 2 (SPSS v20)                   |
|                                                                               |
|  KMO Measure of Sampling Adequacy : 0.824                                     |
|  Bartlett's Test Sig.             : 0.000 (p < 0.001)                         |
|  Total Variance Explained         : 70.131%                                   |
|  Số lượng nhân tố cô đọng         : 7 nhóm nhân tố sạch                       |
+-------------------------------------------------------------------------------+

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

Mô hình hồi quy tuyến tính OLS xác lập mức độ tác động có ý nghĩa thống kê ($p < 0.05$) của các nhân tố đến lòng trung thành của công nhân viên:

  1. Yếu tố Lương và Thu nhập: Có hệ số Beta chuẩn hóa cao nhất, đóng vai trò nền tảng sống còn đối với công nhân sản xuất trực tiếp.
  2. Yếu tố Khen thưởng & Phúc lợi: Đóng góp đáng kể vào động lực thúc đẩy gắn kết lâu dài.
  3. Kiểm định ANOVA nhân khẩu học:
    • Độ tuổi: Sự gắn kết ở nhóm tuổi 26–35 cao hơn rõ rệt nhóm $< 25$ tuổi ($F = 4.12, Sig. < 0.05$).
    • Thâm niên: Khác biệt có ý nghĩa thống kê giữa nhóm thâm niên $> 6$ năm so với nhóm $< 3$ năm ($F = 5.89, Sig. = 0.003$).

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

Nghiên cứu đem lại giá trị học thuật và thực tiễn sâu sắc thông qua các đóng góp vượt trội:

                  SO SÁNH MÔ HÌNH VÀ CẢI TIẾN PHƯƠNG PHÁP
  1. Tinh chỉnh thang đo cho công nhân may mặc miền Trung: Thay vì áp dụng nguyên mẫu các bảng hỏi nước ngoài, đề tài đã loại bỏ 3 biến không phù hợp ngữ cảnh văn hóa nhà xưởng (cohoithangtien1, dklv2, phucloi1), gia tăng độ chuẩn xác của phương sai trích lên $70,131%$.
  2. Làm rõ ranh giới giữa Khen thưởng và Phúc lợi: Chứng minh rằng trong môi trường công nghiệp dệt may, phúc lợi (bảo hiểm, ăn ca, trợ cấp lễ tết) mang tính chất duy trì (Herzberg Hygiene Factor), trong khi khen thưởng năng suất trực tiếp theo sản lượng tạo động lực gắn kết mạnh mẽ (Motivator Factor).
  3. Cơ sở định lượng cho tái cấu trúc nhân sự: Cung cấp cho Ban Giám đốc Công ty May Hòa Thọ luận cứ khoa học để xây dựng quỹ lương thưởng gắn liền KPIs thay vì phương thức cào bằng truyền thống.

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

Kế hoạch ứng dụng (Implementation Roadmap)

Chính sách HRM mới tại Công ty May Hòa Thọ - Đông Hà được thiết kế theo lộ trình 3 giai đoạn:

[Giai đoạn 1: Tháng 1 - 3]  
[Giai đoạn 2: Tháng 4 - 6]  
[Giai đoạn 3: Tháng 7 - 12]

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

+------------------------------------------------------------------------------------+
|                         BẢNG ƯỚC TÍNH HIỆU QUẢ TÀI CHÍNH (ROI)                     |
+------------------------------------+-----------------------------------------------+
| Hạng mục đầu tư                    | Chi phí dự toán / Lợi ích quy đổi             |
+------------------------------------+-----------------------------------------------+
| Cải thiện bữa ăn ca & phúc lợi     | 350.000.000 VNĐ / năm                         |
| Quỹ khen thưởng năng suất chuyền   | 250.000.000 VNĐ / năm                         |
| TỔNG CHI PHÍ ĐẦU TƯ THỰC HIỆN      | 600.000.000 VNĐ                               |
+------------------------------------+-----------------------------------------------+
| Giảm tỷ lệ nghỉ việc (từ 12%->6%)  | Tiết kiệm 450.000.000 VNĐ chi phí tuyển dụng  |
| Giảm thời gian trễ đơn hàng (Final)| Tăng doanh thu đóng góp ~900.000.000 VNĐ/năm  |
| LỢI ÍCH TỔNG HỢP THU VỀ            | 1.350.000.000 VNĐ                             |
+------------------------------------+-----------------------------------------------+
| ROI DỰ KIẾN NĂM ĐẦU TIÊN           | ((1.350M - 600M) / 600M) * 100% = 125%        |
+------------------------------------+-----------------------------------------------+

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

  • Hạn chế mẫu nghiên cứu: Mẫu $n = 190$ tập trung chủ yếu tại Chi nhánh Hòa Thọ Đông Hà (Quảng Trị), chưa bao quát toàn bộ các đơn vị thành viên khác trong Tổng Công ty Dệt May Hòa Thọ (Đà Nẵng, Quảng Nam).
  • Thiên lệch dữ liệu tự báo cáo (Self-report bias): Một số công nhân e ngại khi đánh giá về biến Cấp trên (Captren) dù bảng hỏi đã được bảo mật danh tính.
  • Hướng phát triển tương lai:
    • Áp 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) trên tập mẫu $N > 500$.
    • Bổ sung biến trung gian: Sự thỏa mãn công việc (Job Satisfaction) và Văn hóa tổ chức (Organizational Culture) để khảo sát tương quan gián tiếp đến Lòng trung thành.

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

                   CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI TRỰC TIẾP
Nhóm đối tượng Lợi ích định lượng & Ứng dụng cụ thể
Sinh viên ngành QTKD / Kinh tế Nắm vững quy trình xử lý dữ liệu thực chiến: Kiểm định Cronbach's Alpha $\rightarrow$ EFA 2 bước $\rightarrow$ Hồi quy OLS $\rightarrow$ ANOVA.
Doanh nghiệp Dệt may & Sản xuất Cắt giảm ngay 15-20% chi phí tuyển mộ thay thế nhờ ứng dụng trực tiếp ma trận chính sách Lương - Thưởng - Phúc lợi.
Chuyên viên Nhân sự / HR Analyst Sử dụng bộ thang đo 26 biến đã chuẩn hóa làm công cụ đo lường mức độ gắn kết định kỳ hàng quý (Employee Engagement Survey).
Giảng viên & Nhà nghiên cứu Bổ sung dữ liệu thực nghiệm về thị trường lao động công nghiệp vùng Bắc Trung Bộ vào hệ thống nghiên cứu hành vi tổ chức.

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai bộ phân tích dữ liệu này là gì?

Máy trạm hoặc máy tính cá nhân cấu hình tối thiểu CPU Intel Core i3 / RAM 4GB, cài đặt hệ điều hành Windows 10/11 hoặc Linux, phần mềm IBM SPSS Statistics Version 20.0 trở lên hoặc Python 3.8+ với các thư viện khoa học dữ liệu (pandas, scipy, statsmodels).

2. Tại sao phải thực hiện phân tích EFA lần 2 sau khi đã chạy lần 1?

Trong EFA lần 1, có 3 biến (cohoithangtien1, dklv2, phucloi1) có hệ số tải nhân tố nhỏ hơn 0.3 hoặc tải chéo không đạt giá trị phân biệt. Việc loại bỏ các biến này và thực hiện EFA lần 2 là bắt buộc để đảm bảo cấu trúc dữ liệu sạch, giúp tổng phương sai trích đạt chuẩn $70,131%$.

3. Có thể tích hợp bộ chỉ số này vào phần mềm HRM hiện có của doanh nghiệp không?

Hoàn toàn khả thi. Bộ câu hỏi 26 biến quan sát có thể được cấu hình trực tiếp vào các hệ thống ERP nhân sự (như Odoo, SAP SuccessFactors, Base.vn) dưới dạng khảo sát e-NPS định kỳ, tự động tính toán điểm trung bình từng nhân tố để cảnh báo sớm nguy cơ nghỉ việc.

4. Chi phí triển khai các giải pháp nâng cao lòng trung thành có quá cao đối với các xí nghiệp may nhỏ?

Không. Các giải pháp phi tài chính như: cải thiện phong cách quản lý của tổ trưởng sản xuất (Captren), minh bạch hóa tiêu chí bình xét thi đua (Khenthuong), và nâng cao an toàn vệ sinh lao động (Dklv) hoàn toàn không tốn kém chi phí đầu tư lớn nhưng mang lại mức tăng thỏa mãn tức thì.

5. Bao lâu sau khi áp dụng chính sách mới thì chỉ số lòng trung thành cải thiện?

Theo lộ trình thực nghiệm, các chỉ số thỏa mãn về lương và điều kiện làm việc ghi nhận chuyển biến tích cực sau 1 chu kỳ lương (1 tháng), trong khi chỉ số trung thành gắn bó dài hạn cần từ 6 đến 12 tháng để phản ánh rõ rệt trên tỷ lệ giảm đơn xin thôi việc.

Kết luận

Khóa luận tốt nghiệp của sinh viên Nguyễn Thị Thanh Phong (Đại học Kinh tế Huế) là công trình nghiên cứu ứng dụng công phu, kết hợp chặt chẽ giữa nền tảng học thuyết quản trị nhân sự kinh điển và kỹ thuật phân tích định lượng chuẩn mực trên SPSS v20.0. Đề tài không chỉ bóc tách thành công 7 nhân tố ảnh hưởng then chốt đến lòng trung thành của 2.060 công nhân viên tại Công ty May Hòa Thọ - Đông Hà với phương sai trích $70,131%$, mà còn mở ra định hướng thực thi thiết thực nhằm ổn định nguồn nhân lực, gia tăng năng suất và phát triển bền vững cho ngành dệt may Việt Nam trong thời kỳ hội nhập.