############################################################## ############################################################## ### ### "Offspring development and life-history variation in a ### water flea depends upon clone-specific integration of ### genetic, non-genetic and environmental cues" ### ### Harney, Paterson, & Plaistow ### ### R code 2/6: Plotting univariate reaction norms ### ############################################################## ############################################################## library(Hmisc) library(RColorBrewer) # Package Hmisc for summarize function, # Package RColourBrewer for sensible colours rm(list=ls(all=TRUE)) alldata<-read.table("C://~//Dmagna_ME_alltraits_summary.txt", header=TRUE) alldata=na.omit(alldata) names(alldata) MEcond=factor(alldata$MEcond) #food is coded both as the factor CEfood (A,B,C,D) and as a continuous variable (133,59, 26, 12) CEfood=factor(alldata$CEfood) foodcells=alldata$foodcells clone=factor(alldata$clone) cloMECE=factor(alldata$cloMECE) neon=alldata$bthsz prgr=alldata$grprecoef dugr=alldata$grdurcoef inst=alldata$preinst instjit=jitter(alldata$preinst) clutch1=alldata$c1num cth1jit=jitter(alldata$c1num) c1sz=alldata$c1sz c1ag=alldata$c1a c1aj=jitter(alldata$c1a) allc1<-data.frame(clone,cloMECE,CEfood,foodcells,MEcond,neon,prgr,dugr,inst, instjit,clutch1,cth1jit,c1sz,c1ag,c1aj) allc1=na.omit(allc1) ############################################################## ### ### 1.1 Calculate means and standard errors for ### neonate data using this summary function ### ############################################################## ### the function, used several times in this script ME_summary<-function(x) { xbar<-mean(x,na.rm=T) se<-sd(x,na.rm=T)/sqrt(length(na.omit(x))) n<-length(na.omit(x)) c(Mean=xbar,SE=se,n=n) } ### summarize mean and standard error information neocondclone<-summarize(allc1$neon, llist(allc1$clone,allc1$MEcond),function(x) ME_summary(x)) names <-c("HB5","LB5","HDN3","LDN3","HNs1","LNs1") neocondclone2<-data.frame(names,neocondclone) ### ### 1.2 Plot means and standard errors for neonate data ### ############################################################## HMcol<-brewer.pal(8,"Blues") LMcol<-brewer.pal(8,"Oranges") par(mfrow=c(1,1),mar=c(5,4.5,2,2),xpd=TRUE) x<-neocondclone2$names x2<-c(1,2,3,4,5,6) y <-neocondclone$allc1.neon yplus<-y+neocondclone$SE yminus<-y-neocondclone$SE plot(x2,y,type="n",xlim=c(.75,6.25), ylim=c(0.8,0.94) , ylab="Neonate length (mm)", xlab="Clone" , cex.axis=1,cex.lab=1.2 ,las=1,xaxt='n') par(xpd=TRUE) text(1.5,0.783, "B5", cex=1) text(3.5,0.783, "DKN1-3", cex=1) text(5.5,0.783, "Ness1",cex=1) par(xpd=FALSE) abline(v=2.5,lty=2,col="grey70") abline(v=4.5,lty=2,col="grey70") errbar(x2[neocondclone2$names=="HB5"],y[neocondclone2$names=="HB5"],yplus[neocondclone2$names=="HB5"],yminus[neocondclone2$names=="HB5"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="HB5"],y[neocondclone2$names=="HB5"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) errbar(x2[neocondclone2$names=="LB5"],y[neocondclone2$names=="LB5"],yplus[neocondclone2$names=="LB5"],yminus[neocondclone2$names=="LB5"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="LB5"],y[neocondclone2$names=="LB5"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) errbar(x2[neocondclone2$names=="HDN3"],y[neocondclone2$names=="HDN3"],yplus[neocondclone2$names=="HDN3"],yminus[neocondclone2$names=="HDN3"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="HDN3"],y[neocondclone2$names=="HDN3"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) errbar(x2[neocondclone2$names=="LDN3"],y[neocondclone2$names=="LDN3"],yplus[neocondclone2$names=="LDN3"],yminus[neocondclone2$names=="LDN3"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="LDN3"],y[neocondclone2$names=="LDN3"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) errbar(x2[neocondclone2$names=="HNs1"],y[neocondclone2$names=="HNs1"],yplus[neocondclone2$names=="HNs1"],yminus[neocondclone2$names=="HNs1"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="HNs1"],y[neocondclone2$names=="HNs1"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) errbar(x2[neocondclone2$names=="LNs1"],y[neocondclone2$names=="LNs1"],yplus[neocondclone2$names=="LNs1"],yminus[neocondclone2$names=="LNs1"],cex=0.01,add=TRUE) points(x2[neocondclone2$names=="LNs1"],y[neocondclone2$names=="LNs1"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) legend(0.75,0.94,c("High maternal environment","Low maternal environment"),col=c(HMcol[7],LMcol[7]),pch=c(21,21),cex=c(1,1),pt.lwd=c(2,2),pt.cex=2,pt.bg = c(HMcol[5],LMcol[5]),y.intersp=2,box.lty=1) ############################################################## ### ### 2.1 Calculate means and standard errors for ### Age at maturity and size at maturity ### ############################################################## c1sumsz<-summarize(allc1$c1sz, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) names <-c("HB5A","LB5A","HB5B","LB5B","HB5C","LB5C","HB5D","LB5D" ,"HDN3A","LDN3A","HDN3B","LDN3B","HDN3C","LDN3C","HDN3D","LDN3D" ,"HNs1A","LNs1A","HNs1B","LNs1B","HNs1C","LNs1C","HNs1D","LNs1D") c1sumag<-summarize(allc1$c1ag, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) c1szag2<-data.frame(names,c1sumsz,c1sumag) ### ### 2.2 Plot means and standard errors for age at maturity ### ############################################################## ######### # B5 ######### HMcol<-brewer.pal(8,"Blues") LMcol<-brewer.pal(8,"Oranges") par(mfrow=c(1,3),mar=c(6,6,2,2),xpd=FALSE) ############################### plot(foodcells,c1aj,type="n",xlim=c(0,140),ylim=c(4,17) , ylab="Age (days)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1szag2$allc1.foodcells.1 y <-c1szag2$allc1.c1ag yplus<-y+c1szag2$SE.1 yminus<-y-c1szag2$SE.1 segments(x[c1szag2$names=="HB5A"],y[c1szag2$names=="HB5A"],x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"],x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"],x[c1szag2$names=="HB5D"],y[c1szag2$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LB5A"],y[c1szag2$names=="LB5A"],x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"],x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"],x[c1szag2$names=="LB5D"],y[c1szag2$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],c1aj[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],c1aj[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],c1aj[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],c1aj[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],c1aj[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],c1aj[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],c1aj[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],c1aj[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HB5A"],y[c1szag2$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5D"],y[c1szag2$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5A"],y[c1szag2$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5D"],y[c1szag2$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,c1aj,type="n",xlim=c(0,140),ylim=c(4,17) , ylab="Age (days)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1szag2$names=="HDN3A"],y[c1szag2$names=="HDN3A"],x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"],x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"],x[c1szag2$names=="HDN3D"],y[c1szag2$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LDN3A"],y[c1szag2$names=="LDN3A"],x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"],x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"],x[c1szag2$names=="LDN3D"],y[c1szag2$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],c1aj[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],c1aj[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],c1aj[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],c1aj[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],c1aj[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],c1aj[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],c1aj[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],c1aj[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HDN3A"],y[c1szag2$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3D"],y[c1szag2$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3A"],y[c1szag2$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3D"],y[c1szag2$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,c1aj,type="n",xlim=c(0,140),ylim=c(4,17) , ylab="Age (days)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1szag2$names=="HNs1A"],y[c1szag2$names=="HNs1A"],x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"],x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"],x[c1szag2$names=="HNs1D"],y[c1szag2$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LNs1A"],y[c1szag2$names=="LNs1A"],x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"],x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"],x[c1szag2$names=="LNs1D"],y[c1szag2$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],c1aj[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],c1aj[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],c1aj[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],c1aj[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],c1aj[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],c1aj[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],c1aj[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],c1aj[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HNs1A"],y[c1szag2$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1D"],y[c1szag2$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1A"],y[c1szag2$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1D"],y[c1szag2$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ### ### 2.3 Plot means and standard errors for size at maturity ### ############################################################## ######### # B5 ######### plot(foodcells,c1sz,type="n",xlim=c(0,140),ylim=c(1.9,3.1) , ylab="Size (mm)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1szag2$allc1.foodcells.1 y <-c1szag2$allc1.c1sz yplus<-y+c1szag2$SE yminus<-y-c1szag2$SE segments(x[c1szag2$names=="HB5A"],y[c1szag2$names=="HB5A"],x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"],x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"],x[c1szag2$names=="HB5D"],y[c1szag2$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LB5A"],y[c1szag2$names=="LB5A"],x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"],x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"],x[c1szag2$names=="LB5D"],y[c1szag2$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],c1sz[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],c1sz[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],c1sz[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],c1sz[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],c1sz[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],c1sz[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],c1sz[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],c1sz[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HB5A"],y[c1szag2$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5B"],y[c1szag2$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5C"],y[c1szag2$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HB5D"],y[c1szag2$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5A"],y[c1szag2$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5B"],y[c1szag2$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5C"],y[c1szag2$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LB5D"],y[c1szag2$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,c1sz,type="n",xlim=c(0,140),ylim=c(1.9,3.1) , ylab="Size (mm)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1szag2$names=="HDN3A"],y[c1szag2$names=="HDN3A"],x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"],x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"],x[c1szag2$names=="HDN3D"],y[c1szag2$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LDN3A"],y[c1szag2$names=="LDN3A"],x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"],x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"],x[c1szag2$names=="LDN3D"],y[c1szag2$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],c1sz[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],c1sz[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],c1sz[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],c1sz[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],c1sz[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],c1sz[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],c1sz[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],c1sz[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HDN3A"],y[c1szag2$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3B"],y[c1szag2$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3C"],y[c1szag2$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HDN3D"],y[c1szag2$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3A"],y[c1szag2$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3B"],y[c1szag2$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3C"],y[c1szag2$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LDN3D"],y[c1szag2$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,c1sz,type="n",xlim=c(0,140),ylim=c(1.9,3.1) , ylab="Size (mm)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1szag2$names=="HNs1A"],y[c1szag2$names=="HNs1A"],x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"],x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"],x[c1szag2$names=="HNs1D"],y[c1szag2$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1szag2$names=="LNs1A"],y[c1szag2$names=="LNs1A"],x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"],x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"],x[c1szag2$names=="LNs1D"],y[c1szag2$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],c1sz[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],c1sz[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],c1sz[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],c1sz[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],c1sz[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],c1sz[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],c1sz[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],c1sz[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1szag2$names=="HNs1A"],y[c1szag2$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1B"],y[c1szag2$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1C"],y[c1szag2$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="HNs1D"],y[c1szag2$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1A"],y[c1szag2$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1B"],y[c1szag2$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1C"],y[c1szag2$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1szag2$names=="LNs1D"],y[c1szag2$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ############################################################## ### ### 3.1 Calculate means and standard errors for ### pre-maturation growth and during-maturation growth ### ############################################################## c1sumpg<-summarize(allc1$prgr, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) names <-c("HB5A","LB5A","HB5B","LB5B","HB5C","LB5C","HB5D","LB5D" ,"HDN3A","LDN3A","HDN3B","LDN3B","HDN3C","LDN3C","HDN3D","LDN3D" ,"HNs1A","LNs1A","HNs1B","LNs1B","HNs1C","LNs1C","HNs1D","LNs1D") c1sumdg<-summarize(allc1$dugr, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) c1prdu<-data.frame(names,c1sumpg,c1sumdg) ### ### 3.2 Plot means and standard errors for pre-mat. growth ### ############################################################## ######### # B5 ######### plot(foodcells,prgr,type="n",xlim=c(0,140),ylim=c(0.09,0.4) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1prdu$allc1.foodcells.1 y <-c1prdu$allc1.prgr yplus<-y+c1prdu$SE yminus<-y-c1prdu$SE segments(x[c1prdu$names=="HB5A"],y[c1prdu$names=="HB5A"],x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"],x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"],x[c1prdu$names=="HB5D"],y[c1prdu$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LB5A"],y[c1prdu$names=="LB5A"],x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"],x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"],x[c1prdu$names=="LB5D"],y[c1prdu$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],prgr[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],prgr[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],prgr[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],prgr[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],prgr[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],prgr[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],prgr[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],prgr[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HB5A"],y[c1prdu$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5D"],y[c1prdu$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5A"],y[c1prdu$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5D"],y[c1prdu$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,prgr,type="n",xlim=c(0,140),ylim=c(0.09,0.4) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1prdu$names=="HDN3A"],y[c1prdu$names=="HDN3A"],x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"],x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"],x[c1prdu$names=="HDN3D"],y[c1prdu$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LDN3A"],y[c1prdu$names=="LDN3A"],x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"],x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"],x[c1prdu$names=="LDN3D"],y[c1prdu$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],prgr[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],prgr[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],prgr[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],prgr[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],prgr[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],prgr[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],prgr[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],prgr[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HDN3A"],y[c1prdu$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3D"],y[c1prdu$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3A"],y[c1prdu$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3D"],y[c1prdu$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,prgr,type="n",xlim=c(0,140),ylim=c(0.09,0.4) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1prdu$names=="HNs1A"],y[c1prdu$names=="HNs1A"],x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"],x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"],x[c1prdu$names=="HNs1D"],y[c1prdu$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LNs1A"],y[c1prdu$names=="LNs1A"],x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"],x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"],x[c1prdu$names=="LNs1D"],y[c1prdu$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],prgr[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],prgr[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],prgr[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],prgr[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],prgr[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],prgr[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],prgr[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],prgr[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HNs1A"],y[c1prdu$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1D"],y[c1prdu$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1A"],y[c1prdu$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1D"],y[c1prdu$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ### ### 3.3 Plot means and standard errors for during-mat. growth ### ############################################################## ######### # B5 ######### plot(foodcells,dugr,type="n",xlim=c(0,140) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1prdu$allc1.foodcells.1 y <-c1prdu$allc1.dugr yplus<-y+c1prdu$SE.1 yminus<-y-c1prdu$SE.1 segments(x[c1prdu$names=="HB5A"],y[c1prdu$names=="HB5A"],x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"],x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"],x[c1prdu$names=="HB5D"],y[c1prdu$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LB5A"],y[c1prdu$names=="LB5A"],x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"],x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"],x[c1prdu$names=="LB5D"],y[c1prdu$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],dugr[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],dugr[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],dugr[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],dugr[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],dugr[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],dugr[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],dugr[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],dugr[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HB5A"],y[c1prdu$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5B"],y[c1prdu$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5C"],y[c1prdu$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HB5D"],y[c1prdu$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5A"],y[c1prdu$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5B"],y[c1prdu$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5C"],y[c1prdu$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LB5D"],y[c1prdu$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,dugr,type="n",xlim=c(0,140) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1prdu$names=="HDN3A"],y[c1prdu$names=="HDN3A"],x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"],x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"],x[c1prdu$names=="HDN3D"],y[c1prdu$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LDN3A"],y[c1prdu$names=="LDN3A"],x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"],x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"],x[c1prdu$names=="LDN3D"],y[c1prdu$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],dugr[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],dugr[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],dugr[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],dugr[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],dugr[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],dugr[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],dugr[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],dugr[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HDN3A"],y[c1prdu$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3B"],y[c1prdu$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3C"],y[c1prdu$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HDN3D"],y[c1prdu$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3A"],y[c1prdu$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3B"],y[c1prdu$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3C"],y[c1prdu$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LDN3D"],y[c1prdu$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,dugr,type="n",xlim=c(0,140) , ylab="Growth rate (mm/day)" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1prdu$names=="HNs1A"],y[c1prdu$names=="HNs1A"],x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"],x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"],x[c1prdu$names=="HNs1D"],y[c1prdu$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1prdu$names=="LNs1A"],y[c1prdu$names=="LNs1A"],x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"],x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"],x[c1prdu$names=="LNs1D"],y[c1prdu$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],dugr[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],dugr[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],dugr[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],dugr[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],dugr[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],dugr[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],dugr[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],dugr[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1prdu$names=="HNs1A"],y[c1prdu$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1B"],y[c1prdu$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1C"],y[c1prdu$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="HNs1D"],y[c1prdu$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1A"],y[c1prdu$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1B"],y[c1prdu$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1C"],y[c1prdu$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1prdu$names=="LNs1D"],y[c1prdu$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ############################################################## ### ### 4.1 Calculate means and standard errors for ### No. pre-maturation instars and No. clutch 1 offspring ### ############################################################## c1suminstr<-summarize(allc1$inst, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) names <-c("HB5A","LB5A","HB5B","LB5B","HB5C","LB5C","HB5D","LB5D" ,"HDN3A","LDN3A","HDN3B","LDN3B","HDN3C","LDN3C","HDN3D","LDN3D" ,"HNs1A","LNs1A","HNs1B","LNs1B","HNs1C","LNs1C","HNs1D","LNs1D") c1sumclu<-summarize(allc1$clutch1, llist(allc1$clone,allc1$CEfood,allc1$MEcond,allc1$foodcells),function(x) ME_summary(x)) c1InsClu<-data.frame(names,c1suminstr,c1sumclu) ### ### 4.2 Plot means and standard errors for No. pre-mat. instars ### ############################################################### ######### # B5 ######### plot(foodcells,instjit,type="n",xlim=c(0,140)#,ylim=c(1.9,3.1) , ylab="Pre-maturation instars" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")), cex.axis=1,cex.lab=1.2 , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1InsClu$allc1.foodcells y <-c1InsClu$allc1.inst yplus<-y+c1InsClu$SE yminus<-y-c1InsClu$SE segments(x[c1InsClu$names=="HB5A"],y[c1InsClu$names=="HB5A"],x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"],x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"],x[c1InsClu$names=="HB5D"],y[c1InsClu$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LB5A"],y[c1InsClu$names=="LB5A"],x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"],x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"],x[c1InsClu$names=="LB5D"],y[c1InsClu$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],instjit[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],instjit[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],instjit[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],instjit[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],instjit[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],instjit[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],instjit[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],instjit[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HB5A"],y[c1InsClu$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5D"],y[c1InsClu$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5A"],y[c1InsClu$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5D"],y[c1InsClu$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,inst,type="n",xlim=c(0,140) , ylab="Pre-maturation instars" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1InsClu$names=="HDN3A"],y[c1InsClu$names=="HDN3A"],x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"],x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"],x[c1InsClu$names=="HDN3D"],y[c1InsClu$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LDN3A"],y[c1InsClu$names=="LDN3A"],x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"],x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"],x[c1InsClu$names=="LDN3D"],y[c1InsClu$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],instjit[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],instjit[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],instjit[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],instjit[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],instjit[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],instjit[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],instjit[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],instjit[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HDN3A"],y[c1InsClu$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3D"],y[c1InsClu$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3A"],y[c1InsClu$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3D"],y[c1InsClu$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,inst,type="n",xlim=c(0,140) , ylab="Pre-maturation instars" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1InsClu$names=="HNs1A"],y[c1InsClu$names=="HNs1A"],x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"],x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"],x[c1InsClu$names=="HNs1D"],y[c1InsClu$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LNs1A"],y[c1InsClu$names=="LNs1A"],x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"],x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"],x[c1InsClu$names=="LNs1D"],y[c1InsClu$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],instjit[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],instjit[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],instjit[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],instjit[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],instjit[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],instjit[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],instjit[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],instjit[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HNs1A"],y[c1InsClu$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1D"],y[c1InsClu$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1A"],y[c1InsClu$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1D"],y[c1InsClu$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ### ### 4.3 Plot means and standard errors for No. cluctch 1 offspring ### ################################################################### ######### # B5 ######### plot(foodcells,cth1jit,type="n",xlim=c(0,140)#,ylim=c(1.9,3.1) , ylab="Number of offspring" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")), cex.axis=1,cex.lab=1.2 , cex.axis=1.3,cex.lab=1.3,las=1) x <-c1InsClu$allc1.foodcells y <-c1InsClu$allc1.clutch1 yplus<-y+c1InsClu$SE.1 yminus<-y-c1InsClu$SE.1 segments(x[c1InsClu$names=="HB5A"],y[c1InsClu$names=="HB5A"],x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"],x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"],x[c1InsClu$names=="HB5D"],y[c1InsClu$names=="HB5D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LB5A"],y[c1InsClu$names=="LB5A"],x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"],x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"],x[c1InsClu$names=="LB5D"],y[c1InsClu$names=="LB5D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="B5HFCA"],cth1jit[allc1$cloMECE=="B5HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCB"],cth1jit[allc1$cloMECE=="B5HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCC"],cth1jit[allc1$cloMECE=="B5HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5HFCD"],cth1jit[allc1$cloMECE=="B5HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCA"],cth1jit[allc1$cloMECE=="B5LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCB"],cth1jit[allc1$cloMECE=="B5LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCC"],cth1jit[allc1$cloMECE=="B5LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="B5LFCD"],cth1jit[allc1$cloMECE=="B5LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HB5A"],y[c1InsClu$names=="HB5A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5B"],y[c1InsClu$names=="HB5B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5C"],y[c1InsClu$names=="HB5C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HB5D"],y[c1InsClu$names=="HB5D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5A"],y[c1InsClu$names=="LB5A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5B"],y[c1InsClu$names=="LB5B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5C"],y[c1InsClu$names=="LB5C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LB5D"],y[c1InsClu$names=="LB5D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # DKN1-3 ######### plot(foodcells,cth1jit,type="n",xlim=c(0,140) , ylab="Number of offspring" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1InsClu$names=="HDN3A"],y[c1InsClu$names=="HDN3A"],x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"],x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"],x[c1InsClu$names=="HDN3D"],y[c1InsClu$names=="HDN3D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LDN3A"],y[c1InsClu$names=="LDN3A"],x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"],x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"],x[c1InsClu$names=="LDN3D"],y[c1InsClu$names=="LDN3D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="DN3HFCA"],cth1jit[allc1$cloMECE=="DN3HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCB"],cth1jit[allc1$cloMECE=="DN3HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCC"],cth1jit[allc1$cloMECE=="DN3HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3HFCD"],cth1jit[allc1$cloMECE=="DN3HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCA"],cth1jit[allc1$cloMECE=="DN3LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCB"],cth1jit[allc1$cloMECE=="DN3LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCC"],cth1jit[allc1$cloMECE=="DN3LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="DN3LFCD"],cth1jit[allc1$cloMECE=="DN3LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HDN3A"],y[c1InsClu$names=="HDN3A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3B"],y[c1InsClu$names=="HDN3B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3C"],y[c1InsClu$names=="HDN3C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HDN3D"],y[c1InsClu$names=="HDN3D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3A"],y[c1InsClu$names=="LDN3A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3B"],y[c1InsClu$names=="LDN3B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3C"],y[c1InsClu$names=="LDN3C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LDN3D"],y[c1InsClu$names=="LDN3D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2) ######### # Ness1 ######### plot(foodcells,cth1jit,type="n",xlim=c(0,140) , ylab="Number of offspring" , xlab=expression(paste("Concentration of food (cells/",mu,"l)")) , cex.axis=1.3,cex.lab=1.3,las=1) segments(x[c1InsClu$names=="HNs1A"],y[c1InsClu$names=="HNs1A"],x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"], col = HMcol[6], lwd = 2) segments(x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"],x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"], col = HMcol[5], lwd = 2) segments(x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"],x[c1InsClu$names=="HNs1D"],y[c1InsClu$names=="HNs1D"], col = HMcol[4], lwd = 2) segments(x[c1InsClu$names=="LNs1A"],y[c1InsClu$names=="LNs1A"],x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"], col = LMcol[6], lwd = 2) segments(x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"],x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"], col = LMcol[5], lwd = 2) segments(x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"],x[c1InsClu$names=="LNs1D"],y[c1InsClu$names=="LNs1D"], col = LMcol[4], lwd = 2) points(foodcells[allc1$cloMECE=="Ns1HFCA"],cth1jit[allc1$cloMECE=="Ns1HFCA"],col=HMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCB"],cth1jit[allc1$cloMECE=="Ns1HFCB"],col=HMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCC"],cth1jit[allc1$cloMECE=="Ns1HFCC"],col=HMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1HFCD"],cth1jit[allc1$cloMECE=="Ns1HFCD"],col=HMcol[4],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCA"],cth1jit[allc1$cloMECE=="Ns1LFCA"],col=LMcol[7],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCB"],cth1jit[allc1$cloMECE=="Ns1LFCB"],col=LMcol[6],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCC"],cth1jit[allc1$cloMECE=="Ns1LFCC"],col=LMcol[5],pch=1,cex=1) points(foodcells[allc1$cloMECE=="Ns1LFCD"],cth1jit[allc1$cloMECE=="Ns1LFCD"],col=LMcol[4],pch=1,cex=1) points(x[c1InsClu$names=="HNs1A"],y[c1InsClu$names=="HNs1A"],col=HMcol[7], bg=HMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1B"],y[c1InsClu$names=="HNs1B"],col=HMcol[6], bg=HMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1C"],y[c1InsClu$names=="HNs1C"],col=HMcol[5], bg=HMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="HNs1D"],y[c1InsClu$names=="HNs1D"],col=HMcol[4], bg=HMcol[2], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1A"],y[c1InsClu$names=="LNs1A"],col=LMcol[7], bg=LMcol[5], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1B"],y[c1InsClu$names=="LNs1B"],col=LMcol[6], bg=LMcol[4], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1C"],y[c1InsClu$names=="LNs1C"],col=LMcol[5], bg=LMcol[3], pch=21,cex=2,lwd=2) points(x[c1InsClu$names=="LNs1D"],y[c1InsClu$names=="LNs1D"],col=LMcol[4], bg=LMcol[2], pch=21,cex=2,lwd=2)