Mercurial > repos > vipints > deseq_hts
view deseq-hts_1.0/src/genes_cell2struct.m @ 0:94a108763d9e draft
deseq-hts version 1.0 wraps the DESeq 1.6.0
author | vipints |
---|---|
date | Wed, 09 May 2012 20:43:47 -0400 |
parents | |
children |
line wrap: on
line source
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