DBH_SIMULATIONS/R_script at master · LynchLab/DBH_SIMULATIONS · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
A=400
B=400
t<-read.csv("Data.csv", sep=" ")
jpeg(filename="figureS1.jpg",  width=A, height=B, units="px", pointsize=12)
plot(t$cor, t$err, ylab=expression(paste("Correlation between ", N[e], " and ", mu) ), xlab=expression(paste("% error in measurement of ", mu) ) )
dev.off()
jpeg(filename="figureS2.jpg",  width=A, height=B, units="px", pointsize=12)
plot(t$cor, t$err, ylab=expression(paste("Correlation between ", N[e], " and ", mu) ), xlab=expression(paste("Correlation between errors in measuring ", mu[bp], " and ", mu[in]) ) )
dev.off()
jpeg(filename="figureS3.jpg",  width=A, height=B, units="px", pointsize=12)
plot(t$cor, t$err, ylab=expression(paste("Effective genome size (", G[e], ") ) ), xlab=expression(paste("Silent site diversity (", pi, ") ) ) )
dev.off()