SimpleToast Data analysis - respiration

RScript written by Erinne Stirling 2018

R version 3.5.1 (2018-07-02) -- "Feather Spray"
Copyright (C) 2018 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin15.6.0 (64-bit)

#Create the environment

library(readr)
library(dplyr)
library(ggplot2)
library(agricolae)


###Prepare the Data
TOAST<- read_csv("CSV/CSV toasted_all.csv", col_names=TRUE)%>%#'SimpleToast chemical data.csv'
  group_by(SOIL,LIT,TIME)%>%
  summarise(MBN_m = mean(MBN, na.rm = TRUE),
            MBP_m = mean(MBP, na.rm = TRUE),
            AVN_m = mean(AVN, na.rm = TRUE),
            AVP_m = mean(AVP, na.rm = TRUE))

#Stats
TOAST_0<-TOAST%>%
  filter(TIME==0)

TOAST_1<-TOAST%>%
  filter(TIME==1)

TOAST_2<-TOAST%>%
  filter(TIME==2)

anova(lm(MBN_m~LIT, data=TOAST_1))#ns
HSD.test(aov(MBN_m~LIT, data=TOAST_1), c("LIT"), console=TRUE, alpha=0.05)

anova(lm(MBN_m~LIT, data=TOAST_2))#ns
HSD.test(aov(MBN_m~LIT, data=TOAST_2), c("LIT"), console=TRUE, alpha=0.05)


anova(lm(MBN_m~LIT, data=TOAST_1))#ns
HSD.test(aov(MBN_m~LIT, data=TOAST_1), c("LIT"), console=TRUE, alpha=0.05)

anova(lm(MBN_m~LIT, data=TOAST_2))#ns
HSD.test(aov(MBN_m~LIT, data=TOAST_2), c("LIT"), console=TRUE, alpha=0.05)


anova(lm(MBP_m~LIT, data=TOAST_1))#ns
HSD.test(aov(MBP_m~LIT, data=TOAST_1), c("LIT"), console=TRUE, alpha=0.05)

anova(lm(MBP_m~LIT, data=TOAST_2))#ns
HSD.test(aov(MBP_m~LIT, data=TOAST_2), c("LIT"), console=TRUE, alpha=0.05)


anova(lm(AVP_m~LIT, data=TOAST_1))#ns
HSD.test(aov(AVP_m~LIT, data=TOAST_1), c("LIT"), console=TRUE, alpha=0.05)

anova(lm(AVP_m~LIT, data=TOAST_2))#ns
HSD.test(aov(AVP_m~LIT, data=TOAST_2), c("LIT"), console=TRUE, alpha=0.05)


anova(lm(AVN_m~LIT, data=TOAST_1))#***
HSD.test(aov(AVN_m~LIT, data=TOAST_1), c("LIT"), console=TRUE, alpha=0.05)

anova(lm(AVN_m~LIT, data=TOAST_2))#***
HSD.test(aov(AVN_m~LIT, data=TOAST_2), c("LIT"), console=TRUE, alpha=0.05)