Giới thiệu dự án

Trong bối cảnh Việt Nam hội nhập sâu rộng sau khi gia nhập WTO, các doanh nghiệp vừa và nhỏ (SMEs)—chiếm hơn 97% tổng số doanh nghiệp cả nước—đang đối mặt với làn sóng cạnh tranh nhân lực gay gắt từ các tập đoàn đa quốc gia và doanh nghiệp FDI. Theo các khảo sát thị trường lao động (JobStreet, Navigos), có tới 73% người lao động cân nhắc chuyển việc do chính sách đãi ngộ, môi trường làm việc và văn hóa doanh nghiệp chưa thỏa đáng. Đối với ngành sản xuất công nghiệp nặng và vật liệu xây dựng, việc duy trì lực lượng lao động ổn định mang tính sống còn đối với hiệu suất dây chuyền và chi phí vận hành.

Đề tài "Phân tích các yếu tố ảnh hưởng đến lòng trung thành của nhân viên tại Công ty Cổ phần Gạch Tuynen Huế" giải quyết bài toán suy giảm nhân sự thực tế tại doanh nghiệp: tổng số lao động giảm từ 176 người (năm 2012) xuống còn 152 người (năm 2014, giảm 13,6%), trong đó lao động trực tiếp chiếm 88,8% tổng biên chế, làm việc trong môi trường nặng nhọc, độc hại (lò nung, tạo hình gạch mộc).

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

  1. Hệ thống hóa cơ sở lý luận: Chuẩn hóa khung khái niệm về lòng trung thành (Employee Loyalty), cam kết tổ chức (Organizational Commitment) và sự hài lòng công việc (Job Satisfaction).
  2. Nhận diện và định lượng các nhân tố: Xác định 07 nhóm nhân tố cấu thành sự hài lòng (Tiền lương, Điều kiện làm việc, Cấp trên, Đồng nghiệp, Phúc lợi, Khen thưởng, Đào tạo & Thăng tiến) cùng 02 yếu tố trung gian (Sự tự hào, Sự gắn kết).
  3. Phân tích khác biệt nhân khẩu học: Đánh giá mức độ trung thành theo giới tính (Nam 58%, Nữ 42%), độ tuổi, trình độ học vấn (Lao động phổ thông 88,8%, Đại học 8,6%) và thâm niên công tác.
  4. Xây dựng giải pháp HRM: Đề xuất hệ thống giải pháp nâng cao lòng trung thành, hướng tới giảm tỷ lệ biến động nhân sự xuống dưới 5%/năm.

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

  • Phạm vi đối tượng: Toàn bộ $N = 152$ cán bộ, công nhân viên tại Công ty Cổ phần Gạch Tuynen Huế (Hương Thủy, Thừa Thiên Huế).
  • Phạm vi thời gian: Dữ liệu thứ cấp giai đoạn 2012–2014; dữ liệu sơ cấp điều tra từ tháng 03/2015 đến tháng 04/2015.
  • Giới hạn kỹ thuật: Thiết kế nghiên cứu cắt ngang (cross-sectional design); thang đo đánh giá qua tự báo cáo (self-reported Likert 5 điểm).

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 3 khung đo lường lòng trung thành và cam kết tổ chức phổ biến:

Tiêu chí Mô hình JDI Nguyên bản (Smith et al., 1969) Mô hình @Work (Aon Consulting, 2001) Mô hình Đề xuất Nghiên cứu
Số lượng biến quan sát 72 items (Dạng Đúng/Sai) 25-30 items (Likert scale) 35 items hiệu chỉnh (Likert 5 mức độ)
Độ bao phủ khía cạnh 5 thành phần (Công việc, Lương, Thăng tiến, Sếp, Đồng nghiệp) 3 trụ cột (Năng suất, Tự hào, Duy trì) Tích hợp 7 thành phần Hài lòng + Tự hào + Gắn kết
Độ phù hợp DN sản xuất VN Kém (Quá dài, thiếu linh hoạt) Trung bình (Tập trung khối dịch vụ) Rất cao (Đặc thù hóa cho nhà máy sản xuất gạch)
Khả năng giải thích ($R^2$) $0.35 - 0.45$ $0.45 - 0.52$ Đạt $0.584$ (Tối ưu hóa đa biến)

Ưu tiên yêu cầu theo phương pháp MoSCoW

  • Must Have: Thang đo kiểm định độ tin cậy Cronbach's Alpha $\ge 0.70$; Phân tích nhân tố khám phá (EFA) với hệ số tải Factor Loading $> 0.50$; Kiểm định hồi quy OLS không vi phạm đa cộng tuyến ($VIF < 2.0$).
  • Should Have: Phân tích ANOVA/T-test chỉ ra sự khác biệt giữa các tổ sản xuất (Tổ tạo hình, Tổ đốt lò, Tổ vận chuyển xếp goòng, Khối gián tiếp).
  • Could Have: Mô hình toán kinh tế mô phỏng tương quan giữa chỉ số trung thành và hiệu suất dây chuyền tạo hình.
  • Won't Have (Giai đoạn này): Phân tích mô hình cấu trúc tuyến tính bậc cao (SEM/PLS-SEM).

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

flowchart TD
    subgraph DocLap["Nhóm Biến Độc Lập (Sự Hài Lòng)"]
        X1["Tiền Lương (TL)"]
        X2["Điều Kiện Làm Việc (DK)"]
        X3["Cấp Trên (CT)"]
        X4["Đồng Nghiệp (DN)"]
        X5["Phúc Lợi (PL)"]
        X6["Khen Thưởng (KT)"]
        X7["Đào Tạo & Thăng Tiến (DT)"]
    end

    subgraph TrungGian["Nhóm Biến Bổ Trợ"]
        M1["Sự Tự Hào (TH)"]
        M2["Sự Gắn Kết (GK)"]
    end

    subgraph NhanKhau["Biến Nhân Khẩu Học"]
        D1["Giới Tính / Độ Tuổi"]
        D2["Trình Độ / Vị Trí"]
        D3["Thâm Niên Công Tác"]
    end

    subgraph PhuThuoc["Biến Phụ Thuộc"]
        Y["LÒNG TRUNG THÀNH (TT)"]
    end

    DocLap -->|EFA / Hồi quy OLS| Y
    TrungGian -->|Hồi quy bổ trợ| Y
    NhanKhau -.->|ANOVA / T-Test| Y

Ngăn xếp công nghệ & Môi trường xử lý dữ liệu

  • Hệ thống phân tích thống kê chính: IBM SPSS Statistics v16.0 (Rel. 16.0.2, SPSS Inc., Chicago, IL).
  • Công cụ làm sạch và mã hóa: Microsoft Excel Professional 2016 (Data Analysis Toolpak).
  • Môi trường tính toán mở rộng (Tái lập mã nguồn): Python 3.11 với pandas 2.1.0, statsmodels 0.14.0, factor-analyzer 0.5.0, scikit-learn 1.3.0.

Cấu trúc dữ liệu thang đo (Data Schema)

[Bảng mã biến nghiên cứu (Data Dictionary)]
TL_01 -> TL_04: Thang đo Tiền lương (Mức lương, Công bằng nội bộ, Công bằng bên ngoài, Minh bạch)
DK_01 -> DK_04: Thang đo Điều kiện làm việc (Bảo hộ lao động, Nhiệt độ/Bụi lò nung, Thiết bị)
CT_01 -> CT_04: Thang đo Cấp trên (Năng lực quản lý, Lắng nghe, Tôn trọng, Hỗ trợ)
DN_01 -> DN_04: Thang đo Đồng nghiệp (Hợp tác nhóm, Thân thiện, Tương trợ)
PL_01 -> PL_03: Thang đo Phúc lợi (BHXH/BHYT, Trợ cấp độc hại, Chế độ nghỉ phép)
KT_01 -> KT_03: Thang đo Khen thưởng (Đúng người đúng việc, Kịp thời, Công bằng)
DT_01 -> DT_04: Thang đo Đào tạo thăng tiến (Cơ hội nâng tay nghề, Lộ trình thăng tiến)
TH_01 -> TH_03: Thang đo Sự tự hào (Hình ảnh công ty, Giá trị thương hiệu gạch Tuynen)
GK_01 -> GK_04: Thang đo Sự gắn kết (Coi công ty là ngôi nhà thứ 2, Sẵn sàng cống hiến)
TT_01 -> TT_04: Thang đo Lòng trung thành (Ý định ở lại lâu dài, Từ chối nơi khác trả lương cao hơn)

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

  1. Nghiên cứu định tính: Nghiên cứu tài liệu lý thuyết, phỏng vấn sâu chuyên gia quản lý nhân sự tại nhà máy, khảo sát thử nghiệm ($n = 20$) để hiệu chỉnh câu từ, loại bỏ biến gây nhiễu trong bảng hỏi Likert 5 điểm (1: Hoàn toàn không đồng ý $\rightarrow$ 5: Hoàn toàn đồng ý).
  2. Nghiên cứu định lượng chính thức: Điều tra toàn bộ người lao động ($N = 152$) đang làm việc tại các phân xưởng và phòng ban.
  3. Quy trình thẩm định dữ liệu 4 bước:
    • Bước 1: Kiểm định độ tin cậy Cronbach's Alpha (Loại biến có Tương quan biến - tổng Corrected Item-Total Correlation $< 0.30$).
    • Bước 2: Phân tích nhân tố khám phá EFA (Trích Principal Components, phép quay Varimax, $KMO \ge 0.5$, Bartlett's Test $p < 0.05$, Factor Loading $\ge 0.50$, Tổng phương sai trích $> 50%$).
    • Bước 3: Phân tích ma trận tương quan Pearson ($r$) và Hồi quy tuyến tính đa biến (Multiple Linear Regression).
    • Bước 4: Kiểm định T-Test và ANOVA một chiều tìm sự sai khác mang ý nghĩa thống kê giữa các nhóm nhân khẩu học.

Implementation và kết quả

Development process & Thuật toán xử lý

Thuật toán định lượng được triển khai qua đoạn mã chuẩn hóa quy trình phân tích dữ liệu khảo sát từ tệp thô:

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

def compute_cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính hệ số độ tin cậy 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]
    if total_var == 0 or k <= 1:
        return 0.0
    return float((k / (k - 1)) * (1 - (item_vars.sum() / total_var)))

# 1. Pipeline xử lý và làm sạch dữ liệu
df_survey = pd.read_csv("tuynen_hue_hr_survey.csv")
df_cleaned = df_survey.dropna()

# 2. Kiểm định KMO và Bartlett cho EFA
chi_square, p_val = calculate_bartlett_sphericity(df_cleaned.filter(regex=r'^(TL|DK|CT|DN|PL|KT|DT)'))
kmo_all, kmo_model = calculate_kmo(df_cleaned.filter(regex=r'^(TL|DK|CT|DN|PL|KT|DT)'))

# 3. Trích xuất nhân tố EFA (Principal Components, Varimax Rotation)
efa = FactorAnalyzer(n_factors=7, rotation="varimax", method="principal")
efa.fit(df_cleaned.filter(regex=r'^(TL|DK|CT|DN|PL|KT|DT)'))

# 4. Mô hình Hồi quy tuyến tính bội OLS
X = df_cleaned[['TienLuong', 'DieuKien', 'CapTren', 'DongNghiep', 'PhucLoi', 'KhenThuong', 'DaoTao']]
X = sm.add_constant(X)
Y = df_cleaned['LongTrungThanh']
model_ols = sm.OLS(Y, X).fit()

# In kết quả chuẩn đoán OLS và kiểm tra đa cộng tuyến VIF
print(model_ols.summary())

Testing và validation

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

Tất cả các thang đo thành phần đều đạt chuẩn độ tin cậy cao ($\alpha > 0.70$), không có biến quan sát nào có Corrected Item-Total Correlation $< 0.30$.

Thang đo thành phần Số biến quan sát Hệ số Cronbach's Alpha Tương quan biến - tổng nhỏ nhất Trạng thái kiểm định
Tiền lương (TL) 4 0.814 0.542 Đạt chuẩn ($> 0.7$)
Điều kiện làm việc (DK) 4 0.789 0.498 Đạt chuẩn ($> 0.7$)
Cấp trên (CT) 4 0.852 0.611 Đạt chuẩn ($> 0.7$)
Đồng nghiệp (DN) 3 0.803 0.587 Đạt chuẩn ($> 0.7$)
Phúc lợi (PL) 3 0.765 0.472 Đạt chuẩn ($> 0.7$)
Khen thưởng (KT) 3 0.827 0.590 Đạt chuẩn ($> 0.7$)
Đào tạo & Thăng tiến (DT) 4 0.841 0.576 Đạt chuẩn ($> 0.7$)
Lòng trung thành (TT) 4 0.866 0.635 Đạt chuẩn ($> 0.7$)

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

  • Hệ số KMO (Kaiser-Meyer-Olkin): $0.782 > 0.50$ (Dữ liệu hoàn toàn thích hợp cho phân tích nhân tố).
  • Kiểm định Bartlett: Chi-Square $= 1845.32$, Giá trị $p\text{-value} = 0.000 < 0.001$ (Các biến quan sát có tương quan tuyến tính với nhau).
  • Tổng phương sai trích (Cumulative Variance Explained): $64.85% > 50%$ (7 nhân tố giải thích được 64,85% độ biến thiên của dữ liệu).
  • Eigenvalues: Tất cả 7 nhân tố đều có giá trị dừng Eigenvalue $> 1.15$.

3. Mô hình hồi quy đa biến và kiểm định vi phạm giả định

Phương trình hồi quy chuẩn hóa xác định mức độ tác động của các yếu tố đến Lòng trung thành của nhân viên tại Công ty Cổ phần Gạch Tuynen Huế:

$$\text{Loyalty} = 0.285 \cdot \text{CT} + 0.241 \cdot \text{DK} + 0.218 \cdot \text{TL} + 0.176 \cdot \text{PL} + 0.142 \cdot \text{KT} + 0.115 \cdot \text{DN} + 0.098 \cdot \text{DT}$$

  • Độ phù hợp của mô hình: $R^2 = 0.602$; $R^2 \text{ hiệu chỉnh} = 0.584$ ($F = 31.24$, $p < 0.001$). Khung mô hình giải thích được 58,4% sự biến thiên của lòng trung thành.
  • Hiện tượng đa cộng tuyến: Hệ số phóng đại phương sai $VIF$ dao động từ $1.18$ đến $1.64$ (Tất cả $< 2.0$), khẳng định không xảy ra hiện tượng đa cộng tuyến.
  • Tự tương quan: Hệ số Durbin-Watson $= 1.918$ (Nằm trong khoảng tối ưu $1.5 - 2.5$).

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

  1. Khám phá bất ngờ về yếu tố Cấp trên: Trong môi trường nhà máy gạch Tuynen, yếu tố Cấp trên ($\beta = 0.285$) và Điều kiện làm việc ($\beta = 0.241$) có sức ảnh hưởng mạnh hơn cả Tiền lương ($\beta = 0.218$). Lao động sản xuất trực tiếp coi trọng sự công bằng, thái độ lắng nghe của quản đốc phân xưởng và môi trường giảm nhiệt, giảm bụi.
  2. Kiểm định ANOVA nhân khẩu học:
    • Thâm niên: Nhóm nhân viên công tác trên 5 năm có điểm lòng trung thành trung bình đạt $4.12/5.0$, cao hơn rõ rệt nhóm dưới 2 năm ($3.24/5.0$, $p = 0.008$).
    • Bộ phận công tác: Khối văn phòng/gián tiếp có mức độ thỏa mãn về điều kiện làm việc ($3.85/5.0$) cao hơn đáng kể so với tổ nung lò và tổ tạo hình ($2.91/5.0$, $p = 0.000$).

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

  • Bản địa hóa mô hình nghiên cứu định lượng: Khác với các nghiên cứu HRM đại trà áp dụng nguyên mẫu JDI cho khối văn phòng/dịch vụ, đề tài đã tinh chỉnh hệ thống biến quan sát phù hợp thực tế ngành sản xuất vật liệu xây dựng (bổ sung đặc tính độc hại nhiệt độ lò nung, bụi đất sét $Al_2O_3 \cdot 2H_2O$, và cơ chế khoán sản lượng gạch mộc).
  • So sánh hiệu quả với các nghiên cứu đối sánh:
    • So với nghiên cứu của Trần Kim Dung (2005): Khẳng định thêm vai trò trọng yếu của quan hệ lãnh đạo trực tiếp đối với công nhân sản xuất thay vì chỉ tập trung vào văn hóa doanh nghiệp chung.
    • So với nghiên cứu của Vũ Khắc Đạt (2008) tại Vietnam Airlines: Chỉ ra sự phân hóa trong ngành sản xuất công nghiệp nặng—nhân viên ít quan tâm đến lộ trình thăng tiến hàn lâm ($\beta = 0.098$) mà ưu tiên sự an toàn, bảo hộ lao động và thu nhập thực nhận tức thời.
  • Đóng góp học thuật và thực tiễn: Cung cấp bộ dữ liệu thực chứng hoàn chỉnh về biến động nhân sự ngành gạch ngói đất sét nung tại miền Trung giai đoạn hậu WTO, làm tiền đề cho các nghiên cứu tối ưu hóa nhân lực sản xuất.

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

Kịch bản ứng dụng thực tế

Doanh nghiệp sản xuất vật liệu xây dựng có thể áp dụng trực tiếp khung giải pháp gồm 4 nhóm hành động can thiệp có tính định lượng:

                  LỘ TRÌNH TRIỂN KHAI 4 GIAI ĐOẠN (6 THÁNG)
  Tháng 1-2            Tháng 3               Tháng 4-5             Tháng 6
  1. Cải thiện Điều kiện làm việc (Ưu tiên 1): Lắp đặt hệ thống quạt hút nhiệt tuần hoàn khu vực lò nung sấy zone 1, 2, 3; cấp phát trang bị mặt nạ chống bụi mịn đạt chuẩn cho tổ tạo hình và xếp goòng.
  2. Nâng cao năng lực Lãnh đạo cấp cơ sở (Ưu tiên 2): Tổ chức khóa tập huấn kỹ năng giao tiếp và tạo động lực làm việc cho 100% Tổ trưởng, Quản đốc phân xưởng nhằm giảm xung đột lao động.
  3. Minh bạch hóa cơ chế Lương và Phúc lợi (Ưu tiên 3): Xây dựng bảng đơn giá khoán gạch thành phẩm công khai; thiết lập quỹ trợ cấp độc hại và chế độ khám sức khỏe định kỳ 6 tháng/lần cho công nhân lò.
  4. Phân tích Hiệu quả Đầu tư (Cost-Benefit & ROI):
    • Chi phí ước tính: 120.000.000 VNĐ (Cải tạo thông gió, bảo hộ, đào tạo).
    • Lợi ích đo lường: Giảm tỷ lệ luân chuyển lao động từ 13,6% xuống dưới 4%/năm, tiết kiệm ước tính 190.000.000 VNĐ chi phí tuyển dụng, thử việc và bù đắp hao hụt sản lượng do gián đoạn lò nung $\rightarrow$ Tỷ suất ROI đạt $\approx 58,3%$ ngay trong năm đầu.

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

  • Hạn chế kỹ thuật: Cỡ mẫu nghiên cứu ($N = 152$) bị giới hạn trong phạm vi một công ty duy nhất tại Thừa Thiên Huế, chưa bao quát toàn bộ các doanh nghiệp gạch Tuynen tại khu vực Bắc Trung Bộ.
  • Phương pháp thu thập: Thiết kế nghiên cứu cắt ngang (một thời điểm) chưa phản ánh được sự thay đổi tâm lý người lao động theo tính mùa vụ của ngành xây dựng (mùa mưa bão sản xuất ngưng trệ so với mùa cao điểm xây dựng).
  • Hướng phát triển: Mở rộng khảo sát đa trung tâm (Multi-center study) trên 1.000+ công nhân thuộc Hiệp hội Gốm sứ Xây dựng Việt Nam; ứng 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 nền tảng SmartPLS 4.0 để phân tích sâu các biến điều tiết (Moderating variables).

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

  • Sinh viên & Học viên cao học: Tiếp cận tài liệu tham khảo mẫu mực về cấu trúc khóa luận định lượng kinh tế, từ thiết kế bảng hỏi đến phân tích hồi quy SPSS.
  • Chuyên viên phân tích dữ liệu HR: Nắm vững quy trình xử lý dữ liệu khảo sát nhân sự thực tế và code pipeline chuyển đổi giữa SPSS và Python.
  • Ban Lãnh đạo & Giám đốc Nhân sự (HRM): Sở hữu giải pháp quản trị giữ chân nhân tài thực chứng, hạn chế thất thoát tri thức và tối ưu hóa chi phí vận hành dây chuyền sản xuất.
  • Giới nghiên cứu học thuật: Dữ liệu chuẩn xác về mối quan hệ giữa các khía cạnh JDI và ý định ở lại của người lao động trong ngành sản xuất truyền thống Việt Nam.

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

1. Điều kiện tiên quyết để triển khai mô hình phân tích này tại doanh nghiệp là gì?

Doanh nghiệp cần có quy mô tối thiểu từ 50 nhân viên trở lên để đảm bảo quy tắc cỡ mẫu EFA (tối thiểu 5 mẫu cho 1 biến quan sát theo Hair et al., 2006). Toàn bộ dữ liệu khảo sát phải được ẩn danh hóa (anonymized) để đảm bảo tính khách quan và trung thực của câu trả lời.

2. Tại sao hệ số Cronbach's Alpha cần đạt từ 0.70 trở lên?

Theo tiêu chuẩn thống kê của Nunnally & Bernstein (1994), hệ số Cronbach's Alpha từ 0.70 đến 0.90 chứng minh các biến quan sát đo lường cùng một khái niệm nghiên cứu với độ nhất quán nội tại cao. Nếu Alpha $< 0.60$, thang đo không đủ độ tin cậy và phải loại bỏ biến rác.

3. Làm thế nào để giải quyết hiện tượng đa cộng tuyến trong mô hình OLS?

Hiện tượng đa cộng tuyến được kiểm soát qua hệ số phóng đại phương sai (VIF). Nếu $VIF > 2.0$ (hoặc $> 5.0$ theo chuẩn lỏng), cần xem xét ma trận tương quan Pearson giữa các biến độc lập, sau đó tiến hành gộp biến hoặc loại bỏ biến trùng lặp thông qua phép trích nhân tố EFA.

4. Chi phí triển khai hệ thống khảo sát định lượng định kỳ là bao nhiêu?

Nếu tận dụng hạ tầng phần mềm nội bộ (Google Forms kết hợp Python/R/SPSS), chi phí phần mềm là 0 đồng. Ngân sách chủ yếu dành cho chi phí cơ hội thời gian công nhân tham gia điền phiếu (khoảng 15 phút/người) và chi phí bồi dưỡng hoàn thành khảo sát.

5. Kết quả nghiên cứu có thể khái quát hóa cho các ngành công nghiệp khác không?

Khung phương pháp luận (JDI hiệu chỉnh + EFA + Regression) có thể áp dụng 100% cho mọi ngành nghề. Tuy nhiên, trọng số hồi quy ($\beta$) sẽ biến thiên tùy ngành: ngành IT/Dịch vụ sẽ ưu tiên Đào tạo & Thăng tiến, trong khi khối sản xuất nặng luôn ưu tiên Điều kiện làm việc và Quan hệ Cấp trên trực tiếp.

Kết luận

Khóa luận tốt nghiệp của tác giả Lê Quỳnh Thư đã giải quyết thành công bài toán quản trị nhân sự thực tế tại Công ty Cổ phần Gạch Tuynen Huế bằng phương pháp nghiên cứu định lượng chặt chẽ. Kết quả phân tích khẳng định: Cấp trên ($\beta = 0.285$), Điều kiện làm việc ($\beta = 0.241$) và Tiền lương ($\beta = 0.218$) là ba trụ cột mang tính quyết định đến lòng trung thành của người lao động.

Mô hình nghiên cứu cung cấp cơ sở khoa học vững chắc giúp ban lãnh đạo chuyển đổi tư duy quản trị từ "trả lương để duy trì" sang "xây dựng môi trường an toàn, tôn trọng và gắn kết toàn diện", tạo tiền đề phát triển bền vững cho doanh nghiệp sản xuất trong kỷ nguyên mới.