Mercurial > repos > vipints > rdiff
comparison rDiff/src/genes_cell2struct.m @ 0:0f80a5141704
version 0.3 uploaded
| author | vipints |
|---|---|
| date | Thu, 14 Feb 2013 23:38:36 -0500 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| -1:000000000000 | 0:0f80a5141704 |
|---|---|
| 1 function genes_cell2struct(anno_fname) | |
| 2 % GENES_CELL2STRUCT Converts genes stored as a cell to struct. | |
| 3 % | |
| 4 % genes_cell2struct(anno_fname) | |
| 5 % | |
| 6 % -- input -- | |
| 7 % anno_fname: name of file where genes as cell are stored | |
| 8 % | |
| 9 % -- output -- | |
| 10 % genes as a struct | |
| 11 % | |
| 12 % | |
| 13 % This program is free software; you can redistribute it and/or modify | |
| 14 % it under the terms of the GNU General Public License as published by | |
| 15 % the Free Software Foundation; either version 3 of the License, or | |
| 16 % (at your option) any later version. | |
| 17 % | |
| 18 % Written (W) 2009-2011 Regina Bohnert, Gunnar Raetsch | |
| 19 % Copyright (C) 2009-2011 Max Planck Society | |
| 20 % | |
| 21 | |
| 22 load(anno_fname, 'genes'); | |
| 23 if iscell(genes) | |
| 24 genes_cell = genes; | |
| 25 clear genes; | |
| 26 for g = 1:length(genes_cell), | |
| 27 gene = genes_cell{g}; | |
| 28 for e = 1:length(gene.exons) | |
| 29 gene.exons{e} = double(gene.exons{e}); | |
| 30 end | |
| 31 gene.exons = reshape(gene.exons, 1, length(gene.exons)); | |
| 32 gene.id = double(gene.id); | |
| 33 gene.start = double(gene.start); | |
| 34 gene.stop = double(gene.stop); | |
| 35 genes(g) = gene; | |
| 36 end | |
| 37 save(anno_fname, 'genes'); | |
| 38 end |