====== mismatch distribution ====== * file: lines_Mismatch.r * document page: 1 \\ **code:** read.table("D:/Heidi/Master/R_Daten/Mismatch/Mismatch_mt4.txt")-> mismatch4 mismatch4 attach(mismatch4) V1 -> Diff V2 -> Observed V3 -> Low.bound V4 -> Up.bound V5 -> ModelFreq. max(Up.bound) -> max1 plot(Diff, Observed, type="l", ylim=c(0,max1), xlab="differences", ylab="number", main="Mismatch distribution") lines(Diff, Low.bound, lty=2) lines(Diff, Up.bound, lty=2) legend("topleft", expression (Observed, CI[0.05]), lty=1:2, bty="n") detach(mismatch4)