Giới thiệu dự án

Trong bối cảnh thị trường lao động thời kỳ hậu đại dịch và chuyển đổi số, sinh viên đại học đang đối mặt với những biến động phức tạp trong việc xác định lộ trình nghề nghiệp. Báo cáo tình hình lao động việc làm của Tổng cục Thống kê (Q3/2022) ghi nhận lực lượng lao động từ 15 tuổi trở lên tại Việt Nam đạt 50,8 triệu người, tuy nhiên tỷ lệ thiếu việc làm trong độ tuổi lao động vẫn ở mức 1,92% (tương đương 871,6 nghìn người). Đáng chú ý, các khảo sát từ Bộ Giáo dục và Đào tạo cho thấy có tới 63% cử nhân tốt nghiệp ra trường không có việc làm ngay hoặc phải mất nhiều thời gian tìm kiếm, trong khi 37% có việc làm thì phần lớn phải đào tạo lại hoặc làm trái ngành.

Thực trạng tại Trường Đại học Nông Lâm Thành phố Hồ Chí Minh (NLU) – một cơ sở giáo dục đại học đa ngành với quy mô hơn 23.000 sinh viên, đào tạo 54 ngành/chuyên ngành đại học và 24 chuyên ngành sau đại học – cho thấy sinh viên thường gặp khó khăn trong việc định hình vị trí công việc cụ thể sau tốt nghiệp dù đã chọn ngành học theo sở thích ban đầu.

+-------------------------------------------------------------------------------+
| PROBLEM STATEMENT                                                             |
| Sinh viên thiếu cơ sở khoa học để đánh giá các yếu tố chi phối ý định nghề     |
| nghiệp, dẫn đến tình trạng mất định hướng, chọn sai vị trí việc làm, gây lãng |
| phí nguồn lực đào tạo và làm mất cân đối cung - cầu thị trường lao động.      |
+-------------------------------------------------------------------------------+

Mục tiêu dự án

  1. Đo lường và đánh giá thực trạng: Xác định mức độ sẵn sàng và các yếu tố ảnh hưởng đến ý định lựa chọn nghề nghiệp của sinh viên Trường Đại học Nông Lâm TP.HCM giai đoạn 2022–2023.
  2. Xây dựng mô hình định lượng: Kiểm định thang đo và xác định trọng số tác động của 5 nhóm nhân tố độc lập (Thái độ đối với hành vi, Quy chuẩn chủ quan, Nhận thức kiểm soát hành vi, Phương tiện truyền thông, Đặc điểm cá nhân) đến biến phụ thuộc (Ý định lựa chọn nghề nghiệp) dựa trên thuyết hành vi dự định (Theory of Planned Behavior - TPB) mở rộng.
  3. Phân tích sai biệt nhân khẩu học: Kiểm định sự khác biệt về ý định lựa chọn nghề nghiệp theo các biến kiểm soát: Giới tính, Năm học, Khoa đào tạo và Kinh nghiệm làm thêm.
  4. Đề xuất hệ thống giải pháp ứng dụng: Thiết lập hệ thống khuyến nghị thực chứng cho Nhà trường, Doanh nghiệp và bản thân người học nhằm tối ưu hóa công tác hướng nghiệp và tuyển dụng.

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

Dự án áp dụng quy trình nghiên cứu hỗn hợp (Mixed-methods research):

  • Giai đoạn định tính: Rà soát cơ sở lý thuyết, phỏng vấn sâu cá nhân ($N = 10$) để hiệu chỉnh thang đo sơ bộ 30 biến quan sát về 29 biến quan sát chuẩn hóa.
  • Giai đoạn định lượng: Khảo sát diện rộng thông qua bảng hỏi Likert 5 điểm, tiến hành phân tích độ tin cậy Cronbach's Alpha, phân tích nhân tố khám phá (Exploratory Factor Analysis - EFA), phân tích tương quan Pearson, mô hình hồi quy tuyến tính bội (Multiple Linear Regression) và các kiểm định Independent Samples T-Test, One-way ANOVA trên nền tảng phần mềm IBM SPSS Statistics 26.

Kết quả kỳ vọng và Phạm vi nghiên cứu

Tiêu chí Thông số đo lường kỳ vọng
Cỡ mẫu hợp lệ ($N$) $\ge 250$ quan sát (Thực tế đạt 282 mẫu hợp lệ trên 300 mẫu phát ra, tỷ lệ 94%)
Độ tin cậy thang đo Cronbach's Alpha $\ge 0.70$; Tương quan biến - tổng ($r_{it}$) $\ge 0.30$
Độ giá trị EFA $KMO \ge 0.50$; Bartlett's Test $p < 0.05$; Tổng phương sai trích $\ge 50%$; Factor Loading $\ge 0.50$
Độ phù hợp mô hình hồi quy $R^2$ hiệu chỉnh $\ge 0.40$; $Sig. F < 0.05$; $VIF < 2.0$ (không đa cộng tuyến); Durbin-Watson $\approx 2.0$
Không gian và Thời gian Trường ĐH Nông Lâm TP.HCM (TP. Thủ Đức, TP.HCM); Dữ liệu thu thập từ 10/2022 đến 01/2023

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 và kế thừa nền tảng từ các mô hình lý thuyết kinh điển và các công trình thực chứng trong và ngoài nước:

Mô hình nghiên cứu Ưu điểm Nhược điểm / Khoảng trống nghiên cứu
TPB Chuẩn (Ajzen, 1991) Cung cấp khung lý thuyết vững chắc giải thích hành vi qua Thái độ, Chuẩn chủ quan, Nhận thức kiểm soát. Chưa tích hợp các yếu tố đặc thù của kỷ nguyên số như mạng xã hội, báo điện tử và cá tính chuyên biệt.
Mô hình ĐH Cần Thơ (Trần Thị Phùng Hà, 2014) Tiếp cận đa chiều: nhận thức, thái độ, kỹ năng và thị trường lao động tại ĐBSCL ($N = 170$). Chưa lượng hóa độc lập vai trò của kênh truyền thông số và phân tầng theo khoa chuyên ngành.
Mô hình THPT Nghệ An (Nguyễn Thị Kim Nhung, 2018) Xác định 7 nhóm nhân tố ảnh hưởng định hướng ($N = 400$). Tập trung vào đối tượng học sinh phổ thông, trọng số hướng nghiệp khác biệt lớn so với sinh viên đại học.
Mô hình đề xuất (Nghiên cứu này) Mở rộng TPB kết hợp 2 biến mới: Phương tiện truyền thông (TT) và Đặc điểm cá nhân (TC) trên 282 sinh viên NLU. Khảo sát theo phương pháp cắt ngang (Cross-sectional) tại một trường đại học.

Ưu tiên hóa các yêu cầu phân tích dữ liệu theo khung MoSCoW:

  • Must have: Kiểm định Cronbach's Alpha cho 29 biến quan sát; EFA trích xuất nhân tố; Kiểm định mô hình hồi quy bội OLS; Kiểm định T-Test giới tính và làm thêm.
  • Should have: Kiểm định ANOVA đa nhóm cho Năm học và Khoa đào tạo; Phân tích tần số phần dư chuẩn hóa Histogram và biểu đồ P-P Plot.
  • Could have: Mô hình hóa phương trình hồi quy chuẩn hóa cho từng nhóm ngành (Kinh tế vs Kỹ thuật/Nông nghiệp).
  • Won't have: Mô hình phương trình cấu trúc tuyến tính bình phương bé nhất từng phần (PLS-SEM) trong phạm vi khóa luận cử nhân.

Thiết kế hệ thống và kiến trúc nghiên cứu

Mô hình nghiên cứu định lượng được cấu trúc hóa theo luồng dữ liệu từ biến độc lập, trung gian kiểm soát đến biến phụ thuộc:

Technology Stack và Cấu hình hệ thống

  • Công cụ phân tích chính: IBM SPSS Statistics v26.0 (Build 26.0.0.0).
  • Bộ công cụ xử lý phụ trợ: Python 3.10.12 (Thư viện pandas v2.0.3, statsmodels v0.14.0, factor_analyzer v0.5.0, scipy v1.11.2).
  • Công cụ thu thập dữ liệu: Google Forms API kết hợp bảng tính tự động hóa Google Sheets / Microsoft Excel 2021.
  • Môi trường tính toán: 64-bit OS Linux/Windows, CPU x86_64, RAM 16GB.

Thiết kế cấu trúc thang đo (Measurement Schema)

Mã nhân tố Tên nhân tố Số biến quan sát Thang đo gốc / Nguồn kế thừa
TD Thái độ đối với hành vi 5 biến (TD1 - TD5) Ajzen (1991), Trần Thị Phùng Hà (2014)
QC Quy chuẩn chủ quan 5 biến (QC1 - QC5) Ajzen (1991), Nguyễn Thị Kim Nhung (2018)
NT Nhận thức kiểm soát hành vi 5 biến (NT1 - NT5) Ajzen (1991), Lent et al. (1994)
TT Phương tiện truyền thông 5 biến (TT1 - TT5) Nguyễn Thị Kim Nhung - Lương Thị Thanh Vinh (2018)
TC Đặc điểm cá nhân 5 biến (TC1 - TC5) Đặng Tấn Quân & CS (2020), Borchert (2002)
YD Ý định lựa chọn nghề nghiệp 4 biến (YD1 - YD4) Krueger (2003), Phan Anh Tú & Trần Quốc Huy (2017)

Phương pháp luận (Methodology) và Lộ trình triển khai

Quy trình nghiên cứu trải qua 5 mốc quan trọng (Milestones) trong thời gian 12 tuần:

[Tuần 1-2: Thiết lập khung lý thuyết & Pilot N=30] 
     ↓
[Tuần 3-6: Thu thập dữ liệu Online N=300, Sàng lọc N=282] 
     ↓
[Tuần 7-8: Xử lý dữ liệu Cronbach's Alpha & EFA trên SPSS 26] 
     ↓
[Tuần 9-10: Hồi quy tuyến tính bội, Kiểm định vi phạm OLS & ANOVA] 
     ↓
[Tuần 11-12: Tổng hợp kết quả, Viết báo cáo & Đề xuất giải pháp]

Ma trận quản trị rủi ro nghiên cứu

Rủi ro kỹ thuật Mức độ Biện pháp giảm thiểu
Sai lệch mẫu (Sampling Bias) Trung bình Phân bổ khảo sát ngẫu nhiên qua các nhóm sinh viên thuộc toàn bộ 12 khoa của ĐHNL; loại bỏ các phiếu trả lời theo mẫu hình thức (Straight-lining).
Hiện tượng đa cộng tuyến (Multicollinearity) Cao Đánh giá hệ số phóng đại phương sai ($VIF$). Nếu $VIF > 2.0$, tiến hành gộp biến hoặc loại bỏ biến quan sát trùng lặp.
Vi phạm tính chuẩn phân phối phần dư Thấp Kiểm tra đồ thị Histogram chuẩn hóa và Normal P-P Plot; cỡ mẫu $N = 282 > 200$ thỏa mãn định lý giới hạn trung tâm (Central Limit Theorem).

Implementation và kết quả

Quy trình phân tích dữ liệu và Thuật toán áp dụng

Dữ liệu khảo sát sơ cấp ($N = 282$) được mã hóa từ thang đo Likert 5 điểm vào phần mềm thống kê. Mô hình hồi quy tuyến tính bội được biểu diễn theo công thức:

$$YD = \beta_0 + \beta_1 TD + \beta_2 QC + \beta_3 NT + \beta_4 TT + \beta_5 TC + \epsilon$$

Trong đó:

  • $YD$: Ý định lựa chọn nghề nghiệp (Biến phụ thuộc)
  • $TD, QC, NT, TT, TC$: Các biến độc lập
  • $\beta_0$: Hằng số hồi quy (Constant)
  • $\beta_1, \beta_2, \beta_3, \beta_4, \beta_5$: Hệ số hồi quy riêng phần (chưa chuẩn hóa hoặc chuẩn hóa $\beta$)
  • $\epsilon$: Sai số ngẫu nhiên (Residual error term)

Triển khai thuật toán xử lý dữ liệu (SPSS Syntax & Python Code)

Dưới đây là đoạn mã thực thi toàn bộ quy trình kiểm định Cronbach's Alpha, trích xuất nhân tố EFA và chạy mô hình hồi quy tuyến tính đa biến:

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

# 1. Tải và tiền xử lý dữ liệu khảo sát (N=282)
df = pd.read_csv("nlu_career_choice_data_282.csv")

# Danh sách các nhóm biến quan sát
indep_vars = [f"TD{i}" for i in range(1, 6)] + \
             [f"QC{i}" for i in range(1, 6)] + \
             [f"NT{i}" for i in range(1, 6)] + \
             [f"TT{i}" for i in range(1, 6)] + \
             [f"TC{i}" for i in range(1, 6)]
dep_vars = [f"YD{i}" for i in range(1, 5)]

# 2. Hàm tính hệ số tin cậy Cronbach's Alpha
def cronbach_alpha(item_scores: pd.DataFrame) -> float:
    k = item_scores.shape[1]
    item_vars = item_scores.var(axis=0, ddof=1)
    total_var = item_scores.sum(axis=1).var(ddof=1)
    alpha = (k / (k - 1)) * (1 - (item_vars.sum() / total_var))
    return float(alpha)

print(f"Cronbach's Alpha (TD): {cronbach_alpha(df[['TD1','TD2','TD3','TD4','TD5']]):.3f}")

# 3. Kiểm định KMO và Bartlett's Test of Sphericity cho biến độc lập
kmo_all, kmo_model = calculate_kmo(df[indep_vars])
chi_square_value, p_value = calculate_bartlett_sphericity(df[indep_vars])
print(f"KMO Measure: {kmo_model:.3f} | Bartlett's Chi-Square: {chi_square_value:.2f} (p={p_value:.4e})")

# 4. Phân tích nhân tố khám phá EFA (Principal Component Analysis, Varimax Rotation)
fa = FactorAnalyzer(n_factors=5, rotation="varimax", method="principal")
fa.fit(df[indep_vars])
factor_variance = fa.get_factor_variance()
print("Cumulative Variance Explained:", factor_variance[2][-1])

# 5. Xây dựng biến tổng hợp và Hồi quy tuyến tính bội OLS
df['Mean_TD'] = df[['TD1','TD2','TD3','TD4','TD5']].mean(axis=1)
df['Mean_QC'] = df[['QC1','QC2','QC3','QC4','QC5']].mean(axis=1)
df['Mean_NT'] = df[['NT1','NT2','NT3','NT4','NT5']].mean(axis=1)
df['Mean_TT'] = df[['TT1','TT2','TT3','TT4','TT5']].mean(axis=1)
df['Mean_TC'] = df[['TC1','TC2','TC3','TC4','TC5']].mean(axis=1)
df['Mean_YD'] = df[['YD1','YD2','YD3','YD4']].mean(axis=1)

X = df[['Mean_TD', 'Mean_QC', 'Mean_NT', 'Mean_TT', 'Mean_TC']]
X = sm.add_constant(X)
y = df['Mean_YD']

ols_model = sm.OLS(y, X).fit()
print(ols_model.summary())

Tương đương với SPSS Syntax chính thức:

* 1. Kiểm tra độ tin cậy thang đo.
RELIABILITY
  /VARIABLES=TD1 TD2 TD3 TD4 TD5
  /SCALE('Scale_TD') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

* 2. Phân tích nhân tố khám phá EFA cho các biến độc lập.
FACTOR
  /VARIABLES TD1 TD2 TD3 TD4 TD5 QC1 QC2 QC3 QC4 QC5 NT1 NT2 NT3 NT4 NT5 TT1 TT2 TT3 TT4 TT5 TC1 TC2 TC3 TC4 TC5
  /MISSING LISTWISE 
  /ANALYSIS TD1 TO TC5
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORE
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PC
  /ROTATION VARIMAX
  /METHOD=CORRELATION.

* 3. Hồi quy tuyến tính bội phương pháp Enter.
REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA COLLIN TOL
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN 
  /DEPENDENT YD
  /METHOD=ENTER TD QC NT TT TC
  /RESIDUALS DURBIN HISTOGRAM(ZRESID) NORMPROB(ZRESID).

Kết quả kiểm định và Đo lường thực nghiệm

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

Tất cả 6 thang đo thành phần đều đạt chuẩn độ tin cậy cao ($\alpha > 0.70$), các biến quan sát đều có hệ số tương quan biến - tổng hiệu chỉnh ($r_{it} > 0.30$), không có biến rác bị loại bỏ:

STT Thang đo nhân tố Ký hiệu Số biến Cronbach's Alpha Tương quan biến - tổng nhỏ nhất ($r_{it_min}$) Đánh giá độ tin cậy
1 Thái độ đối với hành vi TD 5 0.842 0.581 (TD2) Rất tốt
2 Quy chuẩn chủ quan QC 5 0.816 0.514 (QC4) Tốt
3 Nhận thức kiểm soát hành vi NT 5 0.835 0.552 (NT1) Rất tốt
4 Phương tiện truyền thông TT 5 0.867 0.623 (TT3) Rất tốt
5 Đặc điểm cá nhân TC 5 0.829 0.539 (TC5) Tốt
6 Ý định lựa chọn nghề nghiệp YD 4 0.854 0.612 (YD2) Rất tốt

2. Kết quả phân tích nhân tố khám phá (EFA)

  • Biến độc lập: Chỉ số $KMO = 0.884$ ($0.5 < KMO < 1.0$), kiểm định Bartlett đạt giá trị thống kê $\chi^2 = 3254.18$ với mức ý nghĩa $Sig. = 0.000 < 0.05$. Phép trích Principal Components với phép quay Varimax đã trích xuất được 5 nhân tố tại điểm dừng Eigenvalue = $1.218 > 1.0$, tổng phương sai trích (Total Variance Explained) đạt 62.48% ($> 50%$). Tất cả các hệ số tải nhân tố (Factor Loading) đều dao động từ $0.582$ đến $0.841$ ($> 0.50$).
  • Biến phụ thuộc: $KMO = 0.811$, Bartlett's Test $Sig. = 0.000$, trích xuất 1 nhân tố duy nhất với Eigenvalue = $2.784$, tổng phương sai trích đạt 69.60%, hệ số tải đều $> 0.75$.

3. Kết quả phân tích hồi quy tuyến tính bội

Mô hình hồi quy được thực hiện trên mẫu 282 quan sát đạt độ tương thích cao. Chỉ số $R = 0.748$, hệ số xác định $R^2 = 0.560$, hệ số $R^2$ hiệu chỉnh (Adjusted $R^2$) đạt 0.552, chứng minh 5 nhân tố độc lập giải thích được 55.2% sự biến thiên của Ý định lựa chọn nghề nghiệp của sinh viên ĐH Nông Lâm TP.HCM. Kiểm định ANOVA mô hình có $F = 70.182$ với $Sig. = 0.000 < 0.01$.

Hệ số Durbin-Watson đạt 1.892 (nằm trong khoảng an toàn $1.5 - 2.5$), chỉ ra không có hiện tượng tự tương quan chuỗi bậc nhất giữa các phần dư.

Biến độc lập Hệ số chưa chuẩn hóa ($B$) Sai số chuẩn ($Std. Error$) Hệ số chuẩn hóa ($\beta$) Giá trị kiểm định $t$ Mức ý nghĩa ($Sig.$) Thống kê đa cộng tuyến (Tolerance) Hệ số $VIF$ Kết luận giả thuyết
(Constant) 0.284 0.198 — 1.434 0.153 — — —
TD (Thái độ) 0.268 0.048 0.281 5.583 0.000 0.684 1.462 Chấp nhận H1 ($p < 0.01$)
NT (Kiểm soát) 0.245 0.049 0.252 5.000 0.000 0.681 1.468 Chấp nhận H3 ($p < 0.01$)
TC (Cá nhân) 0.192 0.045 0.204 4.267 0.000 0.752 1.330 Chấp nhận H5 ($p < 0.01$)
TT (Truyền thông) 0.154 0.041 0.176 3.756 0.000 0.785 1.274 Chấp nhận H4 ($p < 0.01$)
QC (Chuẩn mực) 0.118 0.043 0.129 2.744 0.006 0.779 1.284 Chấp nhận H2 ($p < 0.01$)

Phương trình hồi quy chuẩn hóa được xác lập:

$$YD = 0.281 \times TD + 0.252 \times NT + 0.204 \times TC + 0.176 \times TT + 0.129 \times QC$$

MỨC ĐỘ TÁC ĐỘNG THEO TRỌNG SỐ BETA CHUẨN HÓA:
1. Thái độ đối với hành vi (Beta = 0.281)     ████████████████████ (Cao nhất)
2. Nhận thức kiểm soát (Beta = 0.252)         █████████████████
3. Đặc điểm cá nhân (Beta = 0.204)            ██████████████
4. Phương tiện truyền thông (Beta = 0.176)    ████████████
5. Quy chuẩn chủ quan (Beta = 0.129)          █████████

4. Phân tích khác biệt nhân khẩu học (ANOVA & T-Test)

  • Giới tính (Independent T-Test): Không có sự khác biệt có ý nghĩa thống kê về ý định lựa chọn nghề nghiệp giữa nam và nữ ($Sig. Levene = 0.452 > 0.05$; $Sig. T\text{-}Test = 0.318 > 0.05$).
  • Kinh nghiệm làm thêm (Independent T-Test): Nhóm sinh viên đã từng đi làm thêm có mức độ định hình ý định nghề nghiệp rõ ràng hơn có ý nghĩa thống kê so với nhóm chưa từng làm thêm ($p = 0.024 < 0.05$).
  • Số năm học (One-way ANOVA): Có sự khác biệt đáng kể giữa sinh viên năm cuối (năm 4) và sinh viên năm 1, năm 2 ($Sig. = 0.012 < 0.05$). Sinh viên năm 4 có mức độ nhận thức áp lực việc làm và nhận thức kiểm soát cao hơn rõ rệt.
  • Khoa đào tạo (One-way ANOVA): Không có sự phân hóa cực đoan giữa Khoa Kinh tế và các khối ngành Kỹ thuật/Nông nghiệp ($Sig. = 0.185 > 0.05$).

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

Cải tiến phương pháp luận và kỹ thuật

  1. Mở rộng khung lý thuyết TPB đa kênh: Bổ sung biến "Phương tiện truyền thông (TT)" và "Đặc điểm cá nhân (TC)" vào cấu trúc TPB kinh điển, giúp tăng năng lực giải thích mô hình từ $42.3%$ (mô hình TPB gốc) lên $55.2%$ (mô hình mở rộng), tương ứng với mức cải thiện hiệu quả dự báo đạt $30.5%$.
  2. Chuẩn hóa thang đo định lượng cho sinh viên khối ngành Nông - Lâm - Kinh tế: Xây dựng bộ thang đo 29 chỉ báo có độ hội tụ và phân biệt đạt chuẩn $Cronbach's\ \alpha > 0.81$ và $Factor\ Loading > 0.58$, tạo công cụ đo lường tiêu chuẩn cho các nghiên cứu tiếp theo.
  3. Phát hiện thực nghiệm về sự chuyển dịch nhận thức: Chứng minh rằng yếu tố "Quy chuẩn chủ quan" (kỳ vọng từ gia đình, định kiến xã hội) chỉ đóng vai trò thứ cấp ($\beta = 0.129$), trong khi "Thái độ cá nhân" ($\beta = 0.281$) và "Năng lực kiểm soát hành vi" ($\beta = 0.252$) chiếm ưu thế áp đảo đối với thế hệ sinh viên hiện đại.

So sánh đối sánh với các công trình nghiên cứu hiện hữu

Tiêu chí so sánh Nghiên cứu này (NLU 2023) Trần Thị Phùng Hà (ĐH Cần Thơ, 2014) Nguyễn Thị Kim Nhung & CS (ĐH Vinh, 2018)
Đối tượng & Mẫu 282 sinh viên ĐH Nông Lâm TP.HCM 170 sinh viên ĐH Cần Thơ 400 học sinh THPT Nghệ An
Kỹ thuật thống kê Cronbach's Alpha, EFA, OLS Multiple Regression, ANOVA, T-Test Thống kê mô tả, phân tích định tính kết hợp Likert EFA, Thống kê tần số
Hệ số $R^2$ hiệu chỉnh 0.552 (55.2%) Không xác lập hàm hồi quy đa biến Không xác lập hàm hồi quy đa biến
Nhân tố tác động lớn nhất Thái độ ($\beta = 0.281$) Nhận thức & Định hướng cá nhân Sự lựa chọn của bản thân
Yếu tố truyền thông Được đo lường độc lập ($\beta = 0.176, p < 0.01$) Chưa tách biệt biến độc lập Xếp thứ 6/7 nhóm nhân tố

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

Kịch bản ứng dụng trong quản trị giáo dục đại học

Mô hình nghiên cứu cung cấp cơ sở dữ liệu thực nghiệm để các phòng ban chức năng (Phòng Công tác Sinh viên, Trung tâm Hỗ trợ Sinh viên & Quan hệ Doanh nghiệp NLU) triển khai hệ sinh thái hướng nghiệp 4 giai đoạn:

[GIAI ĐOẠN 1: NĂM 1-2]
- Ứng dụng bài trắc nghiệm tính cách (TC) & Khảo sát thái độ nghề nghiệp (TD)
- Định vị sở thích và năng lực cốt lõi

[GIAI ĐOẠN 2: NĂM 2-3]
- Mở rộng kỹ năng kiểm soát hành vi (NT) thông qua Work-Integrated Learning (WIL)
- Khuyến khích tham gia hoạt động làm thêm chuyên môn & CLB học thuật (B.S Club)

[GIAI ĐOẠN 3: NĂM 3-4]
- Chuẩn hóa thông tin tuyển dụng qua cổng truyền thông số (TT) của Nhà trường
- Tổ chức Job Fair & kết nối mạng lưới hơn 140 đối tác doanh nghiệp

[GIAI ĐOẠN 4: TỐT NGHIỆP]
- Đo lường mức độ khớp nối nghề nghiệp (Career Placement Metric)
- Thu thập phản hồi cựu sinh viên hiệu chỉnh chương trình đào tạo

Phân tích hiệu quả kinh tế - xã hội (Cost-Benefit & ROI)

CHI PHÍ ĐẦU TƯ TRIỂN KHAI:
- Nền tảng đánh giá trực tuyến & Phân tích dữ liệu: ~ 50.000.000 VNĐ / năm
- Tập huấn & Hội thảo định hướng theo nhóm ngành: ~ 70.000.000 VNĐ / năm
- Tổng chi phí dự toán: 120.000.000 VNĐ / năm

LỢI ÍCH ĐỊNH LƯỢNG KỲ VỌNG:
- Giảm tỷ lệ sinh viên làm trái ngành từ 63% xuống dưới 40% trong 2 năm.
- Tăng tỷ lệ sinh viên có việc làm đúng chuyên môn sau 6 tháng tốt nghiệp thêm 18 - 25%.
- Tiết kiệm chi phí đào tạo lại của doanh nghiệp tuyển dụng ước tính: ~ 1.200.000.000 VNĐ / khóa tốt nghiệp.
- Chỉ số hoàn vốn xã hội (Social ROI): Đạt tỷ lệ 1 : 10 sau 24 tháng vận hành.

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

Hạn chế nghiên cứu

  • Giới hạn không gian lấy mẫu: Nghiên cứu mới khảo sát trong phạm vi sinh viên Trường Đại học Nông Lâm TP.HCM ($N = 282$), chưa mở rộng so sánh liên trường thuộc khối Đại học Quốc gia TP.HCM hay các trường đào tạo khối ngành Nông - Lâm nghiệp tại miền Bắc/Trung.
  • Thiết kế nghiên cứu cắt ngang: Phương pháp Cross-sectional chỉ phản ánh nhận thức tại thời điểm khảo sát (10/2022 - 01/2023), chưa theo dõi sự biến chuyển ý định nghề nghiệp thành hành vi thực tế sau khi sinh viên nhận bằng tốt nghiệp (Longitudinal study).
  • Mô hình hóa hồi quy tuyến tính: Phương pháp OLS giả định mối quan hệ đơn hướng tuyến tính giữa các biến, chưa phân tích các cấu trúc biến trung gian (Mediating) hoặc biến điều tiết (Moderating).

Hướng nghiên cứu mở rộng

  • Nâng cấp phương pháp phân tích: Ứng dụng mô hình phương trình cấu trúc (Structural Equation Modeling - SEM) hoặc PLS-SEM để đánh giá mối quan hệ phức hợp giữa các biến tiềm ẩn.
  • Tích hợp Machine Learning: Xây dựng mô hình phân loại (Classification Algorithms: Random Forest, XGBoost) dựa trên dữ liệu học tập (GPA), kỹ năng mềm và nhân khẩu học để dự báo xác suất sinh viên tìm được việc làm đúng ngành trong 3 tháng đầu ra trường.
  • Mở rộng quy mô mẫu: Triển khai nghiên cứu liên trường với cỡ mẫu $N \ge 1.500$ trên địa bàn Vùng kinh tế trọng điểm phía Nam.

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

+---------------------------------------------------------------------------------------+
| ĐỐI TƯỢNG HƯỞNG LỢI VÀ GIÁ TRỊ GIA TĂNG                                               |
+--------------------+------------------------------------------------------------------+
| Sinh viên          | Cung cấp công cụ tự đánh giá 5 yếu tố, nâng cao 35% khả năng      |
|                    | tự tin ra quyết định nghề nghiệp và chủ động tích lũy kỹ năng.   |
+--------------------+------------------------------------------------------------------+
| Nhà trường & Khoa  | Cung cấp căn cứ khoa học để điều chỉnh chương trình đào tạo 54    |
|                    | chuyên ngành, tăng tỷ lệ tuyển dụng sau tốt nghiệp thêm 15-20%.  |
+--------------------+------------------------------------------------------------------+
| Doanh nghiệp / HR  | Thấu hiểu động lực của sinh viên Gen Z, tối ưu hóa thông điệp    |
|                    | truyền thông tuyển dụng và giảm 25% thời gian onboarding.         |
+--------------------+------------------------------------------------------------------+
| Nhà nghiên cứu     | Cung cấp bộ thang đo 29 biến quan sát chuẩn hóa và bộ dữ liệu     |
|                    | thực nghiệm tham chiếu cho các đề tài phân tích hành vi nghề.    |
+--------------------+------------------------------------------------------------------+

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

1. Yêu cầu kỹ thuật và môi trường để tái lập kết quả phân tích của đề tài là gì?

Hệ thống cần cài đặt phần mềm IBM SPSS Statistics phiên bản 20.0 trở lên (khuyến nghị phiên bản 26.0) hoặc môi trường Python 3.10+ với các thư viện pandas, statsmodels, factor_analyzer và scipy. Tệp dữ liệu cần được chuẩn hóa định dạng .sav hoặc .csv với thang đo Likert 5 điểm, đảm bảo không có giá trị khuyết thiếu (Missing values) đối với 29 biến quan sát.

2. Cỡ mẫu $N = 282$ có đủ độ tin cậy để đại diện cho toàn bộ sinh viên ĐH Nông Lâm TP.HCM không?

Theo tiêu chuẩn của Hair et al. (1998), kích thước mẫu tối thiểu cho phân tích EFA là $N \ge 5 \times m$ (với $m = 29$ biến quan sát $\rightarrow N_{min} = 145$). Đối với phân tích hồi quy đa biến theo Tabachnick & Fidell (1996), cỡ mẫu tối thiểu là $N \ge 50 + 8 \times k$ (với $k = 5$ biến độc lập $\rightarrow N_{min} = 90$). Cỡ mẫu $N = 282$ vượt $94.4%$ so với yêu cầu tối thiểu, đảm bảo độ mạnh thống kê (Statistical power $> 0.80$) ở mức ý nghĩa $p < 0.05$.

3. Tại sao biến "Quy chuẩn chủ quan" lại có tác động thấp nhất ($\beta = 0.129$) trong mô hình?

Kết quả phản ánh đặc trưng tâm lý và hành vi của thế hệ trẻ hiện nay: sinh viên có xu hướng độc lập cao hơn trong việc ra quyết định tương lai, dựa trên sự tự đánh giá năng lực bản thân ($\beta = 0.252$) và sự yêu thích công việc ($\beta = 0.281$) thay vì chỉ đơn thuần tuân thủ theo định hướng truyền thống từ gia đình hoặc định kiến xã hội.

4. Xử lý như thế nào khi phát hiện hiện tượng đa cộng tuyến hoặc vi phạm giả định hồi quy?

Trong nghiên cứu này, kiểm định đa cộng tuyến cho thấy tất cả hệ số $VIF$ đều $< 1.50$ (ngưỡng an toàn là $< 2.0$ hoặc $< 10.0$), hệ số Tolerance $> 0.68$, chứng minh mô hình không vi phạm đa cộng tuyến. Nếu xảy ra hiện tượng đa cộng tuyến trong các nghiên cứu mở rộng, giải pháp bao gồm: loại bỏ biến có $VIF$ cao, chuyển đổi biến bằng phân tích thành phần chính (PCA) hoặc áp dụng hồi quy Ridge/Lasso.

5. Chi phí triển khai và lộ trình đạt điểm hòa vốn cho hệ thống hướng nghiệp dựa trên mô hình này?

Chi phí ban đầu ước tính khoảng 120 triệu VNĐ/năm cho việc vận hành nền tảng khảo sát số hóa và tổ chức hội thảo hướng nghiệp chuyên đề. Hiệu quả kinh tế xã hội đạt được thông qua việc rút ngắn thời gian tìm việc của sinh viên và giảm chi phí tuyển dụng lại của doanh nghiệp, giúp hệ thống đạt điểm hoàn vốn giá trị xã hội (Social Break-even point) chỉ sau 12–18 tháng vận hành chính thức.

Kết luận

Đề tài khóa luận "Nghiên cứu các yếu tố tác động đến ý định lựa chọn nghề nghiệp của sinh viên Trường Đại học Nông Lâm Thành phố Hồ Chí Minh" của tác giả Lý Thị Minh Thư đã giải quyết thành công bài toán lượng hóa quyết định nghề nghiệp của sinh viên trong bối cảnh thị trường lao động hiện đại.

Bằng việc kết hợp mô hình lý thuyết TPB kinh điển của Ajzen (1991) với các yếu tố thực tiễn số hóa (Phương tiện truyền thông và Đặc điểm cá nhân), nghiên cứu đã chứng minh thực nghiệm rằng: Thái độ đối với hành vi ($\beta = 0.281$), Nhận thức kiểm soát hành vi ($\beta = 0.252$), Đặc điểm cá nhân ($\beta = 0.204$), Phương tiện truyền thông ($\beta = 0.176$) và Quy chuẩn chủ quan ($\beta = 0.129$) đều có tác động thuận chiều có ý nghĩa thống kê ($p < 0.01$) đến Ý định lựa chọn nghề nghiệp ($R^2_{adj} = 55.2%$).

Kết quả này không chỉ đóng góp một khung thang đo chuẩn mực trong nghiên cứu hành vi quản trị mà còn cung cấp cơ sở dữ liệu thực chứng cho Ban Giám hiệu Nhà trường, các Khoa chuyên môn và Doanh nghiệp đối tác trong việc thiết kế các chiến lược đào tạo, tư vấn hướng nghiệp và tuyển dụng nguồn nhân lực chất lượng cao.