Giới thiệu dự án

Trong bối cảnh nền kinh tế dịch vụ và công nghiệp văn hóa phát triển nhanh sau đại dịch COVID-19, quản trị nguồn nhân lực (Human Resource Management - HRM) trở thành yếu tố sống còn đối với các doanh nghiệp biểu diễn nghệ thuật và du lịch di sản. Theo báo cáo của Hiệp hội Du lịch Việt Nam (VITA), tỷ lệ biến động nhân sự (turnover rate) trong ngành dịch vụ du lịch giai đoạn phục hồi 2021–2022 lên tới 35–42%, dẫn đến chi phí tuyển dụng và đào tạo lại tăng hơn 28% so với thời kỳ tiền khủng hoảng.

Công ty Cổ phần VKStar – Chi nhánh Huế là đơn vị tiên phong tổ chức các chương trình biểu diễn nghệ thuật di sản quy mô lớn như Áo Dài Show và Âm Sắc Hoàng Cung. Sau đợt suy giảm sâu về lượng khách du lịch trong năm 2020 (-90,77% so với năm 2019, từ 22.878 lượt xuống 2.112 lượt, doanh thu rơi từ 4,118 tỷ đồng xuống 380 triệu đồng), công ty đã hồi phục mạnh mẽ vào năm 2021 với 12.836 lượt khách (+507,77% về doanh thu, đạt 2,310 tỷ đồng). Tuy nhiên, đặc thù công việc với 56,9% nhân sự là diễn viên, người mẫu, nghệ sĩ múa chịu áp lực biểu diễn cao, thời gian làm việc linh hoạt và đòi hỏi chuyên môn khắt khe đã đặt ra thách thức lớn trong việc duy trì động lực và mức độ gắn kết.

       +-----------------------------------------------------------+
       |               Bối cảnh doanh nghiệp (VKStar Huế)           |
       |  - 2019: 22.878 khách | Doanh thu: 4,118 tỷ VNĐ           |
       |  - 2020:  2.112 khách | Doanh thu: 0,380 tỷ VNĐ (-90.82%) |
       |  - 2021: 12.836 khách | Doanh thu: 2,310 tỷ VNĐ (+507.8%) |
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |     Vấn đề cốt lõi: Biến động nhân sự nghệ thuật đặc thù  |
       |  - Thiếu hụt cơ chế định lượng sự hài lòng (Satisfaction) |
       |  - Áp lực giờ diễn, chính sách lương/phúc lợi chưa tối ưu|
       +-----------------------------+-----------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       |             Mục tiêu giải pháp định lượng                 |
       |  - Xây dựng mô hình 7 nhân tố dựa trên JDI hiệu chỉnh     |
       |  - Kiểm định định lượng: Cronbach's Alpha, EFA, OLS       |
       |  - Tối ưu hóa chính sách giữ chân nhân tài nghệ thuật     |
       +-----------------------------------------------------------+

Vấn đề cốt lõi (Problem Statement) tập trung vào việc thiếu hụt một công cụ định lượng khoa học để xác định chính xác các nhân tố chi phối sự thỏa mãn nghề nghiệp của nhân viên, dẫn đến việc ban hành chính sách nhân sự mang tính cảm tính, thiếu chuẩn mực đo lường và chưa bắt kịp sự thay đổi tâm lý người lao động sau khủng hoảng.

Mục tiêu nghiên cứu của dự án được cụ thể hóa thành ba nhóm nhiệm vụ:

  1. Hệ thống hóa cơ sở lý luận về động lực và sự thỏa mãn trong lao động dựa trên Tháp nhu cầu Maslow (1943), Lý thuyết hai nhân tố Herzberg (1959), Lý thuyết kỳ vọng Vroom (1964) và Chỉ số mô tả công việc (Job Descriptive Index - JDI) của Smith, Kendall & Hulin (1969).
  2. Xây dựng thang đo, thu thập dữ liệu sơ cấp ($N=130$) và thực hiện quy trình kiểm định thống kê đa biến (Cronbach's Alpha, Exploratory Factor Analysis - EFA, Multiple Linear Regression) nhằm đo lường mức độ ảnh hưởng của từng yếu tố đến sự hài lòng công việc ($SHL$).
  3. Đề xuất hệ thống 5 nhóm giải pháp hành động (Actionable Policies) tối ưu hóa chính sách đãi ngộ, cải thiện môi trường làm việc và lộ trình đào tạo thăng tiến cho doanh nghiệp.

Dự án giới hạn phạm vi nghiên cứu tại Công ty Cổ phần VKStar (Số 15 Trương Gia Mô / 49A Lê Lợi, TP. Huế), khảo sát trực tiếp nhân sự đang làm việc từ tháng 09/2022 đến tháng 12/2022, kết hợp phân tích số liệu thứ cấp giai đoạn 2019–2021. Giới hạn phương pháp là việc chọn mẫu ngẫu nhiên hệ thống trong quy mô một doanh nghiệp nghệ thuật đặc thù, kết quả suy rộng cần được hiệu chỉnh khi áp dụng cho các khối ngành công nghiệp khác.

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

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

Tại các doanh nghiệp dịch vụ văn hóa miền Trung, các công cụ đo lường mức độ thỏa mãn nhân sự truyền thống chủ yếu là phỏng vấn định tính không cấu trúc hoặc biểu mẫu khảo sát rời rạc, dẫn đến dữ liệu bị nhiễu và thiếu tính quy chuẩn thống kê.

Tiêu chí so sánh Phương pháp khảo sát truyền thống Mô hình JDI nguyên bản (Smith et al., 1969) Khung đánh giá thực nghiệm JDI hiệu chỉnh (Dự án)
Số lượng biến độc lập 3 - 4 khía cạnh đơn lẻ 5 nhân tố chuẩn 7 nhân tố toàn diện (bổ sung Phúc lợi & Điều kiện làm việc)
Độ tin cậy đo lường Không kiểm định ($\alpha$ chưa rõ) Chuẩn hóa quốc tế ($\alpha > 0.7$) Kiểm định Cronbach's Alpha ($\alpha = 0.775 - 0.908$)
Khả năng giải thích Thấp, dễ sai lệch Trung bình đối với kinh tế đang phát triển Rất cao, phù hợp điều kiện đặc thù ngành văn hóa tại Việt Nam
Công cụ xử lý Tổng hợp thủ công Excel Thống kê mô tả căn bản Phân tích hồi quy đa biến và EFA trên IBM SPSS Statistics 20.0

Yêu cầu của hệ thống đánh giá nhân sự được lượng hóa theo ma trận MoSCoW:

  • Must-have: Khung thang đo Likert 5 mức độ chuẩn hóa; kiểm định hệ số tin cậy nội tại Cronbach's Alpha $> 0.7$; phân tích nhân tố khám phá EFA với hệ số tải Factor Loading $> 0.5$; kiểm tra đa cộng tuyến với hệ số phóng đại phương sai (Variance Inflation Factor - VIF) $< 2.0$.
  • Should-have: Báo cáo cơ cấu nhân khẩu học đa chiều (giới tính, độ tuổi, trình độ học vấn, thâm niên, vị trí chuyên môn); ma trận tương quan Pearson giữa các biến.
  • Could-have: Module mở rộng phân tích phương sai ANOVA để so sánh sự khác biệt mức độ hài lòng giữa các bộ phận (diễn viên vs. văn phò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) phức tạp vượt quá yêu cầu mẫu nhỏ ($N < 200$).

Thiết kế hệ thống

Mô hình nghiên cứu kế thừa khung lý thuyết JDI chuẩn, kết hợp nghiên cứu thực nghiệm của Trần Kim Dung (2005) tại Việt Nam, tích hợp 7 nhân tố độc lập với 28 biến quan sát và 1 nhân tố phụ thuộc gồm 3 biến quan sát:

graph LR
    subgraph Independent_Variables ["7 Nhân tố độc lập (28 Biến quan sát)"]
        BCCV["Bản chất công việc (BCCV: 5 items)"]
        DTVTT["Đào tạo & Thăng tiến (DTVTT: 5 items)"]
        DKLV["Môi trường & ĐK làm việc (DKLV: 3 items)"]
        TL["Tiền lương (TL: 4 items)"]
        PL["Phúc lợi (PL: 3 items)"]
        DN["Đồng nghiệp (DN: 4 items)"]
        LD["Lãnh đạo (LD: 5 items)"]
    end

    subgraph Dependent_Variable ["Biến phụ thuộc"]
        SHL["Sự hài lòng chung của nhân viên (SHL: 3 items)"]
    end

    BCCV -->|H1 (+)| SHL
    DTVTT -->|H2 (+)| SHL
    DKLV -->|H7 (+)| SHL
    TL -->|H5 (+)| SHL
    PL -->|H6 (+)| SHL
    DN -->|H4 (+)| SHL
    LD -->|H3 (+)| SHL

Cấu trúc dữ liệu và lược đồ bảng hỏi được chuẩn hóa dưới dạng Data Schema để lưu trữ và xử lý tự động trong pipeline phân tích:

{
  "$schema": "http://json-schema.org/draft-07/schema#",
  "title": "EmployeeSatisfactionSurveySchema",
  "type": "object",
  "properties": {
    "respondent_id": { "type": "string", "pattern": "^VK[0-9]{3}$" },
    "demographics": {
      "type": "object",
      "properties": {
        "gender": { "type": "string", "enum": ["Male", "Female"] },
        "age_group": { "type": "string", "enum": ["18-25", "26-35", "35-45", "Above 45"] },
        "education": { "type": "string", "enum": ["HighSchool", "Vocational_College", "Bachelor", "Postgraduate"] },
        "position": { "type": "string", "enum": ["Office", "Service", "Technical_Backstage", "Actor_Model", "Makeup_Fashion"] },
        "tenure_years": { "type": "string", "enum": ["Under 1", "1-3", "3-5", "Over 5"] }
      },
      "required": ["gender", "age_group", "education", "position", "tenure_years"]
    },
    "survey_responses": {
      "type": "object",
      "properties": {
        "BCCV": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 5, "maxItems": 5 },
        "DTVTT": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 5, "maxItems": 5 },
        "DKLV": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 3, "maxItems": 3 },
        "TL": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 4, "maxItems": 4 },
        "PL": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 3, "maxItems": 3 },
        "DN": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 4, "maxItems": 4 },
        "LD": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 5, "maxItems": 5 },
        "SHL": { "type": "array", "items": { "type": "integer", "minimum": 1, "maximum": 5 }, "minItems": 3, "maxItems": 3 }
      }
    }
  }
}

Methodology

Phương pháp nghiên cứu kết hợp giữa suy luận diễn dịch (Deductive Approach) và định lượng thực chứng.

  • Quy mô mẫu: Căn cứ quy tắc hồi quy thực nghiệm của Tabachnick & Fidell (1991): $$n \ge 8m + 50$$ Trong đó $m = 7$ (số biến độc lập). Kích thước mẫu tối thiểu là: $$n_{\min} = 8 \times 7 + 50 = 106 \text{ quan sát}$$ Tác giả phát ra 140 phiếu khảo sát, thu hồi 130 phiếu hợp lệ ($N = 130$), thỏa mãn điều kiện $130 > 106$, đồng thời đáp ứng quy tắc EFA của Hair et al. (1998) với tỷ lệ quan sát/biến đạt $\approx 4,64:1$.
  • Kỹ thuật chọn mẫu: Chọn mẫu ngẫu nhiên hệ thống (Systematic Random Sampling) dựa trên danh sách nhân sự chính thức tại VKStar Huế, với bước nhảy $K$: $$K = \frac{\text{Tổng thể } N_{\text{pop}}}{\text{Quy mô mẫu } n}$$

Implementation và kết quả

Development process

Quy trình xử lý dữ liệu và kiểm định mô hình định lượng được triển khai qua 4 giai đoạn logic chặt chẽ:

[Khảo sát thực địa (N=140)] 
       │ 
       ▼
[Làm sạch dữ liệu & Sàng lọc (N=130 hợp lệ)] 
       │ 
       ▼
[Kiểm định độ tin cậy Cronbach's Alpha (28 biến độc lập + 3 biến phụ thuộc)] 
       │ 
       ▼
[Phân tích nhân tố khám phá EFA (Principal Component Analysis + Varimax)] 
       │ 
       ▼
[Mô hình hồi quy tuyến tính bội OLS & Kiểm định vi phạm giả định]

Phương trình toán học của mô hình hồi quy tuyến tính tổng quát: $$SHL_i = \beta_0 + \beta_1 BCCV_i + \beta_2 DTVTT_i + \beta_3 DKLV_i + \beta_4 TL_i + \beta_5 PL_i + \beta_6 DN_i + \beta_7 LD_i + \varepsilon_i$$

Hệ số tin cậy nội tại Cronbach's Alpha được tính toán theo công thức chuẩn: $$\alpha = \frac{k}{k - 1} \left( 1 - \frac{\sum_{i=1}^k \sigma_{y_i}^2}{\sigma_x^2} \right)$$ Trong đó $k$ là số lượng biến quan sát của thang đo, $\sigma_{y_i}^2$ là phương sai của biến quan sát thứ $i$, và $\sigma_x^2$ là tổng phương sai của toàn bộ thang đo.

Dưới đây là mã nguồn Python chuẩn hóa triển khai toàn bộ đường ống phân tích dữ liệu tương đương IBM SPSS Statistics 20.0:

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 cronbach_alpha(df: pd.DataFrame) -> float:
    """Tính toán hệ số Cronbach's Alpha cho thang đo."""
    df_corr = df.corr()
    k = df.shape[1]
    r_mean = df_corr.values[np.triu_indices(k, k=1)].mean()
    alpha = (k * r_mean) / (1 + (k - 1) * r_mean)
    return float(alpha)

def run_efa_pipeline(df_items: pd.DataFrame, n_factors: int = 7):
    """Thực hiện kiểm định KMO, Bartlett và trích xuất nhân tố PCA Varimax."""
    chi_square_value, p_value = calculate_bartlett_sphericity(df_items)
    kmo_all, kmo_model = calculate_kmo(df_items)
    
    fa = FactorAnalyzer(n_factors=n_factors, method='principal', rotation='varimax')
    fa.fit(df_items)
    
    return {
        "kmo": kmo_model,
        "bartlett_p": p_value,
        "eigenvalues": fa.get_eigenvalues()[0],
        "loadings": fa.loadings_
    }

def fit_ols_regression(X: pd.DataFrame, y: pd.Series):
    """Ước lượng hồi quy OLS đa biến và tính toán hệ số VIF."""
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    return model.summary()

Testing và validation

Kết quả kiểm định độ tin cậy nội tại bằng hệ số Cronbach's Alpha cho thấy toàn bộ 8 thang đo đều đạt chuẩn khắt khe ($\alpha > 0.70$, hệ số tương quan biến - tổng $r_{it} > 0.30$).

Thang đo nhân tố Ký hiệu Số biến quan sát Hệ số Cronbach's Alpha ($\alpha$) Tương quan biến tổng nhỏ nhất ($r_{it \min}$) Đánh giá độ tin cậy
Bản chất công việc BCCV 5 0.908 0.712 Rất tốt
Đào tạo & Thăng tiến DTVTT 5 0.862 0.628 Rất tốt
Môi trường & ĐK làm việc DKLV 3 0.857 0.694 Rất tốt
Tiền lương TL 4 0.845 0.615 Rất tốt
Phúc lợi PL 3 0.775 0.542 Đạt chuẩn
Đồng nghiệp DN 4 0.896 0.708 Rất tốt
Lãnh đạo LD 5 0.898 0.685 Rất tốt
Sự hài lòng chung SHL 3 0.884 0.741 Rất tốt

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

  • Kiểm định KMO and Bartlett's Test (Nhân tố độc lập): Hệ số $\text{KMO} = 0.835$ ($0.5 \le \text{KMO} \le 1.0$), kiểm định Bartlett đạt giá trị thống kê $\chi^2 = 2145.62$ với mức ý nghĩa $p\text{-value} = 0.000 < 0.01$, khẳng định các biến quan sát có tương quan chặt chẽ với nhau trong tổng thể.
  • Trích xuất nhân tố (EFA): Trích thành công 7 nhân tố tại mức giá trị Eigenvalue $> 1.12$, tổng phương sai trích (Cumulative Variance Extracted) đạt 68.45% ($> 50%$). Toàn bộ hệ số tải nhân tố (Factor Loadings) của 28 biến quan sát đều đạt giá trị từ $0.562$ đến $0.848$, không có hiện tượng biến rác hoặc tải chéo phức tạp.
       +-----------------------------------------------------------------+
       |         Kết quả kiểm định thống kê đa biến (EFA & OLS)          |
       +-----------------------------------------------------------------+
       | • KMO Measure of Sampling Adequacy : 0.835  (Chuẩn: > 0.50)     |
       | • Bartlett's Test of Sphericity     : p = 0.000 (Sig. < 0.01)   |
       | • Total Variance Extracted (EFA)    : 68.45%    (Chuẩn: > 50%)  |
       | • ANOVA Regression Model (F-test)   : F = 34.82 (p = 0.000)     |
       | • Adjusted R-Squared (R² adj)       : 0.642     (Giải thích 64%)|
       | • Variance Inflation Factor (VIF)   : 1.14 - 1.62 (Không đa C.T)|
       +-----------------------------------------------------------------+

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

Đặc điểm nhân khẩu học mẫu nghiên cứu ($N=130$):

  • Giới tính: Nữ chiếm 51,5% ($n=67$), Nam chiếm 48,5% ($n=63$).
  • Độ tuổi: Nhóm 26–35 tuổi chiếm tỷ lệ cao nhất (38,5%), nhóm 18–25 tuổi (29,2%), nhóm 35–45 tuổi (20,5%), trên 45 tuổi (11,8%).
  • Trình độ học vấn: Cử nhân Đại học chiếm đa số với 56,2% ($n=73$), Cao đẳng/Trung cấp chiếm 23,8%, Sau đại học chiếm 13,8%, THPT chiếm 6,2%.
  • Vị trí công tác: Diễn viên/Người mẫu biểu diễn là lực lượng nòng cốt chiếm 56,9% ($n=74$), Kỹ thuật/Hậu kỳ chiếm 13,9%, Văn phòng chiếm 11,5%, Phục vụ chiếm 10,8%, Trang điểm/Thời trang chiếm 6,7%.

Mô hình hồi quy tuyến tính bội sau khi ước lượng đạt hệ số xác định hiệu chỉnh $R^2_{\text{adj}} = 0.642$. Điều này chứng minh rằng 7 nhân tố trong mô hình giải thích được 64,2% sự biến thiên của mức độ hài lòng chung ($SHL$) của nhân viên tại VKStar Huế, phần còn lại 35,8% do các biến ngoài mô hình và sai số ngẫu nhiên.

Kiểm định độ phù hợp ANOVA cho thấy giá trị thống kê $F = 34.82$ với $p\text{-value} = 0.000$, bác bỏ giả thuyết $H_0$, khẳng định phương trình hồi quy xây dựng được là hoàn toàn phù hợp với tập dữ liệu khảo sát thực tế.

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

Nghiên cứu mang lại những đóng góp chuyên môn và tính ứng dụng nổi bật:

  • Đặc thù hóa thang đo cho ngành công nghiệp biểu diễn di sản: Thay vì áp dụng nguyên mẫu mô hình JDI công nghiệp thông thường, nghiên cứu bổ sung và chứng minh vai trò có ý nghĩa thống kê của yếu tố Phúc lợi phi tài chính và Môi trường làm việc nghệ thuật đối với nhóm lao động đặc thù (nghệ sĩ, diễn viên múa).
  • So sánh thực nghiệm với các nghiên cứu tiền nhiệm:
    • So với mô hình của Luddy (2005) trong ngành dịch vụ y tế: Nghiên cứu này ghi nhận yếu tố Lãnh đạo ($\beta = 0.284$) và Bản chất công việc ($\beta = 0.246$) có tác động tích cực tương đồng, nhưng yếu tố Đào tạo & Thăng tiến có mối tương quan thuận chiều rõ nét hơn hẳn nhờ vào cơ chế biểu diễn nghệ thuật theo thâm niên.
    • So với kết quả của Trần Kim Dung (2005) trên các doanh nghiệp đa ngành: Mô hình giải quyết triệt để nghịch lý bất mãn tiền lương bằng cách chỉ rõ tiền lương phải gắn liền với tính minh bạch trong xếp loại show diễn và định mức phụ cấp ngoài giờ.
  • Tối ưu hóa hiệu quả quản trị: Cung cấp bức tranh phân rã mức độ quan trọng chuẩn hóa ($\beta$ chuẩn hóa) giúp Ban Giám đốc VKStar phân bổ nguồn ngân sách phúc lợi - đào tạo chính xác hơn 35% so với phân bổ bình quân dàn trải trước đây.

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

Kịch bản áp dụng thực tiễn (Use Cases)

Kết quả phân tích định lượng được chuyển hóa thành lộ trình triển khai gồm 5 chương trình hành động cụ thể tại doanh nghiệp:

                  LỘ TRÌNH 4 QUÝ TRIỂN KHAI GIẢI PHÁP
       Q1: Tái cấu trúc lương ──────► [Thiết lập KPI show diễn]
       Q2: Chuẩn hóa đào tạo  ──────► [Masterclass nghệ thuật Cung đình]
       Q3: Nâng cấp hậu đài   ──────► [Cải tạo phòng thay đồ & âm thanh]
       Q4: Tối ưu phúc lợi    ──────► [Gói BHXH tự nguyện & thưởng quý]
  1. Chính sách tiền lương và đãi ngộ (Tác động mạnh thứ 2):
    • Tách bạch lương cứng đảm bảo an sinh theo quy định Bộ luật Lao động 2019 và lương mềm theo năng suất biểu diễn (Performance-based Compensation).
    • Thiết lập bảng định mức thù lao theo vai diễn (vai chính, vai phụ, múa tập thể, biểu diễn võ thuật).
  2. Lãnh đạo và quản trị gắn kết (Tác động mạnh nhất):
    • Xây dựng cơ chế tiếp nhận phản hồi hai chiều định kỳ hàng tuần thông qua hệ thống họp giao ban kỹ thuật số.
    • Chuẩn hóa kỹ năng quản trị cảm xúc và lãnh đạo hỗ trợ (Supportive Leadership) cho các trưởng bộ phận nghệ thuật.
  3. Cải tiến môi trường và điều kiện làm việc:
    • Nâng cấp hệ thống ánh sáng, thông gió tại khu vực hậu trường sân khấu số 15 Trương Gia Mô và 49A Lê Lợi.
    • Định kỳ kiểm định an toàn sàn diễn và trang thiết bị biểu diễn võ thuật.

Phân tích chi phí - lợi ích (Cost-Benefit Analysis)

  • Chi phí dự toán: Khoảng 180–220 triệu đồng/năm (cho việc chuẩn hóa đào tạo masterclass, mua sắm bảo hộ biểu diễn và phụ cấp thâm niên).
  • Lợi ích định lượng (ROI):
    • Giảm tỷ lệ nhân sự biểu diễn bỏ việc từ 30% xuống dưới 12%/năm.
    • Tiết kiệm chi phí tuyển dụng và đào tạo diễn viên mới: ước tính 340 triệu đồng/năm.
    • Tỷ suất hoàn vốn đầu tư quản trị (HR-ROI) dự kiến đạt 154,5% sau 18 tháng thực hiện.

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

  • Hạn chế kỹ thuật: Kích thước mẫu khảo sát $N=130$ dù đạt chuẩn thống kê hồi quy OLS nhưng chưa đủ lớn để chạy các mô hình phân tích phân nhóm đa cấp (Multigroup SEM) hoặc phân tích hồi quy phi tuyến.
  • Giới hạn không gian: Dữ liệu chỉ thu thập tại Chi nhánh VKStar Huế, chưa bao gồm các đoàn nghệ thuật biểu diễn liên kết tại các tỉnh lân cận (Đà Nẵng, Quảng Nam).
  • Hướng phát triển tương lai:
    • Tích hợp mô hình học máy (Machine Learning Classification như Random Forest, XGBoost) trên hệ thống phần mềm HRM nội bộ để dự báo sớm nguy cơ nghỉ việc (Churn Prediction) của từng nhân viên dựa trên biến động điểm số khảo sát vi mô (Pulse Survey).
    • Mở rộng phân tích sang mô hình chỉ số trung thành eNPS (Employee Net Promoter Score) kết hợp đánh giá văn hóa tổ chức (Organizational Culture Assessment Instrument - OCAI).

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

       +-----------------------------------------------------------+
       |               CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI                |
       +-----------------------------------------------------------+
          │
          ├─► [Sinh viên]: Tài liệu mẫu nghiên cứu định lượng chuẩn mực
          │   (Cronbach's Alpha, EFA, OLS theo quy chuẩn học thuật)
          │
          ├─► [HR Managers]: Khung thang đo 7 nhân tố định lượng hóa
          │   (Giải pháp giảm tỷ lệ nghỉ việc ngành biểu diễn > 15%)
          │
          ├─► [Doanh nghiệp]: Lộ trình tối ưu hóa chi phí quản trị
          │   (Tăng năng suất lao động, ROI nhân sự đạt 154.5%)
          │
          └─► [Nhà nghiên cứu]: Tập dữ liệu & Hệ số hồi quy thực chứng
              (Nền tảng kiểm chứng các mô hình JDI/Herzberg tại VN)
  • Sinh viên & Học viên cao học: Tiếp cận quy trình nghiên cứu khoa học chuẩn chỉ, phương pháp xử lý dữ liệu định lượng trên SPSS từ khâu thiết kế thang đo đến kiểm định giả định hồi quy.
  • Chuyên viên & Giám đốc Nhân sự (HRM): Nhận được bộ công cụ đo lường mức độ hài lòng có sẵn, áp dụng trực tiếp cho các đơn vị kinh doanh dịch vụ giải trí, sự kiện, nhà hàng - khách sạn.
  • Doanh nghiệp dịch vụ văn hóa: Nâng cao năng lực giữ chân nhân tài, ổn định chất lượng nghệ thuật của các show diễn, gián tiếp gia tăng trải nghiệm du khách.
  • Nhà nghiên cứu xã hội học & kinh tế: Có thêm dữ liệu thực chứng về hành vi tổ chức tại khu vực kinh tế trọng điểm miền Trung giai đoạn hậu COVID-19.

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

1. Yêu cầu kỹ thuật và công cụ cần thiết để tái lập nghiên cứu này là gì?

Cần 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 Python 3.9+ với các thư viện pandas, statsmodels, factor_analyzer, scipy). Dữ liệu thô từ bảng hỏi Likert 5 mức độ cần được mã hóa sạch thành ma trận dữ liệu số nguyên ($1 \le x \le 5$).

2. Giới hạn quy mô mẫu tối thiểu để mô hình đạt giá trị thống kê là bao nhiêu?

Theo công thức của Tabachnick & Fidell (1991), với $m=7$ biến độc lập, quy mô mẫu tối thiểu là $n \ge 8(7) + 50 = 106$. Tuy nhiên, để đảm bảo phân tích EFA ổn định với tỷ lệ 5 quan sát trên 1 biến quan sát (Hair et al., 1998), kích thước mẫu khuyến nghị lý tưởng là từ 140 đến 150 quan sát.

3. Làm thế nào để tích hợp khảo sát này vào hệ thống thông tin nhân sự (HRIS) hiện có?

Doanh nghiệp có thể chuyển đổi bảng hỏi 31 biến quan sát thành form biểu mẫu tự động trên cổng thông tin nội bộ (Intranet) hoặc Google Forms, sau đó thiết lập API gửi dữ liệu JSON về cơ sở dữ liệu phân tích tập trung theo chu kỳ 6 tháng một lần.

4. Chi phí triển khai hệ thống đánh giá định lượng này trong thực tế?

Chi phí ban đầu gần như bằng không nếu tận dụng nhân lực nội bộ và các nền tảng phân tích mã nguồn mở (R/Python). Nếu thuê đơn vị tư vấn độc lập, chi phí khảo sát và báo cáo chuyên sâu dao động từ 25–40 triệu đồng cho mỗi chu kỳ đánh giá.

5. Thứ tự ưu tiên can thiệp khi nguồn lực tài chính của doanh nghiệp bị hạn chế?

Nên ưu tiên can thiệp vào các yếu tố có hệ số hồi quy chuẩn hóa $\beta$ cao nhất nhưng ít đòi hỏi ngân sách tài chính lớn: đầu tiên là Cải tiến phong cách Lãnh đạo và sự công nhận ($\beta = 0.284$), tiếp đến là Minh bạch hóa chính sách phân công Bản chất công việc ($\beta = 0.246$), sau đó mới đến các chính sách điều chỉnh Tiền lương và Phúc lợi.

Kết luận

Đồ án đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, thiết lập một khung đánh giá khoa học và thực chứng về các yếu tố chi phối sự hài lòng của người lao động tại Công ty Cổ phần VKStar Huế. Bằng việc áp dụng các kỹ thuật thống kê đa biến nâng cao (Cronbach's Alpha $> 0.77$, EFA trích $68.45%$ phương sai, OLS $R^2_{\text{adj}} = 0.642$), nghiên cứu đã định lượng hóa thành công tầm quan trọng của 7 nhóm yếu tố cốt lõi. Kết quả này không chỉ giải quyết bài toán cấp bách về ổn định nhân sự nghệ thuật hậu COVID-19 cho VKStar mà còn cung cấp tài liệu tham khảo giá trị cho công tác quản trị nguồn nhân lực trong ngành công nghiệp văn hóa - du lịch Việt Nam. Các nhà quản lý được khuyến nghị sớm áp dụng lộ trình 5 nhóm giải pháp đồng bộ để xây dựng môi trường làm việc hạnh phúc, bền vững và hiệu suất cao.