In-class_EX04

pacman::p_load(tidyverse, ggstatsplot,tidyr,dplyr,rstantools,readxl, performance, parameters)

exam <- read_csv("data/Exam_data.csv")

set.seed(1234)

gghistostats(
  data = exam,
  x = ENGLISH,
  type = "bayes",
  #type = "parametric","bayes",
  test.value = 60,
  bin.args = list(color = "black",
                  fill = "grey50",
                  alpha = 0.7),
  normal.curve = FALSE,
  normal.curve.args = list(linewidth = 2),
  xlab = "English scores"
)

ggdotplotstats(
  data = exam,
  x = ENGLISH,
  y = CLASS,
  title = "",
  xlab = ""
  
)

exam_long <- exam %>%
  pivot_longer(
  cols = ENGLISH:SCIENCE,
  names_to = "SUBJECT",
  values_to = "SCORES") %>%
 filter(CLASS == "3A")
ggwithinstats(data = filter(exam_long,
                            SUBJECT %in% c("ENGLISH", "SCIENCE")),
              x = "SUBJECT",
              y = "SCORES",
              type = "p")

ggscatterstats(data = exam,
               x = MATHS,
               y = ENGLISH,
               marginal = TRUE,
               label.var = ID,
               label.expression = ENGLISH > 90 & MATHS > 90,
  
  )

ggscatterstats(data = exam,
               x = MATHS,
               y = ENGLISH,
               marginal = TRUE,
               label.var = ID,
               label.expression = ENGLISH > 90 & MATHS > 90,
  
)

data_cars <- read_xls("data/ToyotaCorolla.xls","data")