###Set Working Directory (fill in path to data in parenthesis)
setwd("~/Desktop/Research/")
###Read in data in .csv file, and assign name of "DO_Stats2" (left of arrow)
DO_Stats2 <- read.csv(DO_Stats2, header=TRUE, sep = ",", dec = ".")
###Set-up linear model of data named DO_Stats2.lm, 
  ## "DO_Stats2$..." specifies data within chart DO_Stats2 and coulumn $... where ... is a header
 ## "DO_Stats2$Push_ups~DO_Stats2$DO_levels" indicates test of push-up influenced by DO_levels
  ## "+ DO_Stats2$Temperature_C + DO_Stats2$Group" indicates additon of tempeature and Group(order of experiments) as covariates.
DO_Stats2.lm = lm(DO_Stats2$Push_ups~DO_Stats2$DO_levels + DO_Stats2$Temperature_C + DO_Stats2$Group)
####Perform anova on model
group3.anova = Anova(DO_Stats2.lm)
####Display anova Results
DO_Stats2.anova