00_迷你小案例.R

# A mini but complete case with only 50 lines of code

# Data Import -------------------------------------------------------------
library(readr)
cjb_url <- "data/cjb.csv"
cjb <- read_csv(cjb_url, locale = locale(encoding = "CP936"))
View(cjb) # View the data

# Data Exploration --------------------------------------------------------
library(tidyverse)
cjb %>%
    dplyr::select(sx, wlfk) %>%
    ggplot(aes(x = wlfk,
               y = sx,
               fill = wlfk)) +
    geom_boxplot(width = 0.5)

# Data Preparation --------------------------------------------------------
as_five_grade_scores <- function(x) {
    cut(
        x,
        breaks = c(0, seq(60, 100, by = 10)),
        include.lowest = TRUE,
        right = FALSE,
        ordered_result = TRUE,
        labels = c("不及格", "及格", "中", "良", "优")
    )
}
cjb <- cjb %>%
    mutate(zcj = rowSums(.[4:12])) %>%
    filter(zcj != 0) %>% #剔除脏数据
    mutate_at(vars(xb, wlfk), factor) %>% #类型转换
    mutate_at(vars(yw:sw), as_five_grade_scores)#数据分箱
View(cjb) #查看转换后的数据

# Model -------------------------------------------------------------------
library(arulesViz)
my_model <- cjb %>%
    select(xb:wlfk) %>%
    apriori(parameter = list(supp = 0.06, conf = 0.8, maxlen=11),
            appearance = list(rhs = paste0("wlfk=", c("文科", "理科"))))

# Visualization -----------------------------------------------------------
inspectDT(my_model) #输出挖掘到的规则，查看其性能指标

plot(my_model, method = "graph") #结果可视化

# The End ^-^ -------------------------------------------------------------