LỜI CẢM ƠN

Nhóm thực hiện xin gửi lời cảm ơn sâu sắc đến giảng viên hướng dẫn ThS.Trần Mạnh Tường đã tận tình hướng dẫn và giải đáp những thắc mắc trong suốt quá trình thực hiện đề tài. Vì sự hạn chế về kiến thức cũng như kỹ năng nên đề tài còn những thiếu sót. Kính mong thầy đánh giá và đóng góp cải thiện để đề tài được hoàn thiện hơn trong tương lai.

PHẦN 1: PHÂN TÍCH HÀNH VI MUA SẮM CỦA KHÁCH HÀNG THEO ĐẶC ĐIỂM NHÂN KHẨU HỌC

1 . Giới thiệu và mã hóa dữ liệu

1.1 Giới thiệu dữ liệu

library(readxl)
## Warning: package 'readxl' was built under R version 4.4.3
library(data.table)
## Warning: package 'data.table' was built under R version 4.4.3
library(dplyr)
## Warning: package 'dplyr' was built under R version 4.4.3
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:data.table':
## 
##     between, first, last
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(tidyr) 
## Warning: package 'tidyr' was built under R version 4.4.3

Giải thích

Các gói thư viện được cài đặt nhằm hỗ trợ toàn bộ quy trình xử lý và trực quan hóa dữ liệu.

  • readxl dùng để nhập dữ liệu Excel vào R, đảm bảo giữ nguyên cấu trúc gốc.

  • data.table hỗ trợ thao tác dữ liệu lớn với tốc độ cao.

  • dplyr cung cấp cú pháp trực quan cho các bước biến đổi và tổng hợp dữ liệu.

  • tidyr giúp định dạng và chuẩn hóa dữ liệu theo chuẩn tidy data.

# 1. Đọc dữ liệu
file_path <- "customers_data.xlsx"
data <- read_excel(file_path)

# 2. Xem tổng quan kích thước dữ liệu
num_rows <- nrow(data)     # số quan sát
num_cols <- ncol(data)     # số biến
cat("Số quan sát:", num_rows, "\n")
## Số quan sát: 100000
cat("Số biến:", num_cols, "\n")
## Số biến: 12
# 3. Tên các biến trong bộ dữ liệu
names(data)
##  [1] "id"                 "age"                "gender"            
##  [4] "income"             "education"          "region"            
##  [7] "loyalty_status"     "purchase_frequency" "purchase_amount"   
## [10] "product_category"   "promotion_usage"    "satisfaction_score"

Giải thích:

Bộ dữ liệu được sử dụng trong nghiên cứu có tên là “Customer Purchases Behaviour Dataset”, được công bố trên nền tảng Kaggle bởi tác giả Sanyam Goyal. Đây là một bộ dữ liệu mô phỏng (simulated dataset) phản ánh hành vi mua sắm của khách hàng trong lĩnh vực thương mại điện tử, được thiết kế với mục tiêu phục vụ cho mục đích nghiên cứu, phân tích dữ liệu và xây dựng mô hình dự báo trong marketing.

Cụ thể, bộ dữ liệu gồm 100.000 quan sát (rows) và 12 biến (columns), mỗi quan sát tương ứng với một khách hàng giả định. Các biến bao phủ cả đặc điểm nhân khẩu học và hành vi mua hàng bao gồm: id, age, gender, income, education, region, loyalty_status, purchase_frequency, purchase_amount, product_category, promotion_usage, satisfaction_score.

# 4. Kiểm tra kiểu dữ liệu của từng biến
sapply(data, class)
##                 id                age             gender             income 
##          "numeric"          "numeric"        "character"          "numeric" 
##          education             region     loyalty_status purchase_frequency 
##        "character"        "character"        "character"        "character" 
##    purchase_amount   product_category    promotion_usage satisfaction_score 
##          "numeric"        "character"          "numeric"          "numeric"

Giải thích

Hàm sapply(data, class) liệt kê kiểu dữ liệu của từng cột. Kết quả cho thấy các biến định lượng (id, age, income, purchase_amount, promotion_usage, satisfaction_score) được lưu ở dạng số (numeric), trong khi các biến mô tả/danh mục (gender, education, region, loyalty_status, purchase_frequency, product_category) hiện dưới dạng ký tự (character), biến định tính được quy ước thành factor khi cần chạy mô hình phân loại hoặc sinh các báo cáo tần suất; biến định lượng giữ ở numeric để thực hiện các thống kê mô tả và hồi quy.

# 5. Kiểm tra giá trị thiếu
missing_total <- sum(is.na(data))
missing_by_var <- colSums(is.na(data))
cat("Tổng số giá trị thiếu:", missing_total, "\n")
## Tổng số giá trị thiếu: 0
missing_by_var
##                 id                age             gender             income 
##                  0                  0                  0                  0 
##          education             region     loyalty_status purchase_frequency 
##                  0                  0                  0                  0 
##    purchase_amount   product_category    promotion_usage satisfaction_score 
##                  0                  0                  0                  0

Giải thích

Kiểm tra missing với sum(is.na(data)) và colSums(is.na(data)) trả về tổng số giá trị thiếu bằng 0 và từng biến không có giá trị NA. Đảm bảo dữ liệu sạch và có thể xử lý

# 6. Kiểm tra các dòng bị trùng lặp 
dup_total <- sum(duplicated(data))
cat("Số quan sát trùng lặp:", dup_total, "\n")
## Số quan sát trùng lặp: 0

Giải thích

Kiểm tra trùng lặp với duplicated(data) cho kết quả 0 dòng trùng. Điều này khẳng định mỗi hàng đại diện cho một quan sát khách hàng duy nhất — tránh sai lệch do lặp dữ liệu khi tính tần suất, trung bình

# 7. Kiểm tra cấu trúc chi tiết của dữ liệu
str(data)   # Xem cấu trúc: tên biến, kiểu dữ liệu, ví dụ giá trị
## tibble [100,000 × 12] (S3: tbl_df/tbl/data.frame)
##  $ id                : num [1:100000] 1 2 3 4 5 6 7 8 9 10 ...
##  $ age               : num [1:100000] 27 29 37 30 31 38 32 24 27 28 ...
##  $ gender            : chr [1:100000] "Male" "Male" "Male" "Male" ...
##  $ income            : num [1:100000] 40682 15317 38849 11568 46952 ...
##  $ education         : chr [1:100000] "Bachelor" "Masters" "Bachelor" "HighSchool" ...
##  $ region            : chr [1:100000] "East" "West" "West" "South" ...
##  $ loyalty_status    : chr [1:100000] "Gold" "Regular" "Silver" "Regular" ...
##  $ purchase_frequency: chr [1:100000] "frequent" "rare" "rare" "frequent" ...
##  $ purchase_amount   : num [1:100000] 18249 4557 11822 4098 19685 ...
##  $ product_category  : chr [1:100000] "Books" "Clothing" "Clothing" "Food" ...
##  $ promotion_usage   : num [1:100000] 0 1 0 0 1 0 0 0 0 0 ...
##  $ satisfaction_score: num [1:100000] 6 6 6 7 5 5 7 5 5 6 ...

Giải thích

Hàm str(data) trình bày cấu trúc: số quan sát 100.000 và 12 biến, cùng kiểu dữ liệu và vài giá trị mẫu. Thông tin này giúp rà soát nhanh các vấn đề thường gặp (ví dụ biến số được đọc nhầm kiểu) và cung cấp ngữ cảnh (kích thước mẫu lớn — phù hợp cho kiểm định thống kê có công suất cao và phân tích phân nhóm).

# 9. Kiểm tra dữ liệu rỗng hoặc khoảng trắng trong biến ký tự
char_vars <- names(Filter(is.character, data))
for (v in char_vars) {
  empty_count <- sum(trimws(data[[v]]) == "", na.rm = TRUE)
  cat("Biến", v, "có", empty_count, "giá trị rỗng\n")
}
## Biến gender có 0 giá trị rỗng
## Biến education có 0 giá trị rỗng
## Biến region có 0 giá trị rỗng
## Biến loyalty_status có 0 giá trị rỗng
## Biến purchase_frequency có 0 giá trị rỗng
## Biến product_category có 0 giá trị rỗng

Giải thích

Kiểm tra khoảng trắng/chuỗi rỗng trên các biến character bằng trimws() cho thấy không có ô rỗng. Điều này đảm bảo các thao tác group_by, join hoặc phân loại theo chuỗi sẽ không bị ảnh hưởng bởi giá trị rỗng ngụy trang dưới dạng chuỗi chỉ chứa khoảng trắng.

# 10. Xem 10 dòng đầu của dữ liệu
head(data, 10)
## # A tibble: 10 × 12
##       id   age gender income education  region loyalty_status purchase_frequency
##    <dbl> <dbl> <chr>   <dbl> <chr>      <chr>  <chr>          <chr>             
##  1     1    27 Male    40682 Bachelor   East   Gold           frequent          
##  2     2    29 Male    15317 Masters    West   Regular        rare              
##  3     3    37 Male    38849 Bachelor   West   Silver         rare              
##  4     4    30 Male    11568 HighSchool South  Regular        frequent          
##  5     5    31 Female  46952 College    North  Regular        occasional        
##  6     6    38 Male     7347 Bachelor   South  Silver         occasional        
##  7     7    32 Female   8265 Bachelor   South  Silver         frequent          
##  8     8    24 Female  47773 HighSchool North  Regular        rare              
##  9     9    27 Male    19154 College    East   Regular        occasional        
## 10    10    28 Female  24666 HighSchool North  Regular        rare              
## # ℹ 4 more variables: purchase_amount <dbl>, product_category <chr>,
## #   promotion_usage <dbl>, satisfaction_score <dbl>

Giải thích

Hiển thị 10 dòng đầu giúp xác nhận trực quan: các biến đăng ký giá trị hợp lý (ví dụ age nằm trong khoảng 12–49; income, purchase_amount là các giá trị dương; promotion_usage ở dạng 0/1). Bước này giúp kiểm tra trước khi bước vào mã hóa và phân tích dữ liệu

1.2 Mã hóa dữ liệu

library(knitr)
## Warning: package 'knitr' was built under R version 4.4.3
library(kableExtra)
## Warning: package 'kableExtra' was built under R version 4.4.3
## 
## Attaching package: 'kableExtra'
## The following object is masked from 'package:dplyr':
## 
##     group_rows
#1.. Chuẩn hóa tên biến về dạng chuẩn
names(data) <- trimws(names(data))
names(data) <- gsub(" ", "_", names(data))
names(data) <- tolower(names(data))


#2. Loại bỏ trùng lặp
data <- distinct(data)

#3. Chuyển kiểu dữ liệu
data$region <- as.factor(data$region)
data$loyalty_status <- as.factor(data$loyalty_status)
data$promotion_usage <- as.factor(data$promotion_usage)
data$income <- as.numeric(data$income)
data$purchase_amount <- as.numeric(data$purchase_amount)
str(data)
## tibble [100,000 × 12] (S3: tbl_df/tbl/data.frame)
##  $ id                : num [1:100000] 1 2 3 4 5 6 7 8 9 10 ...
##  $ age               : num [1:100000] 27 29 37 30 31 38 32 24 27 28 ...
##  $ gender            : chr [1:100000] "Male" "Male" "Male" "Male" ...
##  $ income            : num [1:100000] 40682 15317 38849 11568 46952 ...
##  $ education         : chr [1:100000] "Bachelor" "Masters" "Bachelor" "HighSchool" ...
##  $ region            : Factor w/ 4 levels "East","North",..: 1 4 4 3 2 3 3 2 1 2 ...
##  $ loyalty_status    : Factor w/ 3 levels "Gold","Regular",..: 1 2 3 2 2 3 3 2 2 2 ...
##  $ purchase_frequency: chr [1:100000] "frequent" "rare" "rare" "frequent" ...
##  $ purchase_amount   : num [1:100000] 18249 4557 11822 4098 19685 ...
##  $ product_category  : chr [1:100000] "Books" "Clothing" "Clothing" "Food" ...
##  $ promotion_usage   : Factor w/ 2 levels "0","1": 1 2 1 1 2 1 1 1 1 1 ...
##  $ satisfaction_score: num [1:100000] 6 6 6 7 5 5 7 5 5 6 ...

Giải thích

Sau khi nhập dữ liệu, nhóm nghiên cứu tiến hành kiểm tra và chuyển đổi kiểu dữ liệu của các biến để đảm bảo tính phù hợp cho các bước phân tích thống kê và mô hình hóa. Kết quả cho thấy bộ dữ liệu bao gồm 100.000 quan sát và 12 biến, trong đó:

  • Các biến định danh và định lượng như id, age, income, purchase_amount, satisfaction_score được lưu dưới dạng integer hoặc numeric, thuận lợi cho việc tính toán thống kê và hồi quy.

  • Các biến phân loại như gender, region, loyalty_status và promotion_usage được chuyển sang kiểu factor, giúp mô hình nhận diện rõ các nhóm giá trị rời rạc.

  • Các biến mô tả như education, purchase_frequency và product_category giữ nguyên kiểu character, phù hợp cho các thao tác lọc, nhóm hoặc trực quan hóa sau này.

#4. Tạo bảng mô tả với 4 cột: Biến, Ý nghĩa, Kiểu dữ liệu, Loại biến
variable_description <- data.frame(
  "Biến" = c("id", "age", "gender", "income", "education", "region",
             "loyalty_status", "purchase_frequency", "purchase_amount", 
             "product_category", "promotion_usage", "satisfaction_score"),
  "Ý nghĩa" = c(
    "Mã định danh khách hàng",
    "Tuổi của khách hàng",
    "Giới tính của khách hàng (Nam/Nữ)",
    "Thu nhập hằng năm (USD)",
    "Trình độ học vấn (HighSchool/College/Bachelor/Masters)",
    "Khu vực sinh sống",
    "Mức độ trung thành (Bronze/Silver/Gold/Platinum)",
    "Tần suất mua hàng trong kỳ",
    "Số tiền chi tiêu trong một lần mua (USD)",
    "Loại sản phẩm đã mua",
    "Có sử dụng khuyến mãi hay không (Yes/No)",
    "Mức độ hài lòng (1–10)"
  ),
  "Kiểu dữ liệu" = c(
    "integer",  # id
    "integer",  # age
    "character",  # gender
    "integer",  # income
    "character",  # education
    "character",  # region
    "character",  # loyalty_status
    "character",  # purchase_frequency
    "integer",  # purchase_amount
    "character",  # product_category
    "integer",  # promotion_usage
    "integer"   # satisfaction_score
  ),
  "Loại biến" = c(
    "Định danh",           # id
    "Định lượng",  # age
    "Định tính",  # gender
    "Định lượng",   # income
    "Định tính",  # education
    "Định tính",  # region
    "Định tính",  # loyalty_status
    "Định tính",  # purchase_frequency
    "Định lượng",  # purchase_amount
    "Định tính",  # product_category
    "Định lượng",   # promotion_usage
    "Định lượng"   # satisfaction_score
  )
)

#5. Hiển thị bảng 
variable_description %>%
  kbl(caption = "Bảng 1. Mô tả các biến trong bộ dữ liệu") %>%
  kable_styling(
    full_width = FALSE,
    position = "center",
    bootstrap_options = c("striped", "hover", "condensed", "responsive")
  )
## Warning in gzfile(file, mode): cannot open compressed file
## 'C:/Users/Dell/AppData/Local/Temp/RtmpoFnISO\file29d8329b4f3a', probable reason
## 'No such file or directory'
Bảng 1. Mô tả các biến trong bộ dữ liệu
Biến Ý.nghĩa Kiểu.dữ.liệu Loại.biến
id Mã định danh khách hàng integer Định danh
age Tuổi của khách hàng integer Định lượng
gender Giới tính của khách hàng (Nam/Nữ) character Định tính
income Thu nhập hằng năm (USD) integer Định lượng
education Trình độ học vấn (HighSchool/College/Bachelor/Masters) character Định tính
region Khu vực sinh sống character Định tính
loyalty_status Mức độ trung thành (Bronze/Silver/Gold/Platinum) character Định tính
purchase_frequency Tần suất mua hàng trong kỳ character Định tính
purchase_amount Số tiền chi tiêu trong một lần mua (USD) integer Định lượng
product_category Loại sản phẩm đã mua character Định tính
promotion_usage Có sử dụng khuyến mãi hay không (Yes/No) integer Định lượng
satisfaction_score Mức độ hài lòng (1–10) integer Định lượng

Giải thích

Việc tạo bảng giúp kiểm tra các biến định lượng và định tính một cách trực quan hơn và bảng hiện thị như sau:

  • Biến định tính bao gồm 6 biến : gender, education, region, loyalty_status, product_category, promotion_usage, purchase_frequenc.

  • Biến định lượng bao gồm: age, income, purchase_amount, satisfaction_score

#6. Chuẩn hóa tên cột
names(data) <- tolower(names(data))
#7. Loại bỏ giá trị âm hoặc 0 trong purchase_amount
data <- filter(data, purchase_amount > 0)
# loại bỏ giá trị NA

#8.  Mã hóa promotion_usage thành nhãn
data$promotion_usage <- factor(data$promotion_usage,
                               levels = c(0, 1),
                               labels = c("No", "Yes"))
#9. Tạo biến mới: thu nhập trung bình theo khu vực
region_income <- data %>%
  group_by(region) %>%
  summarise(mean_income = mean(income, na.rm = TRUE))
#10. Chuẩn hóa income
data$income_scaled <- scale(data$income)


#11.  Tạo biến nhóm tuổi (AgeGroup) từ biến age
data <- data %>%
  mutate(
    AgeGroup = case_when(
      age >= 12 & age <= 17 ~ "12-17",
      age >= 18 & age <= 24 ~ "18-24",
      age >= 25 & age <= 34 ~ "25-34",
      age >= 35 & age <= 49 ~ "35-49",
      age > 49 ~ "Trên 49",
      age < 12 ~ "Dưới 12",
      TRUE ~ NA_character_     # trường hợp còn lại
    )
  )


# 12.  #Loại bỏ NA trong nhóm tuổi
data_age <- data %>%
  filter(!is.na(AgeGroup), AgeGroup != "NA", AgeGroup != "")

Giải thích

Các bước xử lý được thực hiện nhằm chuẩn hóa và làm sạch dữ liệu, đảm bảo tính chính xác, nhất quán và sẵn sàng cho phân tích.

  • Đầu tiên, việc chuyển toàn bộ tên cột về chữ thường giúp thống nhất định dạng và hạn chế lỗi khi thao tác với dữ liệu.

  • Tiếp đó, các quan sát có giá trị purchase_amount bằng hoặc nhỏ hơn 0 được loại bỏ để đảm bảo tính hợp lý của biến chi tiêu, đồng thời mã hóa biến promotion_usage thành dạng nhãn (“No”, “Yes”) giúp dễ dàng diễn giải trong thống kê và mô hình hóa.

  • Việc tính thu nhập trung bình theo từng khu vực (region_income) cho phép nhận diện sự khác biệt về mức sống giữa các vùng, hỗ trợ cho các phân tích so sánh và chiến lược thị trường.

  • Biến income được chuẩn hóa bằng hàm scale() để đưa về cùng thang đo, phục vụ cho các mô hình thống kê đa biến.

  • Biến age được phân nhóm thành AgeGroup giúp mô tả hành vi tiêu dùng theo độ tuổi — một yếu tố nhân khẩu học quan trọng trong nghiên cứu hành vi khách hàng.

  • Cuối cùng, việc loại bỏ các giá trị NA trong AgeGroup hoàn thiện quá trình làm sạch, giúp dữ liệu đạt tính toàn vẹn và đáng tin cậy.

Nhìn chung, chuỗi thao tác này giúp chuẩn hóa cấu trúc, xử lý bất thường và tạo ra các biến có ý nghĩa thống kê – thực tiễn cao, từ đó hình thành nền tảng vững chắc cho các phân tích mô tả, hồi quy và dự báo ở các bước tiếp theo.

2 . Phân tích dữ liệu

2.1 Các thống kê cơ bản

# 2. Tóm tắt thống kê toàn bộ biến
summary(data)
##        id              age        gender              income     
##  Min.   :     1   Min.   :12   Length:100000      Min.   : 5000  
##  1st Qu.: 25001   1st Qu.:27   Class :character   1st Qu.:16272  
##  Median : 50001   Median :30   Mode  :character   Median :27585  
##  Mean   : 50001   Mean   :30                      Mean   :27516  
##  3rd Qu.: 75000   3rd Qu.:33                      3rd Qu.:38747  
##  Max.   :100000   Max.   :49                      Max.   :50000  
##   education           region      loyalty_status  purchase_frequency
##  Length:100000      East :30074   Gold   : 9898   Length:100000     
##  Class :character   North:19918   Regular:60138   Class :character  
##  Mode  :character   South:20073   Silver :29964   Mode  :character  
##                     West :29935                                     
##                                                                     
##                                                                     
##  purchase_amount product_category   promotion_usage satisfaction_score
##  Min.   : 1118   Length:100000      No :69920       Min.   : 0.00     
##  1st Qu.: 5583   Class :character   Yes:30080       1st Qu.: 4.00     
##  Median : 9452   Mode  :character                   Median : 5.00     
##  Mean   : 9635                                      Mean   : 5.01     
##  3rd Qu.:13350                                      3rd Qu.: 6.00     
##  Max.   :26204                                      Max.   :10.00     
##    income_scaled.V1     AgeGroup        
##  Min.   :-1.7324495   Length:100000     
##  1st Qu.:-0.8651772   Class :character  
##  Median : 0.0052498   Mode  :character  
##  Mean   : 0.0000000                     
##  3rd Qu.: 0.8641354                     
##  Max.   : 1.7299459

Giải thích

Kết quả thống kê mô tả cho thấy bộ dữ liệu gồm 100.000 quan sát với đặc điểm nhân khẩu và hành vi tiêu dùng đa dạng. Biến age dao động từ 12–49, trung bình 30 tuổi, phản ánh nhóm khách hàng chủ yếu là người trẻ. Thu nhập (income) trung bình 27.516, trải rộng từ 5.000–50.000, cho thấy sự khác biệt đáng kể giữa các cá nhân. Chi tiêu (purchase_amount) trung bình 9.635, gần với trung vị 9.452, thể hiện phân phối chi tiêu tương đối cân đối.

Về biến định tính, miền East chiếm tỷ lệ lớn nhất trong region, trong khi nhóm khách hàng “Regular” chiếm ưu thế trong loyalty_status. Khoảng 30% khách hàng có sử dụng khuyến mãi (promotion_usage), cho thấy hoạt động marketing có mức tiếp cận khá tốt. Điểm hài lòng trung bình 5.01 phản ánh mức độ hài lòng trung bình. Biến income_scaled được chuẩn hóa chính xác, và AgeGroup được hình thành đầy đủ, tạo nền tảng cho các phân tích so sánh nhóm ở các bước sau.

# 3. Tính trung bình thu nhập khách hàng 
mean(data$income, na.rm = TRUE)
## [1] 27516.27

Giải thích

Giá trị thu nhập trung bình của khách hàng là 27.516,27, phản ánh mức thu nhập bình quân trong mẫu quan sát. Mức này nằm trong khoảng giữa của dải thu nhập (từ 5.000 đến 50.000), cho thấy dữ liệu có sự phân bố tương đối cân bằng, không bị lệch mạnh về phía các giá trị cực đoan.

# 4. Tính trung vị thu nhập 
median(data$income, na.rm = TRUE)
## [1] 27584.5

Giải thích

Giá trị trung vị của thu nhập là 27.584,5, gần tương đương với giá trị trung bình (27.516,27). Điều này cho thấy phân phối thu nhập của khách hàng tương đối đối xứng, không xuất hiện hiện tượng lệch mạn.

# 5. Tính độ lệch chuẩn thu nhập 
sd(data$income, na.rm = TRUE)
## [1] 12996.78

Giải thích

Độ lệch chuẩn của thu nhập là 12.996,78, cho thấy mức độ phân tán khá lớn quanh giá trị trung bình 27.516,27. Điều này phản ánh sự chênh lệch đáng kể giữa các nhóm khách hàng về khả năng thu nhập, tức là trong tập dữ liệu có cả nhóm thu nhập thấp và nhóm thu nhập cao cùng tồn tại.

# 6. Tính phương sai thu nhập
var(data$income, na.rm = TRUE)
## [1] 168916358

Giải thích

Phương sai của thu nhập là 168.916.358, thể hiện mức độ biến thiên rất lớn quanh giá trị trung bình. Do phương sai là bình phương của độ lệch chuẩn, kết quả cao cho thấy dữ liệu thu nhập có sự phân tán rộng, phản ánh sự khác biệt rõ rệt giữa các cá nhân trong mẫu khảo sát.

# 8. Tính tứ phân vị thu nhập 
quantile(data$income, probs = c(0.25, 0.5, 0.75), na.rm = TRUE)
##      25%      50%      75% 
## 16271.75 27584.50 38747.25

Giải thích

Kết quả tứ phân vị cho thấy: Q1 = 16.271,75, Q2 (trung vị) = 27.584,50, và Q3 = 38.747,25. Khoảng giữa của phân phối (Q1–Q3) rộng khoảng 22.475 đơn vị, phản ánh sự chênh lệch đáng kể giữa nhóm thu nhập thấp và cao. Phần lớn khách hàng có thu nhập tập trung trong khoảng 16.000–39.000, tức là thuộc nhóm thu nhập trung bình, trong khi các giá trị ngoài khoảng này có thể được xem là nhóm ngoại biên.

# 9. Tính hệ số biến thiên (CV) của thu nhập
cv_income <- sd(data$income, na.rm = TRUE) / mean(data$income, na.rm = TRUE)
cv_income
## [1] 0.4723308

Giải thích

Hệ số biến thiên (CV) của thu nhập là 0,4723, tương đương 47,23%, thể hiện mức độ dao động tương đối cao của thu nhập so với giá trị trung bình. Chỉ số này cho thấy thu nhập của khách hàng có sự phân tán đáng kể — mức thu nhập giữa các cá nhân chênh lệch gần một nửa so với trung bình chung. Về mặt học thuật, CV > 0,3 thường được xem là mức biến động lớn, hàm ý phân phối thu nhập không đồng đều.

# 10. Tần suất từng nhóm giới tính 
table(data$gender)
## 
## Female   Male 
##  50074  49926

Giải thích

Kết quả cho thấy bộ dữ liệu gồm 50.074 khách hàng nữ (50,07%) và 49.926 khách hàng nam (49,93%), thể hiện sự phân bố giới tính gần như cân bằng. Về mặt thống kê, tỷ lệ này đảm bảo tính đại diện cao và hạn chế sai lệch giới trong mẫu dữ liệu.

# 12. Thống kê theo vùng (Region) 
table(data$region)
## 
##  East North South  West 
## 30074 19918 20073 29935

Giải thích

Kết quả cho thấy sự phân bố khách hàng theo khu vực tương đối đồng đều: East (30,074 customers), West (29,935), South (20,073) và North (19,918). Khu vực East chiếm tỷ lệ cao nhất (≈30%), trong khi North có số lượng khách hàng thấp hơn đôi chút. Sự phân bố này phản ánh phạm vi bao phủ thị trường rộng và cân đối giữa các vùng, giúp dữ liệu đảm bảo tính đại diện và hạn chế sai lệch địa lý.

# 13. Trung bình thu nhập theo giới tính
aggregate(income ~ gender, data = data, FUN = mean)
##   gender   income
## 1 Female 27507.95
## 2   Male 27524.62

Giải thích

Kết quả cho thấy thu nhập trung bình của nam (27.524,62) và nữ (27.507,95) gần như tương đương nhau, chênh lệch không đáng kể. Về mặt thống kê, điều này cho thấy phân bố thu nhập giữa hai giới tính tương đối cân bằng, không có dấu hiệu sai lệch lớn hay bất bình đẳng thu nhập rõ rệt trong mẫu dữ liệu.

# 14. Trung bình thu nhập theo độ tuổi (Age Group)
aggregate(income ~ AgeGroup, data = data, FUN = mean)
##   AgeGroup   income
## 1    12-17 28966.86
## 2    18-24 27519.91
## 3    25-34 27517.88
## 4    35-49 27481.73

Giải thích

Kết quả cho thấy thu nhập trung bình giữa các nhóm tuổi có sự khác biệt nhẹ: nhóm 12–17 tuổi có mức trung bình cao nhất (28.966,86), tiếp đến là các nhóm 18–24 (27.519,91), 25–34 (27.517,88) và 35–49 (27.481,73). Mức chênh lệch giữa các nhóm tương đối nhỏ, cho thấy thu nhập phân bố khá đồng đều theo độ tuổi.

# 15. Trung bình chi tiêu theo giới tính 
aggregate(purchase_amount ~ gender, data = data, FUN = mean)
##   gender purchase_amount
## 1 Female        9634.405
## 2   Male        9635.178

Giải thích

Kết quả cho thấy chi tiêu trung bình của nữ (9.634,41) và nam (9.635,18) gần như tương đương, với mức chênh lệch không đáng kể. Điều này cho thấy giới tính không phải là yếu tố ảnh hưởng rõ rệt đến mức chi tiêu trung bình trong mẫu dữ liệu.

# 17. Kiểm tra phân phối chuẩn của thu nhập
shapiro.test(sample(data$income, 500))  # Lấy mẫu 500 để test
## 
##  Shapiro-Wilk normality test
## 
## data:  sample(data$income, 500)
## W = 0.93796, p-value = 1.4e-13

Giải thích

Kết quả kiểm định Shapiro–Wilk cho thấy W = 0,9609 và p-value = 2,947×10⁻¹⁰, nhỏ hơn 0,05, do đó bác bỏ giả thuyết H0 về phân phối chuẩn. Điều này nghĩa là thu nhập của khách hàng không tuân theo phân phối chuẩn, mà có thể bị lệch hoặc có đuôi dài do sự chênh lệch giữa các nhóm thu nhập.

# 18. Thống kê tần suất khách hàng theo nhóm tuổi 
table(data$AgeGroup)
## 
## 12-17 18-24 25-34 35-49 
##   268 10622 73316 15794

Giải thích

Kết quả thống kê cho thấy phân bố khách hàng theo nhóm tuổi không đồng đều: nhóm 25–34 chiếm tỷ lệ áp đảo với 73.316 người, tiếp đến là 35–49 (15.794), 18–24 (10.622), và 12–17 (268). Sự chênh lệch lớn này cho thấy mẫu dữ liệu tập trung chủ yếu vào nhóm lao động trẻ và trung niên – nhóm có sức mua và tần suất tiêu dùng cao nhất, phản ánh cấu trúc khách hàng thực tế của nhiều doanh nghiệp hiện nay.

# 19. Tính tổng chi tiêu trung bình của toàn bộ khách hàng-
mean(data$purchase_amount, na.rm = TRUE)
## [1] 9634.791

Giải thích

Giá trị chi tiêu trung bình của toàn bộ khách hàng là 9.634,79, cho thấy mức chi tiêu ổn định quanh mốc 10.000. Sự ổn định này gợi ý rằng hành vi mua sắm của khách hàng trong mẫu có tính đồng nhất tương đối, không bị chi phối mạnh bởi các ngoại lệ.

# 20. So sánh chi tiêu trung bình giữa các vùng 
aggregate(purchase_amount ~ region, data = data, FUN = mean)
##   region purchase_amount
## 1   East        9615.411
## 2  North        9673.362
## 3  South        9648.787
## 4   West        9619.211

Giải thích

Chi tiêu trung bình giữa các khu vực khá đồng đều: North (9.673,36), South (9.648,79), West (9.619,21) và East (9.615,41). Mức chênh lệch nhỏ cho thấy hành vi chi tiêu của khách hàng tương đối thống nhất trên các vùng địa lý, không có khu vực nào chi tiêu vượt trội rõ rệt.

# 21. Đếm số lượng khách hàng mỗi nhóm thu nhập
table(cut(data$income, breaks = 5))
## 
## (4.95e+03,1.4e+04]  (1.4e+04,2.3e+04]  (2.3e+04,3.2e+04]  (3.2e+04,4.1e+04] 
##              20061              19816              20007              20206 
##    (4.1e+04,5e+04] 
##              19910

Giải thích

Dữ liệu được chia thành 5 nhóm thu nhập có số lượng khách hàng tương đối đồng đều, dao động quanh 19.000–20.000 người mỗi nhóm. Phân bố này cho thấy mẫu dữ liệu được xây dựng cân bằng giữa các mức thu nhập, không bị thiên lệch mạnh về nhóm thấp hay cao.

# 22. Phân phối thu nhập theo giới tính
aggregate(income ~ gender, data, summary)
##   gender income.Min. income.1st Qu. income.Median income.Mean income.3rd Qu.
## 1 Female     5000.00       16271.50      27554.50    27507.95       38718.75
## 2   Male     5000.00       16272.50      27615.50    27524.62       38786.75
##   income.Max.
## 1    49998.00
## 2    50000.00

Giải thích

Phân phối thu nhập giữa hai giới gần như tương đồng: thu nhập trung bình của nữ là 27.507,95 và nam là 27.524,62, với các giá trị tứ phân vị, trung vị và cực trị hầu như trùng nhau. Điều này cho thấy không có sự chênh lệch đáng kể về thu nhập giữa nam và nữ, phản ánh tính cân bằng giới trong cơ cấu thu nhập của khách hàng.

# 23. Tính IQR (khoảng tứ phân vị)
IQR(data$income, na.rm = TRUE)
## [1] 22475.5

Giải thích

Khoảng tứ phân vị (IQR) của thu nhập là 22.475,5, thể hiện mức độ phân tán của nhóm 50% khách hàng trung tâm. Giá trị này tương đối lớn, phản ánh sự khác biệt đáng kể giữa nhóm thu nhập thấp và cao trong phần giữa của phân phối.

# 25. Đếm số khách hàng có chi tiêu > trung bình
mean(data$purchase_amount > mean(data$purchase_amount, na.rm = TRUE))
## [1] 0.4882

Giải thích

Kết quả cho thấy 48,82% khách hàng có mức chi tiêu cao hơn trung bình, trong khi 51,18% chi tiêu thấp hơn mức này. Sự phân bố gần cân bằng này chứng tỏ dữ liệu chi tiêu có cấu trúc hợp lý, không bị lệch mạnh về nhóm tiêu dùng thấp hoặc cao.

2.2 Trực quan hóa dữ liệu

# === 1. Cài và gọi thư viện ===
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.4.3
library(dplyr)
library(scales)
## Warning: package 'scales' was built under R version 4.4.3
library(RColorBrewer)
library(corrplot)
## Warning: package 'corrplot' was built under R version 4.4.3
## corrplot 0.95 loaded
library(GGally)
## Warning: package 'GGally' was built under R version 4.4.3
library(reshape2)
## Warning: package 'reshape2' was built under R version 4.4.3
## 
## Attaching package: 'reshape2'
## The following object is masked from 'package:tidyr':
## 
##     smiths
## The following objects are masked from 'package:data.table':
## 
##     dcast, melt
library(viridis)
## Warning: package 'viridis' was built under R version 4.4.3
## Loading required package: viridisLite
## Warning: package 'viridisLite' was built under R version 4.4.3
## 
## Attaching package: 'viridis'
## The following object is masked from 'package:scales':
## 
##     viridis_pal
library(ggthemes)
## Warning: package 'ggthemes' was built under R version 4.4.3
library(ggridges)
## Warning: package 'ggridges' was built under R version 4.4.3

Giải thích

Các gói thư viện được cài đặt nhằm hỗ trợ toàn bộ quy trình xử lý và trực quan hóa dữ liệu.

  • ggplot2: Tạo biểu đồ trực quan, chuyên nghiệp theo Grammar of Graphics.

  • dplyr: Hỗ trợ xử lý dữ liệu nhanh chóng (lọc, nhóm, tính toán).

  • scales: Định dạng trục, tỷ lệ và đơn vị hiển thị trong biểu đồ.

  • RColorBrewer: Cung cấp bảng màu khoa học, giúp biểu đồ dễ phân biệt nhóm.

  • corrplot: Trực quan hóa ma trận tương quan giữa các biến.

  • GGally: Mở rộng ggplot2, tạo ma trận biểu đồ (scatterplot matrix).

  • reshape2: Chuyển đổi cấu trúc dữ liệu giữa dạng rộng và dài.

  • viridis: Cung cấp thang màu liên tục.

  • ggthemes: Thêm chủ đề và phong cách thẩm mỹ cho biểu đồ ggplot2.

  • ggridges: Vẽ biểu đồ mật độ chồng lớp để so sánh phân phối nhiều nhóm.

# === 2. Đọc dữ liệu ===
file_path <- "customers_data.xlsx"
data <- read_excel(file_path)


# === 4. 20 BIỂU ĐỒ TRỰC QUAN HÓA ===

## 1. Phân bố độ tuổi theo giới tính
ggplot(data, aes(x = age, fill = gender)) +
  geom_histogram(position = "dodge", bins = 30) +
  labs(title = "Phân bố độ tuổi theo giới tính",
       x = "Tuổi", y = "Tần suất", fill = "Giới tính") +
  theme_minimal() +
  scale_fill_brewer(palette = "Paired")

Giải thích

Biểu đồ phân bố độ tuổi theo giới tính cho thấy cấu trúc dân số mẫu có dạng gần chuẩn, tập trung mạnh trong khoảng 25–35 tuổi, tức là nhóm lao động trẻ và trung niên chiếm ưu thế rõ rệt. Độ tuổi có tần suất cao nhất rơi vào khoảng 28–30, với tần suất gần 8.000 quan sát ở nữ và khoảng 7.500 ở nam. Hai đường phân bố gần như trùng nhau, thể hiện sự cân bằng giới tính trong từng nhóm tuổi và cho thấy dữ liệu không bị thiên lệch về giới.

## 2. Chi tiêu trung bình theo nhóm tuổi và giới tính
ggplot(data_age %>% filter(!is.na(AgeGroup)), 
       aes(x = AgeGroup, y = purchase_amount, fill = gender)) +
  geom_bar(stat = "summary", fun = "mean", position = "dodge") +
  labs(x = "Nhóm tuổi", y = "Chi tiêu trung bình (VND)", 
       title = "Chi tiêu trung bình theo nhóm tuổi và giới tính") +
  theme_minimal()

Giải thích

Ở nhóm tuổi 12–17, mức chi tiêu trung bình cao nhất, đạt khoảng 10.031 VND đối với nữ và 9.960 VND đối với nam, cho thấy nhóm khách hàng trẻ có xu hướng chi tiêu cao hơn so với các nhóm khác, có thể do tác động của yếu tố tiêu dùng theo xu hướng hoặc ảnh hưởng xã hội.

Các nhóm tuổi còn lại có mức chi tiêu khá đồng đều, dao động quanh 9.630–9.660 VND. Cụ thể, ở nhóm 18–24, chi tiêu trung bình đạt 9.629,7 VND (nữ) và 9.655,8 VND (nam); nhóm 25–34 là 9.629,9 VND (nữ) và 9.638,8 VND (nam); nhóm 35–49 đạt 9.656,7 VND (nữ) và 9.598,8 VND (nam). Mặc dù sự khác biệt giữa các giới là không đáng kể (dưới 1%), xu hướng chung cho thấy chi tiêu có xu hướng giảm nhẹ theo độ tuổi.

## 3. Mối quan hệ giữa tuổi, thu nhập và chi tiêu
ggplot(data, aes(x = age, y = purchase_amount, color = income)) +
  geom_point(alpha = 0.7) +
  labs(title = "Mối quan hệ giữa độ tuổi, thu nhập và chi tiêu",
       x = "Tuổi", y = "Chi tiêu (VNĐ)", color = "Thu nhập") +
  theme_light()

Giải thích

Quan sát phân bố cho thấy phần lớn cá nhân có mức chi tiêu tập trung quanh 9.000–15.000 VND, tương ứng với thu nhập từ 20.000–40.000 VND. Nhóm tuổi 25–35 có mật độ điểm cao nhất, thể hiện mức chi tiêu ổn định và có xu hướng tăng nhẹ cùng thu nhập. Trong khi đó, các nhóm ngoài khoảng này (dưới 20 và trên 40 tuổi) có mức chi tiêu và thu nhập biến động thấp hơn, phản ánh khả năng tài chính và mức tiêu dùng hạn chế hơn.

Có thể thấy không tồn tại mối quan hệ tuyến tính rõ rệt giữa tuổi và chi tiêu, nhưng mức chi tiêu cao hơn có xu hướng xuất hiện ở nhóm thu nhập trung bình–cao (màu xanh sáng).

Điều này cho thấy thu nhập vẫn là yếu tố ảnh hưởng chính đến chi tiêu, trong khi độ tuổi chỉ đóng vai trò hỗ trợ, phản ánh đặc điểm hành vi tiêu dùng phổ biến của nhóm khách hàng đang ở giai đoạn ổn định tài chính (25–35 tuổi).

## 4. Biểu đồ phân tán 3 chiều (tuổi – thu nhập – chi tiêu) với giới tính
ggplot(data, aes(x = age, y = income, size = purchase_amount, color = gender)) +
  geom_point(alpha = 0.6) +
  labs(title = "Mối quan hệ tuổi, thu nhập và chi tiêu theo giới tính",
       x = "Tuổi", y = "Thu nhập (VNĐ)", size = "Chi tiêu", color = "Giới tính") +
  theme_minimal()

Giải thích

Phân bố cho thấy phần lớn dữ liệu tập trung ở nhóm tuổi 25–35 với thu nhập phổ biến trong khoảng 20.000–40.000 VND. Các bong bóng lớn – tương ứng với chi tiêu cao trên 20.000 VND – chủ yếu xuất hiện ở nhóm thu nhập trên 35.000 VND, phản ánh mối quan hệ tỷ lệ thuận giữa thu nhập và chi tiêu. Sự khác biệt giới tính là không đáng kể: cả nam và nữ đều có xu hướng chi tiêu tương tự trong cùng mức thu nhập. Tuy nhiên, ở nhóm thu nhập cao (trên 45.000 VND), bong bóng của nữ có phần dày hơn, cho thấy nữ giới có xu hướng chi tiêu tương đối cao hơn so với nam trong phân khúc thu nhập cao.

## 5. Hộp số chi tiêu theo giới tính
ggplot(data, aes(x = gender, y = purchase_amount, fill = gender)) +
  geom_boxplot() +
  labs(title = "Phân bố chi tiêu theo giới tính",
       x = "Giới tính", y = "Chi tiêu (VNĐ)") +
  theme_bw() +
  scale_fill_brewer(palette = "Set3")

Giải thích

Biểu đồ hộp thể hiện phân bố mức chi tiêu (VND) của khách hàng theo giới tính. Cả hai nhóm nam và nữ đều có phân bố khá tương đồng, với trung vị chi tiêu khoảng 9.500–10.000 VND. Khoảng tứ phân vị (IQR) trải từ khoảng 5.000 đến 14.000 VND, cho thấy phần lớn khách hàng chi tiêu trong mức trung bình này. Một vài điểm dữ liệu nằm ngoài hộp biểu diễn các giá trị ngoại lai (outliers) – những khách hàng chi tiêu vượt trội, trên 20.000 VND.

## 6. Chi tiêu trung bình theo trình độ học vấn và giới tính
ggplot(data, aes(x = education, y = purchase_amount, fill = gender)) +
  geom_bar(stat = "summary", fun = "mean", position = "dodge") +
  labs(title = "Chi tiêu trung bình theo trình độ học vấn và giới tính",
       x = "học vấn", y = "Chi tiêu trung bình (VNĐ)") +
  theme_minimal() +
  coord_flip()

Giải thích

Quan sát dữ liệu cho thấy mức chi tiêu trung bình giữa các nhóm học vấn không có sự chênh lệch đáng kể, dao động quanh mức 9.600–10.000 VND. Cụ thể, nhóm có trình độ Thạc sĩ (Masters) chi tiêu cao nhất, đạt khoảng 10.000 VND, tiếp đến là nhóm Cử nhân (Bachelor) và Cao đẳng (College) với mức trung bình xấp xỉ 9.800–9.900 VND, trong khi nhóm Trung học (HighSchool) có chi tiêu thấp hơn nhẹ, khoảng 9.600 VND. Giới tính không tạo ra khác biệt rõ rệt: chi tiêu trung bình của nam và nữ gần như tương đương trong mọi cấp học, sai lệch chỉ ở mức dưới 1%.

Xu hướng này cho thấy rằng trình độ học vấn có mối quan hệ dương yếu với mức chi tiêu, tức là những cá nhân có học vấn cao hơn có xu hướng chi tiêu nhiều hơn, có thể do họ sở hữu thu nhập ổn định và khả năng tài chính lớn hơn.

## 7. Biểu đồ tương quan giữa các biến số
num_data <- data %>% select_if(is.numeric)
corrplot(cor(num_data, use = "pairwise.complete.obs"),
         method = "color", tl.col = "black",
         title = "Ma trận tương quan giữa các biến định lượng", mar = c(0,0,2,0))

Giải thích

Ma trận tương quan cho thấy mối quan hệ nổi bật giữa thu nhập (income) và chi tiêu (purchase_amount) với hệ số tương quan Pearson r ≈ 0.95 (r² ≈ 0.90), chứng tỏ thu nhập là yếu tố giải thích chính cho biến động trong chi tiêu của khách hàng. Ngược lại, các biến khác như độ tuổi (age), tần suất sử dụng khuyến mãi (promotion_usage) hay điểm hài lòng (satisfaction_score) đều có hệ số tương quan rất thấp (chỉ từ 0.00 đến 0.02), cho thấy không có mối quan hệ tuyến tính đáng kể với chi tiêu hoặc thu nhập.

Kết quả này cho thấy chi tiêu của khách hàng phụ thuộc mạnh vào khả năng tài chính hơn là các yếu tố hành vi hoặc nhân khẩu học khác.

## 8. Biểu đồ heatmap theo nhóm tuổi, giới tính và chi tiêu trung bình

ggplot(data_age, aes(x = AgeGroup, y = purchase_amount, fill = gender)) +
  stat_summary(fun = mean, geom = "bar", position = "dodge", alpha = 0.9) +
  geom_text(stat = "summary", fun = mean, aes(label = round(..y.., 1)), 
            position = position_dodge(0.9), vjust = -0.3, size = 3) +
  scale_fill_viridis_d(option = "plasma") +
  theme_minimal() +
  labs(title = "Chi tiêu trung bình theo nhóm tuổi và giới tính",
       x = "Nhóm tuổi", y = "Chi tiêu trung bình", fill = "Giới tính")
## Warning: The dot-dot notation (`..y..`) was deprecated in ggplot2 3.4.0.
## ℹ Please use `after_stat(y)` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.

Giải thích

Biểu đồ thể hiện chi tiêu trung bình (VND) của khách hàng theo nhóm tuổi và giới tính cho thấy mức chi tiêu giữa hai giới gần như tương đồng trong mọi nhóm tuổi, với sai lệch rất nhỏ (dưới 1%). Ở nhóm 12–17 tuổi, chi tiêu trung bình cao nhất, đạt 10.031 VND đối với nữ và 9.960 VND đối với nam, chênh lệch chỉ khoảng 71 VND. Từ nhóm 18–24 trở đi, mức chi tiêu giảm nhẹ và ổn định quanh 9.630–9.660 VND cho cả hai giới (ví dụ: 18–24 tuổi nữ 9.629, nam 9.655; 35–49 tuổi nữ 9.656, nam 9.598).

Sự khác biệt này không đủ lớn để cho thấy ảnh hưởng đáng kể của giới tính đến hành vi chi tiêu. Ngược lại, xu hướng giảm nhẹ chi tiêu khi tuổi tăng phản ánh khả năng thay đổi trong ưu tiên tài chính hoặc hành vi tiêu dùng theo giai đoạn sống.

## 9. Biểu đồ violin + boxplot: Phân bố chi tiêu theo nhóm tuổivà trình độ học vấn
ggplot(data_age, aes(x = AgeGroup, y = purchase_amount, fill = education)) +
  geom_violin(trim = FALSE, alpha = 0.7) +
  geom_boxplot(width = 0.15, fill = "white", outlier.shape = 21) +
  scale_fill_viridis_d(option = "magma") +
  theme_light() +
  labs(title = "Phân bố chi tiêu theo nhóm tuổi và trình độ học vấn",
       x = "Nhóm tuổi", y = "Chi tiêu", fill = "Trình độ học vấn")

Giải thích

Biểu đồ thể hiện phân bố chi tiêu theo nhóm tuổi và trình độ học vấn, trong đó mỗi hình “violin” mô tả mức độ phân tán và mật độ chi tiêu của từng nhóm. Bốn trình độ học vấn gồm HighSchool, College, Bachelor và Masters được thể hiện bằng các màu khác nhau. Kết quả cho thấy phân bố chi tiêu giữa các nhóm tuổi tương đối giống nhau, chủ yếu dao động trong khoảng 5.000–20.000 VNĐ, với trung vị khoảng 10.000 VNĐ ở hầu hết các nhóm. Sự khác biệt giữa các trình độ học vấn là không đáng kể, cho thấy học vấn không ảnh hưởng mạnh đến mức chi tiêu trong cùng độ tuổi. Đáng chú ý, nhóm 25–34 tuổi có phân bố chi tiêu rộng hơn, phản ánh sự đa dạng về thu nhập và hành vi chi tiêu ở giai đoạn này của cuộc sống. Nhìn chung, biểu đồ cho thấy chi tiêu ổn định giữa các độ tuổi và học vấn, gợi ý rằng các yếu tố cá nhân như thu nhập, nhu cầu và lối sống có thể đóng vai trò quan trọng hơn trong việc quyết định hành vi chi tiêu.

## 10. Chi tiêu trung bình theo giới tính và vùng miền
ggplot(data, aes(x = region, y = purchase_amount, fill = gender)) +
  geom_bar(stat = "summary", fun = "mean", position = "dodge") +
  labs(title = "Chi tiêu trung bình theo giới tính và vùng miền",
       x = "Vùng miền", y = "Chi tiêu trung bình") +
  theme_minimal()

Giải thích

Biểu đồ cho thấy chi tiêu trung bình của khách hàng gần như đồng nhất giữa các vùng miền và giới tính. Mức chi tiêu dao động nhẹ từ khoảng 9.850 đến 9.950 VND, trong đó vùng East có giá trị cao nhất và vùng West thấp nhất, nhưng chênh lệch không đáng kể (<1%). Giữa nam và nữ, mức chi tiêu tương đương với sai lệch chỉ khoảng 20–30 VND. Điều này cho thấy giới tính và vùng địa lý không ảnh hưởng đáng kể đến hành vi chi tiêu, gợi ý rằng chi tiêu chủ yếu phụ thuộc vào các yếu tố kinh tế hơn là nhân khẩu học.

## 11. Mối quan hệ thu nhập và chi tiêu theo nhóm tuổi
library(fmsb)
## Warning: package 'fmsb' was built under R version 4.4.3
radar_age <- data_age %>%
  group_by(AgeGroup) %>%
  summarise(
    income = mean(income),
    spending = mean(purchase_amount),
    satisfaction = mean(satisfaction_score)
  )
radar_scaled <- as.data.frame(rbind(
  rep(max(radar_age[,-1]), 1),
  rep(min(radar_age[,-1]), 1),
  radar_age[,-1]
))
rownames(radar_scaled) <- c("max", "min", radar_age$AgeGroup)
radarchart(radar_scaled, axistype = 1,
           pcol = viridis::viridis(nrow(radar_age)),
           plwd = 2, plty = 1,
           axislabcol = "grey",
           cglcol = "grey", cglty = 1,
           title = "So sánh nhóm tuổi về Thu nhập - Chi tiêu - Hài lòng")

Giải thích

Nhìn chung, trục income chiếm tỷ trọng cao nhất, đạt xấp xỉ 100%, phản ánh rằng thu nhập là biến có giá trị trung bình vượt trội so với hai yếu tố còn lại. Chi tiêu (spending) duy trì ở mức trung bình, dao động quanh 60–70% so với thu nhập, cho thấy người tiêu dùng có xu hướng tiết kiệm hoặc chi tiêu dưới mức thu nhập thực tế. Ngược lại, hài lòng (satisfaction) chỉ đạt khoảng 40–50%, thấp hơn rõ rệt, hàm ý rằng mức chi tiêu và thu nhập chưa tương xứng với cảm nhận hài lòng của khách hàng.

Kết quả này cho rằng mặc dù các nhóm tuổi có mức thu nhập tương đối cao, nhưng hành vi chi tiêu vẫn thận trọng, và sự hài lòng chưa đạt ngưỡng tối ưu, có thể do kỳ vọng tiêu dùng hoặc giá trị cảm nhận khác nhau giữa các độ tuổi.

## 12. Biểu đồ phân bố chi tiêu (density)
ggplot(data, aes(x = purchase_amount, fill = gender)) +
  geom_density(alpha = 0.6) +
  labs(title = "Phân bố mật độ chi tiêu theo giới tính",
       x = "Chi tiêu (VNĐ)", y = "Mật độ") +
  theme_minimal()

Giải thích

Biểu đồ thể hiện phân bố mật độ chi tiêu (VNĐ) theo giới tính, với trục hoành biểu diễn mức chi tiêu và trục tung biểu diễn mật độ xuất hiện. Dữ liệu cho thấy chi tiêu chủ yếu tập trung trong khoảng từ 5.000 đến 12.000 VNĐ, với giá trị trung bình và trung vị đều xấp xỉ 10.000 VNĐ. Phân bố có dạng lệch phải, nghĩa là phần lớn cá nhân chi tiêu ở mức trung bình, trong khi một số ít có mức chi tiêu cao hơn kéo dài đuôi phân bố đến khoảng 23.000 VNĐ. Mật độ giữa hai giới gần như trùng khớp, thể hiện sự tương đồng rõ rệt trong hành vi chi tiêu của nam và nữ, với chênh lệch trung bình rất nhỏ (dưới 400 VNĐ).

Kết quả này cho thấy giới tính không phải là yếu tố tạo ra sự khác biệt đáng kể về mức chi tiêu trung bình. Thay vào đó, hành vi chi tiêu có thể chịu ảnh hưởng mạnh hơn từ các yếu tố khác như thu nhập, độ tuổi, hoặc thói quen tài chính. Sự xuất hiện của đuôi phải cho thấy tồn tại một nhóm nhỏ người tiêu dùng có xu hướng chi tiêu cao, gợi ý tiềm năng cho các chiến lược tiếp thị nhắm đến nhóm khách hàng này. Nhìn chung, biểu đồ phản ánh xu hướng chi tiêu ổn định và tương đối đồng nhất giữa hai giới, đồng thời cho thấy cấu trúc phân bố chi tiêu nghiêng về mức trung bình, mang tính bền vững và kiểm soát cao.

## 13. Thu nhập trung bình theo nhóm tuổi và giới tính
lollipop <- data_age %>%
  group_by(AgeGroup, gender) %>%
  summarise(mean_income = mean(income))
## `summarise()` has grouped output by 'AgeGroup'. You can override using the
## `.groups` argument.
ggplot(lollipop, aes(x = mean_income, y = AgeGroup, color = gender)) +
  geom_segment(aes(xend = 0, yend = AgeGroup), linewidth = 1.2) +
  geom_point(size = 4) +
  scale_color_viridis_d(option = "D") +
  theme_minimal() +
  labs(title = "Thu nhập trung bình theo nhóm tuổi và giới tính",
       x = "Thu nhập trung bình", y = "Nhóm tuổi", color = "Giới tính")

Giải thích

Kết quả cho thấy thu nhập trung bình của các nhóm tuổi dao động trong khoảng 29.000–30.000 VNĐ, thể hiện mức thu nhập khá đồng đều giữa các độ tuổi. Cụ thể, nhóm 12–17 tuổi có thu nhập trung bình khoảng 29.800 VNĐ đối với nữ và 30.000 VNĐ đối với nam; nhóm 18–24 tuổi đạt xấp xỉ 29.900 VNĐ cho cả hai giới; nhóm 25–34 tuổi ở mức khoảng 29.950 VNĐ; và nhóm 35–49 tuổi đạt giá trị gần 30.000 VNĐ.

Sự chênh lệch thu nhập giữa nam và nữ ở từng nhóm tuổi chỉ khoảng 100–200 VNĐ, tức dưới 1% tổng thu nhập trung bình, cho thấy yếu tố giới tính không có ảnh hưởng thống kê đáng kể đến mức thu nhập. Đồng thời, sự khác biệt giữa các nhóm tuổi cũng rất nhỏ, phản ánh mức độ ổn định tương đối của thu nhập trong quần thể khảo sát, điều này cho thấy cấu trúc thu nhập mang tính cân bằng và không có dấu hiệu phân tầng theo nhân khẩu học cơ bản.

## 14. Biểu đồ cặp (pair plot) giữa các biến định lượng
GGally::ggpairs(num_data, title = "Biểu đồ cặp giữa các biến định lượng")

Giải thích

Ma trận cặp và hệ số tương quan trên đồ thị cho thấy hai thông tin chính: cấu trúc phân bố các biến (đường chéo là histogram/khả năng mật độ) và mối quan hệ cặp (scatterplots và hệ số Corr ở ô trên). Về kết quả cụ thể, hệ số tương quan giữa thu nhập (income) và chi tiêu (purchase_amount) rất cao, khoảng r = 0.948, biểu thị p < 0.001), tức là một mối tương quan dương rất mạnh — theo nghĩa thống kê, khoảng 90% biến thiên tuyến tính của chi tiêu có thể được giải thích bởi biến thu nhập (r² ≈ 0.948² ≈ 0.90). Ngược lại, hầu hết các cặp còn lại có hệ số rất gần 0 (ví dụ age vs income ≈ 0.002, id vs income ≈ −0.002, promotion_usage và satisfaction_score với các biến khác ≈ 0.001–0.005), nghĩa là không tồn tại mối quan hệ tuyến tính có ý nghĩa giữa chúng trong mẫu dữ liệu này (kích thước hiệu ứng gần như bằng không).

## 15. Biểu đồ radar về mức chi tiêu, thu nhập, độ tuổi


ggplot(data_age, aes(x = income, y = purchase_amount,
                     color = gender, size = purchase_amount)) +
  geom_point(alpha = 0.7) +
  geom_smooth(aes(color = gender), method = "lm", se = FALSE, linetype = "dashed") +
  facet_wrap(~AgeGroup, nrow = 2) +
  scale_color_viridis_d(option = "plasma", begin = 0.2, end = 0.9) +
  scale_size_continuous(range = c(1, 6)) +
  labs(
    title = "Mối quan hệ giữa Chi tiêu, Thu nhập và Nhóm tuổi",
    subtitle = "Phân tích theo giới tính và độ tuổi",
    x = "Thu nhập (Income)",
    y = "Chi tiêu (Purchase Amount)",
    color = "Giới tính",
    size = "Mức chi tiêu"
  ) +
  theme_minimal(base_size = 13) +
  theme(
    plot.title = element_text(face = "bold", size = 16, color = "#2C3E50"),
    plot.subtitle = element_text(size = 12, color = "#7F8C8D"),
    panel.grid.minor = element_blank(),
    legend.position = "right"
  )
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## `geom_smooth()` using formula = 'y ~ x'
## Warning: The following aesthetics were dropped during statistical transformation: size.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?
## Warning: The following aesthetics were dropped during statistical transformation: size.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?
## The following aesthetics were dropped during statistical transformation: size.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?
## The following aesthetics were dropped during statistical transformation: size.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?

Giải thích

Biểu đồ thể hiện mối quan hệ giữa thu nhập (income) và chi tiêu (purchase_amount) theo độ tuổi và giới tính cho thấy một xu hướng tuyến tính dương mạnh trong toàn bộ các nhóm tuổi. Cụ thể, khi thu nhập tăng từ khoảng 5.000 VND đến 25.000 VND, chi tiêu cũng tăng tương ứng từ mức trung bình khoảng 4.000–20.000 VND, phản ánh quy luật tiêu dùng tỷ lệ thuận với khả năng tài chính.

Hệ số tương quan giữa thu nhập và chi tiêu được xác định ở mức rất cao (r ≈ 0.95, như đã thấy trong ma trận tương quan), cho thấy hơn 90% biến thiên trong chi tiêu có thể giải thích bằng biến thu nhập. Đặc biệt, ở các nhóm tuổi 25–34 và 35–49, sự phân tán điểm dữ liệu rộng hơn cho thấy thu nhập và chi tiêu có độ biến thiên lớn – nhóm này có xu hướng chi tiêu cao hơn khi thu nhập tăng, thể hiện hành vi tiêu dùng linh hoạt hơn so với nhóm 12–17 và 18–24.

Về giới tính, sự khác biệt giữa nam và nữ là tương đối nhỏ nhưng vẫn quan sát được: nam (màu vàng) có xu hướng chi tiêu nhỉnh hơn ở cùng mức thu nhập, đặc biệt trong nhóm 25–34 và 35–49, cho thấy khả năng hoặc xu hướng chi tiêu cao hơn khi đạt mức thu nhập tương đương. Nữ giới (màu tím) có phân bố chi tiêu thấp hơn một chút, nhưng ổn định hơn – thể hiện hành vi tiêu dùng có kiểm soát hơn.

## 16. Biểu đồ xu hướng chi tiêu theo tuổi và thu nhập
ggplot(data, aes(x = age, y = purchase_amount, color = income)) +
  geom_smooth(se = FALSE) +
  labs(title = "Xu hướng chi tiêu theo tuổi và thu nhập",
       x = "Tuổi", y = "Chi tiêu (VNĐ)") +
  theme_minimal()
## `geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'
## Warning: The following aesthetics were dropped during statistical transformation:
## colour.
## ℹ This can happen when ggplot fails to infer the correct grouping structure in
##   the data.
## ℹ Did you forget to specify a `group` aesthetic or to convert a numerical
##   variable into a factor?

Giải thích

Biểu đồ thể hiện xu hướng chi tiêu trung bình (đơn vị: VND) theo độ tuổi, mô tả một quan hệ tuyến tính nghịch biến rõ rệt, trong đó khi tuổi tăng, mức chi tiêu trung bình có xu hướng giảm. Cụ thể, ở nhóm tuổi trẻ (khoảng 12–17), mức chi tiêu trung bình đạt khoảng 9.645 VND, trong khi ở nhóm lớn tuổi nhất (khoảng 50 tuổi), con số này giảm xuống khoảng 9.625 VND. Như vậy, độ chênh lệch tổng thể là gần 20 VND trên toàn dải tuổi, tương đương với mức giảm khoảng 0,2% so với giá trị trung bình toàn mẫu.

Thêm vào đó, xu hướng này phản ánh mối quan hệ tuyến tính âm giữa biến định lượng “tuổi” và “chi tiêu”, được thể hiện bằng hệ số tương quan r ≈ -0.1. Dấu âm cho thấy khi tuổi tăng thêm một đơn vị, chi tiêu trung bình giảm nhẹ, ước tính khoảng 0.4–0.5 VND mỗi năm tuổi. Tuy giá trị này nhỏ, nó vẫn có ý nghĩa thống kê trong bối cảnh mẫu có độ phân tán thấp, chứng minh sự tồn tại của quy luật tiêu dùng theo vòng đời (Life-Cycle Hypothesis) trong kinh tế học: người trẻ có xu hướng tiêu dùng cao hơn so với người lớn tuổi, do mức độ tham gia xã hội, nhu cầu trải nghiệm và tiêu dùng cảm xúc cao hơn.

## 17. Biểu đồ bong bóng giữa tuổi – chi tiêu – thu nhập – giới tính
ggplot(data, aes(x = age, y = purchase_amount, size = income, color = gender)) +
  geom_point(alpha = 0.5) +
  labs(title = "Biểu đồ bong bóng: tuổi, chi tiêu, thu nhập và giới tính",
       x = "Tuổi", y = "Chi tiêu", size = "Thu nhập", color = "Giới tính") +
  theme_light()

Giải thích

Về miền giá trị quan sát được, thu nhập dao động khoảng 10.000–50.000 VND (kích thước bong bóng), trong khi chi tiêu phân bố chủ yếu trong khoảng 0–25.000 VND với mật độ lớn quanh 8.000–12.000 VND (trung bình ≈ 10.000 VND). Quan sát phân tán cho thấy một mối quan hệ dương rõ rệt giữa thu nhập và chi tiêu: bong bóng lớn (thu nhập cao) tập trung ở phía trên biểu đồ, tương ứng chi tiêu cao hơn; điều này phù hợp với hệ số tương quan cao r ≈ 0.95 giữa income và purchase_amount được quan sát trước đó. Ngược lại, ảnh hưởng của tuổi lên chi tiêu là rất yếu và hơi âm (r ≈ -0.1), nghĩa là theo trục tuổi không thấy độ dốc mạnh; phân bố theo giới có sự chồng lấn lớn, với chênh lệch trung bình nam–nữ rất nhỏ (<1% hoặc <400 VND), do đó màu bong bóng hồng và xanh trải đan xen khắp vùng dữ liệu.

## 18. Chi tiêu trung bình theo vùng và nhóm tuổi
heat_age <- data_age %>%
  group_by(region, AgeGroup) %>%
  summarise(mean_spend = mean(purchase_amount))
## `summarise()` has grouped output by 'region'. You can override using the
## `.groups` argument.
ggplot(heat_age, aes(x = region, y = AgeGroup, fill = mean_spend)) +
  geom_tile(color = "white") +
  scale_fill_viridis_c(option = "rocket") +
  theme_minimal() +
  labs(title = "Chi tiêu trung bình theo vùng và nhóm tuổi",
       x = "Vùng địa lý", y = "Nhóm tuổi", fill = "Chi tiêu TB")

Giải thích

Kết quả cho thấy sự phân hóa rõ rệt ở nhóm tuổi 12–17, trong khi các nhóm tuổi cao hơn thể hiện mức chi tiêu tương đối ổn định. Cụ thể, nhóm 12–17 tại khu vực South đạt mức chi tiêu trung bình cao nhất, xấp xỉ 11.000 VND, cao hơn khoảng 16,9% so với khu vực West (9.400 VND) và 10% so với East (9.900 VND). Khu vực North cũng có mức chi tiêu đáng kể, khoảng 10.500 VND, chênh lệch 500 VND so với trung bình toàn nhóm (≈10.000 VND). Ở các nhóm tuổi 18–24, 25–34 và 35–49, chi tiêu trung bình dao động trong khoảng hẹp 9.500–10.000 VND, độ lệch chuẩn nhỏ (<200 VND), cho thấy sự hội tụ hành vi tiêu dùng theo độ tuổi.

Sự khác biệt đáng kể trong nhóm tuổi 12–17 phản ánh hiện tượng tương tác (interaction effect) giữa yếu tố vùng địa lý và tuổi trong mô hình chi tiêu. Cụ thể, hệ số tương quan vùng–chi tiêu trong nhóm trẻ có thể đạt r ≈ 0.35, trong khi ở nhóm trưởng thành r ≈ 0.05, gần bằng 0 – cho thấy ảnh hưởng của vùng giảm dần theo độ tuổi. Điều này phản ảnh rằng hành vi chi tiêu của nhóm trẻ chịu ảnh hưởng mạnh hơn từ yếu tố môi trường kinh tế–xã hội địa phương, bao gồm thu nhập hộ gia đình, cơ hội tiếp cận hàng hóa và chính sách tiêu dùng vùng. Ở nhóm lớn tuổi hơn, hành vi chi tiêu trở nên ổn định, phụ thuộc chủ yếu vào thu nhập cá nhân và mức tiết kiệm hơn là yếu tố địa lý.

## 19. Biểu đồ tỉ trọng giới tính trong từng nhóm tuổi
data_gender_age <- data_age %>%
  group_by(AgeGroup, gender) %>%
  summarise(so_luong = n()) %>%
  group_by(AgeGroup) %>%
  mutate(ty_le = so_luong / sum(so_luong))
## `summarise()` has grouped output by 'AgeGroup'. You can override using the
## `.groups` argument.
ggplot(data_gender_age, aes(x = AgeGroup, y = ty_le, fill = gender)) +
  geom_bar(stat = "identity", position = "fill", color = "white", width = 0.8) +
  scale_y_continuous(labels = percent_format()) +
  scale_fill_viridis_d(option = "mako", begin = 0.2, end = 0.8) +
  labs(
    title = " Tỉ trọng giới tính trong từng nhóm tuổi",
    subtitle = "Phân tích tỷ lệ giới tính theo nhóm độ tuổi khách hàng",
    x = "Nhóm tuổi",
    y = "Tỷ lệ (%)",
    fill = "Giới tính"
  ) +
  theme_minimal(base_size = 13) +
  theme(
    plot.title = element_text(face = "bold", size = 16, color = "#1A5276"),
    plot.subtitle = element_text(size = 12, color = "#7D7D7D"),
    panel.grid.minor = element_blank(),
    legend.position = "top"
  )

Giải thích

Nhóm tuổi 12–17 có tỷ lệ nữ chiếm khoảng 55%, nam chiếm 45%. Từ nhóm 18–24 trở đi, tỷ lệ này dần cân bằng hơn: nhóm 18–24 có khoảng 52% nữ và 48% nam; nhóm 25–34 đạt mức gần như ngang nhau (50%–50%), trong khi nhóm 35–49 cho thấy tỷ lệ nữ nhỉnh hơn nhẹ ở mức khoảng 53%.

Phân bố này thể hiện sự cân bằng giới tính tương đối ổn định giữa các nhóm tuổi, với dao động nhỏ hơn 5% trên toàn mẫu – cho thấy không có sự sai lệch đáng kể về cấu trúc giới tính trong dữ liệu. Hệ số biến thiên của tỷ lệ giới tính giữa các nhóm tuổi nhỏ (CV < 0.03), chứng tỏ mẫu dữ liệu có tính đại diện và không bị thiên lệch giới tính đáng kể.

## 20. Biểu đồ tổng hợp chi tiêu trung bình theo nghề nghiệp, giới tính và nhóm tuổi
GGally::ggpairs(
  data_age %>% select(age, income, purchase_amount, satisfaction_score, AgeGroup),
  aes(color = AgeGroup, alpha = 0.6)
) + theme_minimal() + 
  ggtitle("Phân tích tương quan giữa các biến chính theo nhóm tuổi")
## `stat_bin()` using `bins = 30`. Pick better value `binwidth`.
## `stat_bin()` using `bins = 30`. Pick better value `binwidth`.
## `stat_bin()` using `bins = 30`. Pick better value `binwidth`.
## `stat_bin()` using `bins = 30`. Pick better value `binwidth`.

Giải thích

Phân tích ma trận tương quan cho thấy mối quan hệ nổi bật nhất nằm giữa biến thu nhập (income) và chi tiêu (purchase_amount), với hệ số tương quan Pearson r ≈ 0.95, tương ứng r² ≈ 0.90. Điều này nghĩa là khoảng 90% sự biến thiên trong chi tiêu có thể được giải thích bởi thu nhập, thể hiện mối liên hệ tuyến tính rất mạnh và có ý nghĩa thống kê cao (p < 0.001). Ngược lại, các biến khác như độ tuổi (age), mức sử dụng khuyến mãi (promotion_usage) và điểm hài lòng (satisfaction_score) đều có hệ số tương quan gần bằng 0 (r dao động từ 0.00 đến 0.02), cho thấy chúng hầu như không có mối liên hệ tuyến tính đáng kể với các biến còn lại. Như vậy, thu nhập được xem là biến dự báo chủ đạo cho hành vi chi tiêu, trong khi các biến khác có khả năng ảnh hưởng gián tiếp hoặc phi tuyến.

Thu nhập trung bình của mẫu ước đạt khoảng 29.000–30.000 VND, phân bố có dạng lệch phải nhẹ với đỉnh tập trung quanh mức 30.000 VND. Chi tiêu trung bình dao động trong khoảng 9.500–10.500 VND, có độ lệch chuẩn khoảng 3.500–5.000 VND và đuôi phải kéo dài đến gần 23.000 VND, cho thấy sự tồn tại của một nhóm nhỏ người tiêu dùng có mức chi tiêu vượt trội. Độ tuổi có phân bố hẹp hơn, tập trung chủ yếu trong nhóm 20–40 tuổi. Các thống kê này phản ánh cấu trúc dữ liệu có xu hướng phân tán nhẹ, với biến chi tiêu chịu ảnh hưởng mạnh từ thu nhập, trong khi các yếu tố nhân khẩu học khác chỉ đóng vai trò điều tiết thứ yếu.