• 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)