fst_matrix
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
fst_matrix [2007/12/21 13:56] – heidi | fst_matrix [2008/07/22 13:31] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 114: | Line 114: | ||
box() | box() | ||
</ | </ | ||
+ | |||
===== read data from a XML file ===== | ===== read data from a XML file ===== | ||
Line 206: | Line 207: | ||
axis(1, at = c(1: | axis(1, at = c(1: | ||
axis(2, at = c(1:b), labels=c(b: | axis(2, at = c(1:b), labels=c(b: | ||
+ | box() | ||
+ | </ | ||
+ | |||
+ | ===== move legend ===== | ||
+ | * file: matrix_PairwiseFst.r | ||
+ | * document page: 22 | ||
+ | |||
+ | \\ | ||
+ | **code:** | ||
+ | < | ||
+ | #----read data-------------------------------------------------------------- | ||
+ | read.table(" | ||
+ | header=TRUE, | ||
+ | as.matrix.data.frame(Data) -> Matrix | ||
+ | |||
+ | a <- ncol(Matrix) | ||
+ | b <- nrow(Matrix) | ||
+ | |||
+ | x <- c(1:a) | ||
+ | y <- c(1:b) | ||
+ | |||
+ | #----draw plot-------------------------------------------------------------- | ||
+ | library(fields) | ||
+ | ColorRamp <- rgb( seq(1, | ||
+ | seq(1, | ||
+ | seq(1, | ||
+ | |||
+ | #----Mirror matrix (left-right)---- | ||
+ | mirror.matrix <- function(x) { | ||
+ | xx <- as.data.frame(x); | ||
+ | xx <- rev(xx); | ||
+ | xx <- as.matrix(xx); | ||
+ | xx; | ||
+ | } | ||
+ | |||
+ | #----Rotate matrix 270 clockworks---- | ||
+ | rotate270.matrix <- function(x) { | ||
+ | mirror.matrix(t(x)) | ||
+ | } | ||
+ | |||
+ | Matrix <- rotate270.matrix(Matrix) | ||
+ | |||
+ | #----draw matrix plot---- | ||
+ | image.plot(x, | ||
+ | xlab=" | ||
+ | cex=1.0, side=4, line=2), axes = FALSE) | ||
+ | contour(Matrix, | ||
+ | axis(1, at = c(1:a)) | ||
+ | axis(2, at = c(1:b), labels=c(b: | ||
box() | box() | ||
</ | </ |
fst_matrix.1198241781.txt.gz · Last modified: 2008/07/22 13:30 (external edit)