----------------------------------------------------------------------------- #Make figure with stepwise increases of delta15N isotopes<-read.csv("bigislandiso.csv",header=TRUE) isotopee<-subset(isotopes,site!="PuuMak") isotope<-subset(isotopee,controlled=="yes") isotope$plot<-as.factor(as.vector(isotope$plot)) par(mar=c(6,6,4,12)+0.1) plot(delta15N~plot,data=isotope) #likewise for delta13C plot(delta13C~plot,data=isotope) -----------------------------------------------------------------------------