# load data - I separated the two data types dR <- read.csv('FST_NeutSNPs.csv', header=T); dT <- read.csv('FST_ColorTraits.csv', header=T); quartz.options(height=5, width=6); plot.new(); par(mar=c(3,5,2,3)); plot.window(ylim=c(-0.1, 0.75), xlim=c(0,10)); axis(1, at=seq(-0.4, 10, by=0.1), labels=NA, lwd.ticks=0, lwd=1.5); axis(2, at=seq(-0.15, 0.75, by=0.01), labels=NA, lwd.ticks=0, lwd=1.5); axis(2, at=seq(0, 0.75, by=0.25), las=1, lwd=0, lwd.ticks=1); # ddRAD points(x=seq(0,10, length.out=nrow(dR)), y=dR$fst, col='grey85', pch=21, cex=0.5); lines(x=c(0,10), y=c(dT[1,2], dT[1,2]), col='grey50', lwd=2); # black lines(x=c(1,1), y=c(dT[4,3], dT[4,4]), col='black', lwd=1.7); points(x=1, y=dT[4,2], col='black', cex=1.4, pch=21, bg='black'); lines(x=c(0.9,1.1), y=c(dT[4,3], dT[4,3]), col='black', lwd=1.7); lines(x=c(0.9,1.1), y=c(dT[4,4], dT[4,4]), col='black', lwd=1.7); # red lines(x=c(3,3), y=c(dT[5,3], dT[5,4]), col='black', lwd=1.9); points(x=3, y=dT[5,2], col='black', cex=1.4, pch=21, bg='black'); lines(x=c(2.9,3.1), y=c(dT[5,3], dT[5,3]), col='black', lwd=1.7); lines(x=c(2.9,3.1), y=c(dT[5,4], dT[5,4]), col='black', lwd=1.7); # cap lines(x=c(7,7), y=c(dT[6,3], dT[6,4]), col='black', lwd=1.9); points(x=7, y=dT[6,2], col='black', cex=1.4, pch=21, bg='black'); lines(x=c(6.9,7.1), y=c(dT[6,3], dT[6,3]), col='black', lwd=1.7); lines(x=c(6.9,7.1), y=c(dT[6,4], dT[6,4]), col='black', lwd=1.7); # collar lines(x=c(9,9), y=c(dT[7,3], dT[7,4]), col='black', lwd=1.9); points(x=9, y=dT[7,2], col='black', cex=1.4, pch=21, bg='black'); lines(x=c(8.9,9.1), y=c(dT[7,3], dT[7,3]), col='black', lwd=1.7); lines(x=c(8.9,9.1), y=c(dT[7,4], dT[7,4]), col='black', lwd=1.7); mtext('Mimetic Traits', side=3, line=-1, at=2, cex=1.3); mtext('Non-mimetic Traits', side=3, line=-1, at=8, cex=1.3); mtext(bquote(italic(F)['ST']), side=2, line=2.7, cex=1.3) mtext('Cross-\nbands', side=1, at=1, line=1.5, cex=1.15) mtext('Red\nStripe', side=1, at=3, line=1.5, cex=1.15) mtext('Black\nCap', side=1, at=7, line=1.5, cex=1.15) mtext('Nuchal\nBand', side=1, at=9, line=1.5, cex=1.15) mtext('Neutral\nSNPs', side=4, line=-0.7, las=1, at=0.06, col='grey30', cex=1.15) mtext('A', side=3, at=-2.2, cex=2, font=2, line=0.1)