SEP-8-301.Rd

# Author: Pawar, Date: Sep 8 2021, Purpose: To perform ANOVA test on sample dataset

library(dplyr) 
PATH <- "https://raw.githubusercontent.com/guru99-edu/R-Programming/master/poisons.csv" 
df <- read.csv(PATH) %>% select(-X) %>% mutate(poison = factor(poison, ordered = TRUE)) 
glimpse(df)
levels(df$poison)
df %>% group_by(poison) %>% summarise( count_poison = n(), mean_time = mean(time, na.rm = TRUE), sd_time = sd(time, na.rm = TRUE) )
library(ggplot2)
ggplot(df, aes(x = poison, y = time, fill = poison)) + geom_boxplot() + geom_jitter(shape = 15, color = "steelblue", position = position_jitter(0.21)) + theme_classic()
anova_one_way <- aov(time~poison, data = df)
summary(anova_one_way)

###IF THE F VALUE IS MORE THAN THE CRTICAL VALUE THEN THERE IS A SIGNIFICANT DIFFERENCE BETWEEN POISON AND TIME############ 

# If you want to perform a pairwise compariosn then you do following

TukeyHSD(anova_one_way)