Examples of BedReader

Examples of fork.lib.bio.anno.genomic.BedReader

35
36
37
38
39
40
41
42
43
44
45
   
    File[] fs = d.listFiles();
    for( int i=0; i<fs.length; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("bed")){
            GenomicRegionsBuilder peaks= new BedReader(f).getGenomicRegionsBuilder();
            SplicedBedReader sr= new SplicedBedReader(anno);

            File od= new File(d+"/locate/"+ FileName.getBaseName(anno)+"/"+FileName.getBaseName(f));
            SplicedRegionsLocator ll= new SplicedRegionsLocator(peaks, sr);
            ll.writeToDir(od, "locate");
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

44
45
46
47
48
49
50
51
52
53
54
    //File anno= new File(dir+"/anno/refseq_hg19_exon.bed");
    File anno= new File(dir+"/anno/ensGene_hg19_exon.bed");
   
    File out= new File(d+"/locate_"+FileName.getBaseName(anno)+".bed");
   
    GenomicRegionsBuilder peaks= new BedReader(f).getGenomicRegionsBuilder();
    SplicedBedReader sr= new SplicedBedReader(anno);
   
    SplicedRegionsLocator ll= new SplicedRegionsLocator(peaks, sr);
    ll.writeToFile(out);
   
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

52
53
54
55
56
57
58
59
60
61
62
    File dir= Dirs.getFile("dir");
    File od= new File(dir+"/data/human_magnus/h2az/norm/peak/novel/scoreChip");
    File cd= new File(dir+"\\data\\human_magnus\\h2az\\norm");
    File rf= new File(dir+"/data\\human_magnus\\h2az\\norm\\peak\\novel/merge.bed");
   
    GenomicRegionsBuilder gb= new BedReader(rf).getGenomicRegionsBuilder();
   
    File[] fs= cd.listFiles();
    for( int i=0; i<fs.length; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("wig")){
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

55
56
57
58
59
60
61
62
63
64
65
    //File d= new File(dir+"/data\\human_magnus\\h2az\\norm\\peak/merge");
    File d= new File(dir+"/data\\human_magnus\\h2az/macs/merge");
    File od= new File(d+"/scoreNorm");
    File rf= new File(d+"/merge.bed");
   
    GenomicRegionsBuilder gb= new BedReader(rf).getGenomicRegionsBuilder();
   
    File[] fs= cd.listFiles();
    for( int i=0; i<fs.length; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("wig")){
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

29
30
31
32
33
34
35
36
37
38
39
   
public void addMergeFile(File f)throws Exception {
    String cntag= FileName.getBaseName(f);
    HashMap<String, String> hma= new HashMap<>();
    HashMap<String, String> hmb= new HashMap<>();
    GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
    Iterator<GenomicRegion> it= gb.iterator();
    while(it.hasNext()){
        GenomicRegion gr= it.next();
        String u= gr.chr+":"+(int)gr.low+"-"+(int)gr.high;
        hma.put( u, gr.getID());
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

33
34
35
36
37
38
39
40
41
42
43
    File[] fs= d.listFiles();
    for( int i=0; i<fs.length; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("bed") && !f.getName().equals(out.getName()) ){
            System.out.println(f);
            GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
            mer.mergeNoSort(gb);
        }
        mer.sort();
        new BedExporter(mer.getMergedGenomicRegionsBuilder()).writeToFile(out);
    }
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

29
30
31
32
33
34
35
36
37
38
39
    double sthr= 0.6;
    double zthr= 80;
    File dir= Dirs.getFile("dir");
   
   
    GenomicRegionsBuilder gbl= new BedReader(new File(dir+"/anno/bed/len/len_150_off_0.bed")).getGenomicRegionsBuilder();
    ArrayList<LandscapeScoringResultEntry> rsz= new LandscapeScoring(
            new BedGraphReader(new File(dir+"/data/chip-seq/norm/norm_wig_chip_742-wt_r1_sacCer1.wig")).getLandscapeBuilder(),
            gbl
            ).getResultEntries();
   
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

45
46
47
48
49
50
51
        String fn= FileName.getBaseName(f);
        if(fn.indexOf("chr")==0){
            chrGet.put(fn, new GetSequenceFasta(f));;
        }
    }
    idReg= new BedReader(bedf).getIDToRegionMap();
}
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

44
45
46
47
48
49
50
51
52
53
54
    od.mkdirs();
    File[] fs= ScoreChipGene.files("sacCer1-xu");
    for(int i=0; i<fs.length ; i++){
        File f=fs[i];
        System.out.println(f);
        ls= new LandscapeScoring(lb, new BedReader(f).getGenomicRegionsBuilder()){
            protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                return new double[]{Landscape2DBuffer.getArea(bufs, gr)};
            }
        };
        File out= new File(od+"/"+FileName.getBaseName(f)+".bed");
View Full Code Here

Examples of fork.lib.bio.anno.genomic.BedReader

31
32
33
34
35
36
37
38
39
40
41
    File ld= new File(d+"/lineup");
    //File regf= new File(ld+"/chr3_cn3300_cz3.bed");
    File regf= new File(ld+"/chr3_cn3300_cnz3.bed");
   
   
    GenomicRegionsBuilder gb= new BedReader(regf).getGenomicRegionsBuilder();
    GenomicRegionsModifier gm= new GenomicRegionsModifier(gb);
    gb= gm.getThreePrimeEnds(300, 100);
   
    File od= new File(ld+"/"+FileName.getBaseName(regf));
    ArrayList<File> fs= new ArrayList<>();
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.