# HG changeset patch # User greg # Date 1515095064 18000 # Node ID cf0a62d084dcd30a48be2530f497b3d395227875 # Parent dcd4be40706353937fcf3881f1fd167ce04a056d Uploaded diff -r dcd4be407063 -r cf0a62d084dc create_heatmap.R --- a/create_heatmap.R Thu Jan 04 12:31:46 2018 -0500 +++ b/create_heatmap.R Thu Jan 04 14:44:24 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=""); @@ -74,9 +74,9 @@ if(regexpr("ctcf", tolower(histone_marks[i])) > 0) { histone_mark_color[i,] = c(200, 0, 250); } - state_color = get_state_color(data_matrix, histone_mark_color)[,]; + state_color = get_state_color(data_matrix, histone_mark_color)[,2]; } - rect(rep(p+0.2, num_rows), 1:num_rows-0.8, rep(p+0.8, num_rows), 1:num_rows-0.2, col=state_color[,2]); + rect(rep(p+0.2, num_rows), 1:num_rows-0.8, rep(p+0.8, num_rows), 1:num_rows-0.2, col=state_color); palette(defpalette); if (!is.null(output_file_name)) { dev.off();