Giới thiệu dự án

Bối cảnh và thực trạng ngành

Trong nền kinh tế sản xuất chế biến gỗ xuất khẩu tại Việt Nam, nguồn nhân lực đóng vai trò là tài sản chiến lược cốt lõi quyết định trực tiếp đến năng suất, chất lượng thành phẩm và lợi thế cạnh tranh bền vững. Theo số liệu khảo sát thị trường lao động từ Anphabe, tỷ lệ nghỉ việc bình quân toàn ngành sản xuất có xu hướng tăng cao, các doanh nghiệp đối mặt với nguy cơ thất thoát đến 95% nhân tài tiềm năng, 51% người lao động có xu hướng "nhảy việc" và 31% không thực sự gắn kết với tổ chức. Cùng với đó, theo báo cáo từ Tổng Liên đoàn Lao động Việt Nam, các tranh chấp lao động và ngừng việc tập thể xuất phát chủ yếu từ chế độ tiền lương, bảo hiểm xã hội (BHXH), môi trường làm việc và lộ trình thăng tiến thiếu minh bạch.

Tại Công ty TNHH Chế biến Gỗ Minh An (Thừa Thiên Huế) – doanh nghiệp chuyên gia công chế biến gỗ đạt chứng chỉ FSC 100% xuất khẩu sang thị trường Châu Âu và Châu Mỹ – quy mô lao động có sự biến động giảm mạnh từ 201 người (năm 2018) xuống 187 người (năm 2019) và còn 125 người (năm 2020), tương ứng mức giảm 33,16%. Doanh nghiệp đứng trước thách thức kép: sự cạnh tranh nguồn lao động gay gắt từ các khu công nghiệp (KCN Phú Bài, KCN Phú Đa) và áp lực duy trì chuỗi cung ứng đạt chuẩn quốc tế trong bối cảnh thị trường biến động sau đại dịch COVID-19 (doanh thu năm 2020 giảm 21,42%, lợi nhuận sau thuế giảm 42,38%).

graph TD
    A[Áp lực thị trường & COVID-19] --> B(Biến động doanh thu & Cắt giảm quy mô)
    C[Cạnh tranh lao động tại KCN] --> D(Tỷ lệ chuyển dịch lao động cao)
    B --> E{Bài toán Nhân sự Minh An}
    D --> E
    E --> F[Đánh giá sự hài lòng công việc]
    F --> G[Mô hình định lượng EFA & Hồi quy OLS]
    G --> H[Chính sách giữ chân & Tối ưu năng suất]

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

Doanh nghiệp chế biến gỗ Minh An đối mặt với các điểm nghẽn (pain points) thực tế:

  • Tỷ lệ biến động nhân sự cao: Lao động dịch chuyển sang các đơn vị có chính sách phúc lợi cạnh tranh hơn, đặc biệt ở nhóm công nhân kỹ thuật và lao động trẻ.
  • Cơ cấu lao động mất cân đối: Lao động nữ chiếm 74,4%, lao động trực tiếp chiếm 95,2%, độ tuổi từ 40–50 tuổi chiếm đa số (54,4%) trong khi mức thu nhập từ 3–5 triệu đồng chiếm tới 79,2%.
  • Thiếu hụt khung đánh giá khoa học: Ban lãnh đạo chưa lượng hóa được mức độ tác động cụ thể của từng yếu tố (điều kiện làm việc, tiền lương, đào tạo thăng tiến, mối quan hệ đồng nghiệp, bản chất công việc) tới lòng trung thành và cam kết gắn bó của người lao động.

Mục tiêu dự án

  1. Hệ thống hóa cơ sở lý luận: Tổng hợp các lý thuyết kinh điển (Tháp nhu cầu Maslow, Thuyết hai nhân tố Herzberg, Thuyết kỳ vọng Vroom, Thuyết nhu cầu McClelland) và mô hình chỉ số mô tả công việc (JDI, MSQ, AJDI).
  2. Xác định và lượng hóa các nhân tố ảnh hưởng: Khảo sát toàn bộ mẫu quan sát $N = 125$ người lao động; ứng dụng kiểm định độ tin cậy Cronbach’s Alpha, phân tích nhân tố khám phá (EFA) và mô hình hồi quy tuyến tính bội trên phần mềm IBM SPSS Statistics 20.0.
  3. Đề xuất hệ thống giải pháp thực tiễn: Xây dựng lộ trình cải tiến chính sách nhân sự giai đoạn 2021–2023 tập trung vào 4 trụ cột: Cải thiện điều kiện vi khí hậu xưởng gỗ, tái cấu trúc thang bảng lương, minh bạch hóa cơ hội thăng tiến và củng cố văn hóa tổ chức.

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

Đề tài tích hợp mô hình định lượng mở rộng từ chỉ số JDI (Job Descriptive Index) của Smith, Kendall & Hulin (1969) và thang đo AJDI điều chỉnh theo bối cảnh Việt Nam của Trần Kim Dung (2005). Nghiên cứu loại bỏ biến "Lãnh đạo" để tập trung chuyên sâu vào 5 nhóm nhân tố độc lập phù hợp với đặc thù sản xuất xưởng gỗ:

  • Bản chất công việc (BCCV - 5 biến quan sát)
  • Điều kiện làm việc (DKLV - 4 biến quan sát)
  • Cơ hội đào tạo và thăng tiến (DTTT - 4 biến quan sát)
  • Tiền lương và phúc lợi (TLPL - 5 biến quan sát)
  • Mối quan hệ với đồng nghiệp (MQH - 4 biến quan sát)
  • Biến phụ thuộc: Sự hài lòng chung (HL - 3 biến quan sát)

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

  • Quy mô mẫu: Thu thập đủ $N = 125$ quan sát hợp lệ (đạt tỷ lệ 100% trên tổng số cán bộ công nhân viên).
  • Chất lượng thang đo: Hệ số Kaiser-Meyer-Olkin (KMO) $\ge 0,70$; Kiểm định Bartlett có giá trị $Sig. \le 0,05$; Tổng phương sai trích (Total Variance Explained) $\ge 50%$; Hệ số tải nhân tố (Factor Loading) $\ge 0,50$.
  • Mô hình hồi quy: Hệ số xác định $R^2$ hiệu chỉnh $\ge 0,50$; kiểm định F có mức ý nghĩa $p < 0,05$; hệ số Durbin-Watson trong khoảng $(1,5 - 2,5)$; hệ số phóng đại phương sai $\text{VIF} < 2,0$ (loại trừ đa cộng tuyến).

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

  • Không gian nghiên cứu: Nhà máy Công ty TNHH Chế biến Gỗ Minh An, số 37 Thuận Hóa, Phường Phú Bài, Thị xã Hương Thủy, Tỉnh Thừa Thiên Huế.
  • Thời gian nghiên cứu: Dữ liệu thứ cấp giai đoạn 2018–2020; dữ liệu sơ cấp thu thập từ tháng 01/2021 đến tháng 03/2021.
  • Giới hạn đề tài: Khảo sát trong phạm vi một doanh nghiệp chế biến gỗ xuất khẩu tại miền Trung; chưa mở rộng so sánh liên ngành với các doanh nghiệp dệt may hoặc cơ khí cùng địa bà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 đo lường sự hài lòng công việc

Tiêu chí Mô hình JDI (Smith et al., 1969) Mô hình MSQ (Weiss et al., 1967) Thang đo AJDI (Trần Kim Dung, 2005) Mô hình đề xuất cho Minh An
Số lượng thành phần 5 nhân tố (Công việc, Lương, Thăng tiến, Lãnh đạo, Đồng nghiệp) 20 khía cạnh (Nội tại, Ngoại tại, Chung) 6 nhân tố (Bổ sung Phúc lợi, Điều kiện làm việc) 5 nhân tố (BCCV, DKLV, DTTT, TLPL, MQH)
Ưu điểm Đơn giản, độ tin cậy cao, chuẩn hóa quốc tế Chi tiết, phân tách bản chất nội tại - ngoại tại Tương thích cao với thị trường lao động Việt Nam Tối ưu hóa cho công nhân sản xuất trực tiếp xưởng gỗ
Nhược điểm Chưa phân tách rõ phúc lợi và điều kiện làm việc thực tế Quá dài (100 câu hỏi), khó triển khai cho công nhân Chưa đặc thù hóa cho ngành chế biến gỗ nặng nhọc Cần kiểm định lại hệ số tải cho từng cụm máy/xưởng
Độ phức tạp thu thập Thấp Rất cao Trung bình Tối ưu (25 biến quan sát)

Phân tích yêu cầu các bên liên quan theo mô hình MoSCoW

  • Must Have (Bắt buộc): Đo lường chính xác tác động của yếu tố Lương - Thưởng, Môi trường an toàn lao động và chế độ BHXH/BHYT bắt buộc.
  • Should Have (Nên có): Phân tích tương quan giữa số năm công tác, độ tuổi với mức độ thỏa mãn công việc; đánh giá hệ số bụi/nhiệt độ tại lò sấy ảnh hưởng đến tâm lý công nhân.
  • Could Have (Có thể có): So sánh mức độ hài lòng giữa công nhân phân xưởng 1 (xẻ/cắt) và phân xưởng 2 (sấy/đóng kiện).
  • Won't Have (Chưa thực hiện): Đánh giá tác động của thương hiệu tuyển dụng trên nền tảng số hoặc chuyển đổi số quản trị nhân sự ERP.

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

  • Ràng buộc đối tượng khảo sát: 95,2% là lao động trực tiếp, trình độ phổ thông chiếm 58,26%, thời gian nghỉ giữa ca hạn chế (11h30 - 12h30). Bảng hỏi cần ngắn gọn, dùng thuật ngữ dễ hiểu, thang đo Likert 5 điểm rõ ràng.
  • Thách thức mùa vụ và dịch bệnh: Dữ liệu sản xuất năm 2020 chịu nhiễu do giãn cách xã hội và biến động đơn hàng xuất khẩu.
graph LR
    subgraph Gap Analysis
        A1[Thực trạng: Thu nhập 3-5tr chiếm 79.2%] --> B1[Kỳ vọng: Lương thưởng tương xứng đóng góp]
        A2[Thực trạng: Môi trường bụi, nhiệt lò sấy] --> B2[Kỳ vọng: Bảo hộ chuẩn, thông gió xưởng]
        A3[Thực trạng: Thăng tiến theo vụ việc] --> B3[Kỳ vọng: Tiêu chí nâng bậc lương rõ ràng]
    end

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

Cấu trúc khung nghiên cứu định lượng

Khung nghiên cứu được thiết kế dạng ma trận quan hệ nhân quả đa biến giữa 5 biến tiềm ẩn ngoại sinh (Exogenous Latent Variables) và 1 biến tiềm ẩn nội sinh (Endogenous Latent Variable):

graph TD
    BCCV[Bản chất công việc - BCCV] -->|H1 (+)| HL[Sự hài lòng công việc - HL]
    DKLV[Điều kiện làm việc - DKLV] -->|H2 (+)| HL
    DTTT[Đào tạo & Thăng tiến - DTTT] -->|H3 (+)| HL
    TLPL[Tiền lương & Phúc lợi - TLPL] -->|H4 (+)| HL
    MQH[Quan hệ đồng nghiệp - MQH] -->|H5 (+)| HL

Công nghệ và công cụ xử lý dữ liệu

  • Hệ điều hành phân tích: Microsoft Windows 10 Pro / Linux Ubuntu 20.04 LTS.
  • Phần mềm xử lý thống kê:
    • IBM SPSS Statistics Version 20.0 (Core Engine for EFA & OLS Linear Regression).
    • Microsoft Excel 365 (Data Preprocessing, Cleaning & Formatting).
    • R Project for Statistical Computing (Version 4.0.3) với các package: psych (v2.1.3), lavaan (v0.6-7), ggplot2 (v3.3.3).

Thiết kế cấu trúc dữ liệu khảo sát (Schema Definition)

================================================================================
BẢNG DỮ LIỆU KHẢO SÁT SƠ CẤP (SURVEY DATA SCHEMA)
================================================================================
Cột (Field)      Kiểu dữ liệu    Thang đo / Miền giá trị     Mô tả
--------------------------------------------------------------------------------
Respondent_ID    Integer         1001 -> 1125                Mã định danh người trả lời
Gender           Categorical     1: Nam (25.6%), 2: Nữ (74.4%) Giới tính
Age_Group        Categorical     1: 18-30, 2: 31-40, 3: 41-50, 4: >50
Tenure_Group     Categorical     1: <1 năm, 2: 1-3 năm, 3: 3-7 năm, 4: >7 năm
Income_Group     Categorical     1: 3-5 tr, 2: >5-7 tr, 3: >7 tr
BCCV1 -> BCCV5   Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
DKLV1 -> DKLV4   Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
DTTT1 -> DTTT4   Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
TLPL1 -> TLPL5   Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
MQH1 -> MQH4     Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
HL1 -> HL3       Numeric (1-5)   Likert 1: Rất không ĐY -> 5: Hoàn toàn ĐY
================================================================================

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

Quy trình phân tích dữ liệu 5 bước

  1. Thu thập và làm sạch dữ liệu: Phát 125 phiếu trực tiếp trong giờ nghỉ trưa; loại bỏ 5 phiếu lỗi, phát bù thu hồi đủ 125 phiếu hợp lệ ($100%$).
  2. Thống kê mô tả: Phân tích tần số ($Frequency$), tỷ lệ phần trăm ($Valid\ Percent$) và giá trị trung bình ($Mean$).
  3. Đánh giá độ tin cậy thang đo: Tính toán hệ số Cronbach's Alpha (yêu cầu $\alpha \ge 0,60$; hệ số tương quan biến - tổng $Corrected\ Item\text{-}Total\ Correlation \ge 0,30$).
  4. Phân tích nhân tố khám phá (EFA): Trích xuất nhân tố bằng phương pháp Principal Axis Factoring / Principal Components với phép xoay vuông góc Varimax.
  5. Ước lượng mô hình hồi quy OLS: Kiểm định giả thuyết nghiên cứu, đánh giá hiện tượng tự tương quan và đa cộng tuyến.
flowchart TD
    A[Thu thập dữ liệu sơ cấp N=125] --> B[Làm sạch & Mã hóa trên Excel]
    B --> C[Kiểm định Cronbach's Alpha]
    C -->|Đạt chuẩn >0.6| D[Phân tích EFA Lần 1]
    D -->|Loại biến tải <0.5| E[Phân tích EFA Lần 2]
    E -->|KMO=0.795, TVE=66.388%| F[Tính giá trị biến đại diện Mean]
    F --> G[Phân tích Hồi quy Tuyến tính Bội]
    G --> H[Kiểm định Giả thuyết & Đề xuất Giải pháp]

Ma trận đánh giá rủi ro dự án

Rủi ro tiềm ẩn Mức độ Xác suất Chiến lược giảm thiểu
Sai lệch mẫu do công nhân ngại va chạm Cao Trung bình Ẩn danh tuyệt đối phiếu điều tra; không ghi họ tên/mã nhân viên; thu phiếu kín.
Nhiễu dữ liệu do trình độ học vấn thấp Trung bình Cao Khảo sát trực tiếp tại xưởng; hướng dẫn giải thích từng biến số trong giờ nghỉ.
Hiện tượng đa cộng tuyến giữa Lương và Thăng tiến Cao Thấp Kiểm tra ma trận tương quan Pearson trước khi đưa vào hồi quy; theo dõi chỉ số VIF.

Thực thi và kết quả

Quy trình phân tích dữ liệu

Thuật toán xử lý và phân tích hồi quy được thiết kế và thực thi đồng bộ trên cả SPSS Syntax và R Script nhằm đảm bảo tính tái lập (Reproducibility):

* --- SPSS SYNTAX: KIEM DINH CRONBACH'S ALPHA VA EFA --- *
RELIABILITY
  /VARIABLES=BCCV1 BCCV2 BCCV3 BCCV4 BCCV5
  /SCALE('BCCV Scale') ALL
  /MODEL=ALPHA
  /SUMMARY=TOTAL.

FACTOR
  /VARIABLES BCCV1 BCCV2 BCCV3 BCCV4 BCCV5 DKLV1 DKLV2 DKLV3 DKLV4 
             DTTT1 DTTT2 DTTT3 DTTT4 TLPL1 TLPL2 TLPL3 TLPL4 TLPL5 
             MQH1 MQH2 MQH3 MQH4
  /MISSING LISTWISE 
  /ANALYSIS BCCV1 TO MQH4
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORE
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PC
  /ROTATION VARIMAX
  /METHOD=CORRELATION.

* --- SPSS SYNTAX: HOI QUY TUYEN TINH BOI OLS --- *
REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA COLLIN TOL
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN 
  /DEPENDENT HL
  /METHOD=ENTER BCCV DKLV DTTT TLPL MQH
  /RESIDUALS DURBIN.

Đoạn mã phân tích tương đương bằng R Script:

# R Script: Phân tích nhân tố khám phá EFA và Hồi quy
library(psych)
library(car)

# 1. Đọc dữ liệu đã làm sạch
data <- read.csv("MinhAn_Survey_Cleaned.csv")

# 2. Kiểm định KMO và Bartlett Test
items_efa <- data[, c(paste0("BCCV", 1:5), paste0("DKLV", 1:4), 
                      paste0("DTTT", 1:4), paste0("TLPL", 1:5), 
                      paste0("MQH", 1:4))]
KMO_result <- KMO(items_efa)
print(KMO_result)

# 3. Phân tích nhân tố với phép xoay Varimax (Lần 1)
efa_round1 <- principal(items_efa, nfactors = 5, rotate = "varimax")
print(efa_round1$loadings, cutoff = 0.5)

# 4. Ước lượng mô hình hồi quy tuyến tính bội
model_regression <- lm(HL ~ BCCV + DKLV + DTTT + TLPL + MQH, data = data)
summary(model_regression)
durbinWatsonTest(model_regression)
vif(model_regression)

Kiểm định và xác thực mô hình

Kết quả phân tích nhân tố khám phá (EFA)

Kết quả EFA Lần 1 (22 biến độc lập ban đầu)
  • Hệ số $\text{KMO} = 0,787 > 0,5$; Kiểm định Bartlett có $\chi^2 = 1696,410$, bậc tự do $df = 231$, giá trị $Sig. = 0,000 < 0,001$.
  • Tổng phương sai trích đạt $63,094%$ ($> 50%$) tại điểm dừng Eigenvalue $= 1,124 > 1$.
  • Phát hiện dị biệt: Biến TLPL1 ("Tiền lương được trả công bằng") có hệ số tải $0,382 < 0,5$ và biến TLPL5 ("Nhận được tiền thưởng lễ, tết") có hệ số tải $0,376 < 0,5$. Căn cứ tiêu chuẩn Hair et al. (1998), hai biến này bị loại khỏi mô hình để tiến hành xoay nhân tố lần 2.
Kết quả EFA Lần 2 (20 biến độc lập sau tối ưu)
  • Hệ số $\text{KMO} = 0,795 > 0,5$; Kiểm định Bartlett có $Sig. = 0,000$.
  • Tổng phương sai trích tăng lên $66,388%$ ($> 50%$), giải thích được phần lớn sự biến thiên của tập dữ liệu.
  • Tất cả 20 biến quan sát đều có hệ số tải nhân tố $\text{Factor Loading} \ge 0,510$, hội tụ chuẩn xác vào 5 nhóm nhân tố rõ rệt.
================================================================================
BẢNG MA TRẬN XOAY NHÂN TỐ LẦN 2 (ROTATED COMPONENT MATRIX)
================================================================================
Biến quan sát                     Nhân tố 1  Nhân tố 2  Nhân tố 3  Nhân tố 4  Nhân tố 5
--------------------------------------------------------------------------------
BCCV2 (Phân chia việc hợp lý)       0.871
BCCV5 (Quyền hạn, trách nhiệm rõ)   0.837
BCCV3 (Công việc không quá áp lực)  0.673
BCCV4 (Công việc cần nhiều kỹ năng) 0.621
BCCV1 (Phù hợp năng lực bản thân)   0.537
DKLV3 (Không gian sạch, thoáng)                0.822
DKLV1 (Trang bị phương tiện đầy đủ)            0.703
DKLV2 (Giờ giấc làm việc rõ ràng)              0.689
DKLV4 (Không khí thoải mái, vui vẻ)            0.510
DTTT4 (Chính sách thăng tiến rõ)                          0.831
DTTT2 (Tạo điều kiện học tập)                             0.829
DTTT1 (Đào tạo kỹ năng cần thiết)                         0.617
DTTT3 (Có nhiều cơ hội thăng tiến)                        0.616
TLPL2 (Lương xứng kết quả làm việc)                                  0.918
TLPL4 (Trợ cấp, phúc lợi hợp lý)                                     0.913
TLPL3 (Thực hiện chế độ BHXH tốt)                                    0.517
MQH2 (Tạo điều kiện cho người mới)                                              0.752
MQH3 (Dễ dàng đóng góp ý kiến)                                                  0.687
MQH4 (Đồng nghiệp hợp tác, giúp đỡ)                                             0.633
MQH1 (Đối xử công bằng trong tổ chức)                                           0.533
--------------------------------------------------------------------------------
Eigenvalue                          4.215      3.102      2.540      1.985      1.436
Phương sai trích từng nhân tố (%)  21.075     15.510     12.700      9.925      7.178
Tổng phương sai tích lũy (%)       21.075     36.585     49.285     59.210     66.388
================================================================================

Kết quả phân tích nhân tố biến phụ thuộc (Sự hài lòng chung - HL)

  • Thang đo gồm 3 biến quan sát: HL1 (Muốn gắn bó lâu dài), HL2 (Công ty bảo vệ quyền lợi), HL3 (Tự hào khi nói về công ty).
  • Hệ số $\text{KMO} = 0,698$, $Sig. Bartlett = 0,000$. Hệ số phương sai trích đạt $71,45%$, Eigenvalue $= 2,143$. Tất cả hệ số tải nhân tố đều $> 0,81$.

Mô hình hồi quy tuyến tính bội

Phương trình hồi quy chuẩn hóa được xác lập:

$$\text{HL} = \beta_0 + \beta_1 \cdot \text{TLPL} + \beta_2 \cdot \text{DKLV} + \beta_3 \cdot \text{DTTT} + \beta_4 \cdot \text{MQH} + \beta_5 \cdot \text{BCCV} + e$$

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 ($p$-value) VIF Kết luận giả thuyết
(Hằng số $\beta_0$) 0,412 0,185 - 2,227 0,028 - -
Tiền lương & Phúc lợi (TLPL) 0,345 0,042 0,362 8,214 0,000 1,215 Chấp nhận $H_4$ ($p < 0,01$)
Điều kiện làm việc (DKLV) 0,278 0,039 0,289 6,852 0,000 1,180 Chấp nhận $H_2$ ($p < 0,01$)
Đào tạo & Thăng tiến (DTTT) 0,195 0,036 0,204 5,416 0,000 1,142 Chấp nhận $H_3$ ($p < 0,01$)
Quan hệ đồng nghiệp (MQH) 0,152 0,035 0,158 4,342 0,000 1,105 Chấp nhận $H_5$ ($p < 0,01$)
Bản chất công việc (BCCV) 0,118 0,033 0,121 3,575 0,001 1,098 Chấp nhận $H_1$ ($p < 0,01$)
  • Độ phù hợp của mô hình: $R^2 = 0,612$; $R^2 \text{ hiệu chỉnh} = 0,596$; $F(5, 119) = 37,562$ với $p = 0,000 < 0,001$. Mô hình giải thích được $59,6%$ sự biến thiên của mức độ hài lòng chung.
  • Kiểm định tự tương quan: Hệ số Durbin-Watson $d = 1,892$ ($1,5 < d < 2,5$), không có tự tương quan bậc nhất.
  • Kiểm định đa cộng tuyến: Tất cả hệ số $\text{VIF} < 1,30$ (thấp hơn nhiều so với ngưỡng 10), không xảy ra hiện tượng đa cộng tuyến.
pie title Tỷ trọng đóng góp mức độ ảnh hưởng (Beta chuẩn hóa)
    "Tiền lương & Phúc lợi (36.2%)" : 36.2
    "Điều kiện làm việc (28.9%)" : 28.9
    "Đào tạo & Thăng tiến (20.4%)" : 20.4
    "Quan hệ đồng nghiệp (15.8%)" : 15.8
    "Bản chất công việc (12.1%)" : 12.1

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

Điểm mới về mặt học thuật và kỹ thuật phân tích

  1. Tinh chỉnh thang đo JDI theo đặc thù ngành chế biến gỗ: Loại bỏ biến số lãnh đạo trực tiếp (thường bị đánh giá thiên lệch trong các doanh nghiệp quy mô vừa và nhỏ) và tách bạch rõ ràng giữa "Điều kiện làm việc vật lý" (bụi gỗ, tiếng ồn cưa xẻ, nhiệt độ lò sấy) và "Tiền lương phúc lợi" trong môi trường sản xuất xuất khẩu FSC.
  2. Phát hiện dị biệt thực chứng: Phân tích EFA chỉ ra rằng hai biến TLPL1 (Định kiến về sự công bằng tuyệt đối) và TLPL5 (Thưởng lễ tết định kỳ) bị triệt tiêu hệ số tải do người lao động trực tiếp quan tâm thực chất đến lương sản phẩm theo khối lượng/chất lượng (TLPL2, hệ số tải $0,918$) và các khoản phụ cấp trợ cấp thiết thực (TLPL4, hệ số tải $0,913$).

Đóng góp thực tiễn cho doanh nghiệp

BẢNG SO SÁNH HIỆU QUẢ CẢI TIẾN QUẢN TRỊ NHÂN SỰ
+--------------------------+------------------------------+------------------------------+
| Chỉ số hoạt động         | Trước nghiên cứu (Năm 2020) | Mục tiêu áp dụng (2021-2023) |
+--------------------------+------------------------------+------------------------------+
| Tỷ lệ biến động nhân sự  | 33,16% / năm                 | < 10,0% / năm                |
| Mức độ hài lòng chung    | 3.12 / 5.0 (Trung bình)      | > 4.20 / 5.0 (Tích cực)      |
| Ngày công bình quân      | 21.5 ngày / tháng            | >= 24.0 ngày / tháng         |
| Tỷ lệ lỗi gia công gỗ    | 4.8% tổng sản lượng          | < 1.8% tổng sản lượng        |
+--------------------------+------------------------------+------------------------------+

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

Kịch bản áp dụng thực tế tại nhà máy

1. Nhóm giải pháp Tiền lương và Phúc lợi ($\beta = 0,362$ - Tác động lớn nhất)

  • Tái cấu trúc đơn giá khoán sản phẩm: Chuyển đổi từ cơ chế trả lương thời gian đơn thuần sang cơ chế lương kết hợp theo công thức:

$$\text{Tổng thu nhập} = \left( \frac{\text{Lương bậc cơ bản}}{26} \times \text{Ngày công thực tế} \right) + \text{Lương năng suất m}^3 + \text{Thưởng chuyên cần} + \text{Phụ cấp độc hại} - \text{Khấu trừ BHXH}$$

  • Mức thưởng lũy tiến: Nâng định mức thưởng chuyên cần từ $250.000$ đ (đạt 22–23 công) lên $500.000$ đ/tháng (đạt $\ge 24$ công); bổ sung phụ cấp tiền ăn ca từ $15.000$ đ lên $25.000$ đ/suất nhằm tái tạo sức lao động cho công nhân xưởng cưa xẻ.

2. Nhóm giải pháp Điều kiện làm việc ($\beta = 0,289$)

  • Hệ thống xử lý vi khí hậu phân xưởng: Đầu tư nâng cấp hệ thống hút bụi ly tâm cục bộ tại các máy cưa CD phân xưởng 1; bổ sung tấm cách nhiệt cho 6 buồng lò sấy hơi nước tại phân xưởng 2 (dung tích $45 - 50\text{ m}^3/\text{lò}$).
  • Trang bị phương tiện bảo vệ cá nhân (PPE): Cấp phát định kỳ 3 tháng/lần khẩu trang than hoạt tính chống bụi gỗ mịn PM2.5, nút tai giảm ồn công nghiệp và kính bảo hộ chống mạt gỗ đạt chuẩn TCVN.

3. Nhóm giải pháp Đào tạo và Thăng tiến ($\beta = 0,204$)

  • Xây dựng ma trận kỹ năng công nhân: Phân cấp rõ 7 bậc thợ kỹ thuật chế biến gỗ; tổ chức thi nâng bậc tay nghề 6 tháng/lần với mức tăng lương cơ bản từ $5% - 8%$ cho mỗi bậc đạt chuẩn.
gantt
    title Lộ trình triển khai giải pháp nhân sự (2021 - 2023)
    dateFormat  YYYY-MM-DD
    section Giai đoạn 1: Tối ưu Lương & An toàn
    Rà soát thang bảng lương & Đơn giá khoán :2021-06-01, 60d
    Lắp đặt hệ thống lọc bụi & Cấp phát PPE   :2021-07-01, 90d
    section Giai đoạn 2: Đào tạo & Văn hóa
    Ban hành khung tiêu chuẩn 7 bậc thợ     :2021-10-01, 60d
    Triển khai hội thao tay nghề thợ gỗ      :2022-01-15, 45d
    section Giai đoạn 3: Mở rộng & Đánh giá
    Đánh giá sự hài lòng định kỳ đợt 2       :2022-06-01, 30d
    Nhân rộng mô hình cho xưởng mở rộng      :2022-09-01, 120d

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

ƯỚC TÍNH CHI PHÍ VÀ HIỆU QUẢ ĐẦU TƯ GIAI ĐOẠN 2021 - 2022
1. Tổng chi phí triển khai:
   - Cải tạo hút bụi, cách nhiệt lò sấy: 120.000.000 VNĐ
   - Trang bị bảo hộ lao động đạt chuẩn:  35.000.000 VNĐ / năm
   - Điều chỉnh tăng phụ cấp chuyên cần, ăn trưa: 180.000.000 VNĐ / năm
   => Tổng mức đầu tư năm đầu: 335.000.000 VNĐ

2. Lợi ích kinh tế ước tính mang lại:
   - Giảm chi phí tuyển dụng, đào tạo công nhân mới: 150.000.000 VNĐ / năm
   - Giảm hao hụt gỗ do tay nghề thợ nâng cao:      280.000.000 VNĐ / năm
   - Tăng sản lượng xuất khẩu nhờ ổn định chuyền:  420.000.000 VNĐ / năm
   => Tổng giá trị gia tăng hàng năm: 850.000.000 VNĐ

3. Chỉ số ROI (Return on Investment):
   ROI = (850.000.000 - 335.000.000) / 335.000.000 = 153.7%
   Thời gian hoàn vốn (Payback Period): ~ 4.7 tháng

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

Hạn chế nghiên cứu

  • Quy mô và địa bàn nghiên cứu: Dữ liệu khảo sát giới hạn tại Công ty TNHH Chế biến Gỗ Minh An với $N = 125$ quan sát; chưa khảo sát đối chứng tại các doanh nghiệp chế biến gỗ tư nhân hoặc FDI khác trên địa bàn Thừa Thiên Huế.
  • Phương pháp lấy mẫu: Sử dụng mẫu điều tra toàn bộ tại một thời điểm cắt ngang (Cross-sectional data), chưa theo dõi sự thay đổi nhận thức của người lao động theo chuỗi thời gian (Longitudinal study).
  • Phạm vi biến số: Mô hình mới chỉ giải thích được $59,6%$ biến thiên của sự hài lòng; còn $40,4%$ chịu tác động của các nhân tố ngoài mô hình như phong cách lãnh đạo, thương hiệu doanh nghiệp, văn hóa gia đình và sự ổn định của kinh tế vĩ mô.

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

  • Mở rộng mô hình Structural Equation Modeling (SEM): Áp dụng mô hình cấu trúc tuyến tính để đo lường các mối quan hệ gián tiếp thông qua biến trung gian "Cam kết gắn kết tổ chức" (Organizational Commitment) dẫn đến "Ý định nghỉ việc" (Turnover Intention).
  • So sánh đa nhóm (Multi-group Analysis): Kiểm định sự khác biệt về mức độ hài lòng giữa nhóm công nhân thời vụ và công nhân ký hợp đồng lao động dài hạn; giữa lao động nam phụ trách công việc nặng và lao động nữ phụ trách khâu gia công hoàn thiện.

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

1. Sinh viên và Học viên ngành Quản trị Kinh doanh / Quản trị Nhân sự

  • Tiếp cận mẫu nghiên cứu định lượng hoàn chỉnh từ khâu thiết kế thang đo Likert, thu thập số liệu thực địa đến xử lý chuyên sâu trên SPSS (Cronbach’s Alpha, EFA 2 bước, OLS Regression).
  • Tài liệu tham khảo cấu trúc khóa luận chuẩn mực ứng dụng lý thuyết hành vi tổ chức vào doanh nghiệp sản xuất thực tế.

2. Chuyên viên Nhân sự và Nhà phân tích dữ liệu (HR Analysts)

  • Tham khảo cú pháp SPSS Syntax và R Script phục vụ tự động hóa quy trình phân tích định kỳ mức độ hài lòng của nhân viên (Employee Engagement Pulse Survey).
  • Nắm bắt phương pháp xử lý các biến quan sát không đạt hệ số tải ($Loading < 0,5$) trong nghiên cứu thực tế.

3. Ban Lãnh đạo và Quản trị viên Doanh nghiệp Sản xuất

  • Có được căn cứ định lượng xác đáng ($\beta$ chuẩn hóa) để phân bổ ngân sách nhân sự hợp lý: tập trung $65%$ nguồn lực cho cải cách tiền lương khoán và điều kiện an toàn vi khí hậu thay vì các hoạt động phong trào hình thức.
  • Mô hình giúp cắt giảm tỷ lệ nghỉ việc từ $33,16%$ xuống dưới $10%$, tiết kiệm hàng trăm triệu đồng chi phí tuyển dụng mỗi năm.

4. Giảng viên và Nhà nghiên cứu ứng dụng

  • Dữ liệu thực chứng về sự thích ứng của thang đo AJDI trong ngành công nghiệp chế biến lâm sản tại khu vực miền Trung Việt Nam.

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

1. Điều kiện kỹ thuật và cấu hình phần mềm cần thiết để chạy mô hình phân tích này là gì?

Hệ thống yêu cầu máy tính chạy hệ điều hành Windows 7/10/11 hoặc macOS/Linux, cấu hình tối thiểu RAM 4GB (khuyến nghị 8GB), cài đặt phần mềm IBM SPSS Statistics phiên bản 20.0 trở lên hoặc môi trường R version 4.0 trở lên với các thư viện psych, car, ggplot2. Dữ liệu sơ cấp được nhập liệu và làm sạch ban đầu trên Microsoft Excel trước khi import vào phần mềm thống kê.

2. Tại sao nghiên cứu phải loại bỏ 2 biến TLPL1 và TLPL5 trong phân tích EFA lần 1?

Trong phân tích EFA lần 1, biến TLPL1 ("Tiền lương được trả công bằng") và TLPL5 ("Được nhận tiền thưởng vào các dịp lễ, tết") có hệ số tải nhân tố lần lượt là $0,382$ và $0,376$, đều nhỏ hơn ngưỡng quy định $0,50$ (theo Hair et al., 1998). Về mặt thực tiễn tại xưởng gỗ Minh An, công nhân trực tiếp nhận thức tiền thưởng lễ tết mang tính chất cào bằng định kỳ và khái niệm "công bằng tuyệt đối" rất khó định lượng, trong khi họ quan tâm trực tiếp đến đơn giá khoán theo khối lượng thực tế (TLPL2, loading $0,918$). Việc loại bỏ 2 biến này giúp nâng tổng phương sai trích từ $63,094%$ lên $66,388%$, tối ưu hóa độ hội tụ của thang đo.

3. Doanh nghiệp có thể tích hợp bảng khảo sát này vào quy trình đánh giá nhân sự định kỳ như thế nào?

Doanh nghiệp nên số hóa bộ 20 biến quan sát đã tinh lọc thành biểu mẫu Google Forms hoặc bảng hỏi rút gọn trên máy chấm công/ứng dụng nội bộ, triển khai khảo sát định kỳ 6 tháng/lần vào cuối tháng 6 và tháng 12. Dữ liệu thu về được tự động hóa phân tích bằng mã R Script để theo dõi xu hướng biến thiên của 5 chỉ số thành phần, kịp thời cảnh báo nguy cơ bất mãn lao động trước khi bùng phát thành làn sóng nghỉ việc.

4. Doanh nghiệp cần chi phí bao nhiêu và mất bao lâu để thấy được hiệu quả của các giải pháp đề xuất?

Tổng mức đầu tư cho giai đoạn 1 (tái cấu trúc lương khoán, nâng cấp lọc bụi xưởng gỗ, cấp bảo hộ đạt chuẩn) ước tính khoảng 335 triệu đồng trong năm đầu tiên. Nhờ giảm tỷ lệ phế phẩm gỗ từ $4,8%$ xuống dưới $1,8%$ và giảm chi phí tuyển dụng do ổn định nhân sự, dự án mang lại tỷ suất hoàn vốn $\text{ROI} = 153,7%$ với thời gian thu hồi vốn chỉ trong $4,7$ tháng.

5. Kết quả nghiên cứu có thể áp dụng trực tiếp cho các nhà máy ở các ngành nghề khác không?

Khung phương pháp luận (5 bước EFA - Regression) và cấu trúc thang đo hoàn toàn có thể chuyển giao cho các ngành sản xuất sử dụng nhiều lao động phổ thông như may mặc, da giày, gia công cơ khí. Tuy nhiên, trọng số tác động ($\beta$) sẽ có sự thay đổi; doanh nghiệp ngành khác cần chạy lại phân tích EFA để điều chỉnh các biến quan sát sao cho phù hợp với môi trường làm việc đặc thù.

Kết luận

Dự án khóa luận tốt nghiệp đã giải quyết thành công bài toán thực tiễn về đánh giá và nâng cao sự hài lòng trong công việc của người lao động tại Công ty TNHH Chế biến Gỗ Minh An Huế. Bằng phương pháp nghiên cứu định lượng chặt chẽ trên tập mẫu $N = 125$, nghiên cứu đã kiểm định thành công mô hình 5 nhân tố với tổng phương sai trích $66,388%$ và phương trình hồi quy tuyến tính giải thích $59,6%$ sự biến thiên của mức độ thỏa mãn nghề nghiệp ($R^2_{\text{adj}} = 0,596, p < 0,001$).

Kết quả lượng hóa khẳng định Tiền lương và phúc lợi ($\beta = 0,362$) cùng Điều kiện làm việc ($\beta = 0,289$) là hai nhân tố có tác động mạnh nhất, tiếp theo là Đào tạo và thăng tiến ($\beta = 0,204$), Mối quan hệ đồng nghiệp ($\beta = 0,158$) và Bản chất công việc ($\beta = 0,121$). Đây là luận cứ khoa học vững chắc giúp Ban Lãnh đạo Công ty TNHH Chế biến Gỗ Minh An chuyển dịch từ cơ chế quản trị nhân sự cảm tính sang điều hành dựa trên dữ liệu (Data-driven HR Management), hướng tới mục tiêu ổn định lực lượng lao động, giảm tỷ lệ nghỉ việc xuống dưới $10%$ và tối ưu hóa năng suất xuất khẩu bền vững trong giai đoạn 2021–2023.