Fire influences on needle decomposition: thermal tipping point in Pinus radiata needle carbon chemistry and soil nitrogen transformations 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: Packages library(readr)#for reading data library(dplyr)#for manipulating data library(ggplot2)#for plotting library(agricolae)#for stastics ##########----------########## MBN ##########----------########## MBN<-read.csv("2018 Scorchin n_chem.csv")%>% mutate(MBN=MBN/0.54) test_0.00<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.00) anova(lm(MBN~ID, data=test_0.00))#*** HSD.test(aov(MBN~ID, data=test_0.00), c("ID"), console=TRUE, alpha=0.05) test_0.25<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.25) anova(lm(MBN~ID, data=test_0.25))#*** HSD.test(aov(MBN~ID, data=test_0.25), c("ID"), console=TRUE, alpha=0.05) test_0.5<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.5) anova(lm(MBN~ID, data=test_0.5))#ns HSD.test(aov(MBN~ID, data=test_0.5), c("ID"), console=TRUE, alpha=0.05) test_0.75<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.75) anova(lm(MBN~ID, data=test_0.75))#*** HSD.test(aov(MBN~ID, data=test_0.75), c("ID"), console=TRUE, alpha=0.05) test_1<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==1) anova(lm(MBN~ID, data=test_1))#*** HSD.test(aov(MBN~ID, data=test_1), c("ID"), console=TRUE, alpha=0.05) test_2<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==2) anova(lm(MBN~ID, data=test_2))#*** HSD.test(aov(MBN~ID, data=test_2), c("ID"), console=TRUE, alpha=0.05) test_3<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==3) anova(lm(MBN~ID, data=test_3))#*** HSD.test(aov(MBN~ID, data=test_3), c("ID"), console=TRUE, alpha=0.05) test_4<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==4) anova(lm(MBN~ID, data=test_4))#*** HSD.test(aov(MBN~ID, data=test_4), c("ID"), console=TRUE, alpha=0.05) test_6<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==6) anova(lm(MBN~ID, data=test_6))#*** HSD.test(aov(MBN~ID, data=test_6), c("ID"), console=TRUE, alpha=0.05) test_8<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==8) anova(lm(MBN~ID, data=test_8))#*** HSD.test(aov(MBN~ID, data=test_8), c("ID"), console=TRUE, alpha=0.05) test_11<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==11) anova(lm(MBN~ID, data=test_11))#*** HSD.test(aov(MBN~ID, data=test_11), c("ID"), console=TRUE, alpha=0.05) test_14<-MBN%>% select(TIME, MBN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==14) anova(lm(MBN~ID, data=test_14))#ns HSD.test(aov(MBN~ID, data=test_14), c("ID"), console=TRUE, alpha=0.05) ##########----------########## PMN ##########----------########## PMN<-read.csv("2018 Scorchin n_chem.csv") test_0.00<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.00) anova(lm(PMN~ID, data=test_0.00))#*** HSD.test(aov(PMN~ID, data=test_0.00), c("ID"), console=TRUE, alpha=0.05) test_0.25<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.25) anova(lm(PMN~ID, data=test_0.25))#*** HSD.test(aov(PMN~ID, data=test_0.25), c("ID"), console=TRUE, alpha=0.05) test_0.5<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.5) anova(lm(PMN~ID, data=test_0.5))#*** HSD.test(aov(PMN~ID, data=test_0.5), c("ID"), console=TRUE, alpha=0.05) test_0.75<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.75) anova(lm(PMN~ID, data=test_0.75))#*** HSD.test(aov(PMN~ID, data=test_0.75), c("ID"), console=TRUE, alpha=0.05) test_1<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==1) anova(lm(PMN~ID, data=test_1))#*** HSD.test(aov(PMN~ID, data=test_1), c("ID"), console=TRUE, alpha=0.05) test_2<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==2) anova(lm(PMN~ID, data=test_2))#*** HSD.test(aov(PMN~ID, data=test_2), c("ID"), console=TRUE, alpha=0.05) test_3<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==3) anova(lm(PMN~ID, data=test_3))#** HSD.test(aov(PMN~ID, data=test_3), c("ID"), console=TRUE, alpha=0.05) test_4<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==4) anova(lm(PMN~ID, data=test_4))#*** HSD.test(aov(PMN~ID, data=test_4), c("ID"), console=TRUE, alpha=0.05) test_6<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==6) anova(lm(PMN~ID, data=test_6))#*** HSD.test(aov(PMN~ID, data=test_6), c("ID"), console=TRUE, alpha=0.05) test_8<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==8) anova(lm(PMN~ID, data=test_8))#** HSD.test(aov(PMN~ID, data=test_8), c("ID"), console=TRUE, alpha=0.05) test_11<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==11) anova(lm(PMN~ID, data=test_11))#*** HSD.test(aov(PMN~ID, data=test_11), c("ID"), console=TRUE, alpha=0.05) test_14<-PMN%>% select(TIME, PMN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==14) anova(lm(PMN~ID, data=test_14))#*** HSD.test(aov(PMN~ID, data=test_14), c("ID"), console=TRUE, alpha=0.05) ##########----------########## NMIN ##########----------########## AN<-read.csv("2018 Scorchin n_chem.csv")%>% mutate(AN=N03+NH4) test_0.00<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.00) anova(lm(AN~ID, data=test_0.00))#*** HSD.test(aov(AN~ID, data=test_0.00), c("ID"), console=TRUE, alpha=0.05) test_0.25<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.25) anova(lm(AN~ID, data=test_0.25))#*** HSD.test(aov(AN~ID, data=test_0.25), c("ID"), console=TRUE, alpha=0.05) test_0.5<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.5) anova(lm(AN~ID, data=test_0.5))#*** HSD.test(aov(AN~ID, data=test_0.5), c("ID"), console=TRUE, alpha=0.05) test_0.75<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==0.75) anova(lm(AN~ID, data=test_0.75))#ns HSD.test(aov(AN~ID, data=test_0.75), c("ID"), console=TRUE, alpha=0.05) test_1<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==1) anova(lm(AN~ID, data=test_1))#*** HSD.test(aov(AN~ID, data=test_1), c("ID"), console=TRUE, alpha=0.05) test_2<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==2) anova(lm(AN~ID, data=test_2))#*** HSD.test(aov(AN~ID, data=test_2), c("ID"), console=TRUE, alpha=0.05) test_3<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==3) anova(lm(AN~ID, data=test_3))#*** HSD.test(aov(AN~ID, data=test_3), c("ID"), console=TRUE, alpha=0.05) test_4<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==4) anova(lm(AN~ID, data=test_4))#*** HSD.test(aov(AN~ID, data=test_4), c("ID"), console=TRUE, alpha=0.05) test_6<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==6) anova(lm(AN~ID, data=test_6))#*** HSD.test(aov(AN~ID, data=test_6), c("ID"), console=TRUE, alpha=0.05) test_8<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==8) anova(lm(AN~ID, data=test_8))#*** HSD.test(aov(AN~ID, data=test_8), c("ID"), console=TRUE, alpha=0.05) test_11<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==11) anova(lm(AN~ID, data=test_11))#*** HSD.test(aov(AN~ID, data=test_11), c("ID"), console=TRUE, alpha=0.05) test_14<-AN%>% select(TIME, AN, ID)%>% mutate(ID=str_sub(ID,0,3))%>% filter(TIME==14) anova(lm(AN~ID, data=test_14))#*** HSD.test(aov(AN~ID, data=test_14), c("ID"), console=TRUE, alpha=0.05)