####0_Hemipterans per tree#### > g2=glm(Hemipterans_tree~factor(Treatment),family=quasipoisson) > summary(g2) > a=anova(g2,test="F") > a ####1_Parasitoid_density_Fig_2#### lme1=lme(aphytis_tree~factor(day)*treatment,random=~1|factor(tree),data=data) anova(lme1) #Check normality and homocedasticity# Res <- residuals(lme1, type="normalized") Fit <- fitted(lme1) par(mfrow=c(3,2)) plot(Res ~ Fit, xlab="Fitted values", ylab="Residuals", main="Residuals vs. fitted") abline(h=0) boxplot(Res ~ data$treatment, ylab="Residuals", main = "treatment") abline(h=0, lty=3) boxplot(Res ~ data$day, ylab="Residuals", main="day") abline(h=0, lty=3) hist(Res, main="Histogram of residuals", xlab="Residuals") qqnorm(Res) qqline(Res) #post test analysis# mmod.t.test <- glht(lme1, linfct=mcp(treatment="Tukey")) summary(mmod.t.test) confint(glht(lme1, linfct = mcp(treatment = "Tukey"))) TukeyHSD(lme1) ####2_Sugar_content_Figs 3 4### data_morning=subset(data,Time=="morning") data_evening=subset(data,Time=="evening") ####sugar content#### #morning# > lm1=lm(Total_sugars ~factor(Treatment), data=data_morning) > summary(lm1) > a1=anova(lm1) > a1 #evening# > lm2=lm(Total_sugars ~factor(Treatment), data=data_evening) > summary(lm2) > a2=anova(lm2) > a2 ####3_Fecundity_Fig_5#### data_morning=subset(data,Time=="Morning") data_evening=subset(data,Time=="Evening") #morning# > g1=glm(Eggs~factor(Treatment),family=quasipoisson, data=data_morning) > summary(g1) > a1=anova(g1,test="F") > a1 #evening# > g2=glm(Eggs~factor(Treatment),family=quasipoisson, data=data_evening) > summary(g2) > a2=anova(g2,test="F") > a2 ####4_Parasitism_Fig_6#### data=read.table("Parastism2.txt",h=T) attach(data) y=cbind(Parasit,Vivos) g=glm(y~Tratamiento,family=quasibinomial) summary(g) a=anova(g,test="F") a glht2 <-glht(g, mcp(Tratamiento="Tukey")) summary(glht2) g2=glm(Parasit~Tratamiento,family=quasipoisson) summary(g2) a2=anova(g2,test="F") a2 glht3 <-glht(g2, mcp(Tratamiento="Tukey")) summary(glht3) ####5_Migration_S1#### data=read.table("migration.txt",h=T) lm1=lm(aphytis_captures_tree~ *treatment) anova(lme1) ####6_Survivorship_S2#### > y=cbind(alive,dead) > g=glm(y~factor(Treatment),family=quasibinomial) > summary(g) > a=anova(g,test="F") > a