Examples of GenomicRegion

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
    GenomicRegionsBuilder gb= new BedReader(file).getGenomicRegionsBuilder();
    GenomicRegionsBuilder gbc= new BedReader(chr).getGenomicRegionsBuilder();
    HashMap<String, GenomicRegion> cr= new HashMap<>();
    Iterator<GenomicRegion> it= gbc.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        cr.put(gr.chr, gr);
    }
    it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        GenomicRegion chrr= cr.get(gr.chr);
        if(chrr!=null){
            double dl= Math.abs(gr.low - chrr.low);
            double dh= Math.abs(gr.high - chrr.high);
            double d= Math.min(dl, dh);
            if(d<thr){
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

94
95
96
97
98
99
100
101
102
103
104
105
public void writeToFile(File out) throws Exception{
    GenomicRegionsBuilder gbout= new GenomicRegionsBuilder();
    for( int i=0; i<ress.size() ; i++ ){
        LandscapeScoringResultEntry r=ress.get(i);
        GenomicRegion ge= r.gr;
        Object o= ge.attr==null ? ge.toUCSCFormat() : ge.attr;
        ge.attr= new Object[]{o,
            new Double(r.scores[0]), new Double(r.scores[1])
        };
        gbout.add(ge);
    }
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

44
45
46
47
48
49
50
51
52
53
54
55
    while(it.hasNext()){
        allrs.add( it.next() );
    }
    it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        if(allrs.contains(gr)){
            gr.setValue(1);
        }
    }
    new BedExporter(gb).writeToFile(out);
}
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
                Landscape2DBuffer tsr= new Landscape2DBuffer(lbr.getLandscape2DForChromosome(chr));
               
               
                ArrayList<Region> regs= gb.getRegionsForChromosome(chr);
                for( int j=0; j<regs.size() ; j++ ){
                    GenomicRegion gr= (GenomicRegion)regs.get(j);
                    GenomicRegion grlow, grhigh;
                    int low, high;
                    low= (int)gr.low;
                    high= (int)gr.high;
                    grlow= new GenomicRegion(gr.chr, low-side, low+side);
                    grhigh= new GenomicRegion(gr.chr, high-side, high+side);

                   
                    double forl= score(tsf.reloadAndGetList(grlow), grlow);
                    double forh= score(tsf.reloadAndGetList(grhigh), grhigh);
                    double revl= score(tsr.reloadAndGetList(grlow), grlow);
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
        int low= (int)telr.low+sublen;
        int high= (int)telr.high-sublen;
        Region eur= new Region(low, high);
        ArrayList<Region> rs= gb.getRegionsForChromosome(chr);
        for(int j=0; j<rs.size() ; j++){
            GenomicRegion r= (GenomicRegion)rs.get(j);
            if(eur.contains(r)){
                r.setValue(0);
            }else{
                r.setValue(1);
            }
           
            gbout.add((GenomicRegion)r);
        }
    }
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
            GenomicRegionsBuilder gbac= new GenomicRegionsBuilder();

            Iterator<GenomicRegion> ita= gba.iterator();
            Iterator<GenomicRegion> itb= gbb.iterator();
            while(ita.hasNext()){
                GenomicRegion ra= ita.next();
                GenomicRegion rb= itb.next();
                double va= ra.getValue();
                double vb= rb.getValue();
                if(va>z){
                    rb.setID(Double.toString(va)+"_"+Double.toString(vb));
                    rb.setValue(0);
                    if(vb<ac){
                        gbde.add(rb);
                    }else{
                        gbac.add(rb);
                    }
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
            tsf.keepLeft= 500;
            tsr.keepLeft= 500;
           
            ArrayList<Region> regs= gb.getRegionsForChromosome(chr);
            for( int j=0; j<regs.size() ; j++ ){
                GenomicRegion gr= (GenomicRegion)regs.get(j);
                GenomicRegion grlow, grhigh;
                int low, high;
                low= (int)gr.low;
                high= (int)gr.high;
                grlow= new GenomicRegion(gr.chr, low-side, low+side);
                grhigh= new GenomicRegion(gr.chr, high-side, high+side);
               
                double forl= score(tsf.reloadAndGetList(grlow), grlow);
                double forh= score(tsf.reloadAndGetList(grhigh), grhigh);
                double revl= score(tsr.reloadAndGetList(grlow), grlow);
                double revh= score(tsr.reloadAndGetList(grhigh), grhigh);
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

34
35
36
37
38
39
40
41
42
43
44
45
    while((l=br.readLine())!=null){
        String[] ss= l.split("\t");
        String id= ss[0];
        String loc= ss[6];
        if(id.contains("CUFF")){
            GenomicRegion gr= GenomicRegion.parseGenomicRegoinUCSC(loc);
            gr.setID(id);
            gbout.add(gr);
        }
    }
    br.close();
   
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
   
    ArrayList<GenomicRegion> al= gb.getAllGenomicRegions();
    RandomObject<GenomicRegion> rand= new RandomObject( al, null);
    HashSet<GenomicRegion> set= new HashSet<>();
    while(set.size()<n){
        GenomicRegion reg= rand.randomObject();
        if(!set.contains(reg)){
            set.add(reg);
        }
    }
   
    for( int i=0; i<al.size(); i++ ){
        GenomicRegion gr= al.get(i);
        if(set.contains(gr)){
            gbout.add(gr);
        }else{
            xxx.add(gr);
        }
View Full Code Here

Examples of fork.lib.bio.anno.genomic.region.GenomicRegion

75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
            double loc= Double.parseDouble(r[3]);
            String chr= r[1].substring(0,r[1].length()-12).toLowerCase();
            if(ind!=0){
                if(gba.containsChromosome(pc)){
                    if(pc.equals(chr)){
                        GenomicRegion gr= new GenomicRegion(pc, (int)pp, (int)loc-1);
                        ptr.put(pid, gr);
                    }else{
                        double cl= gba.getRegionsForChromosome(pc).get(0).high;
                        GenomicRegion gr= new GenomicRegion(pc, (int)pp, (int)cl-1);
                        ptr.put(pid, gr);
                    }
                }
            }
            if(i==arr.length-1){
                if(gba.containsChromosome(chr)){
                    double cl= gba.getRegionsForChromosome(chr).get(0).high;
                    GenomicRegion gr= new GenomicRegion(chr, (int)loc, (int)cl-1);
                    ptr.put(id, gr);
                }
            }
            pp= loc;
            pc= chr;
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.