hengill.ord <- hengill[order(stats$Temp)]
tiff("Triangular food webs final.tif", 20, 12, units="cm", res=600, compression="lzw")
par(oma=c(0,4,0,0), mar=c(1,1,2,2), mfrow=c(3,5), font=2)
for(i in 1:14)
{
PlotWebByLevel(hengill[[i]],
pch=16,
highlight.nodes=NULL,
x.layout='narrow',
main= paste(CPS(hengill.ord[[i]])$title, " - ", CPS(hengill.ord[[i]])$temp, "°C"),
cex.main=1.25,
font=2,
level="ShortWeightedTrophicLevel",)
}
mtext("Trophic height", side = 2, line = 1, outer = TRUE, font=1, cex=1)
dev.off()
Task 3: Reproduce "Food web properties versus temperature.tif"