- 22nd Dec 2022
- 06:03 am

setwd("~/Assignments/Job 2")

rm(list = ls())

# Load libraries:

library(reshape2)

library(ggpubr)

# Read the data:

data <- read.csv('length (Autosaved).csv')

# Convert the structure for ANOVA analysis:

data <- melt(data, id = c('Month'))

# Two-way interaction plot

interaction.plot(x.factor = data$variable, trace.factor = data$Month,

response = data$value, fun = mean,

type = "b", legend = TRUE,

xlab = "Treatment", ylab = "Month",

pch=c(1,19), col = c("#00AFBB", "#E7B800"))

# Two lines are not parallel thus implying no interaction between Month and treatment factors.

# Perform a two way ANOVA:

anova.result <- aov(value ~ Month + variable,

data = data)

summary(anova.result)

# Thus at 5% level of significance we can conclude that there is a significant statistical difference

# between the lengths in September to February the following year.