rm(list=ls())
library(cheddar)
library(nortest)
setwd("~/Practical2/Practical exercises")
hengill = LoadCollection("Hengill Food Webs")
stats= CollectionCPS(hengill,c(Stream='title',
Temp='temp',
SpeciesRichness='NumberOfNodes',
LinkRichness= 'NumberOfTrophicLinks',
LinkageDensity = 'LinkageDensity',
Connectance='DirectedConnectance',
PropBasal='FractionBasalNodes',
PropInter='FractionIntermediateNodes',
PropTop='FractionTopLevelNodes'))
stats$MeanFCL= c((mean(ShortWeightedTrophicLevel(hengill[[1]]))),
(mean(ShortWeightedTrophicLevel(hengill[[2]]))),
(mean(ShortWeightedTrophicLevel(hengill[[3]]))),
(mean(ShortWeightedTrophicLevel(hengill[[4]]))),
(mean(ShortWeightedTrophicLevel(hengill[[5]]))),
(mean(ShortWeightedTrophicLevel(hengill[[6]]))),
(mean(ShortWeightedTrophicLevel(hengill[[7]]))),
(mean(ShortWeightedTrophicLevel(hengill[[8]]))),
(mean(ShortWeightedTrophicLevel(hengill[[9]]))),
(mean(ShortWeightedTrophicLevel(hengill[[10]]))),
(mean(ShortWeightedTrophicLevel(hengill[[11]]))),
(mean(ShortWeightedTrophicLevel(hengill[[12]]))),
(mean(ShortWeightedTrophicLevel(hengill[[13]]))),
(mean(ShortWeightedTrophicLevel(hengill[[14]])))
)
stats$TempCat[stats$Temp < 10] = "cold"
stats$TempCat[stats$Temp > 10] = "inter"
stats$TempCat[stats$Temp > 20] = "warm"
sort.stats <- stats[order(stats$Temp) , ]
write.csv(sort.stats, file="stats.csv", row.names = FALSE)
Task 2: Reproduce "Triangular food webs.tif"