Examples of RnaSpecificSmoothCost

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
  }

  private static void testRNAShuffleExample() throws IOException {
    RnaSpecificTypeRelatedTreePruning prune2 = new RnaSpecificTypeRelatedTreePruning(
        false, true);
    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();

    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    final RNASpecificTree t = new RNASpecificTree("T");
    t.buildFromViennaFormat(
        "UUUAGCAGUUGGGAGCGCCUGAAGGGAGGUCCUNCGAUCGAAACCA"
            .toCharArray(),
        "((..((........))(((.....)))....(.......)))...."
            .toCharArray());
    final RNASpecificTree s = new RNASpecificTree("S");
    s.buildFromViennaFormat(
        "GAUAGAUGGUCUGGGUGUCGCCAGAUCGGGGUNCAAUUCCCCGUCGCCA"
            .toCharArray(),
        "((..(......)(((.....)))...(((((.......)))))))...."
            .toCharArray());
    Random rand = new Random(123l);
    final RNASpecificTree sR = new RNASpecificTree("S-Random");
//    sR.buildFromViennaFormat(
//        "GGGAGUAUGUACCGAUAGGGUGAUUCCGGGCCCUUCGACNGGCUCCAUC"
//            .toCharArray(),
//        "(((((............((......))(((.))).(.....)))))).."
//            .toCharArray());
    sR.buildFromViennaFormat(
        "GAUAGAUGGUCUGGGUGUCGCCAGAUCGGGGUNCAAUUCCCCGUCGCCA"
            .toCharArray(),
        "((..(......)(((.....)))...(((((.......)))))))...."
            .toCharArray());
    sR.shuffle(rand);
    sR.toFASTAFile("sR.fasta");
    sR.toDotFile("sR.dot", false);
    s.toDotFile("s.dot", false);
    t.toDotFile("t.dot", false);
    prune2.calculateCost(t);
    prune2.calculateCost(s);
    prune2.calculateCost(sR);
   
    smooth.calculateCost(t);
    smooth.calculateCost(s);
    smooth.calculateCost(sR);
   
   
   
    double ans1 = 0,ans2 = 0;
   
View Full Code Here

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
  public static void testRNAExample() {

    RnaSpecificTypeRelatedTreePruning prune = new RnaSpecificTypeRelatedTreePruning(
        false, true);
    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();
    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    RNASpecificTree t = new RNASpecificTree();
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    t.buildFromViennaFormat(
        "NNNUAUAGUUUGAGUUCGAUUGCGCUUCGUAUGUUGCGUCUACGUAAAAACGCUCAGUUUAAAUUAUAACUGCAAAAAAUAAUAACAAUUCUUACGCUUUAGCUGCCUAAUAAGCGCUUAACGUAGAUCCUCCCAGGAUCGUCCAUGUUCUGGAUCUGGGUCCUAAAUUUAGUGGACUUACGCUCAAAGCUUCCACCUGGAGUUGCGAGAAGAGACUAAUCAGGUUAGUCAUUGCUGGGUGCCCUGUCAUACGGCGUUUGCAAUGAUGAAAUUUAAAUAGUAUGAAUAUGAGCGUAGAUAUCCGAGGGGCAAUAUGCUUAGACGCNNN"
            .toCharArray(),
        "...((((.......((((((..(((((........((((..........))))...........................................(.((((.(((.......))).)))).)..(((((..(((((((.......)))))))....)))))....................(((..((((.........))))..)))..................(((((...(((.((((.........))))).))..)))))......................))))).....)).))))......))))............"
            .toCharArray());

    prune.calculateCost(t);
    smooth.calculateCost(t);

    RNASpecificTree s = new RNASpecificTree();
    // s.buildFromViennaFormat("GGACCGAUGGUAGUGUCUUCGGAUGCGAGAGUAGGUC".toCharArray(),
    // ".((((((((((((((((....))))))))))))))))".toCharArray());
    // s.buildFromViennaFormat("GAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAUCCCCUC".toCharArray(),
    // "(((((((......)))))))(((((((......)))))))(((((((......)))))))(((((((.....)))))))".toCharArray());
    // s.buildFromViennaFormat("CCGGgGGAUCACCACGGCGGgGGAUCACCACGG".toCharArray(),
    // "((((.((....))))))(((.((....))))).".toCharArray());
    s.buildFromViennaFormat(
        "GGCGACACGGAUUCCAGUGCAUAUCUUAGUGAUACUCCAGUUAACUCCAUACUUUCCCUGCAAUACGCUAUUCGCCUCAGAUGUAUUUGGGUGGUUGCUCCACUAAAGCCCAGGAAUAUCCAGCCAGUUACAUUUGAGGCCAUUUGGGCUUAAGCGUAUUCCAUGGAAAGUUUUCUCCCCACAUUUCGGAAAUUAAAUUCCGAGCCAGCAAGAAAAUCUUCUCUGUUACAAUUUGACAUGGCUAAAAACUGUACUAAUCAAAAUGAAAAAUGUUUCUCUUGGGCGUAAUCUCAUACAAUGAUUACCCUUAAAGAUCGAACAUUUAAACAAUAAUAUUUGAUAUGAUAUUUUCAAUUUCUAUGCUAUGCCAAAGUGUCUGACAUAAUCAAACAUUUGCACAUUCUUUGACCAAGAAUAGUCAGCAAAUUGUAUUUUCAAUCAAUGCAGACCAUAUGUUCCAGUUUCGGAGAUUUUUUGCUGCCAAACGGAAUACUUAUAAAAACCCACAUUCUAUUUACAUCACUAAGAAGAGCAUUGCAAUCUGUUUAGCC"
            .toCharArray(),
        "((((..(((((((.((((((...((((((((((.................(((((((.((.((((((((....((((((((((((.((((...((((.(((...........))).)).))..)))).))))))))))))............)))))))).)).)))))))........................(((((((...((((((((((.((((.((((..((((((((.((((((((.(((.((.....((((((.....((((((((.(((((((.(((((...........))))))))..))))..))))))))...........))))))..........)).))).)).)))))).))))(.(((((.(((.........((((((((((((.......))))).))..)))))...............)))))))))..))))..))))....)))))))))))))).....)))))))..........................))))))))))...)))))).)))))))..))))"
            .toCharArray());

    prune.calculateCost(s);
    smooth.calculateCost(s);

    TIntArrayList[] alignment = new TIntArrayList[3];
    MatcherFactory matcherFactory = new UnorderedMatcherFactory();
    // MatcherFactory matcherFactory = new OrderedMatcherFactory();
    HSA hsa = new HSA(cost, matcherFactory);
View Full Code Here

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
  @Test
  public void testRNAShuffleExample() {
    RnaSpecificTypeRelatedTreePruning prune = new RnaSpecificTypeRelatedTreePruning(
        false, true);
    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();

    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    final RNASpecificTree t = new RNASpecificTree("T");
    t.buildFromViennaFormat(
        "UUUAGCAGUUGGGAGCGCCUGAAGGGAGGUCCUNCGAUCGAAACCA".toCharArray(),
        "((..((........))(((.....)))....(.......)))....".toCharArray());
    final RNASpecificTree s = new RNASpecificTree("S");
    s.buildFromViennaFormat(
        "GAUAGAUGGUCUGGGUGUCGCCAGAUCGGGGUNCAAUUCCCCGUCGCCA"
            .toCharArray(),
        "((..(......)(((.....)))...(((((.......)))))))...."
            .toCharArray());
    Random rand = new Random(123l);
    final RNASpecificTree sR = new RNASpecificTree("S-Random");
    sR.buildFromViennaFormat(
        "GAUAGAUGGUCUGGGUGUCGCCAGAUCGGGGUNCAAUUCCCCGUCGCCA"
            .toCharArray(),
        "((..(......)(((.....)))...(((((.......)))))))...."
            .toCharArray());
    sR.shuffle(rand);
    prune.calculateCost(t);
    prune.calculateCost(s);
    prune.calculateCost(sR);

    smooth.calculateCost(t);
    smooth.calculateCost(s);
    smooth.calculateCost(sR);

    double ans1 = 0, ans2 = 0;

    // reuse
    HSA hsaEngine2 = new HSA(1, 1, cost, new UnorderedMatcherFactory(),
View Full Code Here

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
    RnaSpecificTypeRelatedTreePruning pruneRooted = new RnaSpecificTypeRelatedTreePruning(
        false, true);
    RnaSpecificTypeRelatedTreePruning pruneUnrooted = new RnaSpecificTypeRelatedTreePruning(
        false, false);

    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();

    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    final RNASpecificTree t = new RNASpecificTree("T");
    t.buildFromViennaFormat(
        "GGGGGGGGuuuAAAAAAAAAAcccUUUUUUUUUUuuuGGGGGGcccUUUUUUuuuCCCCCCCC"
            .toCharArray(),
        "((((((((...((((((((((...))))))))))...((((((...))))))...))))))))"
            .toCharArray());
    final RNASpecificTree s = new RNASpecificTree("S");
    s.buildFromViennaFormat(
        "UUUUUUUUUUuuuGGGGGGcccUUUUUUuuuCCCCCCCCcccGGGGGGGGuuuAAAAAAAAAA"
            .toCharArray(),
        "((((((((((...((((((...))))))...((((((((...))))))))...))))))))))"
            .toCharArray());

    smooth.calculateCost(t);
    smooth.calculateCost(s);

    double ans1 = 0, ans2 = 0;

    // reuse
    HSA hsaEngine = new HSA(1, 1, cost, new OrderedMatcherFactory(), true);
View Full Code Here

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
  @Test
  public void testSelfSim() {
    RnaSpecificTypeRelatedTreePruning prune = new RnaSpecificTypeRelatedTreePruning(
        true, false);
   
    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();

    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    final RNASpecificTree t = new RNASpecificTree("T");
    /*t.buildFromViennaFormat(
        "UUUAGCAGUUGGGAGCGCCUGAAGGGAGGUCCUNCGAUCGAAACCA".toCharArray(),
        "((..((........))(((.....)))....(.......)))....".toCharArray());
        */
   
    t.buildFromViennaFormat(
        "GCCUGAAGGGACGCCUGAAGGGA".toCharArray(),
        "(((.....))).(((.....)))".toCharArray());
   
    prune.calculateCost(t);
    smooth.calculateCost(t);
   
    try {
      t.toDotFile("/home/milon/tmp/1.dot", true);
    } catch (IOException e) {
      // TODO Auto-generated catch block
View Full Code Here

Examples of bgu.bio.adt.rna.costs.RnaSpecificSmoothCost

732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
  public static void testRNAExample() {

    RnaSpecificTypeRelatedTreePruning prune = new RnaSpecificTypeRelatedTreePruning(
        false, true);
    RnaSpecificSmoothCost smooth = new RnaSpecificSmoothCost();
    SequenceAlignment aligner = new AffineGapGlobalSequenceAlignmentNoMatrix(
        10, 10, RnaAlphabet.getInstance(), new AffineGapScoringMatrix(
            "matrix" + File.separator
                + "interval-RIBOSUM85-60.matrix",
            RnaAlphabet.getInstance()));

    RnaSpecificCostFunction cost = new RnaSpecificCostFunction(
        new ScoringMatrix("matrix" + File.separator
            + "bp-RIBOSUM85-60.matrix",
            RNABasePairAlphabet.getInstance()), aligner);

    RNASpecificTree t = new RNASpecificTree();
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    // t.buildFromViennaFormat("GAGCCGUAUGCGAUGAAAGUCGCACGUACGGUUC".toCharArray(),
    // "((((((((.((((((..))))))...))))))))".toCharArray());
    t.buildFromViennaFormat(
        "NNNUAUAGUUUGAGUUCGAUUGCGCUUCGUAUGUUGCGUCUACGUAAAAACGCUCAGUUUAAAUUAUAACUGCAAAAAAUAAUAACAAUUCUUACGCUUUAGCUGCCUAAUAAGCGCUUAACGUAGAUCCUCCCAGGAUCGUCCAUGUUCUGGAUCUGGGUCCUAAAUUUAGUGGACUUACGCUCAAAGCUUCCACCUGGAGUUGCGAGAAGAGACUAAUCAGGUUAGUCAUUGCUGGGUGCCCUGUCAUACGGCGUUUGCAAUGAUGAAAUUUAAAUAGUAUGAAUAUGAGCGUAGAUAUCCGAGGGGCAAUAUGCUUAGACGCNNN"
            .toCharArray(),
        "...((((.......((((((..(((((........((((..........))))...........................................(.((((.(((.......))).)))).)..(((((..(((((((.......)))))))....)))))....................(((..((((.........))))..)))..................(((((...(((.((((.........))))).))..)))))......................))))).....)).))))......))))............"
            .toCharArray());

    prune.calculateCost(t);
    smooth.calculateCost(t);

    RNASpecificTree s = new RNASpecificTree();
    // s.buildFromViennaFormat("GGACCGAUGGUAGUGUCUUCGGAUGCGAGAGUAGGUC".toCharArray(),
    // ".((((((((((((((((....))))))))))))))))".toCharArray());
    // s.buildFromViennaFormat("GAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAAUCCCCUCGAGGGGAUGAAAUCCCCUC".toCharArray(),
    // "(((((((......)))))))(((((((......)))))))(((((((......)))))))(((((((.....)))))))".toCharArray());
    // s.buildFromViennaFormat("CCGGgGGAUCACCACGGCGGgGGAUCACCACGG".toCharArray(),
    // "((((.((....))))))(((.((....))))).".toCharArray());
    s.buildFromViennaFormat(
        "GGCGACACGGAUUCCAGUGCAUAUCUUAGUGAUACUCCAGUUAACUCCAUACUUUCCCUGCAAUACGCUAUUCGCCUCAGAUGUAUUUGGGUGGUUGCUCCACUAAAGCCCAGGAAUAUCCAGCCAGUUACAUUUGAGGCCAUUUGGGCUUAAGCGUAUUCCAUGGAAAGUUUUCUCCCCACAUUUCGGAAAUUAAAUUCCGAGCCAGCAAGAAAAUCUUCUCUGUUACAAUUUGACAUGGCUAAAAACUGUACUAAUCAAAAUGAAAAAUGUUUCUCUUGGGCGUAAUCUCAUACAAUGAUUACCCUUAAAGAUCGAACAUUUAAACAAUAAUAUUUGAUAUGAUAUUUUCAAUUUCUAUGCUAUGCCAAAGUGUCUGACAUAAUCAAACAUUUGCACAUUCUUUGACCAAGAAUAGUCAGCAAAUUGUAUUUUCAAUCAAUGCAGACCAUAUGUUCCAGUUUCGGAGAUUUUUUGCUGCCAAACGGAAUACUUAUAAAAACCCACAUUCUAUUUACAUCACUAAGAAGAGCAUUGCAAUCUGUUUAGCC"
            .toCharArray(),
        "((((..(((((((.((((((...((((((((((.................(((((((.((.((((((((....((((((((((((.((((...((((.(((...........))).)).))..)))).))))))))))))............)))))))).)).)))))))........................(((((((...((((((((((.((((.((((..((((((((.((((((((.(((.((.....((((((.....((((((((.(((((((.(((((...........))))))))..))))..))))))))...........))))))..........)).))).)).)))))).))))(.(((((.(((.........((((((((((((.......))))).))..)))))...............)))))))))..))))..))))....)))))))))))))).....)))))))..........................))))))))))...)))))).)))))))..))))"
            .toCharArray());

    prune.calculateCost(s);
    smooth.calculateCost(s);

    TIntArrayList[] alignment = new TIntArrayList[3];
    MatcherFactory matcherFactory = new UnorderedMatcherFactory();
    // MatcherFactory matcherFactory = new OrderedMatcherFactory();
    HSA hsa = new HSA(cost, matcherFactory);
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.