Mercurial > repos > greg > ideas_genome_tracks
changeset 81:dcd4be407063 draft
Uploaded
| author | greg |
|---|---|
| date | Thu, 04 Jan 2018 12:31:46 -0500 |
| parents | 6d97ef6bf7b7 |
| children | cf0a62d084dc |
| files | create_heatmap.R |
| diffstat | 1 files changed, 3 insertions(+), 3 deletions(-) [+] |
line wrap: on
line diff
--- a/create_heatmap.R Thu Jan 04 12:31:30 2018 -0500 +++ b/create_heatmap.R Thu Jan 04 12:31:46 2018 -0500 @@ -7,7 +7,7 @@ num_columns = dim(data_frame)[2]; num_rows = dim(data_frame)[1]; p = (sqrt(9 + 8 * (num_columns-1)) - 3) / 2; - data_matrix = as.matrix(data_frame[,1+1:p] / data_frame[,1]); + data_matrix = as.matrix(data_frame[,1+1:p]/data_frame[,1]); colnames(data_matrix) = colnames(data_frame)[1+1:p]; histone_marks = colnames(data_matrix); rownames(data_matrix) = paste(1:num_rows-1, " (", round(data_frame[,1]/sum(data_frame[,1])*10000)/100, "%)", sep=""); @@ -26,8 +26,8 @@ plot(NA, NA, xlim=c(0, p+0.7), ylim=c(0, num_rows), xaxt="n", yaxt="n", xlab=NA, ylab=NA, frame.plot=F); axis(1, at=1:p-0.5, labels=colnames(data_matrix), las=2); axis(4, at=1:num_rows-0.5, labels=rownames(data_matrix), las=2); - color = round((t(data_matrix) - min_max_vector[1]) / (min_max_vector[2] - min_max_vector[1]) * 100); - rect(rep(1:p-1, num_rows), rep(1:num_rows-1, each=p), rep(1:p, num_rows), rep(1:num_rows, each=p), col=color); + col = round((t(data_matrix) - min_max_vector[1]) / (min_max_vector[2] - min_max_vector[1]) * 100); + rect(rep(1:p-1, num_rows), rep(1:num_rows-1, each=p), rep(1:p, num_rows), rep(1:num_rows, each=p), col=col); histone_mark_color = t(col2rgb(terrain.colors(ceiling(p))[1:p])); # Specify a color for common feature names like "h3k4me3".