# *************************************************************************************************** # ** R-scrit to analyse the data on respiratory responses to hypoxia ** # ** Part of manuscript: "Thermal limits in native and alien freshwater peracarid Crustacea: ** # ** the role of habitat use and oxygen limitation"by: WILCO C.E.P. VERBERK, ROB S.E.W. LEUVEN, ** # ** GERARD VAN DER VELDE, and FRIEDERIKE GABEL ** # *************************************************************************************************** require(nlme) d<-read.csv2(file="FE_datafile_Ea.csv", head=TRUE, sep=",", dec=",") traits<-data.frame(Species=c("Asellus aquaticus","Crangonyx pseudogracilis","Dikerogammarus villosus","Echinogammarus berilloni","Gammarus fossarum","Gammarus pulex","Gammarus roeselii","Gammarus tigrinus"),Pref=c("standing","standing","running","indifferent","running","indifferent","indifferent","indifferent"), Origin=c("N","E","E","E","N","N","N","E"), Prefnum=c(1,1,3,2,3,2,2,2)) d$Species[which(d$Species=="Gammarus pullex")]<-"Gammarus pulex" d<-merge(d,traits) #sorting class for numeric factors d$high.temp<-as.numeric(levels(d$high.temp))[d$high.temp] d$Test.temperature<-as.numeric(levels(d$Test.temperature))[d$Test.temperature] d$Ea..eV.<-as.numeric(levels(d$Ea..eV.))[d$Ea..eV.] e<-(d[!is.na(d$Ea..eV.),]) #Figure S3 boxplot(Ea..eV.~Pref,data=e) # tabulation for table 3 aggregate(e[, 14], list(e$Species), mean) aggregate(e[, 14], list(e$Species), sd) table(e$Species) #model mix<-lme(Ea..eV.~Pref, random=~1|Source, method="ML",data=e) mix<-lme(Ea..eV.~Origin, random=~1|Source, method="ML",data=e) summary(mix) shapiro.test(e$Ea..eV.)#Ea values are normally distributed leveneTest(Ea..eV.~Pref,data=e)#equal variances across habitat preference categories leveneTest(Ea..eV.~Origin,data=e)#equal variances across geographic origin categories