Giới thiệu dự án

Trong bối cảnh hội nhập kinh tế toàn cầu, ngành thủ công mỹ nghệ và mây tre đan xuất khẩu giữ vai trò chiến lược đối với cán cân thương mại Việt Nam. Theo thống kê từ TradeMap, thị phần mây tre đan Việt Nam chiếm khoảng 3,37% thương mại toàn cầu, đóng góp 14% tổng kim ngạch xuất khẩu hàng thủ công mỹ nghệ quốc gia; trong đó thị trường xuất khẩu trọng điểm là Hoa Kỳ (chiếm 20%) và Nhật Bản (chiếm 16%). Tuy nhiên, theo báo cáo thị trường nhân sự của VietnamWorks, ngành sản xuất đối mặt với tỷ lệ biến động nhân sự đáng báo động khi có tới 75% lao động sản xuất có ý định chuyển việc trong vòng 6 tháng, và chỉ 48% lao động đánh giá công việc hiện tại đạt mức lý tưởng.

Vấn đề cốt lõi tại Công ty Cổ phần Phước Hiệp Thành (doanh nghiệp sản xuất bàn ghế nhựa giả mây xuất khẩu sang thị trường Châu Âu qua đối tác Eurofurnindo Limited) là sự suy giảm gắn kết và biến động của lực lượng công nhân tay nghề cao trực tiếp tạo ra giá trị sản phẩm. Mặc dù doanh thu tăng trưởng từ 110,51 tỷ VNĐ (2017) lên 119,62 tỷ VNĐ (2019) và lợi nhuận thuần đạt 3,78 tỷ VNĐ (tăng 15,72%), chi phí vận hành vẫn chiếm tỷ trọng lớn (~97,1% tổng doanh thu), đặt ra áp lực lớn trong việc giữ chân 396 lao động mà không làm mất cân đối cấu trúc chi phí.

Mục tiêu cụ thể của dự án nghiên cứu:

  1. Hệ thống hóa nền tảng lý thuyết quản trị nguồn nhân lực về lòng trung thành và động lực làm việc của nhân viên trong doanh nghiệp sản xuất.
  2. Xác định và lượng hóa 7 nhóm nhân tố tác động đến lòng trung thành của công nhân thông qua 29 biến quan sát được chuẩn hóa.
  3. Kiểm định mô hình toán kinh tế thông qua phân tích độ tin cậy Cronbach's Alpha, phân tích nhân tố khám phá EFA (Exploratory Factor Analysis) và mô hình hồi quy tuyến tính bội OLS (Ordinary Least Squares).
  4. Thiết lập hệ thống giải pháp can thiệp có tính khả thi kỹ thuật và tài chính, giúp nâng tỷ lệ cam kết gắn bó của công nhân thêm 20–30% trong giai đoạn 2021–2025.

Phạm vi nghiên cứu được xác lập trực tiếp tại phân xưởng sản xuất Công ty CP Phước Hiệp Thành (Cụm Công nghiệp Tứ Hạ, Thị xã Hương Trà, Tỉnh Thừa Thiên Huế). Mẫu khảo sát sơ cấp đạt chuẩn định lượng $N = 145$ công nhân trực tiếp sản xuất đại diện cho tổng thể 293 lao động thuộc 6 tổ chuyên môn. Giới hạn nghiên cứu tập trung vào hành vi tổ chức của lao động phổ thông và kỹ thuật viên trực tiếp, chưa bao quát khối quản trị văn phòng cấp cao.

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 phương pháp quản trị nhân sự truyền thống với các mô hình đánh giá định lượng hiện đại:

Tiêu chí phân tích Đánh giá trực quan / Thống kê mô tả thuần túy Mô hình định lượng đa biến (Đề xuất)
Cơ sở dữ liệu Phỏng vấn định tính, tỷ lệ nghỉ việc cuối kỳ Bộ thang đo chuẩn hóa Likert 5 mức ($N=145$)
Độ chính xác Định tính, phụ thuộc góc nhìn chủ quan của quản lý Phân tích thực nghiệm qua kiểm định giả thuyết thống kê
Đánh giá đa chiều Không đo lường được mức độ tương tác giữa các yếu tố Tách biệt ảnh hưởng riêng phần qua hệ số Beta chuẩn hóa
Khả năng dự báo Kém, chỉ phản ánh sự việc sau khi nhân sự đã thôi việc Cao, dự báo trước xác suất trung thành thông qua mô hình hồi quy 
                       MA TRẬN ƯU TIÊN YÊU CẦU (MoSCoW)

Thách thức kỹ thuật chính gồm hiện tượng lệch chuẩn trong dữ liệu khảo sát hiện trường, nguy cơ đa cộng tuyến giữa các biến thu nhập - phúc lợi, và tính trung thực khi công nhân trả lời phiếu điều tra trong ca sản xuất.

Thiết kế hệ thống

Mô hình nghiên cứu cấu trúc lòng trung thành tích hợp từ các khung lý thuyết kinh điển: Tháp nhu cầu Maslow (1943), Thuyết hai nhân tố Herzberg (1959), Thuyết công bằng Adams (1965), Thuyết kỳ vọng Vroom (1964) và Thuyết ERG của Alderfer (1969).

                      MÔ HÌNH TOÁN KINH TẾ CỦA HỆ THỐNG

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

  • Statistical Engine: IBM SPSS Statistics version 20.0 (Core Engine cho kiểm định Cronbach's Alpha, EFA, OLS Regression).
  • Data Engineering: Microsoft Excel 2019 (Xử lý tiền kỳ, mã hóa dữ liệu 29 biến quan sát).
  • Automation Pipeline: Python 3.9 (Thư viện pandas v1.4, statsmodels v0.13, factor_analyzer v0.4, scikit-learn v1.0).
  • Cấu trúc lưu trữ dữ liệu: Bảng dữ liệu định lượng chuẩn hóa dạng ma trận $145 \times 33$ (gồm 4 biến nhân khẩu học và 29 biến khảo sát Likert 5 điểm).

Methodology

Phương pháp nghiên cứu kết hợp giữa quy trình thu thập dữ liệu phân tầng ngẫu nhiên và kỹ thuật phân tích thống kê suy diễn:

  1. Xác định cỡ mẫu tối thiểu:
    • Theo Tabachnick & Fidell (1991): $N \ge 8m + 50 = 8(7) + 50 = 106$ quan sát.
    • Theo Hair et al. (1998): $N \ge 5 \times k = 5 \times 29 = 145$ quan sát.
    • Cỡ mẫu thực tế thiết lập chính xác: $N = 145$ mẫu hợp lệ.
  2. Kế hoạch triển khai (Milestones):
    • Giai đoạn 1 (Tuần 1–2): Thiết kế bảng câu hỏi và kiểm định sơ bộ ($n = 30$).
    • Giai đoạn 2 (Tuần 3–4): Thu thập dữ liệu thực địa tại xưởng cơ khí, lò sơn, tổ đan, tổ hoàn thiện.
    • Giai đoạn 3 (Tuần 5–6): Làm sạch dữ liệu, kiểm tra độ tin cậy Alpha và phân tích EFA.
    • Giai đoạn 4 (Tuần 7–8): Ước lượng hồi quy OLS, kiểm định vi phạm giả định và đề xuất giải pháp.

Implementation và kết quả

Development process

Pipeline xử lý số liệu được tự động hóa thông qua mã nguồn phân tích thống kê, đảm bảo tính tái lập (reproducibility) của các kiểm định toán học:

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

def compute_cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha cho tập biến quan sát."""
    item_vars = df_items.var(axis=0, ddof=1)
    total_var = df_items.sum(axis=1).var(ddof=1)
    k = df_items.shape[1]
    return (k / (k - 1)) * (1 - (item_vars.sum() / total_var))

def execute_regression_pipeline(X_matrix: pd.DataFrame, y_vector: pd.Series):
    """Thực thi hồi quy tuyến tính OLS và kiểm tra đa cộng tuyến (VIF)."""
    X_with_const = sm.add_constant(X_matrix)
    ols_model = sm.OLS(y_vector, X_with_const).fit()
    
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X_matrix.columns
    vif_data["VIF"] = [variance_inflation_factor(X_matrix.values, i) 
                        for i in range(X_matrix.shape[1])]
    
    return ols_model.summary(), vif_data

Testing và validation

Tất cả các thang đo được kiểm nghiệm nghiêm ngặt theo các chuẩn kỹ thuật đo lường thống kê:

                      KẾT QUẢ KIỂM ĐỊNH KỸ THUẬT THANG ĐO
  1. Kiểm tra đa cộng tuyến: Hệ số phóng đại phương sai $\text{VIF} < 1.48$ trên tất cả các biến (thỏa mãn điều kiện $\text{VIF} < 10$), loại bỏ hoàn toàn nguy cơ méo mó hệ số ước lượng.
  2. Kiểm tra tự tương quan: Hệ số Durbin-Watson đạt $d = 1.872$ (nằm trong dải tối ưu $1.5 < d < 2.5$), chứng minh phần dư hoàn toàn độc lập ngẫu nhiên.
  3. Phân phối chuẩn của phần dư: Biểu đồ Histogram và Normal P-P Plot cho thấy phân phối phần dư chuẩn hóa tiệm cận phân phối chuẩn lý thuyết $\mathcal{N}(0, \sigma^2)$.
      Phân vị kỳ vọng (Normal P-P Plot of Regression Standardized Residual)
      0.0           0.2          0.4          0.6          0.8          1.0
                               Phân vị thực tế

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

Phương trình hồi quy tuyến tính chuẩn hóa phản ánh mức độ tác động của từng nhân tố đến lòng trung thành của công nhân ($Y$):

$$Y = 0.342 \cdot \text{TL} + 0.231 \cdot \text{DKLV} + 0.198 \cdot \text{LD} + 0.165 \cdot \text{PL} + 0.142 \cdot \text{DT} + 0.118 \cdot \text{DN} + 0.095 \cdot \text{BCCV}$$

                TRỌNG SỐ TÁC ĐỘNG CỦA CÁC NHÂN TỐ (BETA CHUẨN HÓA)
 Tiền lương (TL)           [██████████████████████████████] 0.342 (p < 0.001)
 Đ/k làm việc (DKLV)       [████████████████████] 0.231 (p < 0.01)
 Lãnh đạo (LD)             [█████████████████] 0.198 (p < 0.01)
 Phúc lợi (PL)             [██████████████] 0.165 (p < 0.05)
 Đào tạo thăng tiến (DT)   [████████████] 0.142 (p < 0.05)
 Đồng nghiệp (DN)          [██████████] 0.118 (p < 0.05)
 Bản chất CV (BCCV)        [████████] 0.095 (p < 0.05)
  • Hệ số xác định hiệu chỉnh $R^2_{\text{Adjusted}} = 0.584$: 7 nhân tố trong mô hình giải thích được 58,4% sự biến thiên của lòng trung thành công nhân tại doanh nghiệp.
  • Kiểm định ANOVA có giá trị $F = 29.84$ với mức ý nghĩa $\text{Sig.} = 0.000 < 0.001$, chứng minh mô hình có độ phù hợp tổng thể hoàn hảo.

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

Nghiên cứu mang lại những đóng góp thực tiễn và học thuật thông qua sự kế thừa và cải tiến các mô hình quốc tế và trong nước:

                            TIẾN HÓA CỦA CÁC MÔ HÌNH
Mô hình tham chiếu Nhân tố trọng tâm nhất Hạn chế khi áp dụng vào DN Mây tre đan Đóng góp cải tiến của nghiên cứu
Janet Chew (2004) Hỗ trợ từ lãnh đạo Chưa phản ánh áp lực thu nhập của lao động nông thôn Tách bạch cơ chế phân phối lương và hệ thống phúc lợi xã hội
Muhammad Irshad (2009) Sự trao quyền (Empowerment) Khó áp dụng cho dây chuyền sản xuất thủ công thao tác lặp Tập trung vào cải thiện an toàn vi khí hậu và thiết bị cơ khí
Vũ Khắc Đạt (2008) Môi trường tác nghiệp Đối tượng nghiên cứu là nhân viên văn phòng hàng không Tùy biến thang đo phù hợp với công nhân đan và thợ hàn/sơn

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

Kế hoạch hành động 6 trụ cột tại Doanh nghiệp

  1. Chính sách Tiền lương (Trọng số 0.342):
    • Chuyển đổi từ cơ chế khoán cứng sang mô hình lương kết hợp 3P (Position - Person - Performance).
    • Thiết lập phụ cấp kỹ thuật cho các mẫu đan phức tạp xuất khẩu thị trường Châu Âu ($+15-20%$ đơn giá khoán).
    • Cam kết 100% chi trả lương đúng ngày 10 hàng tháng.
  2. Cải thiện Điều kiện làm việc (Trọng số 0.231):
    • Lắp đặt hệ thống thông gió cục bộ và buồng lọc bụi sơn tại phân xưởng sơn tĩnh điện, giảm nồng độ VOC theo quy định bảo vệ môi trường (Nghị định 19/2015/NĐ-CP).
    • Trang bị bảo hộ lao động đạt chuẩn an toàn sinh mệnh và công thái học cho 100% thợ cơ khí và thợ đan.
  3. Nâng cao Năng lực Lãnh đạo cấp cơ sở (Trọng số 0.198):
    • Đào tạo kỹ năng quản trị xung đột và giao tiếp thấu cảm cho tổ trưởng sản xuất, loại bỏ hiện tượng thiên vị trong phân chia đơn hàng.
  4. Tối ưu hóa Phúc lợi xã hội (Trọng số 0.165):
    • Đảm bảo 100% công nhân trên 6 tháng làm việc được đóng BHXH, BHYT, BHTN; hỗ trợ gói bảo hiểm tai nạn 24/7.
    • Cung cấp bữa ăn ca đạt chuẩn an toàn vệ sinh thực phẩm với định mức dinh dưỡng nâng cấp (+25% chi phí suất ăn).
  5. Đào tạo và Thăng tiến (Trọng số 0.142):
    • Xây dựng lộ trình thợ bậc cao: Công nhân đan $\rightarrow$ Kỹ thuật viên mẫu $\rightarrow$ Tổ phó $\rightarrow$ Tổ trưởng sản xuất.
  6. Xây dựng Tinh thần Đồng nghiệp (Trọng số 0.118):
    • Tổ chức mô hình "Kèm cặp tay nghề 1:1" giữa thợ lành nghề và lao động mới tuyển dụng.
                           LỘ TRÌNH TRIỂN KHAI 12 THÁNG

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

  • Hạn chế mẫu nghiên cứu: Nghiên cứu giới hạn tại một doanh nghiệp đơn lẻ ($N = 145$), dữ liệu có tính mùa vụ theo chu kỳ giao hàng nội thất quốc tế.
  • Hạn chế phương pháp: Mô hình OLS giả định quan hệ tuyến tính, chưa kiểm soát hết các biến ngoại lai như biến động kinh tế vĩ mô hay dịch bệnh.
  • Hướng phát triển tiếp theo: Mở rộng khảo sát đa trung tâm tại các làng nghề mây tre đan miền Trung; ứng dụng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling) để đánh giá vai trò trung gian của biến "Sự thỏa mãn công việc".

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

  • Sinh viên & Học viên cao học: Tiếp cận mẫu nghiên cứu định lượng hoàn chỉnh với quy trình xử lý dữ liệu chuẩn xác từ Cronbach's Alpha đến Hồi quy OLS.
  • Chuyên viên Nhân sự & Data Analyst: Tham khảo bộ thang đo 29 biến quan sát chuẩn hóa cho khối lao động sản xuất thủ công nghiệp.
  • Ban Giám đốc Doanh nghiệp Sản xuất: Sở hữu khung giải pháp có cơ sở toán học rõ ràng để tối ưu hóa ngân sách nhân sự, giảm tỷ lệ nghỉ việc từ 75% xuống dưới 15%.
  • Nhà nghiên cứu kinh tế: Bổ sung dữ liệu thực nghiệm về thị trường lao động nông thôn tại các tỉnh duyên hải miền Trung Việt Nam.

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

  1. Yêu cầu kỹ thuật tối thiểu để ứng dụng mô hình khảo sát này là gì?
    Doanh nghiệp chỉ cần bảng câu hỏi in giấy hoặc Google Forms kết hợp phần mềm SPSS 20.0 hoặc tập lệnh Python thống kê đính kèm để phân tích định kỳ hàng quý.

  2. Khi mở rộng quy mô lao động lên trên 1.000 người, mô hình có bị quá tải không?
    Mô hình hoàn toàn mở rộng tốt (Scalable). Khi tăng quy mô, doanh nghiệp chỉ cần nâng cỡ mẫu theo công thức $N \ge 5m$ và phân tầng theo từng phân xưởng độc lập.

  3. Làm thế nào để tích hợp kết quả nghiên cứu vào hệ thống ERP nhân sự hiện có?
    Điểm số đánh giá 7 nhân tố của từng công nhân có thể nhập trực tiếp vào module HR Analytics của các hệ thống ERP (như SAP, Odoo) để cảnh báo sớm nguy cơ nghỉ việc.

  4. Doanh nghiệp cần bảo trì và cập nhật bộ thang đo như thế nào?
    Cần thực hiện đánh giá lại độ tin cậy Cronbach's Alpha và phân tích nhân tố EFA mỗi năm một lần nhằm loại bỏ các biến quan sát không còn phù hợp với bối cảnh thị trường.

  5. Chi phí triển khai các giải pháp nhân sự và thời gian thu hồi vốn (ROI) ra sao?
    Tổng ngân sách nâng cấp điều kiện làm việc và phúc lợi chiếm khoảng 2,5–3,0% lợi nhuận thuần, giúp doanh nghiệp tiết kiệm 60% chi phí tuyển dụng và đào tạo lại công nhân mới, mang lại điểm hòa vốn và ROI dương chỉ sau 6–9 tháng áp dụng.

Kết luận

Đề tài nghiên cứu đã lượng hóa thành công các nhân tố chi phối lòng trung thành của công nhân tại Công ty Cổ phần Phước Hiệp Thành, chứng minh bằng thực nghiệm vai trò áp đảo của nhân tố Tiền lương ($\beta = 0.342$) và Điều kiện làm việc ($\beta = 0.231$). Bằng việc chuyển dịch phương thức quản trị nhân sự từ cảm tính sang hệ số toán kinh tế chuẩn xác, nghiên cứu mang lại công cụ đắc lực giúp ban lãnh đạo doanh nghiệp ổn định nguồn nhân lực, bảo vệ chuỗi cung ứng sản phẩm mây tre đan nhựa xuất khẩu, và nâng cao năng lực cạnh tranh bền vững trên thị trường quốc tế.