# *************************************************************************************************** # ** 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(lme4) require(MuMIn) library(car) d<-read.csv2(file="FE_datafile_Pcrit.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<-merge(d,traits) colnames(d)[9]<- "drop" d$drop<-as.numeric(levels(d$drop))[d$drop] d$Pcrit<-as.numeric(levels(d$Pcrit))[d$Pcrit] #plot S4 plot(100-d$drop~d$Pcrit,pch=16) points(100-drop~Pcrit, data=subset(d,d$Pref=="standing"),col="red",pch=16) points(100-drop~Pcrit, data=subset(d,d$Pref=="running"),col="blue",pch=16) points(c(15.4,30,50.4,26.5),c(63.2,63.7,74.8,58.7), pch=0,cex=2, col=c("red","blue","black","black")) #model m1<-lm(d$Pcrit~Pref* drop,data=d) summary(m1) shapiro.test(m1$res)#residuals are normally distributed m2<-lm(d$Pcrit~drop,data=d) leveneTest(m2$res~Pref,data=d)#equal variances across habitat preference categories