Giới thiệu dự án

Ngành dệt may Việt Nam đóng vai trò then chốt trong nền kinh tế xuất khẩu quốc gia, đóng góp trên 15% tổng kim ngạch xuất khẩu cả nước nhưng đồng thời phải đối mặt với áp lực cạnh tranh gay gắt về năng suất và tỷ lệ chuyển dịch lao động (turnover rate dao động từ 15-25%/năm). Tại Công ty Cổ phần May Xuất khẩu Ngọc Châu (Thừa Thiên Huế) – đơn vị gia công 100% hàng xuất khẩu sang các thị trường khó tính như Hoa Kỳ, Nhật Bản và Hàn Quốc – năng lực cạnh tranh phụ thuộc trực tiếp vào tay nghề và động lực cống hiến của lực lượng lao động trực tiếp.

                                  +-----------------------------------------+
                                  |    Mô hình Động lực Lao động (Herzberg) |
                                  +-----------------------------------------+
                                                       |
         +---------------------------------------------+---------------------------------------------+
         |                                                                                           |
         v                                                                                           v
+-----------------------------------+                                       +-----------------------------------+
|     Nhân tố Duy trì (Hygiene)     |                                       |    Nhân tố Động viên (Motivator)  |
| - Lương, thưởng & phúc lợi (LTPL) |                                       | - Trách nhiệm với tổ chức (TNTC)  |
| - Môi trường làm việc (MTLV)      |                                       | - Bố trí & sắp xếp CV (BTCV)      |
| - Bố trí công việc (BTCV)         |                                       | - Cơ hội thăng tiến (CH)          |
+-----------------------------------+                                       | - Sự công nhận đóng góp (CNDG)    |
                                                                            +-----------------------------------+
                                                       |
                                                       v
                                    +-------------------------------------+
                                    | Động lực làm việc nhân viên (DLLV)  |
                                    | R2 = 0.570 | Adjusted R2 = 0.553   |
                                    +-------------------------------------+

Vấn đề nghiên cứu thực tiễn

Mặc dù doanh thu thuần năm 2017 đạt 17.348 triệu đồng (tăng 4,38% so với 2016) và lợi nhuận sau thuế phục hồi đạt 889 triệu đồng, công ty phải đối mặt với tình trạng biến động quy mô nhân sự liên tục giảm từ 332 người (2015) xuống 326 người (2016) và còn 320 người (2017). Trong đó, lao động nữ chiếm 82%, nhân viên sản xuất trực tiếp chiếm 93,75% và 68,75% là lao động phổ thông chưa qua đào tạo may công nghiệp bài bản. Áp lực về cường độ làm việc, tiếng ồn, nhiệt độ nhà xưởng và chính sách đãi ngộ chưa đồng bộ khiến mức độ gắn kết suy giảm, tạo nguy cơ chuyển dịch lao động sang các doanh nghiệp cạnh tranh trong khu vực như May An Phước hay May Việt Tiến.

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

  1. Hệ thống hóa cơ sở lý thuyết về động lực lao động dựa trên Thuyết hai nhân tố của Herzberg, Thuyết nhu cầu của Maslow và Thuyết ba nhu cầu của McClelland.
  2. Kiểm định mô hình đo lường 7 nhân tố ảnh hưởng đến động lực làm việc: Lương thưởng & phúc lợi (LTPL), Môi trường làm việc (MTLV), Bố trí và sắp xếp công việc (BTCV), Sự hứng thú trong công việc (HT), Cơ hội thăng tiến (CH), Sự công nhận đóng góp (CNDG), và Trách nhiệm cá nhân với tổ chức (TNTC).
  3. Định lượng mức độ tác động của từng nhân tố đến biến phụ thuộc Động lực làm việc (DLLV) bằng phân tích hồi quy OLS (Ordinary Least Squares).
  4. Kiểm định sự khác biệt về động lực làm việc theo các biến nhân khẩu học (giới tính, độ tuổi, trình độ học vấn, mức thu nhập) qua Independent Samples T-Test và One-Way ANOVA.
  5. Đề xuất hệ thống giải pháp quản trị nhân sự thực chứng nhằm tối ưu hóa năng suất và duy trì tỷ lệ gắn kết bền vững.

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

Nghiên cứu kết hợp phương pháp định tính (nghiên cứu tài liệu, hiệu chỉnh bảng hỏi sơ bộ) và định lượng thực chứng (khảo sát trực tiếp mẫu ngẫu nhiên phân tầng $N = 180$, làm sạch dữ liệu, kiểm định độ tin cậy thang đo Cronbach's Alpha, phân tích nhân tố khám phá EFA, tương quan Pearson và mô hình hồi quy tuyến tính đa biến).

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

  • Xây dựng mô hình hồi quy đạt độ thích hợp $R^2 > 0.50$, hệ số Durbin-Watson đạt ngưỡng tối ưu $[1, 3]$ và hiện tượng đa cộng tuyến được loại trừ ($VIF < 3$).
  • Xác lập trọng số $\beta$ chuẩn hóa cho từng nhóm chính sách quản trị, cung cấp bộ công cụ định lượng hỗ trợ Ban Giám đốc ra quyết định nhân sự.

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

  • Không gian: Nghiên cứu khảo sát nội bộ tại Công ty Cổ phần May Xuất khẩu Ngọc Châu (Phú Lộc, Thừa Thiên Huế).
  • Thời gian: Số liệu thứ cấp tài chính - nhân sự thu thập từ 2015 đến 2017; số liệu sơ cấp điều tra từ ngày 01/02/2018 đến 20/04/2018.

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 các mô hình đo lường động lực làm việc thực nghiệm trong và ngoài nước để lựa chọn khung tham chiếu:

Mô hình nghiên cứu Tác giả & Năm Đối tượng khảo sát Ưu điểm kỹ thuật Hạn chế thực nghiệm
Thuyết Hai nhân tố Mở rộng Boeve (2007) Giảng viên y khoa tại Hoa Kỳ Phân tách rõ nhân tố nội tại và ngoại tại Mẫu mang tính hàn lâm, không phản ánh môi trường gia công cơ chế chuyền may
Mô hình Bán lẻ Đa biến Teck-Hong & Waheed (2011) Nhân viên kinh doanh tại Malaysia Khai thác sâu điều kiện làm việc và đãi ngộ tài chính Trọng số điều kiện làm việc chiếm ưu thế tuyệt đối do đặc thù ngành bán lẻ
Mô hình Cán bộ Quản lý Lê Thanh Dũng (2007) Quản lý cấp trung và cao Bao quát 12 biến số tổ chức 3 biến bị loại khi kiểm định, không áp dụng được cho công nhân trực tiếp
Mô hình May Công nghiệp Đề xuất hiện tại (2018) Công nhân may & khối văn phòng Đo lường toàn diện 7 khía cạnh sản xuất và tâm lý Cần hiệu chỉnh thang đo lương thưởng theo thực tế trả lương sản phẩm

Ma trận ưu tiên yêu cầu giải pháp (MoSCoW)

  • Must have: Hệ thống chuẩn hóa mã hóa dữ liệu 36 biến quan sát, kiểm định Cronbach's Alpha $> 0.70$, loại biến rác có tương quan biến - tổng $< 0.30$.
  • Should have: Phân tích EFA trích xuất phương sai $> 50%$, hệ số Factor Loading $> 0.50$, kiểm tra phân phối chuẩn phần dư.
  • Could have: Kiểm định sự khác biệt bình quân nhóm qua ANOVA đồng nhất phương sai (Levene Test).
  • Won't have: Mô hình hóa phi tuyến SEM-PLS (tập trung tối ưu hóa mô hình hồi quy đa biến OLS chuẩn mực trên SPSS).
                      +---------------------------------------------+
                      | Dữ liệu thô từ phiếu khảo sát (N = 180)     |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      | Kiểm định Cronbach's Alpha (Loại LTPL1)     |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      | Phân tích EFA (KMO = 0.868, Bartlett Sig=0) |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      | Tương quan Pearson & Kiểm tra giả định OLS  |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      | Hồi quy OLS Đa biến -> Mô hình hiệu chỉnh   |
                      +---------------------------------------------+

Thiết kế hệ thống

Ngăn xếp công nghệ xử lý dữ liệu (Technology Stack)

  • Công cụ phân tích thống kê: IBM SPSS Statistics version 20.0 (Core Engine for Descriptive, Reliability, EFA, OLS Regression, ANOVA).
  • Hệ thống bảng tính phụ trợ: Microsoft Excel 2016 (Data Entry, Cleansing, Pivot Charting).
  • Môi trường tái lập thuật toán (Code Simulation): Python 3.9 (Thư viện pandas 1.5.3, scikit-learn 1.2.2, statsmodels 0.14.0, factor_analyzer 0.4.1).

Thiết kế cấu trúc thang đo và mã hóa dữ liệu

Dữ liệu được thu thập qua thang đo Likert 5 mức độ (1 = Rất không đồng ý đến 5 = Rất đồng ý), gồm 32 biến độc lập thuộc 7 nhóm và 3 biến phụ thuộc:

Thang đo độc lập (32 items)

Thang đo phụ thuộc (3 items)

Methodology

Quy trình nghiên cứu áp dụng mô hình thác nước chuẩn hóa (Standardized Waterfall Research Workflow) chia thành 5 giai đoạn:

  1. Giai đoạn 1 (Tuần 1-3): Nghiên cứu lý luận, tổng hợp cơ sở dữ liệu thứ cấp giai đoạn 2015-2017 từ Phòng Kế toán và Phòng Nhân sự.
  2. Giai đoạn 2 (Tuần 4-6): Thiết kế bảng câu hỏi sơ bộ, xác định kích thước mẫu theo công thức của Hair et al. (1998): $N \ge 5 \times m = 5 \times 36 = 180$.
  3. Giai đoạn 3 (Tuần 7-9): Tiến hành thu thập dữ liệu sơ cấp theo phương pháp hạn ngạch phân tầng (169 công nhân sản xuất chiếm 93,89% và 11 nhân viên văn phòng chiếm 6,11%).
  4. Giai đoạn 4 (Tuần 10-11): Nhập liệu, làm sạch dữ liệu, xử lý thống kê bằng phần mềm SPSS 20.0.
  5. Giai đoạn 5 (Tuần 12): Thảo luận kết quả hồi quy, hoàn thiện báo cáo và chuyển giao khuyến nghị.

Implementation và kết quả

Development process

Quá trình phân tích dữ liệu thực nghiệm được hiện thực hóa qua quy trình phân tích định lượng chặt chẽ. Dưới đây là code snippet mô phỏng thuật toán đánh giá độ tin cậy thang đo và hồi quy OLS tương thích với cấu trúc của SPSS:

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

def cronbach_alpha(df_items: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha cho thang đo đa biến."""
    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_ols_model(X: pd.DataFrame, y: pd.Series):
    """Thực thi hồi quy OLS và kiểm tra hệ số phóng đại phương sai (VIF)."""
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    
    # Tính VIF cho từng biến độc lập
    vif_data = pd.DataFrame()
    vif_data["Variable"] = X.columns
    vif_data["VIF"] = [
        variance_inflation_factor(X_with_const.values, i + 1)
        for i in range(X.shape[1])
    ]
    
    return model.summary(), vif_data

Testing và validation

1. Kiểm định độ tin cậy thang đo (Reliability Analysis)

Thang đo sơ bộ gồm 36 biến quan sát được kiểm định bằng hệ số Cronbach's Alpha:

  • Biến bị loại: Biến LTPL1 ("Tiền lương đảm bảo chi trả đủ các nhu cầu cơ bản") có tương quan biến - tổng hiệu chỉnh là $-0.008 < 0.30$; nếu loại LTPL1, hệ số Cronbach's Alpha của thang đo LTPL tăng từ $0.851$ lên $0.879$. Biến LTPL1 bị loại bỏ để đảm bảo tính nhất quán nội tại.
  • Kết quả sau hiệu chỉnh:
    • LTPL (6 biến): $\alpha = 0.879$
    • MTLV (8 biến): $\alpha = 0.895$
    • BTCV (5 biến): $\alpha = 0.817$
    • HT (4 biến): $\alpha = 0.849$
    • CH (3 biến): $\alpha = 0.807$
    • CNDG (3 biến): $\alpha = 0.859$
    • TNTC (3 biến): $\alpha = 0.848$
    • DLLV (3 biến phụ thuộc): $\alpha = 0.813$ Tất cả các thang đo đều đạt chuẩn xuất sắc ($\alpha > 0.80$, tương quan biến - tổng $> 0.50$).

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

  • Biến độc lập:
    • Hệ số Kaiser-Meyer-Olkin ($KMO$) = $0.868$ (thỏa mãn $0.5 \le KMO \le 1.0$).
    • Kiểm định Bartlett: $\chi^2 = 3636.056, df = 496, p = 0.000 < 0.001$.
    • 7 nhân tố được trích xuất với Eigenvalues = $1.089 > 1.0$, tổng phương sai trích (Total Variance Explained) đạt $69.403% > 50%$.
    • Ma trận xoay nhân tố Varimax: Tất cả 32 biến quan sát đều có hệ số tải nhân tố (Factor Loading) $\ge 0.519$ (vượt xa ngưỡng tối thiểu $0.50$).
  • Biến phụ thuộc:
    • $KMO = 0.697$, Bartlett's Test $\chi^2 = 187.495 (df=3, p=0.000)$, Eigenvalues = $2.185$, tổng phương sai trích = $72.845%$. Hệ số tải các biến DLLV1, DLLV2, DLLV3 dao động từ $0.815$ đến $0.886$.

3. Kiểm định các giả định mô hình hồi quy tuyến tính

  • Phân phối chuẩn phần dư: Biểu đồ Histogram có $Mean = -5.76 \times 10^{-15} \approx 0$, độ lệch chuẩn $Std. Dev = 0.980 \approx 1$. Đồ thị Normal P-P Plot thể hiện các điểm phân vị bám sát đường chéo kỳ vọng.
  • Tính tuyến tính và phương sai không đổi: Scatterplot giữa phần dư chuẩn hóa và giá trị dự đoán chuẩn hóa phân bố ngẫu nhiên tập trung quanh trục hoành độ 0.
  • Tự tương quan chuỗi bậc nhất: Hệ số Durbin-Watson = $2.128$ (nằm trong khoảng an toàn $[1.0, 3.0]$).
  • Đa cộng tuyến: Tất cả các biến độc lập có hệ số $VIF$ dao động từ $1.346$ đến $1.918$ ($VIF < 3.0$), chứng minh không có hiện tượng đa cộng tuyến.
                           +----------------------------------------+
                           |   Ma trận Phân tích Hồi quy Đa biến   |
                           +----------------------------------------+
                                                |
          +-------------------------------------+-------------------------------------+
          |                                                                           |
          v                                                                           v
+-----------------------------------+                               +-----------------------------------+
|     Mô hình Thống kê Tổng thể     |                               |      Kiểm định Giả thuyết         |
| - R = 0.755                       |                               | - TNTC: Beta=0.282 (p=0.000) [OK] |
| - R2 = 0.570                      |                               | - BTCV: Beta=0.194 (p=0.002) [OK] |
| - Adjusted R2 = 0.553             |                               | - CH:   Beta=0.170 (p=0.010) [OK] |
| - F = 32.598 (p = 0.000)          |                               | - MTLV: Beta=0.169 (p=0.004) [OK] |
| - Durbin-Watson = 2.128           |                               | - CNDG: Beta=0.154 (p=0.019) [OK] |
+-----------------------------------+                               | - LTPL: Beta=0.121 (p=0.038) [OK] |
                                                                    | - HT:   Beta=-0.062 (p=0.309)[LOẠI]|
                                                                    +-----------------------------------+

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

Bảng kết quả ước lượng tham số hồi quy OLS từ tập dữ liệu $N = 180$:

Nhân tố Hệ số chưa chuẩn hóa ($B$) Sai số chuẩn ($SE$) Hệ số chuẩn hóa ($\beta$) Giá trị $t$ Mức ý nghĩa ($Sig.$) Dung nạp ($Tolerance$) $VIF$ Kết luận giả thuyết
(Hằng số) 1.038 0.199 - 5.218 0.000 - - Có ý nghĩa thống kê
TNTC 0.195 0.048 0.282 4.078 0.000 0.522 1.918 Chấp nhận $H_7$
BTCV 0.148 0.047 0.194 3.173 0.002 0.670 1.493 Chấp nhận $H_4$
CH 0.113 0.051 0.170 2.605 0.010 0.589 1.699 Chấp nhận $H_6$
MTLV 0.093 0.032 0.169 2.909 0.004 0.741 1.350 Chấp nhận $H_1$
CNDG 0.123 0.052 0.154 2.374 0.019 0.591 1.693 Chấp nhận $H_5$
LTPL 0.065 0.031 0.121 2.091 0.038 0.743 1.346 Chấp nhận $H_2$
HT -0.048 0.047 -0.062 -1.021 0.309 0.688 1.454 Bác bỏ $H_3$

Phương trình hồi quy tuyến tính chuẩn hóa thu được: $$\text{DLLV} = 0.282 \times \text{TNTC} + 0.194 \times \text{BTCV} + 0.170 \times \text{CH} + 0.169 \times \text{MTLV} + 0.154 \times \text{CNDG} + 0.121 \times \text{LTPL}$$

  • Hệ số xác định hiệu chỉnh $R^2_{adj} = 0.553$: 6 nhân tố giữ lại giải thích được $55.3%$ sự biến động của động lực làm việc ở nhân viên công ty; $44.7%$ còn lại do các yếu tố ngẫu nhiên ngoài mô hình.
  • Biến HT (Hứng thú trong công việc) bị loại bỏ: Do mức ý nghĩa $p = 0.309 > 0.05$. Điều này phản ánh tính chất đặc thù của lao động ngành may gia công: công việc vận hành theo dây chuyền mang tính chất lặp đi lặp lại thao tác đơn điệu, khiến sự hứng thú cá nhân không trở thành điều kiện tiên quyết kích thích động lực làm việc so với ý thức trách nhiệm và chế độ phân công.

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

Nghiên cứu mang lại những đóng góp học thuật và thực tiễn đáng kể cho lĩnh vực quản trị nhân lực ngành dệt may:

  1. Phát hiện nghịch lý nhân tố "Hứng thú công việc": Khác với các nghiên cứu văn phòng của Nguyễn Khắc Hoàn (2010) tại Ngân hàng ACB (nơi hứng thú công việc tác động mạnh), nghiên cứu chỉ ra rằng trong môi trường sản xuất gia công công nghiệp, biến HT không có ý nghĩa thống kê ($p = 0.309$). Sự gắn bó của công nhân chủ yếu được dẫn dắt bởi Trách nhiệm với tổ chức ($\beta = 0.282$) và Bố trí công việc hợp lý ($\beta = 0.194$).
  2. Loại trừ biến đo lường không phù hợp bằng quy chuẩn thực chứng: Việc loại bỏ biến LTPL1 giúp tăng độ tin cậy thang đo lương thưởng từ $\alpha = 0.851$ lên $\alpha = 0.879$, loại trừ sai số đo lường trong môi trường lao động có thu nhập dao động theo đơn giá sản phẩm.
  3. Cung cấp thang đo chuẩn hóa cho ngành dệt may miền Trung: Bộ thang đo 32 biến quan sát sau EFA có tổng phương sai trích $69.403%$ là công cụ tin cậy có thể tái sử dụng cho các nghiên cứu tiếp theo trong khối doanh nghiệp gia công xuất khẩu.

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

Dựa trên kết quả hồi quy chuẩn hóa, hệ thống giải pháp can thiệp quản trị được đề xuất theo mức độ ưu tiên trọng số:

                          +------------------------------------------+
                          |   Lộ trình 4 Giai đoạn Nâng cao Động lực |
                          +------------------------------------------+
                                               |
         +--------------------+--------------------+--------------------+
         |                    |                    |                    |
         v                    v                    v                    v
+--------------------+ +--------------------+ +--------------------+ +--------------------+
|  Phase 1 (Tháng 1) | |  Phase 2 (Tháng 2) | |  Phase 3 (Tháng 3) | |  Phase 4 (Tháng 4) |
| Chuẩn hóa định mức | | Cải tạo vi khí hậu | | Thiết lập lộ trình | | Tối ưu hóa đơn giá|
| chuyền & KCS (BTCV)| | nhà xưởng (MTLV)   | | thăng tiến (CH)    | | lương thưởng (LTPL)|
+--------------------+ +--------------------+ +--------------------+ +--------------------+

1. Nâng cao trách nhiệm cá nhân đối với tổ chức ($\beta = 0.282$)

  • Giải pháp: Xây dựng văn hóa trách nhiệm chủ động bằng cách minh bạch hóa kế hoạch giao hàng xuất khẩu cho từng tổ chuyền may. Tổ chức đối thoại định kỳ giữa Ban Giám đốc và công nhân sản xuất để gắn mục tiêu cá nhân với sự phát triển của nhà máy.
  • Chỉ số đo lường: Tỷ lệ công nhân tham gia họp giao ban chuyền đạt $100%$, giảm $30%$ số lỗi thao tác do thiếu trách nhiệm kiểm tra ban đầu.

2. Tối ưu hóa bố trí và sắp xếp công việc ($\beta = 0.194$)

  • Giải pháp: Ứng dụng kỹ thuật cân bằng chuyền may (Line Balancing) nhằm giảm thiểu thời gian chờ công đoạn (bottlenecks). Bố trí công nhân có tay nghề bậc cao kèm cặp thợ mới vào các công đoạn phức tạp (may cổ áo sơ mi, tra tay, ráp sườn quần tây).
  • Chỉ số đo lường: Hiệu suất cân bằng chuyền tăng $12-18%$, giảm thời gian ứ đọng bán thành phẩm $25%$.

3. Xây dựng lộ trình thăng tiến và đào tạo ($\beta = 0.170$)

  • Giải pháp: Thiết lập tiêu chuẩn nâng bậc thợ rõ ràng từ Bậc 1 đến Bậc 6; xây dựng cơ chế bổ nhiệm Chuyền trưởng, KCS từ chính lực lượng công nhân có tay nghề giỏi.
  • Chỉ số đo lường: $100%$ công nhân được đánh giá tay nghề định kỳ 6 tháng/lần; tối thiểu $15%$ công nhân đạt tiêu chuẩn được xem xét nâng ngạch lương.

4. Cải thiện môi trường làm việc ($\beta = 0.169$)

  • Giải pháp: Đầu tư hệ thống quạt thông gió làm mát áp suất âm, bổ sung màn chắn cách âm tại khu vực máy may tự động và bố trí hệ thống đèn LED chống lóa đạt chuẩn chiếu sáng công nghiệp.
  • Chi phí dự toán: Khoảng 180-250 triệu đồng cho 2 phân xưởng Ngọc Châu 1 và Ngọc Châu 2.

5. Hoàn thiện chính sách công nhận và đãi ngộ ($\beta = 0.154$ và $\beta = 0.121$)

  • Giải pháp: Khen thưởng kịp thời "Công nhân xuất sắc của tháng" với phần thưởng hiện vật/tiền mặt công khai tại bảng tin phân xưởng. Cải tiến đơn giá tiền lương theo sản phẩm lũy tiến vượt định mức và đảm bảo chất lượng bữa ăn trưa đủ dinh dưỡng.

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

Hạn chế học thuật và thực tiễn

  1. Phương pháp chọn mẫu thuận tiện: Dù cỡ mẫu $N = 180$ đạt chuẩn kỹ thuật ($N \ge 5 \times 36$), việc áp dụng chọn mẫu thuận tiện trong từng bộ phận có thể mang thiên lệch chủ quan nhất định.
  2. Giới hạn biến số ngoài mô hình: Mô hình giải thích $55.3%$ biến thiên; còn $44.7%$ chưa được giải thích có thể thuộc về các nhân tố như: Phong cách lãnh đạo của chuyền trưởng, Văn hóa doanh nghiệp đặc thù, hoặc Yếu tố hoàn cảnh gia đình của lao động nữ.
  3. Thiết kế nghiên cứu cắt ngang (Cross-sectional): Dữ liệu thu thập tại một thời điểm (đầu năm 2018), chưa theo dõi được sự thay đổi động lực theo tính chu kỳ của các đơn hàng may xuất khẩu (mùa cao điểm xuất khẩu vào quý 3 và quý 4).

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

  • Triển khai phương pháp nghiên cứu theo chuỗi thời gian (Longitudinal Study) để đánh giá tác động sau khi ban hành các chính sách nhân sự mới.
  • Ứng dụng mô hình cấu trúc tuyến tính SEM (Structural Equation Modeling) để phân tích các mối quan hệ gián tiếp và biến trung gian (như Sự hài lòng công việc, Cam kết gắn bó tổ chức).

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

  • Sinh viên & Học viên cao học chuyên ngành Quản trị Kinh doanh / Quản trị Nhân lực: Tài liệu tham khảo chuẩn mực về quy trình phân tích thực chứng, từ xây dựng bảng hỏi Likert, làm sạch dữ liệu, kiểm định Cronbach's Alpha, EFA đến hồi quy đa biến OLS trên SPSS.
  • Chuyên viên phân tích dữ liệu & Nhà nghiên cứu định lượng: Cung cấp bộ số liệu thực chứng, các chỉ số thống kê đối sánh ($KMO = 0.868$, $R^2_{adj} = 0.553$) và cách thức xử lý loại biến rác (LTPL1, HT).
  • Ban Lãnh đạo & Giám đốc Nhân sự các Doanh nghiệp Dệt may: Bộ khung giải pháp can thiệp quản trị có cơ sở khoa học để tái phân bổ ngân sách nhân sự (ưu tiên nâng cao trách nhiệm tổ chức và cân bằng chuyền hơn là các biện pháp kích thích hứng thú ngắn hạn).
  • Tập thể Người lao động tại Công ty May Ngọc Châu: Hưởng lợi trực tiếp từ các cải tiến về môi trường làm việc, định mức lao động công bằng, cơ hội thăng tiến minh bạch và chế độ phúc lợi nâng cao.

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

1. Tại sao nhân tố "Sự hứng thú trong công việc" (HT) lại bị loại khỏi mô hình hồi quy?

Trong bảng phân tích hồi quy, biến HT có mức ý nghĩa kiểm định $t$ là $p = 0.309 > 0.05$. Điều này chứng minh rằng sự hứng thú không có tác động mang ý nghĩa thống kê đến động lực làm việc của nhân viên May Ngọc Châu. Lý do thực tiễn: Ngành may gia công xuất khẩu có quy trình sản xuất cơ giới hóa cao, thao tác lặp lại theo dây chuyền khép kín, công nhân chú trọng tính chuẩn xác và thu nhập theo sản phẩm hơn là tìm kiếm sự mới lạ trong công việc.

2. Việc loại bỏ biến quan sát LTPL1 dựa trên cơ sở khoa học nào?

Biến LTPL1 có hệ số tương quan biến - tổng hiệu chỉnh là $-0.008$ (nhỏ hơn ngưỡng chấp nhận $0.30$). Khi loại biến này, hệ số tin cậy Cronbach's Alpha của nhóm thang đo Lương thưởng & phúc lợi tăng từ $0.851$ lên $0.879$. Việc loại bỏ này loại trừ biến rác không nhất quán với tổng thể thang đo.

3. Mô hình có gặp hiện tượng đa cộng tuyến hoặc tự tương quan không?

Không. Hệ số phóng đại phương sai của tất cả các biến độc lập đều nằm trong khoảng $1.346 \le VIF \le 1.918$ (nhỏ hơn ngưỡng cảnh báo 3.0), chứng minh không có đa cộng tuyến. Đồng thời, hệ số Durbin-Watson đạt $2.128$ (nằm trong khoảng an toàn từ $1.0$ đến $3.0$), khẳng định mô hình không bị vi phạm tự tương quan chuỗi bậc nhất.

4. Kết quả nghiên cứu có thể mở rộng cho các ngành nghề khác không?

Mô hình có thể mở rộng tốt cho các doanh nghiệp sản xuất gia công thâm dụng lao động (như giày da, lắp ráp linh kiện điện tử). Tuy nhiên, đối với các ngành dịch vụ, tài chính - ngân hàng hoặc công nghệ thông tin đòi hỏi tính sáng tạo cao, thang đo cần được bổ sung và kiểm định lại trọng số của nhân tố HT (Hứng thú công việc).

5. Ban quản lý nên ưu tiên đầu tư vào yếu tố nào trước để đạt hiệu quả chi phí tối ưu?

Căn cứ vào hệ số $\beta$ chuẩn hóa, Ban quản lý nên tập trung ưu tiên hàng đầu vào Trách nhiệm của cá nhân với tổ chức ($\beta = 0.282$) và Bố trí, sắp xếp công việc ($\beta = 0.194$). Đây là hai nhóm giải pháp phi tài chính, ít tốn kém chi phí đầu tư trực tiếp nhưng mang lại tỷ suất cải thiện động lực làm việc cao nhất.

Kết luận

Nghiên cứu "Các yếu tố ảnh hưởng đến động lực làm việc của nhân viên tại Công ty Cổ phần May Xuất khẩu Ngọc Châu" đã kiểm định thành công mô hình động lực lao động trong môi trường may công nghiệp xuất khẩu. Thông qua tập dữ liệu khảo sát thực chứng $N = 180$ và các thuật toán phân tích định lượng chuẩn mực, nghiên cứu đã chứng minh 6 nhân tố tác động cùng chiều có ý nghĩa thống kê đến động lực làm việc, xếp theo thứ tự trọng số giảm dần: Trách nhiệm đối với tổ chức ($\beta = 0.282$), Bố trí và sắp xếp công việc ($\beta = 0.194$), Cơ hội thăng tiến và phát triển ($\beta = 0.170$), Môi trường làm việc ($\beta = 0.169$), Sự công nhận đóng góp của cá nhân ($\beta = 0.154$) và Lương thưởng & phúc lợi ($\beta = 0.121$), giải thích $55.3%$ sự biến thiên của động lực làm việc ($R^2_{adj} = 0.553, F = 32.598, p < 0.001$).

Kết quả nghiên cứu cung cấp cơ sở khoa học xác đáng để Ban Lãnh đạo Công ty May Ngọc Châu chuyển đổi từ phương thức quản trị cảm tính sang quản trị nhân sự dựa trên dữ liệu thực chứng (Data-Driven HR Management), đảm bảo sự tăng trưởng bền vững và nâng cao năng lực cạnh tranh trong chuỗi cung ứng dệt may toàn cầu.