tiff("Trivariate food web noah atkin.tif", 20, 14, units="cm", res=600, compression="lzw")
par(oma=c(2,4,2,0), mar=c(3,1.5,1.5,1.5), mfrow=c(3,5))
for (i in 1:14) {
NvMLinearRegressions(hengill.ord[[i]])
NvMSlope(hengill.ord[[i]])
NvMIntercept(hengill.ord[[i]])
 
slope1<- NvMSlope(hengill.ord[[i]])
slope1 <- round(slope1,digits=2)
slope1
 
PlotNvM(hengill.ord[[i]],
pch=16,
ylim=c(-2,10),
xlim=c(-8,4),
xlab=" ",
main= " ",
highlight.nodes=NULL,
)
title(main = (paste(CPS(hengill.ord[[i]])$title, " - ", CPS(hengill.ord[[i]])$temp, "°C")), line=0.75, cex=1.25)
mtext(paste("         slope =", slope1), line=-1, cex=0.75)
 
model1 <- NvMLinearRegressions(hengill.ord[[i]])
model1
PlotLinearModels(model1)
 
}
mtext(expression(paste("Log"["10"]*" total species abundance (m"^"-2"*")")), side = 2, line = 1, outer = TRUE, font=1, cex=1)
mtext(expression(paste("Log"["10"]*" mean species body mass (mg)"^""*"")), side = 1, line = 0, outer = TRUE, font=1, cex=1)
 
dev.off()

Task 4:  Reproduce "Food web properties versus temperature.tif"