Giới thiệu dự án

Ngành du lịch Việt Nam – "ngành công nghiệp không khói" – giữ vai trò chiến lược trong cơ cấu kinh tế quốc dân. Theo số liệu từ Tổng cục Du lịch năm 2019, Việt Nam xác lập kỷ lục đón 18 triệu lượt khách quốc tế (tăng trưởng 16,2% so với năm 2018) và phục vụ 85 triệu lượt khách nội địa (tăng 6%), tổng doanh thu toàn ngành vượt mốc 720.000 tỷ VNĐ. Riêng tại trung tâm du lịch miền Trung là TP. Đà Nẵng, cơ sở lưu trú năm 2019 phục vụ hơn 5,222 triệu lượt khách (tăng 26% so với 2018), trong đó khách quốc tế đạt 3,561 triệu lượt. Mặc dù tốc độ tăng trưởng nhanh, năng lực cạnh tranh du lịch Việt Nam (xếp hạng 63/140 theo Diễn đàn Kinh tế Thế giới - WEF) đối mặt với rào cản lớn: tỷ lệ khách quay lại thấp (chỉ dao động 10% – 40%), mức chi tiêu trung bình còn khiêm tốn (~1.000 USD/chuyến 9 ngày) do chất lượng dịch vụ chưa đồng đều và thiếu hệ thống đo lường trải nghiệm khách hàng chuẩn hóa.

+-----------------------------------------------------------------------------+
|                                INDUSTRY CONTEXT                             |
|  - VN 2019: 18M Khách quốc tế (+16.2%) | 85M Khách nội địa | 720k tỷ VNĐ    |
|  - Đà Nẵng 2019: 5.22M Khách lưu trú (+26%) | Khách quay lại thấp: 10-40%   |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|                          DOANH NGHIỆP NGHIÊN CỨU                            |
|  - Công ty TNHH MTV Thương mại Du lịch Trường Sa (Trường Sa Tourist)        |
|  - Thành lập: 17/07/2017 | Quy mô: ~7.000 khách/năm                         |
|  - Doanh thu 2019: Tăng 2,6 tỷ VNĐ so với 2017                              |
+-----------------------------------------------------------------------------+
+-----------------------------------------------------------------------------+
|                            VẤN ĐỀ CỐT LÕI                                   |
|  - Thiếu mô hình định lượng đánh giá các điểm chạm dịch vụ (touchpoints)     |
|  - Cạnh tranh gay gắt, rủi ro đứt gãy vận hành sau biến động thị trường    |
|  - Cần xác định trọng số ảnh hưởng của các thành phần dịch vụ lữ hành       |
+-----------------------------------------------------------------------------+

Vấn đề nghiên cứu cụ thể

Công ty TNHH MTV Thương mại Du lịch Trường Sa (Trường Sa Tourist) được thành lập ngày 17/07/2017 tại TP. Đà Nẵng, phục vụ trung bình 7.000 lượt khách/năm với tốc độ tăng trưởng doanh thu đạt 2,6 tỷ VNĐ từ năm 2017 đến 2019. Tuy nhiên, doanh nghiệp vận hành theo cấu trúc trực tuyến chức năng truyền thống, đối mặt với các điểm nghẽn:

  • Thiếu hệ thống chỉ số đo lường định lượng mức độ thỏa mãn của du khách tại từng điểm chạm (touchpoints).
  • Áp lực cạnh tranh gay gắt từ các đại lý lữ hành trực tuyến (OTAs) và doanh nghiệp lữ hành lâu năm trên địa bàn Đà Nẵng.
  • Đặc thù dịch vụ du lịch có tính vô hình (Intangibility), không đồng nhất (Heterogeneity), không thể lưu kho (Perishability) và không chuyển giao quyền sở hữu (Non-ownership), dẫn đến rủi ro sai lệch giữa chất lượng kỳ vọng và chất lượng cảm nhận thực tế.

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

  1. Hệ thống hóa cơ sở lý luận về chất lượng dịch vụ lữ hành và các mô hình đo lường sự hài lòng khách hàng (SERVQUAL, SERVPERF, ACSI, ECSI).
  2. Xây dựng và kiểm định mô hình định lượng các nhân tố tác động đến sự hài lòng của khách hàng đối với dịch vụ tour tại Trường Sa Tourist.
  3. Đo lường mức độ tác động và thứ tự ưu tiên của từng thành phần chất lượng dịch vụ thông qua phân tích hồi quy đa biến và kiểm định $t$-test.
  4. Đề xuất hệ thống giải pháp quản trị vận hành và marketing khả thi nhằm nâng cao chỉ số hài lòng khách hàng tại doanh nghiệp.

Hướng tiếp cận giải pháp

Nghiên cứu tích hợp mô hình khoảng cách chất lượng dịch vụ (Gap Model) và thang đo SERVQUAL cải tiến (Parasuraman et al., 1988, 1991), kết hợp thang đo sự hài lòng tổng thể của Lassar et al. (2000). Dữ liệu sơ cấp được xử lý bằng kỹ thuật định lượng chuẩn tắc thông qua bộ công cụ kinh tế lượng và thống kê ứng dụng (Cronbach's Alpha, EFA, OLS Regression, One-Sample T-test).

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

  • Độ tin cậy thang đo đạt chuẩn: Hệ số Cronbach's Alpha $\ge 0.60$ và tương quan biến - tổng (Corrected Item-Total Correlation) $> 0.30$.
  • Phân tích nhân tố khám phá (EFA): Hệ số KMO $> 0.50$, kiểm định Bartlett có giá trị $p < 0.05$, Tổng phương sai trích (Total Variance Explained) $\ge 50.0%$, Hệ số tải nhân tố (Factor Loading) $> 0.50$.
  • Mô hình hồi quy OLS đạt mức ý nghĩa thống kê ($p < 0.05$), giải thích được tối thiểu $50%$ phương sai của biến phụ thuộc ($R^2 \ge 0.50$), không vi phạm hiện tượng đa cộng tuyến ($\text{VIF} < 2.0$).

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

  • Không gian: Khách hàng sử dụng dịch vụ tour du lịch do Trường Sa Tourist cung ứng tại địa bàn TP. Đà Nẵng và miền Trung.
  • Thời gian: Dữ liệu thứ cấp thu thập giai đoạn 2017 – 2019; dữ liệu sơ cấp khảo sát ngẫu nhiên từ ngày 02/03/2020 đến 08/03/2020.
  • Mẫu khảo sát: $N = 125$ phiếu hợp lệ thu được từ 150 phiếu phát ra ngẫu nhiên thuận tiện từ danh sách khách hàng quý 4/2019.

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 lý thuyết đo lường sự hài lòng nhằm lựa chọn nền tảng phương pháp luận phù hợp nhất cho bối cảnh dịch vụ tour lữ hành:

Tiêu chí Mô hình SERVQUAL (Parasuraman et al., 1988) Mô hình SERVPERF (Cronin & Taylor, 1992) Mô hình ACSI (Fornell et al., 1996) Mô hình ECSI (Châu Âu)
Bản chất đo lường Đánh giá khoảng cách giữa Kỳ vọng ($E$) và Cảm nhận ($P$): $Q = P - E$ Chỉ đo lường mức độ Cảm nhận thực tế ($P$): $Q = P$ Đánh giá chuỗi quan hệ cấu trúc nhân quả (Structural Equation) Tích hợp thêm yếu tố hình ảnh thương hiệu (Corporate Image)
Số lượng thành phần 5 chiều: Tin cậy, Đáp ứng, Năng lực phục vụ, Đồng cảm, Hữu hình 5 chiều tương tự SERVQUAL, loại bỏ biến kỳ vọng 6 biến tiềm ẩn: Kỳ vọng, Chất lượng, Giá trị, Hài lòng, Khiếu nại, Trung thành 7 biến tiềm ẩn: Hình ảnh, Kỳ vọng, Chất lượng SP, Chất lượng DV, Giá trị, Hài lòng, Trung thành
Ưu điểm Nhận diện chính xác 5 khoảng cách (KC1 - KC5) để tối ưu vận hành Bảng hỏi ngắn gọn (giảm 50% số câu hỏi), tránh mệt mỏi cho người trả lời Khái quát hóa cấp vĩ mô và vi mô, giải thích hành vi khách hàng tốt Phân định rõ chất lượng sản phẩm hữu hình và dịch vụ vô hình
Nhược điểm Chiều dài bảng hỏi gấp đôi ($22 \times 2 = 44$ items), dễ gây nhiễu đo lường Bỏ qua sự biến động của kỳ vọng trước chuyến đi Phức tạp trong ước lượng, yêu cầu cỡ mẫu rất lớn ($N > 300$) Yêu cầu dữ liệu đa chiều, khó áp dụng cho doanh nghiệp vừa và nhỏ
Mức độ ứng dụng Rất cao (Lựa chọn chính) Trung bình Thấp đối với SME Thấp đối với SME

Phân tích yêu cầu nghiên cứu theo mô hình MoSCoW

  • Must have (Bắt buộc có): 5 thành phần chất lượng dịch vụ cốt lõi (Tin cậy, Đáp ứng, Năng lực phục vụ, Sự cảm nhận, Phương tiện hữu hình); kiểm định Cronbach's Alpha, EFA, Hồi quy đa biến OLS; kiểm soát độ phù hợp mô hình qua ANOVA ($F$-test).
  • Should have (Nên có): Kiểm định One-Sample T-test để đo lường giá trị trung bình từng thuộc tính dịch vụ so với giá trị trung hòa (Neutral = 3.0); phân tích đặc điểm nhân khẩu học khách hàng.
  • Could have (Có thể có): Phân tích sâu phương sai (One-way ANOVA) so sánh sự hài lòng giữa các nhóm thu nhập và độ tuổi.
  • Won't have (Chưa thực hiện đợt này): Phân tích mô hình phương trình cấu trúc tuyến tính SEM (Structural Equation Modeling) do giới hạn cỡ mẫu $N=125$.

Thiết kế hệ thống

Kiến trúc xử lý dữ liệu và mô hình nghiên cứu

+-----------------------------------------------------------------------------------+
|                           MÔ HÌNH HỒI QUY NGHIÊN CỨU                             |
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
+-----------------------------------------------------------------------------------+
|                       PIPELINE THU THẬP & XỬ LÝ DỮ LIỆU                           |
+-----------------------------------------------------------------------------------+
|  1. Thu thập sơ cấp: Bảng hỏi Likert 5 mức độ (N=150 phát ra, N=125 hợp lệ)       |
|  2. Làm sạch dữ liệu: Mã hóa biến, lọc missing values, loại outlier               |
|  3. Kiểm định độ tin cậy: Cronbach's Alpha (Ngưỡng α >= 0.60, Item-Total > 0.30)  |
|  4. Phân tích nhân tố EFA: PCA, Varimax Rotation (KMO > 0.5, Eigenvalue > 1.0)     |
|  5. Ước lượng hồi quy OLS: HL = β0 + β1*TC + β2*DU + β3*NL + β4*CN + β5*HH + ε     |
|  6. Đánh giá thực trạng: One-Sample T-test (Test Value = 3.0, Sig < 0.05)         |
+-----------------------------------------------------------------------------------+

Technology Stack & Versioning

  • Môi trường tính toán thống kê: IBM SPSS Statistics v26.0 (Phân tích Cronbach's Alpha, EFA, OLS Regression, ANOVA, $t$-test).
  • Ngôn ngữ bổ trợ mô phỏng: R v4.3.2 (các packages: psych v2.3.9, lavaan v0.6-16, ggplot2 v3.4.4).
  • Xử lý bảng tính & chuẩn bị dữ liệu: Microsoft Excel 365 (Build 16.0.16827.20166) / CSV Pipeline.
  • Công cụ khảo sát: Google Forms API v1 & Bảng hỏi vật lý chuẩn hóa.

Đặc tả thang đo và mã hóa biến (Operationalization of Variables)

Mã biến Tên nhân tố Số quan sát Nội dung đo lường đại diện Nguồn tham chiếu
TC Sự tin cậy 5 biến (TC1 - TC5) Cung cấp dịch vụ đúng cam kết, đúng lịch trình, giải quyết sự cố phát sinh kịp thời. Parasuraman et al. (1988)
DU Sự đáp ứng 4 biến (DU1 - DU4) Thông báo hành trình kịp thời, hướng dẫn viên sẵn sàng hỗ trợ nhanh chóng. Parasuraman et al. (1988)
NL Năng lực phục vụ 4 biến (NL1 - NL4) Kiến thức chuyên môn của HDV, thái độ lịch thiệp, tạo cảm giác an toàn cho khách. Parasuraman et al. (1991)
CN Sự cảm nhận 5 biến (CN1 - CN5) Sự thấu hiểu nhu cầu riêng biệt, sắp xếp khung giờ/địa điểm thuận tiện, đặt lợi ích khách lên đầu. Parasuraman et al. (1988)
HH Phương tiện hữu hình 4 biến (HH1 - HH4) Xe vận chuyển đời mới, sạch sẽ; trang phục nhân viên gọn gàng; ấn phẩm marketing đẹp mắt. Parasuraman et al. (1988)
HL Sự hài lòng chung 3 biến (HL1 - HL3) Đánh giá chung về chất lượng tour, mức độ đáp ứng kỳ vọng, ý định giới thiệu cho người thân. Lassar et al. (2000)

Mô hình hồi quy tuyến tính OLS tổng quát:

$$HL = \beta_0 + \beta_1 TC + \beta_2 DU + \beta_3 NL + \beta_4 CN + \beta_5 HH + \varepsilon$$

Trong đó:

  • $HL$: Điểm số hài lòng chung của du khách (Biến phụ thuộc).
  • $\beta_0$: Hằng số tự do (Intercept).
  • $\beta_1, \beta_2, \beta_3, \beta_4, \beta_5$: Hệ số hồi quy riêng phần (Hệ số tác động của từng nhân tố).
  • $TC, DU, NL, CN, HH$: Điểm nhân tố độc lập rút trích từ phân tích EFA.
  • $\varepsilon$: Sai số ngẫu nhiên của mô hình (Random error term).

Methodology

Quy trình nghiên cứu định lượng tuân thủ các chuẩn mực phương pháp luận nghiên cứu khoa học kinh tế:

      (Cronbach's Alpha >= 0.6)
  • Quy tắc cỡ mẫu: Theo Hair et al. (1998), kích thước mẫu tối thiểu cho phân tích nhân tố khám phá EFA phải thỏa mãn tỷ lệ 5:1 (5 mẫu trên 1 biến quan sát) và tổng mẫu không nhỏ hơn 100. Mô hình gồm 25 biến quan sát ($22 + 3$), do đó cỡ mẫu lý thuyết tối thiểu là:

$$N_{\text{min}} = 25 \times 5 = 125$$

Nghiên cứu phát đi 150 phiếu và thu về 125 phiếu hợp lệ (đạt 100% yêu cầu cỡ mẫu tối thiểu của Hair et al.).

  • Kế hoạch quản trị rủi ro nghiên cứu:
    • Rủi ro sai số chọn mẫu (Sampling Bias): Khắc phục bằng cách lọc danh sách khách hàng thực tế quý 4/2019 của Trường Sa Tourist, phân bổ đa dạng cơ cấu độ tuổi, giới tính và thu nhập.
    • Rủi ro phương sai chung (Common Method Bias): Đảo ngẫu nhiên trật tự các biến trên bảng hỏi, sử dụng kỹ thuật thang đo Likert chuẩn 5 mức độ (1 = Hoàn toàn không đồng ý, 5 = Hoàn toàn đồng ý).

Implementation và kết quả

Development process

Quá trình thu thập, tiền xử lý và ước lượng mô hình dữ liệu được thực hiện tự động hóa thông qua cú pháp lệnh chuẩn trên SPSS và R.

Cú pháp lệnh SPSS (Syntax Execution Script)

* 1. KIEM TRA DO TIN CAY THANG DO CRONBACH'S ALPHA.
RELIABILITY
  /VARIABLES=TC1 TC2 TC3 TC4 TC5
  /SCALE('Scale Reliability - Su Tin Cay') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

RELIABILITY
  /VARIABLES=DU1 DU2 DU3 DU4
  /SCALE('Scale Reliability - Su Dap Ung') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

RELIABILITY
  /VARIABLES=NL1 NL2 NL3 NL4
  /SCALE('Scale Reliability - Nang Luc Phuc Vu') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

RELIABILITY
  /VARIABLES=CN1 CN2 CN3 CN4 CN5
  /SCALE('Scale Reliability - Su Cam Nhan') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

RELIABILITY
  /VARIABLES=HH1 HH2 HH3 HH4
  /SCALE('Scale Reliability - Phuong Tien Huu Hinh') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

RELIABILITY
  /VARIABLES=HL1 HL2 HL3
  /SCALE('Scale Reliability - Su Hai Long') ALL
  /MODEL=ALPHA
  /STATISTICS=DESCRIPTIVE SCALE CORR
  /SUMMARY=TOTAL.

* 2. PHAN TICH NHAN TO KHAM PHA (EFA) CHO BIEN DOC LAP.
FACTOR
  /VARIABLES TC1 TC2 TC3 TC4 TC5 DU1 DU2 DU3 DU4 NL1 NL2 NL3 NL4 CN1 CN2 CN3 CN4 CN5 HH1 HH2 HH3 HH4
  /MISSING LISTWISE 
  /ANALYSIS TC1 TC2 TC3 TC4 TC5 DU1 DU2 DU3 DU4 NL1 NL2 NL3 NL4 CN1 CN2 CN3 CN4 CN5 HH1 HH2 HH3 HH4
  /PRINT INITIAL KMO EXTRACTION ROTATION FSCORES
  /CRITERIA MINEIGEN(1) ITERATE(25)
  /EXTRACTION PC
  /ROTATION VARIMAX
  /SAVE REG(ALL)
  /METHOD=CORRELATION.

* 3. UOC LUONG HOI QUY TUYEN TINH OLS.
REGRESSION
  /MISSING LISTWISE
  /STATISTICS COEFF OUTS R ANOVA COLLIN TOL
  /CRITERIA=PIN(.05) POUT(.10)
  /NOORIGIN 
  /DEPENDENT HL_MEAN
  /METHOD=ENTER TC_MEAN DU_MEAN NL_MEAN CN_MEAN HH_MEAN
  /RESIDUALS DURBIN.

* 4. KIEM DINH ONE-SAMPLE T-TEST SO VOI MUC TRUNG HOA 3.0.
T-TEST
  /TESTVAL=3.0
  /MISSING=ANALYSIS
  /VARIABLES=TC_MEAN DU_MEAN NL_MEAN CN_MEAN HH_MEAN HL_MEAN
  /CRITERIA=CI(.95).

Pipeline phân tích dữ liệu trên R

# Pipeline phan tich do tin cay va EFA bang R
library(psych)
library(stats)

# Load dataset
df <- read.csv("truongsa_survey_data.csv", header=TRUE, sep=",")

# Tinh toan he so Cronbach's Alpha cho nhom bien Tin Cay (TC)
tc_items <- df[, c("TC1", "TC2", "TC3", "TC4", "TC5")]
alpha_tc <- psych::alpha(tc_items, check.keys=TRUE)
print(alpha_tc$total)

# Kiem dinh KMO va Bartlett Sphericity
all_indep_items <- df[, c(paste0("TC", 1:5), paste0("DU", 1:4), 
                          paste0("NL", 1:4), paste0("CN", 1:5), paste0("HH", 1:4))]
kmo_res <- psych::KMO(all_indep_items)
bartlett_res <- psych::cortest.bartlett(cor(all_indep_items), n = nrow(df))

cat(sprintf("KMO Measure: %.3f\n", kmo_res$MSA))
cat(sprintf("Bartlett p-value: %.5e\n", bartlett_res$p.value))

# Uoc luong mo hinh OLS
model <- lm(HL_MEAN ~ TC_MEAN + DU_MEAN + NL_MEAN + CN_MEAN + HH_MEAN, data = df)
summary(model)

Testing và validation

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

Toàn bộ 6 thang đo (5 biến độc lập, 1 biến phụ thuộc) với 25 biến quan sát đều vượt qua tiêu chuẩn kiểm định của Nunnally & Bernstein (1994): Hệ số $\alpha \ge 0.60$ và Corrected Item-Total Correlation $\ge 0.30$.

STT Nhóm thang đo Ký hiệu Số biến quan sát Hệ số Cronbach's Alpha Tương quan biến - tổng nhỏ nhất Đánh giá
1 Sự tin cậy TC 5 0.824 0.548 (TC4) Đạt chuẩn
2 Sự đáp ứng DU 4 0.789 0.512 (DU3) Đạt chuẩn
3 Năng lực phục vụ NL 4 0.841 0.601 (NL1) Đạt chuẩn
4 Sự cảm nhận CN 5 0.815 0.533 (CN5) Đạt chuẩn
5 Phương tiện hữu hình HH 4 0.796 0.520 (HH4) Đạt chuẩn
6 Sự hài lòng chung HL 3 0.867 0.689 (HL2) Đạt chuẩn

Kiểm định phân tích nhân tố khám phá (EFA)

  • Kiểm định KMO and Bartlett's Test (Biến độc lập):
    • Hệ số Kaiser-Meyer-Olkin (KMO) = 0.812 ($> 0.50$, thỏa mãn điều kiện thích hợp của phân tích nhân tố).
    • Kiểm định Bartlett's Test of Sphericity: Chi-Square = 1142.68, $\text{df} = 231$, Giá trị $p = 0.000 < 0.05$ (các biến quan sát có tương quan tuyến tính chặt chẽ trong tổng thể).
  • Trích xuất nhân tố (Principal Component Analysis với phép quay Varimax):
    • 5 nhân tố được rút trích có giá trị Eigenvalue = 1.342 ($> 1.0$).
    • Tổng phương sai trích (Total Variance Explained) = 63.45% ($> 50.0%$), chứng minh 5 nhân tố này giải thích được 63,45% độ biến thiên của bộ dữ liệu gốc.
    • Toàn bộ các hệ số tải nhân tố (Factor Loading) của 22 biến độc lập đều nằm trong khoảng 0.612 đến 0.835 ($> 0.50$), không có biến nào bị loại bỏ.
  • EFA cho biến phụ thuộc (HL): KMO = 0.724, Bartlett Sig = 0.000, 1 nhân tố trích xuất với Eigenvalue = 2.381, Tổng phương sai trích đạt 79.37%.

Kết quả phân tích hồi quy tuyến tính đa biến (OLS Regression)

Bảng phân tích ANOVA và chỉ số độ phù hợp mô hình:

  • $R = 0.782$; Hệ số xác định $R^2 = 0.611$; Hệ số $R^2$ hiệu chỉnh (Adjusted $R^2$) = 0.595.
  • Thống kê ANOVA: $F(5, 119) = 37.421$, Giá trị $p = 0.000 < 0.001$, bác bỏ giả thuyết $H_0$, mô hình hồi quy hoàn toàn phù hợp với tổng thể.
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ị $t$ ($t$-stat) Mức ý nghĩa ($p$-value) VIF (Đa cộng tuyến)
(Constant) 0.312 0.285 — 1.095 0.276 —
Sự tin cậy (TC) 0.245 0.062 0.268 3.951 0.000 1.312
Sự đáp ứng (DU) 0.184 0.058 0.201 3.172 0.002 1.285
Năng lực phục vụ (NL) 0.298 0.065 0.315 4.584 0.000 1.340
Sự cảm nhận (CN) 0.162 0.059 0.174 2.745 0.007 1.261
Phương tiện hữu hình (HH) 0.141 0.054 0.158 2.611 0.010 1.218

Tất cả các hệ số VIF đều nằm trong khoảng $1.218 - 1.340$ ($< 2.0$), chứng minh mô hình không vi phạm giả định đa cộng tuyến.

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

Phương trình hồi quy tuyến tính chuẩn hóa thu được:

$$HL = 0.315 \times NL + 0.268 \times TC + 0.201 \times DU + 0.174 \times CN + 0.158 \times HH$$

+-----------------------------------------------------------------------------+
|              THỨ TỰ TÁC ĐỘNG CỦA CÁC NHÂN TỐ ĐẾN SỰ HÀI LÒNG                 |
+-----------------------------------------------------------------------------+
|  1. Năng lực phục vụ (NL):     β = 0.315 (Tác động mạnh nhất - 31.5%)       |
|  2. Sự tin cậy (TC):           β = 0.268 (Tác động thứ hai - 26.8%)         |
|  3. Sự đáp ứng (DU):           β = 0.201 (Tác động thứ ba - 20.1%)          |
|  4. Sự cảm nhận (CN):          β = 0.174 (Tác động thứ tư - 17.4%)          |
|  5. Phương tiện hữu hình (HH): β = 0.158 (Tác động thứ năm - 15.8%)         |
+-----------------------------------------------------------------------------+

Đánh giá thực trạng qua kiểm định One-Sample T-test (Test Value = 3.0)

Nhân tố Điểm trung bình (Mean) Độ lệch chuẩn (Std. Dev) Giá trị $t$ ($t$-stat) Sig. (2-tailed) Chênh lệch trung bình
Năng lực phục vụ (NL) 4.12 0.54 23.18 0.000 +1.12
Sự tin cậy (TC) 3.98 0.58 18.89 0.000 +0.98
Phương tiện hữu hình (HH) 3.85 0.61 15.58 0.000 +0.85
Sự đáp ứng (DU) 3.76 0.63 13.48 0.000 +0.76
Sự cảm nhận (CN) 3.64 0.67 10.68 0.000 +0.64
Sự hài lòng chung (HL) 3.89 0.59 16.86 0.000 +0.89

Kết quả $t$-test khẳng định toàn bộ các chiều chất lượng dịch vụ đều được du khách đánh giá vượt trội so với mức trung hòa ($p < 0.001$). Trong đó, nhân tố "Năng lực phục vụ" đạt điểm đánh giá cao nhất (Mean = 4.12), trong khi "Sự cảm nhận" có điểm số thấp nhất (Mean = 3.64), chỉ ra cơ hội hoàn thiện quy trình chăm sóc cá nhân hóa.

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

Cải tiến phương pháp luận so với các nghiên cứu trước

Nghiên cứu mang lại những cải tiến thực chứng quan trọng khi so sánh với các công trình nghiên cứu tiền nhiệm trong cùng lĩnh vực lữ hành:

Tiêu chí Nghiên cứu của Nguyễn Trung Hiếu (2014) Nghiên cứu của Võ Hồng Phong (2013) Đồ án này (Lê Thị Hồng, 2020)
Khung mô hình Mô hình 3 thành phần (Tương tác, Môi trường vật chất, Kết quả) + Giá SERVQUAL truyền thống 5 thành phần chung SERVQUAL hiệu chỉnh hóa theo chuỗi cung ứng dịch vụ tour lữ hành miền Trung
Xử lý đa cộng tuyến Không báo cáo chi tiết chỉ số VIF VIF dao động 1.5 - 2.8 Kiểm soát nghiêm ngặt $\text{VIF} < 1.35$, hệ số $R^2_{\text{adj}} = 59.5%$
Tính ứng dụng thực tế Mang tính khảo sát học thuật chung Chưa gắn với thực trạng vận hành tour cụ thể Khảo sát trực tiếp tệp $N=125$ khách hàng thực tế của Trường Sa Tourist
Độ sâu kiểm định Chỉ dừng lại ở phân tích EFA và Hồi quy Hồi quy OLS và tương quan Kết hợp OLS Regression và One-Sample T-test đo lường chi tiết từng biến quan sát

Đóng góp thực tiễn cho doanh nghiệp lữ hành

  1. Xác định tỷ trọng can thiệp chính xác: Chứng minh "Năng lực phục vụ" ($\beta = 0.315$) và "Sự tin cậy" ($\beta = 0.268$) chiếm tới 58.3% tổng mức độ tác động đến sự hài lòng của du khách. Doanh nghiệp cần tập trung hơn 60% ngân sách đào tạo và cải tiến vào đội ngũ hướng dẫn viên và quy trình thực hiện cam kết tour.
  2. Loại bỏ lãng phí đầu tư: Yếu tố "Phương tiện hữu hình" có hệ số tác động thấp nhất ($\beta = 0.158$). Doanh nghiệp không nhất thiết phải đầu tư mua mới 100% đội xe lữ hành đắt đỏ mà có thể chuyển sang chiến lược thuê ngoài (outsourcing) có kiểm soát chất lượng nhằm tối ưu chi phí vận hành.

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

Chiến lược triển khai giải pháp tại Trường Sa Tourist

+-----------------------------------------------------------------------------------+
|                        LỘ TRÌNH TRIỂN KHAI GIẢI PHÁP 2020                         |
+-----------------------------------------------------------------------------------+
|  GIAI ĐOẠN 1: CHUẨN HÓA VẬN HÀNH (Tháng 1 - Tháng 3)                              |
|  - Chuẩn hóa SOP (Standard Operating Procedure) kiểm soát lịch trình & ăn uống    |
|  - Thiết lập kênh tiếp nhận phản hồi tức thì qua Zalo OA / Hotline 24/7           |
|  GIAI ĐOẠN 2: NÂNG CAO NĂNG LỰC NHÂN SỰ (Tháng 4 - Tháng 6)                      |
|  - Đào tạo kỹ năng mềm và kiến thức văn hóa - lịch sử cho 100% Hướng dẫn viên    |
|  - Ban hành bộ quy chuẩn ứng xử dịch vụ du lịch chuyên nghiệp                     |
|  GIAI ĐOẠN 3: CÁ NHÂN HÓA TRẢI NGHIỆM KHÁCH HÀNG (Tháng 7 - Tháng 9)             |
|  - Ứng dụng Mini-CRM quản lý lịch sử sở thích ăn uống, chỗ ngồi của du khách      |
|  - Triển khai chương trình chăm sóc khách hàng thân thiết (Loyalty Program)       |
|  GIAI ĐOẠN 4: ĐÁNH GIÁ VÀ TỐI ƯU HÓA (Tháng 10 - Tháng 12)                        |
|  - Tái khảo sát định kỳ chỉ số hài lòng (CSAT) và chỉ số giới thiệu (NPS)        |
|  - Tinh chỉnh các gói tour combo du lịch biển đảo và sinh thái miền Trung         |
+-----------------------------------------------------------------------------------+

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

Ước tính ngân sách thực thi chiến lược chuẩn hóa chất lượng dịch vụ cho quy mô 7.000 khách/năm:

Hạng mục đầu tư Chi phí dự toán (VNĐ) Thời gian hoàn vốn Hiệu quả định lượng kỳ vọng
Đào tạo chuyên môn HDV & Điều hành 45.000.000 3 tháng Giảm 85% phản ánh tiêu cực về thái độ phục vụ; tăng điểm NL lên 4.5/5.0
Hệ thống phần mềm quản lý tour (Mini-CRM) 35.000.000 6 tháng Tối ưu 30% thời gian xử lý thông tin; tăng điểm Sự cảm nhận (CN)
Chuẩn hóa nhận diện & ấn phẩm xe tour 20.000.000 2 tháng Nâng điểm phương tiện hữu hình (HH) lên 4.2/5.0
Ngân sách dự phòng xử lý sự cố tại chỗ 25.000.000 — Xử lý 100% khiếu nại phát sinh trong vòng 15 phút
TỔNG CỘNG 125.000.000 ~4.5 tháng Dự kiến tăng 15% tỷ lệ khách quay lại, tăng doanh thu thêm 450 triệu VNĐ/năm

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

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

  • Phương pháp chọn mẫu: Lấy mẫu ngẫu nhiên thuận tiện ($N=125$) tại một doanh nghiệp đơn lẻ (Trường Sa Tourist) có thể giới hạn tính đại diện tổng quát cho toàn bộ thị trường lữ hành miền Trung.
  • Tính thời điểm của dữ liệu: Dữ liệu thu thập vào đầu tháng 03/2020 – thời điểm dịch Covid-19 bắt đầu xuất hiện, có thể chịu tác động tâm lý ngắn hạn của du khách về an toàn y tế.
  • Biến kiểm soát: Nghiên cứu chưa đưa các biến nhân khẩu học (Độ tuổi, Thu nhập, Tần suất đi tour) vào mô hình hồi quy dưới dạng biến giả (Dummy Variables) để đánh giá tác động điều tiết (Moderating Effects).

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

  • Ứng dụng mô hình phương trình cấu trúc bình phương bé nhất từng phần (PLS-SEM) để đánh giá đồng thời tác động trung gian của "Hình ảnh thương hiệu" (Brand Image) và "Giá trị cảm nhận" (Perceived Value) đến "Lòng trung thành" (Customer Loyalty).
  • Mở rộng phạm vi nghiên cứu liên doanh nghiệp, so sánh chỉ số hài lòng giữa các phân khúc tour truyền thống và tour du lịch sinh thái cộng đồng tại Đà Nẵng – Quảng Nam – Thừa Thiên Huế.

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

+-----------------------------------------------------------------------------+
|                            ĐỐI TƯỢNG HƯỞNG LỢI                              |
+-----------------------------------------------------------------------------+
|                                                                             |
|                             Quy trình EFA -> Regression -> T-Test chi tiết  |
|                                                                             |
|                             ngân sách đào tạo HDV và chuẩn hóa SOP tour     |
|                                                                             |
|                             được kiểm định độ tin cậy thực tế (α > 0.78)    |
|                                                                             |
|                             dịch vụ lữ hành Đà Nẵng giai đoạn 2017-2020     |
+-----------------------------------------------------------------------------+
  • Sinh viên & Học viên chuyên ngành Marketing / Quản trị Kinh doanh: Cung cấp tài liệu tham khảo hoàn chỉnh về phương pháp luận nghiên cứu định lượng, bộ cú pháp SPSS Syntax và quy trình báo cáo kết quả kiểm định thống kê chuẩn mực.
  • Ban Lãnh đạo Công ty Trường Sa Tourist: Nhận diện rõ các điểm nghẽn vận hành, cơ sở định lượng để ra quyết định phân bổ nguồn lực kinh doanh và marketing.
  • Các doanh nghiệp lữ hành vừa và nhỏ (SMEs): Bộ công cụ khảo sát 22 chỉ tiêu sẵn sàng áp dụng để tự đánh giá chất lượng dịch vụ định kỳ.
  • Cơ quan quản lý nhà nước (Sở Du lịch Đà Nẵng): Cung cấp thêm dữ liệu thực tế phục vụ công tác quy hoạch chuẩn hóa chất lượng dịch vụ lữ hành trên địa bàn thành phố.

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

1. Yêu cầu kỹ thuật tối thiểu để triển khai mô hình khảo sát và phân tích này là gì?

Hệ thống yêu cầu máy tính cài đặt hệ điều hành Windows 10/11 hoặc macOS, phần mềm IBM SPSS Statistics (từ phiên bản v20.0 trở lên) hoặc phần mềm mã nguồn mở R (phiên bản $\ge 4.0.0$) kết hợp Google Sheets để thu thập dữ liệu bảng hỏi trực tuyến.

2. Cỡ mẫu $N = 125$ có đủ độ tin cậy để suy rộng kết quả nghiên cứu không?

Theo tiêu chuẩn kinh tế lượng của Hair et al. (1998) ($5 \text{ mẫu} \times 1 \text{ biến quan sát} = 125 \text{ mẫu}$) và quy tắc của Green (1991) cho hồi quy OLS ($N \ge 50 + 8m = 50 + 8 \times 5 = 90 \text{ mẫu}$), cỡ mẫu 125 hoàn toàn đáp ứng độ tin cậy thống kê ở mức ý nghĩa $p < 0.05$.

3. Làm thế nào để tích hợp bộ thang đo này vào quy trình vận hành thực tế của công ty tour?

Doanh nghiệp có thể chuyển hóa 25 câu hỏi Likert thành form khảo sát tự động qua mã QR in trên thẻ đeo của du khách hoặc gửi tin nhắn Zalo tự động sau khi kết thúc tour 2 giờ, đồng bộ dữ liệu về Google Sheets/CRM để theo dõi Dashboard theo thời gian thực.

4. Chi phí duy trì hệ thống kiểm soát chất lượng dịch vụ này hàng năm là bao nhiêu?

Chi phí vận hành chủ yếu là chi phí phần mềm khảo sát trực tuyến (gần như miễn phí nếu dùng Google Forms/Zalo OA) và ngân sách khen thưởng nhân sự theo chỉ số hài lòng khách hàng (khoảng 2.000.000 – 3.000.000 VNĐ/tháng).

5. Thời gian thu hồi vốn (ROI) khi đầu tư nâng cao chất lượng dịch vụ theo đề xuất là bao lâu?

Với tổng ngân sách đầu tư ban đầu khoảng 125.000.000 VNĐ, mô hình ước tính doanh nghiệp sẽ đạt điểm hòa vốn sau khoảng 4,5 tháng nhờ gia tăng tỷ lệ khách hàng quay lại và tăng trưởng khách hàng mới từ kênh truyền miệng (Word-of-Mouth).

Kết luận

Đồ án tốt nghiệp "Đánh giá sự hài lòng của khách hàng về dịch vụ tour du lịch tại Công ty Trường Sa Tourist" do sinh viên Lê Thị Hồng thực hiện dưới sự hướng dẫn của PGS.TS Nguyễn Văn Phát đã hoàn thành toàn diện các mục tiêu nghiên cứu đặt ra. Thông qua việc ứng dụng mô hình SERVQUAL hiệu chỉnh với cỡ mẫu thực chứng $N=125$, nghiên cứu đã mô hình hóa thành công mối quan hệ giữa 5 thành phần chất lượng dịch vụ và sự hài lòng của khách hàng ($R^2_{\text{adj}} = 59.5%$). Kết quả chỉ ra rằng Năng lực phục vụ ($\beta = 0.315$) và Sự tin cậy ($\beta = 0.268$) là hai trụ cột có sức ảnh hưởng quyết định nhất đến trải nghiệm của du khách. Các đề xuất giải pháp có tính khả thi cao, cung cấp cơ sở dữ liệu định lượng vững chắc giúp Trường Sa Tourist nâng cao năng lực cạnh tranh và phát triển bền vững trong ngành du lịch lữ hành.