improvement <- readxl::read_excel("state-data.xlsx", sheet = "improvement")
attach(improvement)
boxplot(t.before, t.after)
## check the normality
shapiro.test(t.before)
shapiro.test(t.after)
## Wilcoxon signed rank test
wilcox.test(t.before, t.after, paired = TRUE)