Package fork.lib.bio.anno.genomic.region

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


   
   
    GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
    Iterator<GenomicRegion> it = gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        Integer len= (int) gr.getRange()+1;
        gr.setValue(len);
    }
   
    new BedExporter(gb).writeToFile(new File(dir+"/anno/bed/"+org+"/prop/"+org+"_length.bed"));
   
}
View Full Code Here


protected void init() throws Exception{
    gbout= new GenomicRegionsBuilder();
    Iterator<GenomicRegion> it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        int low = (int)gr.low;
        int high= (int)gr.high;
        int ind= low;
        while( (ind+len)<=high ){
            gbout.add(new GenomicRegion(gr.chr, ind, ind+len-1));
            ind+= len;
        }
    }
}
View Full Code Here

    BedGraphExporter.writeTitle(bw, tit);
   
    FastaSequenceExtractor se= new FastaSequenceExtractor(gen, gbout);
    SequenceExtractorEntry en;
    while( (en= se.nextEntry())!=null ){
        GenomicRegion gr= en.getGenomicRegion();
        NucleotideSequence seq= NucleotideSequenceParser.parseSequence(en.getSequence());
        double at= seq.atContent();
        bw.write( gr.toBedFormat()+"\t"+ ((double)Math.round(at*100)/100)+"\n" );
    }
   
    bw.close();
}
View Full Code Here

protected void init() throws Exception{
    GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
    FastaSequenceExtractor se= new FastaSequenceExtractor(sf, gb);
    SequenceExtractorEntry en;
    while((en=se.nextEntry())!=null){
        GenomicRegion reg= en.getGenomicRegion();
        String seq= en.getSequence();
        String id= reg.getID();
        idseq.put(id , seq);
        tab.appendEmptyRow(id);
    }
}
View Full Code Here

        String chr= chrs[i];
        Region r= gb.getRegionsForChromosome(chr).get(0);
        double low = r.low;
        while(low+len< r.high){
            Region win= new Region((double)low, (double)low+len-1);
            gbout.add(new GenomicRegion(chr, win));
            low+=len;
        }
    }
}
View Full Code Here

   
   
protected void init() throws Exception {
    Iterator<GenomicRegion> it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        double v= gr.getValue();
        dis.add(v);
    }
}
View Full Code Here

    for( int i=0; i<binn; i++ ){
        gbouts.add(new GenomicRegionsBuilder());
    }
    Iterator<GenomicRegion> it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        double v= gr.getValue() ;
        for( int i=0; i<bs.size()-1; i++ ){
            if( v>=bs.get(i) && v<bs.get(i+1) ){
                gbouts.get(i).add(gr);
                break;
            }
View Full Code Here

TOP

Related Classes of fork.lib.bio.anno.genomic.region.GenomicRegion

Copyright © 2018 www.massapicom. 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.