head(TempSite2a) names(TempSite2a)<-c("Time", "Marsh_Temp", "Mangrove_Temp","Sand_Flat_Temp") if(!require(default.plot)){install.packages("default.plot")} if(!require(plotrix)){install.packages("plotrix")} libary(default.plot) libary(plotly) plot(TempSite2$Sand_Flat_Temp, TempSite2$Marsh_Temp) plot(TempSite2) cor(TempSite2a$Marsh_Temp,TempSite2a$Sand_Flat_Temp) cor(TempSite2a$Mangrove_Temp, TempSite2a$Sand_Flat_Temp) cor(TempSite2a$Mangrove_Temp, TempSite2a$Marsh_Temp) library(plotrix) plot(TempSite2a$Sand_Flat_Temp, TempSite2a$Marsh_Temp, col = "black", pch = 21, bg = "grey", cex = 2, ylab = "Marsh Grass Temperature ˚C", xlab = "Sand Flat Temperature ˚C") plot(TempSite2a$Sand_Flat_Temp, TempSite2a$Mangrove_Temp, col = "black", pch = 21, bg = "grey", cex = 2, ylab = "Mangrove Temperature ˚C", xlab = "Sand Flat Temperature ˚C") head(ConeHeat2) names(ConeHeat2)<-c("Treatment", "Respiration") aov(Respiration~Treatment, data=ConeHeat2) dataaovConeHeat2<-aov(Respiration~Treatment, data=ConeHeat2) dataaovConeHeat2 summary(aov(Respiration~Treatment, data=ConeHeat2)) TukeyHSD(dataaovConeHeat2) boxplot(ConeHeat2) boxplot(ConeHeat1, xlab="Treatment", ylab="Respiration (ppm Oxygen/g bodymass)") head(Cone11.11.3) names(Cone11.11.3)<-c("Sand", "Mangrove", "Marsh") var(Cone11.11.3$Sand_Temp) var(Cone11.11.3$Mangrove_Temp) var(Cone11.11.3$Marsh_Temp) boxplot(Cone11.11.3, main="Temperature by Habitat", face="bold", xlab="Habitat", ylab="Temperature (˚C)") times<-rep(c("Day 1 0000", "0030", "0100", "0130"), 48) head(Cone11.11.6) names(Cone11.11.6)<-c("Time", "Sand", "Mangrove", "Marsh") plot(Cone11.11.6$Marsh~Cone11.11.6$Time, type="l", col="red", lwd=3, axes=FALSE, xlab= NA, ylab=NA, ylim=c(24,50)) mtext(side=2, line=3, "Temperature (˚C)") axis(side=2) #mtext(side=1, line=3, "Time") axis(side=1, at=c(0, 24, 48, 72, 96, 120, 144, 168, 192), labels=c("Night 1", "Day 1", "Night 2", "Day 2", "Night 3", "Day 3", "Night 4", "Day 4", "Night 5"), las=2) lines(Cone11.11.6$Sand, col="blue", lwd=3) lines(Cone11.11.6$Mangrove, col="forestgreen", lwd=3) legend("topleft", c("Marsh Grass", "Mangrove", "Sand"), lty= c(1,1,1), col=c("red", "forestgreen", "blue"), box.lwd=0, y.intersp=1, lwd= 3)