Post date: Dec 10, 2014 3:43:59 PM
Since it's not appropriate to compare thetas among the four taxa, since the number of SNPs in each dataset is very different, I made a new figure with four plots instead of one:
theta_pt<-read.csv("theta_pt.csv",header=T)
theta_up<-read.csv("theta_up.csv",header=T)
theta_lo<-read.csv("theta_lo.csv",header=T)
pdf("fig_theta_9dec14.pdf",width=8,height=5)
par(mfrow=c(2,2))
plot(1:10,theta_pt$E_pt,pch=16,ylim=c(0.2,0.8),ylab=expression(paste(theta)),xlab="Population",axes=F)
title(main="A) Eurycea sp.",adj=0, cex.main = 0.9)
segments(1:10,theta_lo$E_pt,1:10,theta_up$E_pt,lwd=1.2)
axis(1,at=1:10,labels=as.character(theta_pt$Pop), cex.lab=0.8, cex.axis=0.8)
axis(2)
plot(1:10,theta_pt$H_pt,pch=16,ylim=c(0.2,0.8),ylab=expression(paste(theta)),xlab="Population",axes=F)
title(main="B) Heterelmis comalensis",adj=0, cex.main = 0.9)
segments(1:10,theta_lo$H_pt,1:10,theta_up$H_pt,lwd=1.2)
axis(1,at=1:10,labels=as.character(theta_pt$Pop), cex.lab=0.8, cex.axis=0.8)
axis(2)
plot(1:10,theta_pt$Sb_pt,pch=16,ylim=c(0.2,0.8),ylab=expression(paste(theta)),xlab="Population",axes=F)
title(main="C) Stygobromus pecki",adj=0, cex.main = 0.9)
segments(1:10,theta_lo$Sb_pt,1:10,theta_up$Sb_pt,lwd=1.2)
axis(1,at=1:10,labels=as.character(theta_pt$Pop), cex.lab=0.8, cex.axis=0.8)
axis(2)
plot(1:10,theta_pt$Sp_pt,pch=16,ylim=c(0.2,0.8),ylab=expression(paste(theta)),xlab="Population",axes=F)
title(main="D) Stygoparnus comalensis",adj=0, cex.main = 0.9)
segments(1:10,theta_lo$Sp_pt,1:10,theta_up$Sp_pt,lwd=1.2)
axis(1,at=1:10,labels=as.character(theta_pt$Pop), cex.lab=0.8, cex.axis=0.8)
axis(2)
dev.off()