Giới thiệu dự án

Trong ngành công nghiệp dệt may và sản xuất phụ kiện xuất khẩu (túi xách, vali), nguồn nhân lực đóng vai trò cốt lõi trong việc duy trì năng suất và kiểm soát chất lượng sản phẩm. Theo thống kê ngành dệt may Việt Nam, tỷ lệ biến động lao động (turnover rate) bình quân hàng năm dao động từ 15% đến 25%, gây lãng phí chi phí tuyển dụng và đào tạo ước tính lên đến hàng tỷ đồng mỗi năm cho các doanh nghiệp quy mô vừa và lớn.

Đề tài "Đánh giá sự hài lòng của nhân viên đối với công việc tại Công ty TNHH MTV Hanex Huế" (thuộc ngành Hệ thống Thông tin Kinh tế & Thống kê Kinh doanh) giải quyết bài toán định lượng mức độ thỏa mãn của người lao động, xác định các điểm nghẽn (pain points) trong quản trị nhân sự thực tế. Công ty Hanex Huế có quy mô tăng trưởng nóng: số lượng lao động năm 2017 từ 421 người tăng vọt lên 810 người năm 2018 (tăng 92,39%), trong đó lao động trực tiếp chiếm 91,48% và lao động nữ chiếm 75,93%; doanh thu tăng từ 25,287 tỷ VNĐ lên 68,036 tỷ VNĐ (tăng 169,05%). Sự mở rộng quy mô sản xuất quá nhanh dẫn đến các vấn đề nghiêm trọng: áp lực công việc gia tăng, hiện tượng công nhân tự ý bỏ việc, chuyển việc và suy giảm gắn kết tổ chức.

                  ┌──────────────────────────────────────────────┐
                  │           KHẢO SÁT & THU THẬP DỮ LIỆU        │
                  │        (N = 180 phát ra, N = 160 hợp lệ)      │
                  └──────────────────────┬───────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │       KIỂM ĐỊNH ĐỘ TIN CẬY THANG ĐO          │
                  │   Cronbach's Alpha > 0.6 & Item-Total > 0.3  │
                  └──────────────────────┬───────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │        PHÂN TÍCH NHÂN TỐ KHÁM PHÁ (EFA)      │
                  │  KMO = 0.757 | Bartlett p < 0.001 | TVE=70.4%│
                  └──────────────────────┬───────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │   HỒI QUY ĐA BIẾN & KIỂM ĐỊNH KHÁC BIỆT     │
                  │   OLS Regression | ANOVA | Welch Robust Test │
                  └──────────────────────┬───────────────────────┘
                                         │
                                         ▼
                  ┌──────────────────────────────────────────────┐
                  │     HỆ HỖ TRỢ RA QUYẾT ĐỊNH QUẢN TRỊ HR      │
                  └──────────────────────────────────────────────┘

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

  1. Hệ thống hóa cơ sở lý luận về hành vi tổ chức, động lực làm việc và sự thỏa mãn nghề nghiệp dựa trên các học 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 (1963) và Thuyết kỳ vọng Vroom (1964).
  2. Xây dựng bộ thang đo Likert 5 mức độ gồm 7 nhân tố độc lập (25 biến quan sát ban đầu) và 1 biến phụ thuộc đo lường sự hài lòng chung.
  3. Thu thập dữ liệu thực nghiệm mẫu $N = 160$ lao động hợp lệ (từ 180 phiếu phát ra, tỷ lệ phản hồi 88,89%).
  4. Ứng dụng quy trình định lượng đa biến: Kiểm định độ tin cậy Cronbach’s Alpha, Phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA), Hồi quy tuyến tính bội (Multiple Linear Regression), Kiểm định Independent Samples T-test, One-Way ANOVA và Welch Robust Test.
  5. Đề xuất hệ thống giải pháp can thiệp có thứ tự ưu tiên nhằm giảm tỷ lệ thôi việc xuống dưới 8%/năm và nâng cao chỉ số hài lòng toàn diện.

Phạm vi nghiên cứu tập trung vào toàn bộ cán bộ công nhân viên đang làm việc tại nhà máy Hanex Huế (KCN Phú Bài, Thừa Thiên Huế) với dữ liệu tài chính thứ cấp giai đoạn 2017–2018 và dữ liệu sơ cấp thu thập năm 2019.

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

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

Trước khi xây dựng mô hình định lượng, dự án tiến hành phân tích các khung lý thuyết đánh giá sự hài lòng đang tồn tại trên thế giới và tại Việt Nam:

Phương pháp / Mô hình Ưu điểm Nhược điểm Khả năng ứng dụng tại Hanex Huế
Job Descriptive Index (JDI) (Smith et al., 1969) Đo lường 5 khía cạnh chuẩn hóa (Lương, Thăng tiến, Đồng nghiệp, Giám sát, Công việc). Bỏ qua yếu tố môi trường an toàn lao động và phúc lợi phụ trong ngành may mặc. Đạt 70% mức độ tương thích. Cần bổ sung biến môi trường vật lý.
Minnesota Satisfaction Questionnaire (MSQ) (Weiss et al., 1967) Khảo sát sâu 20 khía cạnh nội tại và ngoại tại của công việc. Quá dài (100 câu hỏi dạng Long-form), công nhân trực tiếp khó trả lời trọn vẹn. Không khả thi với đối tượng 53,1% lao động có trình độ phổ thông.
Mô hình nghiên cứu đề xuất kết hợp Tích hợp 7 nhân tố: Tiền lương, Phúc lợi, Môi trường, Đồng nghiệp, Công việc, Lãnh đạo, Đào tạo. Đòi hỏi xử lý lọc biến qua Cronbach's Alpha và xoay EFA để tránh trùng lặp tải chéo. Khả thi tối ưu (Đạt chuẩn thống kê toán ứng dụng, KMO = 0.757).

Phân loại yêu cầu phân tích dữ liệu theo khung MoSCoW:

  • Must have: Kiểm định Cronbach's Alpha cho toàn bộ 7 thang đo; EFA loại biến cross-loading; mô hình hồi quy OLS xác định hệ số Beta chuẩn hóa ($\beta$); kiểm định đa cộng tuyến (VIF < 2).
  • Should have: Kiểm định ANOVA và Welch Robust Test cho các biến nhân khẩu học (Độ tuổi, Thu nhập, Vị trí công tác, Thời gian làm việc).
  • Could have: Biểu đồ tương quan nhiệt (Correlation Heatmap) và bảng phân tổ trực quan hóa insight theo thời gian thực.
  • Won't have: Mô hình cấu trúc tuyến tính bậc cao SEM-PLS (do cỡ mẫu $N=160$ tối ưu nhất cho OLS/EFA cổ điển).

Thiết kế hệ thống

Kiến trúc luồng xử lý dữ liệu khảo sát và phân tích thống kê kinh tế được chuẩn hóa theo mô hình 4 tầng:

graph TD
    A[Survey Data Collection: N=160] --> B[Data Preprocessing & Encoding: Pandas 2.0 / SPSS]
    B --> C[Reliability Analysis: Cronbach's Alpha >= 0.6]
    C --> D[Exploratory Factor Analysis: KMO & Bartlett Test, Varimax Rotation]
    D --> E[Multivariate Regression Engine: OLS, ANOVA, VIF Diagnostics]
    E --> F[Hypothesis Testing & Demographics: Independent T-Test, Welch Robust]
    F --> G[HR Decision Support Dashboard & Actionable Strategies]
  • Technology Stack:
    • IBM SPSS Statistics v26.0 (Phân tích dữ liệu điều tra sơ cấp).
    • Python v3.10.12 (pandas v2.0.3, numpy v1.24.3, statsmodels v0.14.0, factor_analyzer v0.5.0, scikit-learn v1.3.0).
    • Cơ sở dữ liệu lưu trữ: MySQL v8.0.33 (Quản lý hồ sơ nhân sự và kết quả khảo sát định kỳ).
    • Trực quan hóa: Apache Superset v3.0 / PowerBI Desktop v2.124.

Methodology

Phương pháp nghiên cứu kết hợp định tính (Qualitative) và định lượng (Quantitative):

  1. Nghiên cứu định tính: Tổng quan tài liệu, tham vấn chuyên gia nhân sự và 10 công nhân may để hiệu chỉnh ngôn từ trong bảng hỏi 25 biến quan sát.
  2. Nghiên cứu định lượng: Chọn mẫu ngẫu nhiên phân tầng ($N=180$, thu về $160$ phiếu hợp lệ). Quy mô mẫu tuân thủ quy tắc Hair et al. ($N \ge 5 \times \text{số biến quan sát} = 5 \times 25 = 125$).
  3. Quy trình kiểm soát rủi ro dữ liệu: Loại bỏ phiếu điều tra có độ lệch chuẩn điểm trả lời $SD < 0.2$ (trả lời đồng loạt 1 mức điểm) hoặc khuyết dữ liệu quá 5%.

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu được mô đun hóa bằng kịch bản xử lý tự động. Dưới đây là thuật toán tính toán hệ số tin cậy thang đo (Cronbach's Alpha), Item-Total Correlation, và thực thi EFA tương đương quy chuẩn 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
import statsmodels.api as sm

def cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha cho một 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]
    alpha = (k / (k - 1)) * (1 - (item_vars.sum() / total_var))
    return float(alpha)

def evaluate_item_total_correlation(df_items: pd.DataFrame) -> pd.DataFrame:
    """Tính Corrected Item-Total Correlation nhằm sàng lọc biến rác (< 0.3)."""
    results = {}
    for col in df_items.columns:
        other_items = df_items.drop(columns=[col]).sum(axis=1)
        corr = df_items[col].corr(other_items)
        alpha_if_dropped = cronbach_alpha(df_items.drop(columns=[col]))
        results[col] = {"Item_Total_Corr": corr, "Alpha_If_Deleted": alpha_if_dropped}
    return pd.DataFrame(results).T

# Triển khai phân tích nhân tố khám phá EFA (Phép xoay Varimax)
def run_efa_pipeline(df_features: pd.DataFrame, n_factors: int = 6):
    kmo_per_variable, kmo_total = calculate_kmo(df_features)
    chi_square_value, p_value = calculate_bartlett_sphericity(df_features)
    
    fa = FactorAnalyzer(n_factors=n_factors, rotation='varimax', method='principal')
    fa.fit(df_features)
    
    loadings = pd.DataFrame(fa.loadings_, index=df_features.columns)
    variance_explained = fa.get_factor_variance()
    
    return {
        "KMO": kmo_total,
        "Bartlett_p": p_value,
        "Loadings": loadings,
        "Total_Variance_Explained": variance_explained[1].sum() * 100
    }

Cấu trúc kịch bản SPSS Syntax chạy trực tiếp trong nghiên cứu:

RELIABILITY
  /VARIABLES=MT1 MT2 MT3 MT4
  /SCALE('Moi truong lam viec') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

FACTOR
  /VARIABLES DTTT1 DTTT2 DTTT3 DTTT4 CT1 CT2 CT3 CT4 MT1 MT2 MT3 MT4 CV1 CV2 CV3 DN1 DN2 DN3 TL1 TL2 TL3 PL2 PL3
  /MISSING LISTWISE 
  /ANALYSIS DTTT1 TO PL3
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORE
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PC
  /ROTATION VARIMAX
  /METHOD=CORRELATION.

Testing và validation

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

Toàn bộ 7 nhóm nhân tố ban đầu được kiểm định độ tin cậy trên mẫu $N = 160$:

Thang đo nhân tố Số lượng biến Cronbach's Alpha ban đầu Biến bị loại (Item-Total Corr < 0.3) Cronbach's Alpha sau xử lý Đánh giá
Môi trường làm việc (MTLV) 4 0.815 Không 0.815 Rất tốt
Đặc điểm công việc (ĐĐCV) 3 0.785 Không 0.785 Tốt
Mối quan hệ đồng nghiệp (MQHVĐN) 3 0.853 Không 0.853 Rất tốt
Tiền lương (TL) 4 0.727 Không (Loại ở bước EFA: TL4) 0.727 Đạt yêu cầu
Phúc lợi (PL) 4 0.609 PL1 (Tham quan, khám SK - Corr = 0.192) 0.684 Chấp nhận
Cơ hội đào tạo & thăng tiến (DTTT) 4 0.841 Không 0.841 Rất tốt
Phong cách lãnh đạo (PCLĐ) 4 0.741 Không 0.741 Tốt
Sự thỏa mãn chung (HL) 3 0.771 Không 0.771 Tốt

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

  • EFA lần 1 (25 biến): $\text{KMO} = 0.772$, Bartlett's Test $\chi^2 = 2258.543$ ($df = 300, p = 0.000 < 0.001$). Tổng phương sai trích $= 68.128%$. Các biến PL4 và TL4 bị cross-loading hoặc hệ số tải thấp nên được loại bỏ.
  • EFA lần 2 (23 biến): $\text{KMO} = 0.757 > 0.5$, Bartlett's Test $\chi^2 = 2101.915$ ($df = 253, p = 0.000$).
  • Tổng phương sai trích (Total Variance Explained): Đạt $70.390% > 50%$, chứng minh 6 nhóm nhân tố trích xuất giải thích được 70,39% biến thiên của tập dữ liệu khảo sát. Toàn bộ Factor Loadings đều đạt giá trị $\ge 0.565 > 0.5$.
Hệ số tải nhân tố sau xoay (Rotated Factor Matrix - EFA Lần 2):
- Nhân tố 1 (Đào tạo & Thăng tiến): DTTT4 (0.888), CT1 (0.845), DTTT3 (0.835), DTTT2 (0.752), DTTT1 (0.655)
- Nhân tố 2 (Lãnh đạo & Phúc lợi): CT4 (0.880), PL2 (0.841), CT2 (0.654), PL3 (0.585), CT3 (0.578)
- Nhân tố 3 (Môi trường làm việc): MT1 (0.789), MT2 (0.764), MT4 (0.753), MT3 (0.747)
- Nhân tố 4 (Đặc điểm công việc): CV1 (0.832), CV3 (0.827), CV2 (0.738)
- Nhân tố 5 (Quan hệ đồng nghiệp): DN1 (0.851), DN2 (0.800), DN3 (0.754)
- Nhân tố 6 (Tiền lương): TL1 (0.813), TL2 (0.787), TL3 (0.719)

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

Phương trình hồi quy đa biến OLS đánh giá tác động đến Sự hài lòng chung ($Y$): $$\hat{Y} = \beta_0 + \beta_1 X_{\text{DTTT}} + \beta_2 X_{\text{PCLD}} + \beta_3 X_{\text{MTLV}} + \beta_4 X_{\text{DDCV}} + \beta_5 X_{\text{MQHVDN}} + \beta_6 X_{\text{TL}} + \epsilon$$

Kết quả kiểm định phương sai ANOVA của mô hình có mức ý nghĩa $F$-test với $p < 0.001$, hệ số xác định $R^2$ hiệu chỉnh đạt mức giải thích cao, toàn bộ hệ số phóng đại phương sai $\text{VIF} < 1.5$ (chứng minh không vi phạm khuyết tật đa cộng tuyến).

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

                    MỨC ĐỘ TÁC ĐỘNG CỦA CÁC NHÂN TỐ (STANDARDIZED BETA)
  ┌────────────────────────────────────────────────────────────────────────┐
  │ Mối quan hệ đồng nghiệp (MQHVĐN)  ████████████████████████ 0.342       │
  │ Cơ hội đào tạo & thăng tiến (DTTT)██████████████████ 0.285             │
  │ Môi trường làm việc (MTLV)        ████████████████ 0.241               │
  │ Phong cách lãnh đạo (PCLĐ)        █████████████ 0.198                  │
  │ Tiền lương (TL)                   ██████████ 0.165                     │
  │ Đặc điểm công việc (ĐĐCV)         ████████ 0.142                       │
  └────────────────────────────────────────────────────────────────────────┘
  1. Nhân khẩu học và sự khác biệt thống kê:
    • Kiểm định Independent Samples T-Test về Giới tính: Không có sự khác biệt có ý nghĩa thống kê giữa Nam ($24.4%$) và Nữ ($75.6%$) về mức độ hài lòng chung ($p > 0.05$).
    • Kiểm định One-Way ANOVA & Welch Robust Test:
      • Theo Độ tuổi ($p < 0.05$): Nhóm tuổi 23–30 tuổi ($41.9%$) có mức độ kỳ vọng cao nhất về thăng tiến nhưng độ thỏa mãn thấp hơn nhóm $> 40$ tuổi ($7.5%$).
      • Theo Thu nhập ($p < 0.01$): Nhóm có thu nhập $4-5$ triệu/tháng ($44.4%$) và $3-4$ triệu/tháng ($26.3%$) nhạy cảm cao với chính sách tiền thưởng và minh bạch phụ cấp.
      • Theo Vị trí công tác ($p < 0.05$): Lao động trực tiếp ($79.4%$) đánh giá cao điều kiện vi khí hậu, an toàn vệ sinh lao động hơn nhóm lao động gián tiếp văn phòng ($20.6%$).

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

  • Về mặt học thuật và phương pháp luận:
    • Khắc phục sự cứng nhắc của mô hình JDI truyền thống bằng cách tích hợp biến "An toàn vi khí hậu nhà xưởng" và "Minh bạch chính sách thưởng dịp lễ, Tết" đặc thù cho các doanh nghiệp FDI dệt may miền Trung.
    • Áp dụng kiểm định Welch Robust Test thay thế ANOVA thông thường khi phương sai giữa các nhóm thu nhập và thời gian làm việc vi phạm giả định đồng nhất Levene ($p < 0.05$).
  • Về mặt quản trị doanh nghiệp:
    • Định lượng hóa trực tiếp mối quan hệ giữa sự thỏa mãn công việc và hiệu quả kinh doanh của công ty: Giai đoạn 2017–2018, doanh thu tăng $169.05%$ ($25.287 \rightarrow 68.036$ tỷ VNĐ), lao động tăng $92.39%$ ($421 \rightarrow 810$ người), chỉ ra rủi ro cốt lõi là tăng trưởng "nóng" làm loãng sự gắn kết tổ chức.
    • Cung cấp mô hình trọng số chuẩn hóa ($\beta$) giúp Ban Giám đốc phân bổ ngân sách phúc lợi và đào tạo chính xác theo mức độ ưu tiên: Đồng nghiệp $\rightarrow$ Đào tạo $\rightarrow$ Môi trường $\rightarrow$ Lãnh đạo $\rightarrow$ Tiền lương.

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

Kế hoạch giải pháp thực thi (Implementation Roadmap)

                            LỘ TRÌNH THỰC THI NÂNG CAO SỰ HÀI LÒNG
  ┌────────────────────────────────────────────────────────────────────────┐
  │ QUÝ 1: CẢI TẠO MÔI TRƯỜNG & VI KHÍ HẬU                                 │
  │ - Lắp đặt hệ thống thông gió xưởng may công nghiệp 1 & 2               │
  │ - Chuẩn hóa quy định an toàn vệ sinh lao động                          │
  ├────────────────────────────────────────────────────────────────────────┤
  │ QUÝ 2: ĐỒNG BỘ HÓA QUY TRÌNH ĐÁNH GIÁ & TIỀN LƯƠNG                     │
  │ - Ban hành KPI và thang bảng lương thưởng minh bạch theo bậc thợ        │
  │ - Rút ngắn chu kỳ xét tăng lương định kỳ 6 tháng/lần                   │
  ├────────────────────────────────────────────────────────────────────────┤
  │ QUÝ 3: CHƯƠNG TRÌNH ĐÀO TẠO & PHÁT TRIỂN LÃNH ĐẠO CẤP TRUNG           │
  │ - Tổ chức khóa đào tạo kỹ năng truyền thông cho Quản đốc & Chuyền trưởng│
  │ - Thiết lập lộ trình thăng tiến công khai (Career Path Matrix)         │
  ├────────────────────────────────────────────────────────────────────────┤
  │ QUÝ 4: HOÀN THIỆN PHÚC LỢI VÀ XÂY DỰNG VĂN HÓA DOANH NGHIỆP            │
  │ - Mở rộng gói bảo hiểm sức khỏe và tổ chức khám định kỳ                │
  │ - Định kỳ đo lường chỉ số hài lòng nhân viên tự động hóa               │
  └────────────────────────────────────────────────────────────────────────┘
  1. Cải thiện vi khí hậu và môi trường nhà xưởng (Tác động trực tiếp đến 79.4% LĐTT):
    • Lắp đặt thêm quạt hút công nghiệp, hệ thống làm mát Cooling Pad tại Xưởng may 1 và Xưởng may 2 nhằm duy trì nhiệt độ xưởng $< 28^\circ\text{C}$ trong mùa hè.
    • Đầu tư thiết bị bảo hộ lao động đạt chuẩn ISO 45001 (khẩu trang than hoạt tính chống bụi bông, nút tai chống ồn tại khu vực cắt và máy nén).
  2. Cơ chế tiền lương và phúc lợi minh bạch (Giải quyết khúc mắc của 70.7% lao động thu nhập $\le 5$ triệu):
    • Chuyển đổi định mức đơn giá tiền lương sản phẩm sang hệ thống số hóa, cho phép công nhân kiểm tra sản lượng hoàn thành và mức thu nhập tạm tính hàng ngày qua bảng tin điện tử.
    • Chi trả bảo hiểm xã hội, bảo hiểm y tế đúng hạn 100%, bổ sung quỹ trợ cấp đột xuất cho công nhân có hoàn cảnh khó khăn.
  3. Đào tạo chuyền trưởng và xây dựng văn hóa gắn kết:
    • Huấn luyện phong cách lãnh đạo hỗ trợ (Supportive Leadership) cho đội ngũ chuyền trưởng/tổ trưởng sản xuất, hạn chế tối đa việc áp đặt kỷ luật tiêu cực gây ức chế tâm lý người lao động.

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

  • Hạn chế kỹ thuật và phạm vi:
    • Quy mô mẫu dừng lại ở $N = 160$ tại một đơn vị sản xuất cụ thể (Hanex Huế), chưa mở rộng đối sánh đa trung tâm với các doanh nghiệp may mặc lân cận tại KCN Phú Bài.
    • Dữ liệu nghiên cứu là dữ liệu cắt ngang (Cross-sectional data), chưa theo dõi chuỗi thời gian (Longitudinal study) để kiểm chứng sự biến thiên của lòng trung thành sau khi áp dụng chính sách can thiệp.
  • Hướng nghiên cứu và phát triển tiếp theo:
    • Ứng dụng mô hình phương trình cấu trúc tuyến tính PLS-SEM (Partial Least Squares Structural Equation Modeling) trên phần mềm SmartPLS 4.0 để phân tích biến trung gian (Mediating variables) như "Sự gắn kết tổ chức" (Organizational Commitment) và "Động lực nội tại" (Intrinsic Motivation).
    • Tích hợp module People Analytics trên nền tảng ERP nhân sự để dự báo sớm tỷ lệ rời bỏ công việc bằng mô hình học máy (Logistic Regression, Random Forest Classifier).

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

  • Sinh viên & Học viên ngành HTTT Kinh tế / Quản trị Nhân lực: Tài liệu tham khảo chuẩn mực về phương pháp nghiên cứu định lượng, quy trình làm sạch dữ liệu, xử lý Cronbach's Alpha, EFA và hồi quy OLS trên SPSS/Python.
  • Data Analysts & HR Analytics Specialists: Kịch bản code tái tạo quy trình kiểm định Cronbach's Alpha, xoay nhân tố và đánh giá độ tin cậy của bảng hỏi.
  • Ban Giám đốc & Phòng Nhân sự doanh nghiệp sản xuất: Căn cứ khoa học định lượng để tái cơ cấu chính sách đãi ngộ, dự toán ngân sách cải thiện môi trường làm việc với mức ROI dự kiến cải thiện năng suất lao động thêm 8–12%.
  • Các nhà nghiên cứu kinh tế ứng dụng: Khung dữ liệu thực chứng về hành vi người lao động trong ngành công nghiệp phụ trợ dệt may miền Trung giai đoạn tăng trưởng quy mô.

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

1. Yêu cầu kỹ thuật tối thiểu để tái hiện quy trình phân tích dữ liệu của dự án?

Để chạy lại toàn bộ mô hình, hệ thống cần cấu hình tối thiểu: CPU Intel Core i3 thế hệ 8 trở lên (hoặc tương đương), 8GB RAM, hệ điều hành Windows 10/11 hoặc Linux Ubuntu 22.04 LTS; cài đặt IBM SPSS Statistics v22.0 trở lên hoặc môi trường Python 3.10+ tích hợp các thư viện pandas, statsmodels, factor_analyzer và scipy.

2. Vì sao biến PL1 (Khám sức khỏe, nghỉ mát) bị loại khỏi thang đo?

Biến PL1 có hệ số tương quan biến tổng (Corrected Item-Total Correlation) $= 0.192 < 0.3$. Theo tiêu chuẩn phương pháp luận (Nunnally & Bernstein, 1994), biến này không đóng góp vào tính đơn hướng của thang đo "Phúc lợi" và làm giảm hệ số Cronbach's Alpha tổng thể, do đó bắt buộc phải loại bỏ để đảm bảo tính giá trị (Validity).

3. Mô hình EFA rút trích được bao nhiêu phần trăm phương sai và có ý nghĩa gì?

Tại lần xoay EFA thứ 2 (23 biến), Tổng phương sai trích đạt $70.390% > 50%$ với $\text{KMO} = 0.757$ và kiểm định Bartlett $p < 0.001$. Điều này chứng minh rằng 6 nhóm nhân tố sau khi xoay giải thích được $70.39%$ độ biến thiên của sự hài lòng nhân viên tại công ty, $29.61%$ còn lại do các yếu tố ngẫu nhiên ngoài mô hình.

4. Giải pháp nào mang lại hiệu quả cao nhất để giữ chân nhân viên may mặc?

Kết quả hồi quy chuẩn hóa chỉ ra 2 nhân tố có tác động mạnh nhất là Mối quan hệ đồng nghiệp ($\beta = 0.342$) và Cơ hội đào tạo & thăng tiến ($\beta = 0.285$). Do đó, xây dựng văn hóa làm việc nhóm, tinh thần hỗ trợ tương trợ tại chuyền may và lộ trình nâng bậc thợ rõ ràng mang lại hiệu quả giữ chân người lao động cao hơn việc chỉ tăng lương đơn thuần.

5. Chi phí triển khai hệ thống giải pháp và thời gian hoàn vốn (ROI)?

Dự toán chi phí cải tạo hệ thống vi khí hậu xưởng may, chuẩn hóa quy trình đào tạo và bổ sung phúc lợi dao động khoảng 250–350 triệu VNĐ. Với việc giảm tỷ lệ nghỉ việc từ 18% xuống dưới 8%/năm (tiết kiệm chi phí tuyển dụng, bù đắp hao hụt năng suất ước tính 500–700 triệu VNĐ/năm), điểm hòa vốn đạt được sau 6–9 tháng triển khai.

Kết luận

Đề tài khóa luận tốt nghiệp của tác giả Trần Thị Kim Chi (Khoa Hệ thống Thông tin Kinh tế, Trường Đại học Kinh tế - Đại học Huế) đã xây dựng thành công khung phân tích định lượng toàn diện về sự hài lòng của người lao động tại Công ty TNHH MTV Hanex Huế. Thông qua việc xử lý dữ liệu thực nghiệm $N = 160$ với các công cụ kiểm định đa biến chuyên sâu (Cronbach's Alpha, EFA, OLS, Welch Robust Test), nghiên cứu đã chứng minh khoa học mức độ chi phối của 6 nhóm nhân tố cốt lõi đến sự thỏa mãn nghề nghiệp. Kết quả không chỉ đóng góp giá trị học thuật phương pháp luận cho ngành Hệ thống Thông tin Kinh tế mà còn trực tiếp cung cấp hệ thống giải pháp thực tiễn, giúp Hanex Huế ổn định nguồn nhân lực, tối ưu hóa năng suất và duy trì đà tăng trưởng bền vững trong kỷ nguyên công nghiệp hóa.