User Tools

Site Tools


haplotype_distance_between_within_two_populations

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
haplotype_distance_between_within_two_populations [2007/12/21 13:40] heidihaplotype_distance_between_within_two_populations [2008/07/22 13:31] (current) – external edit 127.0.0.1
Line 80: Line 80:
         mtext(side = 2, at =(3*Row/4), line = 2, text = "Population 1", cex=0.8)         mtext(side = 2, at =(3*Row/4), line = 2, text = "Population 1", cex=0.8)
 </code> </code>
 +
  
 ===== mixed data from several populations ===== ===== mixed data from several populations =====
Line 183: Line 184:
  image.plot(x,y,DistanceMatrixTogether, col=ColorRamp, main="haplotype distancematrix between/within populations",   image.plot(x,y,DistanceMatrixTogether, col=ColorRamp, main="haplotype distancematrix between/within populations", 
              xlab="Haplotype", ylab="Haplotype", axes = FALSE)              xlab="Haplotype", ylab="Haplotype", axes = FALSE)
-  contour(DistanceMatrixTogether, add = TRUE) +        contour(DistanceMatrixTogether, add = TRUE) 
-  axis(1, at = c(1:a), labels=c(Pop1,Pop2), cex.axis=0.6) +        axis(1, at = c(1:a), labels=c(Pop1,Pop2), cex.axis=0.6) 
-  axis(2, at = c(1:b), labels=c(Pop2[NROW(Pop2):1],Pop1[NROW(Pop1):1]), cex.axis=0.6)+        axis(2, at = c(1:b), labels=c(Pop2[NROW(Pop2):1],Pop1[NROW(Pop1):1]), cex.axis=0.6)
         box()           box()  
         lines(c(0,NROW(Pop1)+0.5),c(NROW(Pop2)+0.5,NROW(Pop2)+0.5), lwd=2)         lines(c(0,NROW(Pop1)+0.5),c(NROW(Pop2)+0.5,NROW(Pop2)+0.5), lwd=2)
Line 191: Line 192:
         mtext(side = 1, at =(Row/4), line = 2, text = "Population 1", cex=0.8)         mtext(side = 1, at =(Row/4), line = 2, text = "Population 1", cex=0.8)
         mtext(side = 1, at =(3*Row/4), line = 2, text = "Population 2", cex=0.8)         mtext(side = 1, at =(3*Row/4), line = 2, text = "Population 2", cex=0.8)
-  mtext(side = 2, at =(Row/4), line = 2, text = "Population 2", cex=0.8) +        mtext(side = 2, at =(Row/4), line = 2, text = "Population 2", cex=0.8) 
-  mtext(side = 2, at =(3*Row/4), line = 2, text = "Population 1", cex=0.8)+        mtext(side = 2, at =(3*Row/4), line = 2, text = "Population 1", cex=0.8)
 </code> </code>
  
  
  
haplotype_distance_between_within_two_populations.1198240845.txt.gz · Last modified: 2008/07/22 13:30 (external edit)