Mercurial > repos > vipints > rdiff

diff rDiff/src/tools/compute_testing_region.m @ 0:0f80a5141704

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
version 0.3 uploaded
author vipints
date Thu, 14 Feb 2013 23:38:36 -0500
parents
children
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/rDiff/src/tools/compute_testing_region.m	Thu Feb 14 23:38:36 2013 -0500
@@ -0,0 +1,17 @@
+function [genes]=compute_testing_region(CFG,genes)
+
+%Iterate over gene
+
+for i=1:size(genes,2)
+    [SPLICINGEVENTS,SEQUENCE,EXONSEQUENCE]=splicingsequence(genes(i));
+    genes(i).splicingevents=SPLICINGEVENTS;
+    genes(i).sequence=SEQUENCE;
+    genes(i).exonsequence=EXONSEQUENCE;
+    
+    [UNIQUE_NEW_EXONS,GRAPHNODES,ORDER_OF_GRAPHNODE,EIRS_IN_SEQ]=transform_single_end_reads(genes(i),CFG.sequenced_length-CFG.bases_to_clip*2);
+    genes(i).unique_new_exons=UNIQUE_NEW_EXONS;
+    genes(i).graphnodes=GRAPHNODES;
+    genes(i).order_of_graphnodes=ORDER_OF_GRAPHNODE;
+    genes(i).eirs_in_seq=EIRS_IN_SEQ;
+ 
+end