diff --git a/analysis_gene_ontology.m b/analysis_gene_ontology.m index 5ffd00d..2b15110 100644 --- a/analysis_gene_ontology.m +++ b/analysis_gene_ontology.m @@ -11,7 +11,7 @@ % get go annotation files annot_files = dir('data/gene_ontology/go_*.txt'); -go_matrix = zeros(size(ref_go_terms, 1), size(annot_files, 1)); +go_matrix = ones(size(ref_go_terms, 1), size(annot_files, 1)); for i = 1:size(annot_files, 1) file = annot_files(i); fRead = fopen(sprintf('data/gene_ontology/%s',file.name),'r'); @@ -25,7 +25,10 @@ go_matrix(idx_ref, i) = p_values(idx_qry); end +log_matrix = log10(go_matrix); +log_matrix = abs(log_matrix); -go_table = array2table(go_matrix, 'RowNames', ref_go_terms, 'VariableNames', cellfun(@(s) strtok(s, '.') ,{annot_files.name}, 'UniformOutput', false)); -clustergram(go_matrix,'Standardize', 2,'cluster',1, 'ColumnLabels', cellfun(@(s) strtok(s, '.') ,{annot_files.name}, 'UniformOutput', false)); +go_table = array2table(log_matrix, 'RowNames', ref_go_terms, 'VariableNames', cellfun(@(s) strtok(s, '.'), {annot_files.name}, 'UniformOutput', false)); + +clustergram(log_matrix,'Standardize', 2,'cluster',1, 'ColumnLabels', cellfun(@(s) strtok(s, '.'), {annot_files.name}, 'UniformOutput', false)); diff --git a/figures/gene_ontology_clustergram.fig b/figures/gene_ontology_clustergram.fig index 02a011c..533748d 100644 Binary files a/figures/gene_ontology_clustergram.fig and b/figures/gene_ontology_clustergram.fig differ