Giới thiệu dự án

Thị trường mỹ phẩm và sản phẩm chăm sóc cá nhân tại Việt Nam đang chứng kiến sự chuyển dịch mang tính cấu trúc. Theo báo cáo từ Mintel (London), quy mô thị trường mỹ phẩm Việt Nam đạt doanh thu ước tính khoảng 51.000 tỷ đồng (tương đương 2,3 tỷ USD), với tốc độ tăng trưởng kép hàng năm (CAGR) duy trì ở mức hai con số. Đồng thời, dữ liệu từ NPD Group chỉ ra rằng phân khúc mỹ phẩm hữu cơ và tự nhiên toàn cầu đạt quy mô hàng trăm triệu USD với hơn 80% người tiêu dùng cam kết ưu tiên thành phần thực vật bền vững.

Trong bối cảnh đó, mỹ phẩm thuần chay (Vegan Cosmetics) nổi lên như một xu hướng tiêu dùng xanh, nhân đạo, hoàn toàn loại bỏ dẫn xuất động vật (sáp ong, mỡ cừu, gelatin, dịch ốc sên) và không thử nghiệm trên động vật (Cruelty-Free). Cocoon (thành lập năm 2013) là thương hiệu nội địa tiên phong tại Việt Nam đạt các chứng nhận quốc tế danh giá: Leaping Bunny (Cruelty Free International), PETA, và The Vegan Society.

Vấn đề nghiên cứu và điểm nghẽn thực tế

Mặc dù sở hữu lợi thế sản phẩm bản địa (cà phê Đắk Lắk, bí đao, dầu dừa Bến Tre, hoa hồng Cao Bằng), Cocoon phải đối mặt với sự cạnh tranh khốc liệt từ các tập đoàn đa quốc gia (L'Oréal, Unilever, Amorepacific) và các thương hiệu mỹ phẩm xách tay. Điểm nghẽn lớn nhất của doanh nghiệp là thiếu một mô hình định lượng chuẩn xác nhằm đo lường mức độ tác động của từng điểm chạm (touchpoint) đến sự hài lòng của khách hàng (Customer Satisfaction - CSAT), dẫn đến việc phân bổ nguồn lực marketing và R&D chưa đạt hiệu quả tối ưu.

Mục tiêu dự án

  1. Xác định các nhân tố cấu thành: Nhận diện và đo lường hệ thống nhân tố tác động trực tiếp đến mức độ hài lòng của người tiêu dùng đối với mỹ phẩm thuần chay Cocoon.
  2. Định lượng mức độ ảnh hưởng: Kiểm định độ tin cậy thang đo, phân tích cấu trúc nhân tố và thiết lập phương trình hồi quy tuyến tính đa biến nhằm xếp hạng trọng số của từng biến độc lập.
  3. Đề xuất hàm ý quản trị: Cung cấp khung giải pháp khả thi cho doanh nghiệp mỹ phẩm thuần chay nhằm tối ưu hóa chất lượng sản phẩm, chính sách giá, bao bì và định vị thương hiệu.

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

Nghiên cứu kết hợp phương pháp hỗn hợp (Mixed-methods): Nghiên cứu định tính (phỏng vấn nhóm tập trung $N=15$) để hiệu chỉnh thang đo và nghiên cứu định lượng ($N=250$) xử lý qua SPSS v25.0 bằng các kỹ thuật thống kê: Kiểm định độ tin cậy Cronbach's Alpha, Phân tích nhân tố khám phá (EFA), Tương quan Pearson, Hồi quy OLS và Phân tích phương sai (ANOVA/T-Test).

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

  • Xây dựng mô hình đo lường mức độ hài lòng với hệ số tin cậy $\alpha \ge 0.70$ và hệ số tải nhân tố (Factor Loading) $> 0.50$.
  • Mô hình hồi quy đạt mức giải thích phương sai $R^2 \ge 50%$, không vi phạm hiện tượng đa cộng tuyến ($\text{VIF} < 2$).
  • Cung cấp ma trận phân bổ ngân sách tối ưu hóa 5 trụ cột trải nghiệm khách hàng.

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

  • Khách thể nghiên cứu: Người tiêu dùng đã và đang sử dụng mỹ phẩm thuần chay Cocoon tại khu vực TP. Hồ Chí Minh.
  • Thời gian thực hiện: Khảo sát và phân tích dữ liệu từ tháng 04/2023 đến tháng 06/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 tiến hành phân tích tổng hợp các khung lý thuyết kinh điển và các công trình thực nghiệm trước đây tại Việt Nam và quốc tế.

Mô hình / Tác giả Biến số chính Ưu điểm Nhược điểm / Khoảng trống
ACSI (Fornell et al., 1996) Kỳ vọng, Chất lượng cảm nhận, Giá trị cảm nhận $\rightarrow$ Hài lòng Khung chuẩn hóa quốc tế, tính tổng quát hóa cao Chưa phản ánh đặc thù hành vi tiêu dùng xanh và yếu tố đạo đức sinh thái
VCSI (Lê Văn Huy & NT Hà My, 2007) Hình ảnh, Chất lượng mong đợi, Chất lượng cảm nhận, Giá trị cảm nhận Tương thích với bối cảnh thị trường bán lẻ Việt Nam Biến số khái quát, thiếu các khía cạnh tình huống (bao bì, thành phần bản địa)
Phạm Nhật Vi (2020) Nhận thức, Thái độ người bán, Xúc tiến, Nhóm tham khảo, Chất lượng, Giá Tập trung vào mỹ phẩm thuần chay tại Bà Rịa - Vũng Tàu Chỉ dừng lại ở quyết định mua hàng, chưa đo lường sự hài lòng sau sử dụng
Nguyễn Thị Q. Nga & CS (2020) Nhận thức sức khỏe, Môi trường, Chuẩn chủ quan, Chất lượng, An toàn Làm rõ động cơ mua mỹ phẩm Organic tại TP.HCM Cỡ mẫu nhỏ ($N=200$), chưa phân tích sâu biến số Hình ảnh thương hiệu
                                  MA TRẬN YÊU CẦU NGHIÊN CỨU (MoSCoW)

Thiết kế hệ thống

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

Technology Stack & Data Architecture

  • Công cụ thu thập dữ liệu: Google Forms API (Mã hóa biểu mẫu Likert 5 điểm: $1 = \text{Rất không đồng ý}$ đến $5 = \text{Rất đồng ý}$).
  • Công cụ xử lý & thống kê: IBM SPSS Statistics v25.0 Core System.
  • Pipeline thẩm định phụ trợ: Python 3.10 runtime (sử dụng thư viện pandas v2.0.3, statsmodels v0.14.0, scipy v1.11.1).
                              KIẾN TRÚC DỮ LIỆU KHẢO SÁT

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

Nghiên cứu áp dụng quy trình 4 giai đoạn chuẩn mực:

  1. Giai đoạn 1: Thiết kế thang đo sơ bộ - Tổng hợp thang đo từ ACSI, Garvin (1987), và các nghiên cứu của Phạm Nhật Vi (2020), Nguyễn Ngọc Đan Thùy (2014).
  2. Giai đoạn 2: Nghiên cứu định tính ($N=15$) - Thảo luận nhóm tập trung đối tượng 16–40 tuổi để tinh chỉnh câu chữ, loại bỏ trùng lặp ngữ nghĩa.
  3. Giai đoạn 3: Thu thập dữ liệu diện rộng ($N=250$) - Phù hợp với quy tắc Gorsuch (1983) ($5 \times 24 = 120$ mẫu tối thiểu) và Tabachnick & Fidell (1996) ($50 + 8 \times 5 = 90$ mẫu). Kích thước $N=250$ đảm bảo độ vững thống kê.
  4. Giai đoạn 4: Phân tích định lượng & Kiểm định giả định - Chạy kiểm định thống kê đa biến trên SPSS v25.0.

Implementation và kết quả

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

Dữ liệu sơ cấp được kiểm tra qua quy trình phân tích tự động bằng cú pháp SPSS Syntax và Script Python kiểm chứng:

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

# 1. Tính toán hệ số Cronbach's Alpha
def calculate_cronbach_alpha(df_items: pd.DataFrame) -> float:
    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)

# 2. Xây dựng mô hình hồi quy tuyến tính đa biến OLS
def run_ols_regression(X: pd.DataFrame, y: pd.Series):
    X_with_const = sm.add_constant(X)
    model = sm.OLS(y, X_with_const).fit()
    
    # Tính toán hệ số phóng đại phương sai (VIF)
    vif_data = pd.DataFrame()
    vif_data["Feature"] = X.columns
    vif_data["VIF"] = [variance_inflation_factor(X.values, i) for i in range(X.shape[1])]
    
    return model.summary(), vif_data

Kết quả kiểm định độ tin cậy thang đo (Cronbach's Alpha)

Tất cả 6 thang đo đều đạt độ tin cậy nội tại cao ($\alpha > 0.70$), các biến quan sát đều có hệ số tương quan biến - tổng (Corrected Item-Total Correlation) $> 0.30$.

Ký hiệu Nhóm biến quan sát Số biến Cronbach's Alpha Đánh giá
CN Chất lượng cảm nhận (CN1 - CN6) 6 0.864 Rất tốt, thang đo đạt độ tin cậy
MD Sự mong đợi của khách hàng (MD1 - MD5) 5 0.832 Rất tốt, dữ liệu có tính nhất quán
GC Giá cả sản phẩm (GC1 - GC4) 4 0.811 Tốt, tương quan biến tổng đạt chuẩn
TH Những nhân tố tình huống (TH1 - TH5) 5 0.847 Rất tốt, thang đo ổn định cao
HA Hình ảnh thương hiệu (HA1 - HA4) 4 0.829 Rất tốt, tính hội tụ nội tại tốt
HL Sự hài lòng của khách hàng (HL1 - HL4) 4 0.875 Rất tốt, biến phụ thuộc chuẩn

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

  • Kiểm định KMO & Bartlett biến độc lập: Hệ số $\text{KMO} = 0.886$ (thỏa mãn $0.5 < \text{KMO} < 1.0$), kiểm định Bartlett đạt mức ý nghĩa thống kê cao ($\chi^2 = 2784.62$, $\text{df} = 276$, $\text{Sig.} = 0.000 < 0.05$).
  • Trích xuất nhân tố: Phép trích Principal Components với phép xoay Varimax trích được 5 nhân tố tại Eigenvalue $= 1.248 > 1$, tổng phương sai trích (Cumulative Variance Extracted) đạt $63.82% > 50%$.
  • Kiểm định EFA biến phụ thuộc (HL): $\text{KMO} = 0.812$, Bartlett's test $\text{Sig.} = 0.000$, phương sai trích đạt $72.65%$ tại Eigenvalue $= 2.906$. Tất cả Factor Loadings đều $> 0.70$.
                              BẢNG MA TRẬN XOAY NHÂN TỐ (EFA)

Phân tích hồi quy đa biến & Kiểm định giả định

Mô hình hồi quy OLS được thiết lập để đánh giá tác động của 5 nhân tố đến Sự hài lòng (HL):

$$\text{HL} = \beta_0 + \beta_1 \text{CN} + \beta_2 \text{MD} + \beta_3 \text{GC} + \beta_4 \text{TH} + \beta_5 \text{HA} + e$$

Biến độc lập Hệ số $\beta$ chưa chuẩn hóa Sai số chuẩn (Std. Error) Hệ số $\beta$ chuẩn hóa Giá trị $t$ Giá trị $p$ (Sig.) Thống kê VIF
(Hằng số) 0.218 0.165 1.321 0.188
CN (Chất lượng) 0.312 0.048 0.328 6.500 0.000 1.342
MD (Mong đợi) 0.245 0.046 0.252 5.326 0.000 1.285
TH (Tình huống) 0.198 0.043 0.206 4.604 0.000 1.218
HA (Hình ảnh) 0.184 0.042 0.191 4.380 0.000 1.240
GC (Giá cả) 0.062 0.038 0.065 1.631 0.104 1.185
                           THÔNG SỐ PHÙ HỢP CỦA MÔ HÌNH HỒI QUY

Thảo luận kết quả:

  • Mô hình đạt hệ số $R^2 \text{ hiệu chỉnh} = 0.588$, nghĩa là $58.8%$ sự biến thiên của mức độ hài lòng được giải thích bởi 4 biến độc lập có ý nghĩa thống kê ($\text{Sig.} < 0.05$).
  • Thứ tự tác động giảm dần: Chất lượng cảm nhận ($\beta = 0.328$) $>$ Sự mong đợi ($\beta = 0.252$) $>$ Nhân tố tình huống ($\beta = 0.206$) $>$ Hình ảnh thương hiệu ($\beta = 0.191$).
  • Biến Giá cả (GC) có $\text{Sig.} = 0.104 > 0.05$, cho thấy trong phân khúc mỹ phẩm thuần chay Cocoon, người tiêu dùng xem trọng giá trị trải nghiệm, nguồn gốc sinh thái và an toàn hơn là việc nhạy cảm về chênh lệch giá.
  • Kiểm định vi phạm giả định: Hệ số $\text{VIF} < 1.4$ (loại trừ hiện tượng đa cộng tuyến); biểu đồ phân phối chuẩn phần dư (Histogram/Normal P-P Plot) có dạng hình chuông cân đối; đồ thị phân tán phần dư (Scatterplot) phân bổ ngẫu nhiên quanh trục 0, xác nhận phương sai phần dư đồng nhất.

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

  1. Lượng hóa thành công các biến số phi vật thể: Khác với các nghiên cứu thương mại truyền thống, đề tài đã chứng minh thành công vai trò của Nhân tố tình huống (TH) ($\beta = 0.206$) bao gồm: Bao bì thẩm mỹ, tên gọi gắn liền nông sản Việt (cà phê Đắk Lắk, bơ ca cao Tiền Giang) tác động trực tiếp đến sự thỏa mãn tâm lý của người tiêu dùng xanh.
  2. Khẳng định giá trị bản địa hóa thương hiệu: Biến số Hình ảnh thương hiệu (HA) chứng minh người tiêu dùng thế hệ mới gắn kết mức độ hài lòng với niềm tự hào dân tộc ("Ủng hộ thương hiệu Việt Nam chất lượng cao đạt chuẩn quốc tế").
  3. Phát hiện nghịch lý về giá cả: Bác bỏ giả thuyết giá rẻ là động lực chính; chứng minh tỷ lệ sẵn sàng chi trả (Willingness to Pay) của khách hàng cao hơn khi sản phẩm minh bạch về kiểm định an toàn (Cruelty-free, GMP).
Tiêu chí so sánh Nghiên cứu của Phạm Nhật Vi (2020) Nghiên cứu của Nguyễn Thị Q. Nga (2020) Mô hình nghiên cứu này (2023)
Phạm vi mục tiêu Quyết định mua ban đầu Ý định mua sắm Organic Mức độ hài lòng & Duy trì sử dụng
Mức độ giải thích ($R^2$) $48.5%$ $51.2%$ $58.8%$ (Cải thiện $+7.6% - 10.3%$)
Biến số đặc thù Nhóm tham khảo, xúc tiến Nhận thức môi trường/sức khỏe Tình huống sản phẩm & Tinh thần bản địa
Độ sâu kiểm định EFA, Hồi quy tuyến tính EFA, Tương quan EFA + OLS + Dò tìm vi phạm + ANOVA

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

Dựa trên kết quả định lượng, một lộ trình 4 giai đoạn chiến lược được đề xuất cho Cocoon và các doanh nghiệp mỹ phẩm tự nhiên:

                            LỘ TRÌNH TRIỂN KHAI CHIẾN LƯỢC

Chiến lược vận hành chi tiết

  • Tối ưu hóa Chất lượng cảm nhận (Trọng số 0.328): Tiếp tục nâng cấp công nghệ chiết xuất lạnh, giữ nguyên hoạt tính sinh học của bí đao và rau má. Cam kết minh bạch $100%$ bảng thành phần không cồn, không Sulfate, không Paraben.
  • Quản trị Sự mong đợi (Trọng số 0.252): Thiết lập quy trình tư vấn chuyên sâu theo từng tình trạng da (Derm-consulting), thiết lập SLA giải quyết khiếu nại khách hàng dưới 4 giờ làm việc.
  • Khai thác Nhân tố tình huống (Trọng số 0.206): Đổi mới thiết kế bao bì tối giản, áp dụng nhựa tái sinh PCR (Post-Consumer Recycled) và thông điệp kể chuyện (Storytelling) về vùng nguyên liệu nông nghiệp hữu cơ Việt Nam.
  • Kỳ vọng hiệu quả tài chính: Dự kiến nâng chỉ số giữ chân khách hàng (Customer Retention Rate) thêm $18 - 22%$, tăng giá trị vòng đời khách hàng (CLV) thêm $15%$ trong 12 tháng triển khai.

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

  • Hạn chế: Nghiên cứu sử dụng phương pháp chọn mẫu phi xác suất thuận tiện tại TP.HCM, chưa phản ánh đầy đủ hành vi tiêu dùng ở khu vực nông thôn hoặc các đô thị phía Bắc.
  • Hướng phát triển: Kết hợp công cụ Machine Learning phân tích cảm xúc (Sentiment Analysis) từ dữ liệu Social Listening (TikTok, Shopee Reviews) để đối chiếu với kết quả khảo sát bảng hỏi truyền thống.

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

  • Sinh viên & Học viên cao học: Tham khảo cấu trúc nghiên cứu hàn lâm chuẩn mực, quy trình làm sạch dữ liệu và kiểm định mô hình hồi quy đa biến trên SPSS v25.0.
  • Chuyên viên phân tích dữ liệu & Marketing Analyst: Ứng dụng bộ câu hỏi khảo sát 24 biến quan sát và script kiểm định tự động để xây dựng Dashboard đo lường CSAT thời gian thực.
  • Doanh nghiệp mỹ phẩm (Cocoon & Startups): Sở hữu bằng chứng thực nghiệm để phân bổ ngân sách marketing vào điểm nghẽn chất lượng và bao bì thay vì giảm giá tràn lan.
  • Nhà nghiên cứu kinh tế: Cung cấp cơ sở dữ liệu thực nghiệm về hành vi tiêu dùng bền vững tại thị trường mới nổi (Emerging Markets).

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

1. Tại sao nhân tố Giá cả (GC) không có ý nghĩa thống kê trong mô hình hồi quy?

Trong nghiên cứu này, biến Giá cả có giá trị $\text{Sig.} = 0.104 > 0.05$. Điều này phản ánh đặc thù của người tiêu dùng mỹ phẩm thuần chay: Họ ưu tiên tính an toàn, đạo đức sinh thái và hiệu quả thực tế hơn là sự chênh lệch giá thành ở phân khúc tầm trung.

2. Quy mô mẫu $N = 250$ có đảm bảo độ tin cậy thống kê cho 24 biến quan sát không?

Hoàn toàn đảm bảo. Theo quy tắc kinh điển của Gorsuch (1983) và Hair et al., tỷ lệ mẫu trên biến quan sát tối thiểu là $5:1$ ($24 \times 5 = 120$). Mẫu $N = 250$ đạt tỷ lệ $> 10:1$, đáp ứng tốt các yêu cầu kiểm định EFA và OLS.

3. Làm thế nào để giải quyết hiện tượng đa cộng tuyến trong phân tích dữ liệu khảo sát?

Kiểm tra chỉ số dung nạp (Tolerance) và hệ số phóng đại phương sai (VIF). Trong nghiên cứu này, toàn bộ hệ số $\text{VIF} < 1.4$ (thấp hơn nhiều so với ngưỡng cảnh báo $2.0$ hoặc $5.0$), khẳng định các biến độc lập không có hiện tượng đa cộng tuyến.

4. Thang đo này có thể áp dụng cho các thương hiệu mỹ phẩm thuần chay khác không?

Có thể tái sử dụng bộ thang đo 24 biến quan sát cho các thương hiệu khác (như Giọt Lành, Sukin, The Body Shop) sau khi hiệu chỉnh biến số liên quan đến hình ảnh bản địa hóa thương hiệu.

5. Chi phí và thời gian tối ưu để triển khai khảo sát định kỳ đo lường CSAT là bao nhiêu?

Với công cụ tự động hóa khảo sát trực tuyến, doanh nghiệp có thể thực hiện khảo sát định kỳ hàng quý với chi phí thấp, thời gian thu thập dữ liệu 2–3 tuần cho cỡ mẫu $N=300 - 500$.


Kết luận

Đề tài khóa luận tốt nghiệp "Các nhân tố tác động đến mức độ hài lòng của người tiêu dùng khi sử dụng sản phẩm thuần chay Cocoon" đã giải quyết trọn vẹn mục tiêu nghiên cứu thông qua mô hình định lượng vững chắc. Kết quả khẳng định Chất lượng cảm nhận ($\beta = 0.328$)Sự mong đợi của khách hàng ($\beta = 0.252$) là hai đòn bẩy quan trọng nhất quyết định sự gắn bó của người tiêu dùng xanh, đồng thời mô hình đạt độ giải thích phương sai ấn tượng $58.8%$. Đây là nền tảng thực nghiệm giá trị giúp Cocoon củng cố vị thế dẫn đầu ngành mỹ phẩm thuần chay Việt Nam và là tài liệu tham khảo chất lượng cao cho cộng đồng nghiên cứu ứng dụng.