Giới thiệu dự án

Thị trường đồ uống tại Việt Nam đang chứng kiến sự chuyển dịch hành vi tiêu dùng mạnh mẽ. Theo số liệu từ Tổ chức Cà phê Quốc tế (International Coffee Organization - ICO), Việt Nam giữ vững vị thế quốc gia xuất khẩu cà phê lớn thứ hai toàn cầu với sản lượng xấp xỉ 1,650,000 tấn/năm. Thống kê thị trường của Kantar Worldpanel chỉ ra rằng có tới 60% người dân tại TP. Hồ Chí Minh duy trì thói quen uống cà phê khi ra ngoài. Tuy nhiên, thị trường nội địa đối mặt với khủng hoảng niềm tin nghiêm trọng do vấn nạn cà phê pha tạp bắp, đậu nành rang cháy, hóa chất tạo bọt, tạo đắng không rõ nguồn gốc và các vụ vi phạm an toàn thực phẩm nghiêm trọng. Theo Bộ Y tế, các bệnh lý liên quan đến thực phẩm không an toàn tiêu tốn khoảng 0.22% GDP hàng năm.

Khóa luận tốt nghiệp "Nghiên cứu các yếu tố ảnh hưởng đến hành vi tiêu dùng cà phê nguyên chất của khách hàng từ 18 tuổi trở lên tại Thành phố Hồ Chí Minh" được thực hiện bởi sinh viên Nguyễn Mỹ Tiên (Khoa Quản trị Kinh doanh, Trường Đại học Công nghiệp TP. HCM) dưới sự hướng dẫn của ThS. Lê Thị Thanh Trúc nhằm giải quyết bài toán cốt lõi: Xác định các rào cản và động lực thúc đẩy người tiêu dùng chuyển dịch sang sản phẩm cà phê hạt rang xay nguyên chất 100%.

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

  1. Hệ thống hóa cơ sở lý thuyết về hành vi người tiêu dùng và xây dựng khung thang đo chuẩn hóa cho sản phẩm cà phê nguyên chất.
  2. Đo lường mức độ tác động của từng nhân tố đến hành vi tiêu dùng thực tế thông qua các kiểm định định lượng.
  3. Đề xuất ma trận giải pháp quản trị marketing (4P/7P) có tính khả thi cao nhằm định hướng sản xuất và mở rộng thị phần.

Phạm vi nghiên cứu tập trung vào đối tượng khách hàng từ 18 tuổi trở lên đang sinh sống, học tập và làm việc tại khu vực TP. Hồ Chí Minh với cỡ mẫu khảo sát hợp lệ $N = 300$.


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

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

Nghiên cứu tiến hành tổng hợp và đánh giá các mô hình lý thuyết nền tảng cùng các công trình thực nghiệm trước đây để xác định khoảng trống nghiên cứu (Research Gap).

Mô hình / Tác giả Ưu điểm Hạn chế Khả năng áp dụng vào đề tài
Thuyết TRA (Fishbein & Ajzen, 1975) Giải thích rõ mối liên hệ giữa niềm tin, thái độ và chuẩn chủ quan. Chưa tính đến yếu tố kiểm soát nguồn lực, tài chính và rào cản ngoại cảnh. Cung cấp nền tảng xây dựng biến chuẩn chủ quan và thái độ.
Thuyết TPB (Ajzen, 1991) Bổ sung biến "Nhận thức kiểm soát hành vi", tăng độ chuẩn xác dự báo. Mô hình tổng quát, thiếu các thuộc tính đặc thù của ngành F&B/nông sản. Đóng vai trò khung lý thuyết mở rộng để tích hợp các biến độc lập.
Trần Thị Mỹ Thuận (2016) Khảo sát trực tiếp giới trẻ tại TP.HCM về cà phê nguyên chất ($N = 200$). Giới hạn độ tuổi hẹp (18-30 tuổi), chưa phân tích sâu nhóm khách hàng trung niên. Kế thừa thang đo yếu tố cá nhân và ý thức sức khỏe.
Đoàn Phương Trang & CS (2018) Đo lường chi tiết yếu tố xuất xứ, bao bì và chất lượng cảm nhận ($N = 120$). Khảo sát thị trường nước ngoài (Nhật Bản), sai khác về hành vi văn hóa. Kế thừa thang đo thương hiệu và giá trị cảm nhận.

Bảng ưu tiên yêu cầu nghiên cứu (MoSCoW Prioritization)

  • Must-have: Kiểm định độ tin cậy thang đo (Cronbach's Alpha $\ge 0.60$), phân tích nhân tố khám phá (EFA), hồi quy đa biến (Multiple Linear Regression) kiểm tra đa cộng tuyến ($VIF < 2.0$).
  • Should-have: Kiểm định sâu ANOVA và T-test nhằm phát hiện sai biệt theo nhân khẩu học (độ tuổi, giới tính, thu nhập, tần suất sử dụng).
  • Could-have: Tích hợp mô hình đối sánh giá trị trung bình Mean để xếp hạng mức độ ưu tiên của khách hàng.
  • Won't-have (giai đoạn này): Phân tích 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).

Thiết kế hệ thống và mô hình nghiên cứu

Mô hình nghiên cứu đề xuất bao gồm 6 biến độc lập với 22 biến quan sát và 1 biến phụ thuộc gồm 6 biến quan sát:

$$\text{HV} = \beta_0 + \beta_1 \text{CN} + \beta_2 \text{THSP} + \beta_3 \text{GTCN} + \beta_4 \text{SK} + \beta_5 \text{GC} + \beta_6 \text{CT} + \epsilon$$

Ngăn xếp công nghệ và công cụ xử lý (Technology Stack)

  • Thu thập dữ liệu: Google Forms API kết hợp bảng khảo sát giấy chuẩn hóa theo thang đo khoảng cách Likert 5 điểm ($1 = \text{Rất không đồng ý}$ đến $5 = \text{Rất đồng ý}$).
  • Xử lý và làm sạch dữ liệu: Microsoft Excel 2016 MSO.
  • Phân tích thống kê chuyên sâu: IBM SPSS Statistics phiên bản 20.0 (hỗ trợ phân tích Descriptive, Cronbach's Alpha, EFA với phép xoay Varimax, Pearson Correlation, OLS Regression, One-Way ANOVA).
  • Hệ thống hóa quy trình phân tích tự động: Tương thích cú pháp Python (pandas v2.1.0, factor-analyzer v0.5.0, statsmodels v0.14.0).

Phương pháp luận (Methodology)

Nghiên cứu phối hợp chặt chẽ giữa hai pha nghiên cứu:

  1. Nghiên cứu định tính: Nghiên cứu tài liệu thứ cấp, thảo luận nhóm chuyên gia và giảng viên hướng dẫn để điều chỉnh 28 biến quan sát phù hợp với văn hóa tiêu dùng thực tế tại TP.HCM.
  2. Nghiên cứu định lượng:
    • Pha sơ bộ ($N = 50$): Đánh giá độ tin cậy sơ cấp của bảng hỏi.
    • Pha chính thức ($N = 320$, hợp lệ $N = 300$): Lấy mẫu thuận tiện phi xác suất tại các điểm bán (Take-away, cà phê vỉa hè, chuỗi quán cao cấp) và nền tảng trực tuyến.

Implementation và kết quả

Quy trình phân tích dữ liệu và thuật toán xử lý

Dữ liệu sơ cấp sau khi làm sạch được đưa vào quy trình phân tích định lượng chuẩn thông qua cú pháp phân tích hồi quy và trích xuất nhân tố.

import pandas as pd
import numpy as np
from factor_analyzer import FactorAnalyzer
import statsmodels.api as sm

# 1. Load và tiền xử lý bộ dữ liệu 300 mẫu hợp lệ
df = pd.read_csv("coffee_consumer_behavior_hcm_n300.csv")

# 2. Thuật toán tính toán Cronbach's Alpha
def cronbach_alpha(itemscores):
    itemscores = np.asarray(itemscores)
    itemvars = itemscores.var(axis=0, ddof=1)
    tscores = itemscores.sum(axis=1)
    nitems = itemscores.shape[1]
    return (nitems / (nitems - 1)) * (1 - (itemvars.sum() / tscores.var(ddof=1)))

# 3. Trích xuất nhân tố EFA (KMO & Varimax Rotation)
fa = FactorAnalyzer(n_factors=6, rotation="varimax", method="principal")
fa.fit(df[['CN1','CN2','CN3','THSP1','THSP2','THSP3','THSP4',
           'GTCN1','GTCN2','GTCN3','SK1','SK2','SK3','SK4',
           'GC1','GC2','GC3','CT1','CT2','CT3','CT4','CT5']])

# 4. Mô hình Hồi quy tuyến tính đa biến (OLS Regression Pipeline)
X = df[['CN', 'THSP', 'GTCN', 'SK', 'GC', 'CT']]
y = df['HV']
X_const = sm.add_constant(X)
ols_model = sm.OLS(y, X_const).fit()
print(ols_model.summary())

Kiểm định và đánh giá thang đo (Testing & Validation)

  1. Kiểm định Cronbach's Alpha: Toàn bộ 7 nhóm thang đo đều đạt độ nhất quán nội tại cao ($\alpha > 0.70$), hệ số tương quan biến - tổng (Corrected Item-Total Correlation) của 28 biến quan sát đều vượt ngưỡng chuẩn $0.30$, không có biến rác bị loại bỏ.
  2. Phân tích nhân tố khám phá (EFA):
    • Biến độc lập: Hệ số $KMO = 0.842 \in [0.5, 1.0]$, kiểm định Bartlett đạt ý nghĩa thống kê ($Sig. = 0.000 < 0.05$). Trích xuất thành công 6 nhân tố tại $Eigenvalue = 1.218 > 1.0$, tổng phương sai trích (Total Variance Explained) đạt $63.45% > 50%$. Tất cả các trọng số nhân tố (Factor Loading) đều đạt $\ge 0.55$.
    • Biến phụ thuộc: Hệ số $KMO = 0.865$, phương sai trích đạt $61.20%$, đơn nhân tố hội tụ hoàn hảo.
  3. Phân tích tương quan Pearson: Tất cả 6 biến độc lập đều có tương quan tuyến tính có ý nghĩa thống kê với biến phụ thuộc ($r \in [0.312, 0.654], p < 0.01$).

Kết quả đạt được từ mô hình hồi quy đa biến

Mô hình hồi quy đa biến đạt độ tương thích cao với dữ liệu thực tế, giải thích phần lớn sự biến thiên trong hành vi tiêu dùng.

Biến độc lập Mã hóa Hệ số Beta chuẩn hóa ($\beta$) Mức ý nghĩa ($Sig.$) Hệ số phóng đại phương sai ($VIF$) Kết luận giả thuyết
Thương hiệu sản phẩm THSP 0.312 0.000 1.245 Chấp nhận $H_2$
Nhận thức về sức khỏe SK 0.268 0.000 1.310 Chấp nhận $H_4$
Giá trị cảm nhận GTCN 0.195 0.002 1.187 Chấp nhận $H_3$
Yếu tố cá nhân CN 0.154 0.008 1.220 Chấp nhận $H_1$
Giá cả GC 0.138 0.015 1.156 Chấp nhận $H_5$
Chiêu thị CT 0.114 0.034 1.289 Chấp nhận $H_6$
  • Hệ số xác định: $R^2 = 0.584$, $R^2 \text{ hiệu chỉnh} = 0.575$ (Mô hình giải thích được $57.5%$ sự biến thiên của hành vi tiêu dùng cà phê nguyên chất).
  • Kiểm định F (ANOVA): $F = 68.74, Sig. = 0.000$, khẳng định mô hình hồi quy hoàn toàn phù hợp.
  • Kiểm định hiện tượng tự tương quan & đa cộng tuyến: Hệ số Durbin-Watson đạt $d = 1.872$ (nằm trong khoảng $1.5 - 2.5$), tất cả $VIF < 1.4$ (bác bỏ hoàn toàn hiện tượng đa cộng tuyến).

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

  1. Chuẩn hóa bộ thang đo chuyên biệt: Đề tài đã tinh chỉnh và hoàn thiện bộ 28 biến quan sát thích ứng chuẩn xác cho phân khúc cà phê sạch/nguyên chất tại đô thị lớn, khắc phục tính chung chung của các thang đo F&B truyền thống.
  2. Khám phá tác động áp đảo của Thương hiệu & Sức khỏe: Khác với các nghiên cứu cà phê thông thường nơi "Giá cả" đóng vai trò tiên quyết, nghiên cứu chứng minh người tiêu dùng cà phê nguyên chất tại TP.HCM sẵn sàng chi trả mức giá cao hơn (mean biến $HV4 = 3.82/5.0$) nếu sản phẩm chứng minh được tính an toàn thực phẩm và minh bạch xuất xứ.
  3. Phát hiện khác biệt nhân khẩu học: Kiểm định One-Way ANOVA chỉ ra sự khác biệt có ý nghĩa thống kê ($p < 0.05$) về hành vi tiêu dùng giữa các nhóm thu nhập ($> 15$ triệu đồng/tháng có xu hướng chọn thương hiệu có chứng nhận hữu cơ) và nhóm nghề nghiệp (nhân viên văn phòng tiêu thụ với tần suất cao nhất: 4-6 lần/tuần).

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

Dựa trên trọng số hồi quy chuẩn hóa $\beta$, doanh nghiệp F&B cần triển khai ma trận giải pháp theo lộ trình 4 giai đoạn:

  1. Chiến lược Thương hiệu (Ưu tiên số 1 - $\beta = 0.312$): Đầu tư bao bì màng nhôm bảo quản van một chiều, in ấn minh bạch vùng trồng (Cầu Đất, Buôn Ma Thuột), tỷ lệ phối trộn Robusta/Arabica và tích hợp mã QR truy xuất nguồn gốc nông trại.
  2. Chiến lược Định vị Sức khỏe ($\beta = 0.268$): Đẩy mạnh chiến dịch truyền thông "Cà phê mộc - Không phụ gia", công bố kết quả kiểm nghiệm hàm lượng Caffeine và chỉ số an toàn độc tố vi nấm Ochratoxin A.
  3. Chiến lược Giá trị cảm nhận ($\beta = 0.195$): Đào tạo nhân viên pha chế (Barista) tư vấn khách hàng cách phân biệt nước cà phê nâu cánh gián nguyên bản với màu đen đậm đặc của cà phê trộn bắp cháy.

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

  • Hạn chế mẫu nghiên cứu: Phương pháp chọn mẫu thuận tiện phi xác suất với quy mô $N = 300$ tập trung chủ yếu tại khu vực nội thành TP.HCM, chưa bao phủ toàn diện các huyện ngoại thành.
  • Biến ngoại sinh: Chưa tích hợp các yếu tố công nghệ mới như phương thức đặt hàng qua ứng dụng giao đồ ăn (Food Delivery Apps) hay xu hướng tiêu dùng cà phê đặc sản (Specialty Coffee).
  • Hướng phát triển: Các nghiên cứu tiếp theo nên mở rộng quy mô mẫu đa vùng miền ($N \ge 1000$), ứng dụng phương pháp lấy mẫu phân tầng ngẫu nhiên và sử dụng mô hình cấu trúc tuyến tính SEM/PLS-SEM để kiểm định các biến trung gian như lòng trung thành thương hiệu và niềm tin khách hàng.

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

  • Sinh viên & Học viên cao học: Bản mẫu nghiên cứu định lượng chuẩn mực từ thiết kế bảng hỏi, kiểm định EFA đến hồi quy OLS trong lĩnh vực Marketing và Quản trị Kinh doanh.
  • Doanh nghiệp F&B & Cơ sở rang xay: Cung cấp cơ sở dữ liệu thực chứng và trọng số định lượng để tối ưu hóa ngân sách R&D, thiết kế bao bì và phân bổ chi phí truyền thông.
  • Nhà nghiên cứu thị trường: Khung thang đo chuẩn hóa đã được kiểm định độ tin cậy và giá trị hội tụ, sẵn sàng kế thừa cho các nghiên cứu hành vi người tiêu dùng ngành hàng thực phẩm sạch.

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

1. Tại sao nhân tố "Chiêu thị" lại có hệ số Beta thấp nhất trong mô hình ($\beta = 0.114$)?

Khách hàng uống cà phê nguyên chất có mức độ am hiểu và quan tâm sâu sắc đến sức khỏe cùng chất lượng cốt lõi. Do đó, các chương trình quảng cáo rầm rộ hay giảm giá ngắn hạn ít có tính quyết định bằng uy tín thương hiệu và hương vị tự nhiên.

2. Mô hình có thể ứng dụng cho các tỉnh thành khác ngoài TP.HCM không?

Có thể ứng dụng trực tiếp bộ thang đo 28 biến quan sát. Tuy nhiên, cần tái kiểm định độ tin cậy Cronbach's Alpha và EFA do đặc tính văn hóa thưởng thức cà phê giữa các miền (miền Bắc chuộng cà phê đậm/pha phin truyền thống, miền Nam chuộng cà phê đá, tiện lợi) có sự phân hóa nhất định.

3. Cần cấu hình phần mềm tối thiểu nào để tái lập phân tích này?

Cần máy tính cá nhân trang bị phần mềm IBM SPSS Statistics (phiên bản 20.0 trở lên) hoặc môi trường Python 3.8+ cài đặt sẵn các thư viện numpy, pandas, statsmodels, factor_analyzer để chạy lại toàn bộ mô hình hồi quy và kiểm định thống kê.


Kết luận

Đồ án khóa luận tốt nghiệp của tác giả Nguyễn Mỹ Tiên đã hoàn thành xuất sắc các mục tiêu nghiên cứu đặt ra, cung cấp cái nhìn thực nghiệm sâu sắc về thị trường cà phê nguyên chất tại TP.HCM. Thông qua hệ thống phân tích định lượng chặt chẽ trên cỡ mẫu $N = 300$, đề tài khẳng định vai trò then chốt của Thương hiệu sản phẩm ($\beta = 0.312$) và Nhận thức về sức khỏe ($\beta = 0.268$) trong việc định hình hành vi mua sắm. Kết quả này không chỉ đóng góp giá trị học thuật quan trọng cho chuyên ngành Marketing mà còn là kim chỉ nam thực tiễn giúp các doanh nghiệp cà phê Việt Nam xây dựng chuỗi giá trị minh bạch, nâng tầm nông sản sạch và bảo vệ sức khỏe cộng đồng.