setwd()setwd("C:/Users/Shelby/Desktop/")


fruitfly <- read.table("fruitfly3.txt", sep="", header=TRUE)
head(fruitfly)
str(fruitfly)

hist(fruitfly$lifespan)
boxplot(lifespan ~ group, data=fruitfly, main = "Group effect on Lifespan",xlab = "Group", ylab = "Lifespan")

fruitfly.lm <- lm(lifespan ~ group, data=fruitfly)
anova(fruitfly.lm)

fruitfly.res = aov(lifespan ~ group, data=fruitfly)
summary(fruitfly.res)

TukeyHSD(fruitfly.res)