isotopes<-read.csv("bigislandiso.csv")
isotopee<-subset(isotopes,site!="PuuMak")
isotope<-subset(isotopee,controlled=="yes")
summary(isotope)
Nall<-lm(delta15N~site*category,data=isotope)
anova(Nall)
summary(Nall)

Call<-lm(delta13C~site*category,data=isotope)
anova(Call)
summary(Call)