Giới thiệu dự án

Nghiên cứu ứng dụng và định lượng hóa các yếu tố tác động đến mức độ thỏa mãn trong công việc của người lao động đóng vai trò sống còn đối với các doanh nghiệp thâm dụng lao động tại Việt Nam. Theo thống kê của tổ chức Gallup, Việt Nam từng thuộc nhóm quốc gia có tỷ lệ nhân viên gắn kết và hài lòng với công việc thấp tại khu vực Châu Á. Trong bối cảnh ngành sản xuất thủ công mỹ nghệ và vật liệu đan lát sợi nhựa tổng hợp cả nước đóng góp 2,35 tỷ USD kim ngạch xuất khẩu (tăng trưởng bình quân 9,5% - 10%/năm) và giải quyết việc làm cho hơn 2,3 triệu lao động nông thôn, việc duy trì sự ổn định của lực lượng sản xuất trực tiếp là bài toán then chốt.

+---------------------------------------------------------------------------------------------------+
|                                   THỰC TRẠNG NGÀNH & DOANH NGHIỆP                                 |
+------------------------------------+--------------------------------------------------------------+
| Chỉ số ngành Thủ công Mỹ nghệ VN   | 2.556 làng nghề toàn quốc; Xuất khẩu đạt 2,35 tỷ USD        |
| Quy mô lao động nghiên cứu         | 396 nhân sự (88,5% công nhân trực tiếp, 6,5% văn phòng)       |
| Doanh thu Phước Hiệp Thành (2019)  | 119,98 tỷ VNĐ (Tăng trưởng doanh thu 3,31%, LN: 3,78 tỷ VNĐ)  |
| Chi phí vận hành/Doanh thu         | > 96,8% (Chi phí nhân công và nguyên vật liệu chiếm tỷ trọng)|
+------------------------------------+--------------------------------------------------------------+

Vấn đề cốt lõi (Problem Statement) tại Công ty Cổ phần 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ế) là sự gia tăng chi phí vận hành (đạt 116,2 tỷ VNĐ năm 2019) và áp lực giữ chân đội ngũ 396 lao động trực tiếp trong dây chuyền đan lát nhựa thủ công xuất khẩu sang thị trường Châu Âu (Eurofurnindo Limited), Mỹ, Hà Lan. Tỷ lệ biến động lao động phổ thông gây phát sinh chi phí tuyển dụng, suy giảm năng suất và rủi ro trễ tiến độ các đơn hàng xuất khẩu.

Dự án xác định 3 mục tiêu cụ thể:

  1. Hệ thống hóa cơ sở lý luận về hành vi tổ chức và sự thỏa mãn công việc dựa trên mô hình Chỉ số Mô tả Công việc (Job Descriptive Index - JDI), thuyết 2 nhân tố Herzberg và tháp nhu cầu Maslow.
  2. Đo lường, kiểm định định lượng mức độ ảnh hưởng của 6 nhóm nhân tố (Bản chất công việc, Đào tạo - thăng tiến, Lãnh đạo, Đồng nghiệp, Thu nhập, Điều kiện làm việc) tới sự hài lòng của 139 nhân sự được chọn mẫu phân tầng.
  3. Thiết lập hệ thống hàm hồi quy tuyến tính bội ($Y = \beta_0 + \sum \beta_i X_i$) và đề xuất ma trận giải pháp quản trị nhân sự tối ưu hóa chi phí - gia tăng mức độ gắn kết.

Phạm vi nghiên cứu tập trung vào toàn bộ khối văn phòng, đội ngũ vận tải và công nhân phân xưởng cơ khí - lò sơn - đan lát - hoàn thiện tại nhà máy Tứ Hạ, sử dụng dữ liệu thứ cấp giai đoạn 2017 - 2019 và dữ liệu sơ cấp điều tra thực địa.

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

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

Phương pháp đánh giá sự thỏa mãn nhân sự tại các doanh nghiệp sản xuất truyền thống thường dựa vào đánh giá cảm tính hoặc khảo sát định tính rời rạc, dẫn đến sai lệch trong việc phân bổ ngân sách phúc lợi và tiền lương.

Tiêu chí so sánh Khảo sát định tính truyền thống Thang đo Minnesota Satisfaction Questionnaire (MSQ) Mô hình JDI tích hợp định lượng (Đề xuất)
Công cụ xử lý Phiếu đánh giá tự do / Phỏng vấn Thang đo 20 khía cạnh (100 câu hỏi) Thang đo Likert 5 mức độ, xử lý SPSS 26.0
Độ phức tạp Thấp, dễ triển khai Quá dài, gây quá tải cho công nhân Tối ưu hóa 25 biến quan sát cho 6 nhân tố
Độ tin cậy Không thể kiểm định định lượng Cronbach's Alpha $> 0.80$ Cronbach's Alpha $\ge 0.70$, KMO $\ge 0.70$
Chi phí khảo sát Thấp nhưng không hiệu quả Rất cao, thời gian xử lý kéo dài Tối ưu hóa nguồn lực, phân tích hồi quy OLS
Tính ứng dụng Khó quy đổi thành KPI quản trị Phù hợp nghiên cứu hàn lâm Trực tiếp hỗ trợ thiết kế thang bảng lương

Phân tích yêu cầu hệ thống quản trị theo mô hình MoSCoW:

  • Must have: Thang đo chuẩn hóa cho 6 nhân tố độc lập; Kiểm định độ tin cậy Cronbach's Alpha; Phân tích nhân tố khám phá EFA; Kiểm định đa cộng tuyến qua hệ số phóng đại phương sai (Variance Inflation Factor - VIF < 2.0).
  • Should have: Kiểm định sự khác biệt One-Sample T-Test theo từng vị trí việc làm (Văn phòng, Lái xe, Công nhân đan lát).
  • Could have: Tích hợp thuật toán tự động phân tích dữ liệu trên nền tảng Python/R thay thế thao tác thủ công.
  • Won't have: Mô hình hóa phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) trong giai đoạn 1.

Thiết kế hệ thống

Kiến trúc luồng xử lý dữ liệu nghiên cứu định lượng được xây dựng qua pipeline 5 giai đoạn:

graph TD
    A["Dữ liệu thô từ bảng hỏi (N=139)"] --> B["Tiền xử lý & Mã hóa (Excel / SPSS 26)"]
    B --> C["Kiểm định độ tin cậy Thang đo (Cronbach's Alpha >= 0.7)"]
    C -->|Loại biến có Item-Total < 0.3| D["Phân tích Nhân tố Khám phá (EFA - Varimax)"]
    D -->|Kiểm định KMO >= 0.5 & Bartlett p < 0.05| E["Phân tích Tương quan Pearson & Hồi quy Tuyến tính Bội"]
    E --> F["Kiểm định Giả định: VIF < 2.0, Durbin-Watson, ANOVA F-test"]
    F --> G["Kiểm định One-Sample T-Test & Đề xuất Ma trận Quyết định"]

Quy chuẩn ngăn xếp công nghệ (Technology Stack) áp dụng trong đề tài:

  • Xử lý thống kê chuyên sâu: IBM SPSS Statistics v26.0 (Syntax Engine OLS/Factor Analysis).
  • Xử lý bảng tính & Làm sạch dữ liệu: Microsoft Excel 2019 / Office 365 ProPlus.
  • Ngôn ngữ kịch bản tự động hóa: Python 3.10.x với các thư viện pandas 2.1.0, numpy 1.24.0, statsmodels 0.14.0, factor_analyzer 0.5.0, scipy 1.11.0.

Thiết kế từ điển dữ liệu (Data Schema) mã hóa cho 25 biến quan sát thuộc 6 nhóm nhân tố:

[Nhân tố Độc lập]

[Nhân tố Phụ thuộc]

Phương pháp nghiên cứu

Nghiên cứu áp dụng quy trình định lượng đa biến. Công thức xác định kích thước mẫu tuân thủ nghiêm ngặt hai tiêu chuẩn thống kê:

  1. Theo Green (1991) cho hồi quy tuyến tính bội: $N \ge 104 + k = 104 + 6 = 110$ quan sát (với $k=6$ là số biến độc lập).
  2. Theo Hair et al. (2014) cho phân tích EFA: $N \ge 5 \times M' = 5 \times 25 = 125$ quan sát (với $M'=25$ là tổng số biến quan sát).

Với tỷ lệ phản hồi dự kiến đạt 90%, kích thước mẫu tối ưu được xác định là $N = 125 / 0,9 = 139$ mẫu. Phương pháp chọn mẫu phân tầng tỷ lệ (Stratified Proportional Sampling) được triển khai cụ thể:

$$\text{Cơ cấu mẫu} = \begin{cases} \text{Công nhân sản xuất trực tiếp:} & 123 \text{ mẫu } (88,49%) \ \text{Nhân viên văn phòng, kế toán:} & 9 \text{ mẫu } (6,47%) \ \text{Đội ngũ lái xe, vận tải:} & 7 \text{ mẫu } (5,04%) \end{cases}$$

Implementation và kết quả

Quy trình triển khai xử lý dữ liệu

Mã nguồn triển khai toàn bộ pipeline phân tích dữ liệu định lượng (Cronbach's Alpha, EFA, Hồi quy OLS, VIF và One-Sample T-Test) được cấu trúc bằng Python:

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

def calculate_cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha cho thang đo."""
    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 run_efa_pipeline(df_independent: pd.DataFrame):
    """Thực hiện kiểm định KMO, Bartlett và trích xuất nhân tố EFA."""
    chi_square, p_value = calculate_bartlett_sphericity(df_independent)
    kmo_per_variable, kmo_total = calculate_kmo(df_independent)
    
    fa = FactorAnalyzer(n_factors=6, rotation="varimax", method="principal")
    fa.fit(df_independent)
    
    loadings = pd.DataFrame(fa.loadings_, index=df_independent.columns)
    eigenvalues, _ = fa.get_eigenvalues()
    return kmo_total, p_value, eigenvalues, loadings

def fit_ols_regression(X: pd.DataFrame, y: pd.Series):
    """Xây dựng mô hình hồi quy tuyến tính bội OLS và kiểm tra VIF."""
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X.columns
    vif_data["VIF"] = [
        variance_inflation_factor(X.values, i) for i in range(X.shape[1])
    ]
    return model, vif_data

def run_one_sample_ttest(df_factors: pd.DataFrame, test_value: float = 3.0):
    """Kiểm định One-Sample T-Test so sánh giá trị trung bình với mức 3.0."""
    results = {}
    for col in df_factors.columns:
        t_stat, p_val = stats.ttest_1samp(df_factors[col].dropna(), popmean=test_value)
        mean_val = df_factors[col].mean()
        results[col] = {"Mean": mean_val, "T-Statistic": t_stat, "p-value": p_val}
    return pd.DataFrame(results).T

Kiểm định và đánh giá thống kê

Kết quả phân tích thống kê trên bộ dữ liệu $N=139$ đạt các chỉ số kiểm định chuẩn:

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

    • Thang đo Đào tạo & Thăng tiến (DT): $\alpha = 0,842$; Tất cả hệ số tương quan biến tổng (Corrected Item-Total Correlation) $> 0,45$.
    • Thang đo Lãnh đạo (LD): $\alpha = 0,876$; Corrected Item-Total Correlation $> 0,52$.
    • Thang đo Thu nhập (TN): $\alpha = 0,815$; Corrected Item-Total Correlation $> 0,41$.
    • Thang đo Điều kiện làm việc (DK): $\alpha = 0,798$; Corrected Item-Total Correlation $> 0,38$.
    • Thang đo Đồng nghiệp (DN): $\alpha = 0,785$; Corrected Item-Total Correlation $> 0,35$.
    • Thang đo Biến phụ thuộc Sự hài lòng (HL): $\alpha = 0,854$.

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

    • Kiểm định KMO đạt $0,812$ ($> 0,5$), kiểm định Bartlett có mức ý nghĩa $p = 0,000 < 0,05$.
    • Trích xuất 6 nhân tố với chỉ số Eigenvalue $> 1,0$ tại điểm dừng; Tổng phương sai trích (Total Variance Explained) đạt $63,48% (> 50%)$.
    • Hệ số tải nhân tố (Factor Loading) của tất cả 25 biến quan sát đều $> 0,55$, không có biến bị phân tán chéo.
+----------------------------------------------------------------------------------------------------+
|                                KẾT QUẢ PHÂN TÍCH HỒI QUY TUYẾN TÍNH BỘI                            |
+---------------------+-------------------+-------------------+-------------+---------+--------------+
| Biến độc lập        | Hệ số Chưa chuẩn  | Hệ số Chuẩn hóa   | Giá trị t   | Sig.    | Hệ số VIF    |
|                     | hóa (B)           | Beta ($\beta$)    |             | (p-val) | (< 2.0)      |
+---------------------+-------------------+-------------------+-------------+---------+--------------+
| (Hằng số / Const)   | 0,312             | --                | 1,215       | 0,226   | --           |
| Thu nhập (TN)       | 0,285             | 0,314             | 4,512       | 0,000   | 1,245        |
| Lãnh đạo (LD)       | 0,241             | 0,268             | 3,891       | 0,000   | 1,312        |
| Điều kiện LV (DK)   | 0,195             | 0,215             | 3,124       | 0,002   | 1,189        |
| Đào tạo - TT (DT)   | 0,162             | 0,178             | 2,654       | 0,009   | 1,204        |
| Đồng nghiệp (DN)    | 0,134             | 0,145             | 2,145       | 0,034   | 1,156        |
| Bản chất CV (BC)    | 0,098             | 0,108             | 1,985       | 0,049   | 1,112        |
+---------------------+-------------------+-------------------+-------------+---------+--------------+
| R = 0,782 | R² = 0,611 | R² hiệu chỉnh = 0,594 | F(6, 132) = 34,612 (p = 0,000) | Durbin-Watson = 1,876 |
+----------------------------------------------------------------------------------------------------+

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

Hàm hồi quy chuẩn hóa biểu diễn sự tác động của các nhân tố đến sự hài lòng chung:

$$\text{HL} = 0,314 \times \text{TN} + 0,268 \times \text{LD} + 0,215 \times \text{DK} + 0,178 \times \text{DT} + 0,145 \times \text{DN} + 0,108 \times \text{BC}$$

Mô hình giải thích được $59,4%$ ($R^2 \text{ hiệu chỉnh} = 0,594$) sự biến thiên của mức độ thỏa mãn trong công việc. Kết quả kiểm định One-Sample T-Test với giá trị kiểm định $\mu_0 = 3,0$ (Mức trung lập) cho thấy:

  • Nhân tố Thu nhập đạt điểm trung bình thấp nhất ($\bar{X} = 3,18 / 5,0$), thể hiện nhu cầu bức thiết của công nhân về cải thiện đơn giá sản phẩm đan lát.
  • Nhân tố Lãnh đạo ($\bar{X} = 3,65$) và Đồng nghiệp ($\bar{X} = 3,82$) được đánh giá tích cực nhất, khẳng định môi trường gắn kết nội bộ cao.

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

Đổi mới phương pháp luận ứng dụng

Khác với các nghiên cứu học thuật thuần túy chỉ sử dụng mô hình JDI nguyên bản (vốn được thiết kế cho môi trường công nghiệp phương Tây), đề tài đã tinh chỉnh hệ thống biến quan sát phù hợp với đặc thù sản xuất tiểu thủ công nghiệp tại Miền Trung:

  1. Tách biệt rõ ràng chi phí lương cố định và lương khoán sản phẩm trong biến đo lường Thu nhập ($\text{TN}$), phản ánh đúng đặc điểm trả lương theo mét đan/sản phẩm hoàn thiện.
  2. Lồng ghép tiêu chuẩn an toàn hóa chất tại xưởng sơn tĩnh điện và xử lý bụi sợi nhựa tổng hợp vào thang đo Điều kiện làm việc ($\text{DK}$).
+----------------------------------------------------------------------------------------------------+
|                         SO SÁNH MÔ HÌNH NGHIÊN CỨU VỚI CÁC CÔNG TRÌNH TIỀN NHIỆM                   |
+----------------------+-------------------------+-------------------------+-------------------------+
| Đặc tính kỹ thuật    | JDI gốc (Smith, 1969)   | Trần Kim Dung (2005)    | Mô hình Đề tài (2021)   |
+----------------------+-------------------------+-------------------------+-------------------------+
| Đối tượng khảo sát   | Đa ngành tổng quát      | Doanh nghiệp TP.HCM     | Doanh nghiệp đan lát    |
|                      |                         |                         | sợi nhựa xuất khẩu      |
| Cấu trúc nhân tố     | 5 nhóm nhân tố          | 7 nhóm nhân tố          | 6 nhóm nhân tố tối ưu   |
| Số lượng quan sát    | 72 câu hỏi Có/Không     | 34 câu hỏi Likert 5     | 25 câu hỏi Likert 5     |
| Độ phức tạp bảng hỏi | Cao, mất 25 phút        | Trung bình, 15 phút     | Tối ưu, 7-10 phút       |
| Kiểm soát VIF        | Không đề cập            | Có kiểm soát            | VIF < 1,35 (Rất thấp)   |
| Khả năng nhân rộng   | Thấp đối với công nhân  | Trung bình              | Rất cao cho làng nghề   |
+----------------------+-------------------------+-------------------------+-------------------------+

Đóng góp thực tiễn cho ngành

Nghiên cứu mang lại đóng góp cụ thể về mặt định lượng:

  • Chứng minh nhân tố Thu nhập có trọng số tác động lớn nhất ($\beta = 0,314$), vượt qua các yếu tố tinh thần trong việc giữ chân thợ thủ công.
  • Chỉ ra rằng việc cải thiện điều kiện vi khí hậu nhà xưởng (thông gió, chống nóng cho phân xưởng đan lát 123 công nhân) mang lại hiệu quả gia tăng sự hài lòng cao gấp 2 lần so với việc tái cấu trúc bản chất công việc ($\beta_{DK} = 0,215$ vs $\beta_{BC} = 0,108$).

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

Kịch bản ứng dụng trong doanh nghiệp

Mô hình và kết quả phân tích được tích hợp vào hệ thống quản trị nhân sự của Công ty Cổ phần Phước Hiệp Thành qua 3 giải pháp trọng tâm:

  1. Cải tiến chính sách Đơn giá & Thu nhập:

    • Thiết lập bảng định mức đơn giá lũy tiến cho các đơn hàng xuất khẩu có độ khó cao (Ghế sofa đôi, Bộ ghế thư giãn ngoài trời theo tiêu chuẩn Greenguard).
    • Tăng tỷ lệ thưởng chuyên cần thêm 15% đối với các công đoạn có độ phức tạp cao, giúp tăng thu nhập bình quân của công nhân từ mức 4,5 triệu lên 5,8 triệu VNĐ/tháng.

  2. Nâng cấp Điều kiện lao động & Vi khí hậu xưởng sản xuất:

    • Lắp đặt hệ thống quạt thông gió công nghiệp tại phân xưởng đan lát và buồng hút lọc bụi sơn tại bộ phận lò sơn tĩnh điện.
    • Đầu tư trang thiết bị bảo hộ lao động đạt chuẩn kiểm định an toàn môi trường theo Nghị định 19/2015/NĐ-CP.

  3. Chuẩn hóa quy trình Đào tạo & Thăng tiến:

    • Tổ chức mô hình "Thợ bậc cao kèm thợ mới", rút ngắn thời gian làm quen công việc từ 30 ngày xuống 12 ngày.
+----------------------------------------------------------------------------------------------------+
|                                 LỘ TRÌNH TRIỂN KHAI GIẢI PHÁP QUẢN TRỊ                             |
+-----------+-----------------------------------+-----------------------------------+----------------+
| Giai đoạn | Nội dung công việc                | Kết quả đầu ra (Deliverables)     | Thời gian      |
+-----------+-----------------------------------+-----------------------------------+----------------+
| Phase 1   | Khảo sát vi khí hậu & rà soát     | Báo cáo đánh giá mức độ nóng/bụi; | Quý 1/2021     |
|           | đơn giá công đoạn đan lát         | Dự thảo bảng lương khoán mới      | (3 tháng)      |
| Phase 2   | Lắp đặt hệ thống thông gió & ban  | Nghiệm thu 8 quạt hút công nghiệp;| Quý 2-3/2021   |
|           | hành quy chế thưởng năng suất     | Tăng 12% năng suất xuất xưởng     | (6 tháng)      |
| Phase 3   | Đánh giá sau can thiệp (Post-test)| Chỉ số thỏa mãn HL tăng > 18%;    | Quý 4/2021     |
|           | & chuẩn hóa quy trình đào tạo     | Tỷ lệ nghỉ việc giảm xuống < 4%   | (3 tháng)      |
+-----------+-----------------------------------+-----------------------------------+----------------+

Đánh giá hiệu quả kinh tế (Cost-Benefit Analysis)

  • Chi phí triển khai giải pháp: 185 triệu VNĐ (Nâng cấp quạt thông gió, đồ bảo hộ, quỹ đào tạo và hiệu chỉnh phần mềm quản lý chấm công).
  • Lợi ích kinh tế ước tính:
    • Giảm tỷ lệ nghỉ việc từ 8,5% xuống dưới 3,5%/năm $\rightarrow$ Tiết kiệm 120 triệu VNĐ chi phí tuyển dụng và đào tạo thợ mới.
    • Tăng năng suất lao động thêm 6,8% $\rightarrow$ Đóng góp trực tiếp vào mức tăng trưởng lợi nhuận thuần (đạt trên 4,2 tỷ VNĐ trong niên độ tiếp theo).
    • Tỷ suất hoàn vốn đầu tư (ROI) ước tính đạt $148%$ sau 12 tháng vận hành.

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

Dự án ghi nhận các giới hạn kỹ thuật và phương pháp:

  • Kích thước mẫu: Mẫu điều tra $N=139$ dù đạt chuẩn thống kê Green (1991) và Hair (2014) nhưng mới chỉ giới hạn trong một nhà máy tại Thừa Thiên Huế, chưa bao phủ toàn bộ các chi nhánh hợp tác xã liên kết.
  • Phương pháp lấy mẫu: Mẫu phân tầng tỷ lệ vẫn có rủi ro sai lệch chọn lựa ngẫu nhiên do đặc thù ca kíp của công nhân.
  • Phương pháp lượng hóa: Sử dụng mô hình OLS truyền thống chỉ đánh giá được tác động trực tiếp, chưa phân tích được các biến trung gian (Mediating Variables) như Cam kết tổ chức (Organizational Commitment) hay Động lực nội tại (Intrinsic Motivation).

Hướng phát triển nghiên cứu tiếp theo:

  • Ứng dụng mô hình Phương trình Cấu trúc Tuyến tính (Structural Equation Modeling - PLS-SEM) trên SmartPLS 4.0 để phân tích mối quan hệ đa tầng giữa Sự thỏa mãn $\rightarrow$ Cam kết gắn bó $\rightarrow$ Năng suất thực tế.
  • Xây dựng mô hình học máy phân loại (Machine Learning Classification - Random Forest / XGBoost) dự báo sớm nguy cơ nghỉ việc của công nhân dựa trên dữ liệu chấm công và biến động sản lượng hàng tháng.

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

  • Sinh viên & Học viên cao học: Tiếp cận khung phương pháp nghiên cứu định lượng chuẩn tắc từ khâu tính cỡ mẫu, thiết kế bảng hỏi Likert, đến xử lý dữ liệu thực chứng trên SPSS/Python.
  • Lập trình viên & Chuyên viên Phân tích Dữ liệu (Data Analysts): Kế thừa mã nguồn Python tự động hóa toàn bộ quy trình kiểm định Cronbach's Alpha, EFA, OLS Regression và tính toán ma trận VIF.
  • Chủ doanh nghiệp & Nhà quản trị Nhân sự (HR Managers): Nắm bắt bằng chứng số học chính xác về mức độ nhạy cảm của người lao động đối với thu nhập và điều kiện làm việc để hoạch định chính sách giữ chân nhân tài.
  • Các nhà nghiên cứu kinh tế ứng dụng: Tài liệu tham khảo về thực trạng chuỗi cung ứng và quản trị nguồn nhân lực ngành sản xuất thủ công mỹ nghệ sợi nhựa tổng hợp Việt Nam.

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

  1. Điều kiện kỹ thuật để tái lập pipeline phân tích dữ liệu của đề tài là gì?
    Cần cài đặt môi trường Python 3.8+ với các gói thư viện pandas, statsmodels, factor-analyzer, scipy hoặc phần mềm IBM SPSS Statistics phiên bản 20.0 trở lên.

  2. Tại sao nhân tố Bản chất công việc có hệ số Beta thấp nhất ($\beta = 0,108$)?
    Do 88,5% lực lượng lao động là công nhân phổ thông thực hiện các thao tác đan lát thủ công có tính lặp đi lặp lại. Đối với nhóm lao động này, yếu tố thu nhập và điều kiện an toàn, môi trường làm việc có tác động chi phối mạnh mẽ hơn bản thân tính chất phức tạp của công việc.

  3. Mô hình có gặp hiện tượng đa cộng tuyến nghiêm trọng không?
    Không. Tất cả các hệ số phóng đại phương sai (VIF) của 6 biến độc lập đều nằm trong khoảng $1,112 - 1,312$, thấp hơn rất nhiều so với ngưỡng cảnh báo $2,0$ (hoặc $5,0$ theo chuẩn lỏng), chứng minh các biến độc lập không có hiện tượng cộng tuyến vi phạm giả định OLS.

  4. Làm thế nào để đảm bảo tính trung thực của câu trả lời từ công nhân trực tiếp?
    Bảng hỏi được thiết kế ẩn danh hoàn toàn, tiến hành thu thập trực tiếp trong giờ nghỉ giao ca, giải thích rõ mục đích nghiên cứu học thuật không ảnh hưởng đến đánh giá thi đua cá nhân.

  5. Chi phí đầu tư nâng cấp điều kiện làm việc có ảnh hưởng đến biên lợi nhuận của công ty không?
    Với tổng chi phí đầu tư khoảng 185 triệu VNĐ trên nền doanh thu 119,98 tỷ VNĐ (chiếm chưa đầy $0,15%$ chi phí hoạt động), khoản đầu tư này nhanh chóng được bù đắp bởi mức giảm chi phí phế phẩm và chi phí tuyển dụng lao động mới.

Kết luận

Đề tài đã thực hiện thành công việc lượng hóa các yếu tố tác động đến sự hài lòng của người lao động tại Công ty Cổ phần Phước Hiệp Thành thông qua bộ công cụ thống kê đa biến chuẩn tắc. Kết quả chứng minh vai trò dẫn dắt của nhân tố Thu nhập ($\beta = 0,314$) và Lãnh đạo ($\beta = 0,268$), đồng thời xác lập mô hình hồi quy đạt độ tin cậy cao với $R^2 \text{ hiệu chỉnh} = 59,4%$.

Các hàm ý quản trị và lộ trình triển khai 3 giai đoạn mang lại giải pháp thực tiễn giúp ban lãnh đạo công ty nâng cao năng lực cạnh tranh, ổn định nguồn nhân lực phục vụ các đơn hàng xuất khẩu quốc tế, đóng góp mô hình tham chiếu có giá trị cho hơn 2.500 làng nghề và doanh nghiệp sản xuất tiểu thủ công nghiệp trên cả nước.