annotate deseq-hts_1.0/src/genes_cell2struct.m @ 0:94a108763d9e draft

deseq-hts version 1.0 wraps the DESeq 1.6.0

function genes_cell2struct(anno_fname)
% GENES_CELL2STRUCT   Converts genes stored as a cell to struct.
%
%   genes_cell2struct(anno_fname)
%
%   -- input --
%   anno_fname:   name of file where genes as cell are stored
%
%   -- output --
%   genes as a struct

load(anno_fname, 'genes');
if iscell(genes)
  genes_cell = genes;
  clear genes;
  for g = 1:length(genes_cell), 
    gene = genes_cell{g};
    for e = 1:length(gene.exons)
      gene.exons{e} = double(gene.exons{e});
    end    
    gene.exons = reshape(gene.exons, 1, length(gene.exons));
    gene.id = double(gene.id);
    gene.start = double(gene.start);
    gene.stop = double(gene.stop);
    genes(g) = gene;
  end
save(anno_fname, 'genes');
end