tailgating <- read.delim('https://raw.githubusercontent.com/IowaBiostat/data-sets/main/tailgating/tailgating.txt')

summary(tailgating$Distance)

sd(tailgating$Distance)

by(tailgating$Distance, tailgating$Group, summary)

hist(tailgating$Distance)

hist(tailgating$Distance, main = "Histogram of Tailgating Distance", 
     xlab = "Following Distance")
abline(v=mean(tailgating$Distance), 
       col="steelblue4", lwd=2) #lwd makes the line thicker (line width)
abline(v=median(tailgating$Distance), col="orange", 
       lwd=2, lty = 2)#lty makes the line dashed (line type) 
legend(x = "topright", 
       legend = c("Median", "Mean"), 
       col = c("orange", "steelblue4"),
       lwd = 2, lty = c(2,1))

# the mean is the solid line and the median is the dashed line 

hist(tailgating$Distance, breaks = seq(0, 400, 2)) # bin size of 2

hist(tailgating$Distance, breaks = seq(0, 400, 100)) # bin size of 100

hist(tailgating$Distance, breaks = seq(0, 400, 10)) # bin size of 10

#customized labels, solid color & white border   
hist(tailgating$Distance, col= "pink", border="white", breaks = seq(1, 400, 10),
     xlab="Distance",
     ylab="Frequency",
     main = "")

par(mfrow=c(2,2)) # view all four histograms in a 2 by 2 window

hist(tailgating$Distance[tailgating$Group=="ALC"], col= "yellow", breaks = seq(1, 400, 10), 
     main = "", xlab = "ALC")
hist(tailgating$Distance[tailgating$Group=="MDMA"], col= "red", breaks = seq(1, 400, 10), 
     main = "", xlab = "MDMA")
hist(tailgating$Distance[tailgating$Group=="NODRUG"], col= "blue", breaks = seq(1, 400, 10), 
     main = "", xlab = "NoDrug")
hist(tailgating$Distance[tailgating$Group=="THC"], col= "green", breaks = seq(1, 400, 10), 
     main = "", xlab = "THC")

boxplot(tailgating$Distance, 
        ylab = "Distance")

boxplot(tailgating$Distance ~ tailgating$Group, 
        col=rainbow(4),
        xlab = "Group",
        ylab = "Distance")

boxplot(tailgating$Distance ~ tailgating$Group, col= rainbow(4), outline=FALSE)

quantile(tailgating$Distance, 0.30)

library(kableExtra)

tab2 <- data.frame(. = c("Don't Reject", "Reject", "Total"),
                   NullT = character(3),
                   NullF = character(3),
                   Total = character(3))
colnames(tab2) <- c("", "True Null", "False Null", "Total")

kable(tab2) %>%
  kable_styling(bootstrap_options = c("striped", "hover"))

tab3 <- data.frame(. = c("Don't Reject", "Reject", "Total"),
                   NullT = c(630, 70, 700),
                   NullF = c(20, 80, 100),
                   Total = c(650, 150, 800))
colnames(tab3) <- c("", "True Null", "False Null", "Total")

kable(tab3) %>%
  kable_styling(bootstrap_options = c("striped", "hover"))

sitstand <- data.frame("Group" = c("Sitting", "Standing"),
                       Perfect.Responses = c(71, 94),
                       Total.Responses = c(134, 142),
                       Perfect_Responses = c(88, 57),
                       Total_Responses = c(104, 90))
kable(sitstand)%>%
  kable_styling()%>%
  add_header_above(c(" " = 1, "Drinking" = 2, "Eating" = 2))