# FIGURE 1 # Melanin-photo regression from DNA Damage Experiment data<-read.csv("MelaninPhotoRegression_DNAdamageExpt.csv") head(data) model1<-lm(melanin.per.length~photo.grey,data=data) summary(model1) # Multiple R-squared: 0.435 # F-statistic: 35.36 on 1 and 46 DF, p-value: 3.479e-07 # Melanin-photo regression from UV Tolerance Experiment data2<-read.csv("MelaninPhotoRegression_UVtoleranceExpt.csv") model2<-lm(melanin.per.length~photo.grey,data=data2) summary(model2) # F1,57 = 20.42, p = 3.18e-05, R2 = 0.2509 #### FIGURE 1 PLOTTED quartz(width=7.12,height=10.5) par(mfrow=c(2,1)) par(las=1,bty="l",mar=c(1,5,1,1)) col1=rgb(213,94,0,maxColorValue=255) #Vermillion col2=rgb(0,114,178,maxColorValue=255) #Blue # DNA Damage Expt par(mgp=c(3,1,0),cex.axis=1.5,cex.lab=2.0) plot(melanin.per.length~photo.grey,data=data, ylab=expression(paste("Assayed melanin (",mu,"g/mm)")), pch=16,col=col2, cex=1.5, ylim=c(0.2,0.7), xlim=c(45,138), xlab="" # xlab="Mean pixel grey value", ) abline(coef(model1),col=col2,lwd=3,lty=2) text(50,0.7,"(a)",font=3,cex=2.0,xpd=NA) ### UV Tolerance Expt par(mar=c(5,5,1,1)) par(mgp=c(3,1,0),cex.axis=1.5,cex.lab=2.0) plot(melanin.per.length~photo.grey,data=data2, xlab="Mean pixel grey value", ylab=expression(paste("Assayed melanin (",mu,"g/mm)")), pch=16,col=col2, cex=1.5, ylim=c(0.2,0.7), xlim=c(45,138) ) abline(coef(model2),col=col2,lwd=3,lty=2) text(50,0.7,"(b)",font=3,cex=2.0,xpd=NA)