/matlab/get_training_sequences_intergenic_long.m

http://github.com/jcnossen/hmmgenefinder · Objective C · 14 lines · 13 code · 1 blank · 0 comment · 2 complexity · 63f980838867caf8b1e1b9061cb618b9 MD5 · raw file 

    

  1. % Returns sequences for training HMM (used to export data for GHMM)
    2. 

  2. % Input:
    3. 

  3. %       seq - sequence with annotated genes
    4. 

  4. % Output:
    5. 

  5. %       res - cell array with intergnenic regions
    6. 

  6. function [res] = get_training_sequences_intergenic_long(seq)
    7. 

  7.     % Just in case, but better do it when loading the genome
    8. 

  8.     seq.Sequence = upper(seq.Sequence);
    9. 

  9.     
    10. 

  10.     seq = HMM_Intergenic.get_intergenic(seq, true, true);
    11. 

  11.     seq = HMM_Intergenic.select_intergenic(seq, 10, Inf);
    12. 

  12.     seq = HMM.select_genes(seq, true, false, false);
    13. 

  13.     res = HMM.get_training_sequences(seq);
    14. 

  14. end
    15.