Examples of LandscapeScoring

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

43
44
45
46
47
48
49
50
51
52
53
54
    }
   
   
   
public void writeToFile(File out)throws Exception {
    LandscapeScoring ls= new LandscapeScoring(lb, gb){
            protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                return new double[]{Landscape2DBuffer.getArea(bufs, gr) / (gr.getRange()+1)*1000 };
            }
        };
    ls.writeToBedFile(out);
}
View Full Code Here

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

27
28
29
30
31
32
33
34
35
36
37
38
39
40
    File lbf= new File("/home/forksapien/mystudy/phd/files/data/chip-seq/norm/norm_742.wig");
    //File lbf= new File("/home/forksapien/mystudy/phd/files/data/chip-seq/norm/norm_743.wig");
    //File lbf= new File("/home/forksapien/mystudy/phd/files/data/chip-seq/norm/norm_744.wig");
   
   
    LandscapeScoring ls= new LandscapeScoring(
            new BedGraphReader(lbf).getLandscapeBuilder(),
            new BedReader(reg).getGenomicRegionsBuilder()
            );
   
    File out = new File(d+"/score/score_"+ FileName.getBaseName(lbf)+"_"+FileName.getBaseName(reg)+".bed");
    ls.writeToBedFile(out);
   
}
View Full Code Here

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

41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
    GenomicRegionsBuilder gb= new BedReader(rf).getGenomicRegionsBuilder();
   
    LandscapeBuilder lb= new BedGraphReader(new File(datd+"/norm_wig_chip_htz1-anti660_772-wt_sd_20131104-r1_pe1_sacCer1.wig")).getLandscapeBuilder();
   
   
    LandscapeScoring ls= new LandscapeScoring(lb, gb){
            protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                return new double[]{Landscape2DBuffer.getArea(bufs, gr) / ((gr.getRange()+1)/1000) / 100 };
            }
        };
   
   
    ls.writeToBedFile(new File(zd+"/score_"+FileName.getBaseName(rf)+".bed") );

}
View Full Code Here

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

40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
public void writeToDir(File od)throws Exception {
    ArrayList<File> fs= GenerateGenic.getFiles();
    for( int i=0; i<fs.size(); i++ ){
        File f= fs.get(i);
        GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
        LandscapeScoring ls= new LandscapeScoring(lb, gb){
        protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
            return new double[]{Landscape2DBuffer.getArea(bufs, gr) / (gr.getRange()+1) *1000 };
        }

        };
        File out= new File(od+"/"+tag+"_"+FileName.getBaseName(f)+".bed");
        ls.writeToBedFile(out);
    }
}
View Full Code Here

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

39
40
41
42
43
44
45
46
47
48
49
50
    }
   
   
   
public void writeToFile(File out)throws Exception {
    LandscapeScoring ls= new LandscapeScoring(lb, gb){
            protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                return new double[]{Landscape2DBuffer.getArea(bufs, gr)  };
            }
        };
    ls.writeToBedFile(out);
}
View Full Code Here

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

56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
    for( int i=0; i<fs.length ; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("bed")){
            gfs.add(f);
            GenomicRegionsBuilder gb= new BedReader(f).getGenomicRegionsBuilder();
            LandscapeScoring ls= new LandscapeScoring(lb, gb){
                protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                    return new double[]{Landscape2DBuffer.getArea(bufs, gr), gr.getRange()+1 };
                }
            };
            Distribution dis= LandscapeScoringResultEntry.toDistribution(ls.getResultEntries(), 0);
            lens.add(LandscapeScoringResultEntry.toDistribution(ls.getResultEntries(), 1).sum());
            double s= dis.sum();
            vs.add(s); sum+=s;
            //LandscapeScoringResultEntry.writeToFile(ls.getResultEntries(), new File(out.getParentFile()+"/file_"+f.getName()) );
        }
    }
View Full Code Here

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

46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
    for( int i=0; i<fs.length; i++ ){
        File f= fs[i];
        if(FileName.getExt(f).equals("wig")){
            System.out.println(f);
            LandscapeBuilder lb= new BedGraphReader(f).getLandscapeBuilder();
            LandscapeScoring ls= new LandscapeScoring(lb, gb){
                protected double[] score(ArrayList<Region> bufs, GenomicRegion gr) throws RegionException{
                    return new double[]{Landscape2DBuffer.getArea(bufs, gr)  / (gr.getRange()+1)
                    };
                }

            };
            ls.writeToBedFile(new File(od+"/score-lambda_"+FileName.getBaseName(f)+".bed"));
        }
    }
   
   
   
View Full Code Here

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

35
36
37
38
39
40
41
42
public void writeToFile(File out)throws Exception {
    File dir= Dirs.getFile("dir");
    File anno= new File(dir+"/anno/bed/sacCer1-xu/sacCer1-xu_gene_five_0.150.bed");
    GenomicRegionsBuilder gb= new BedReader(anno).getGenomicRegionsBuilder();
    LandscapeScoring ls= new LandscapeScoring(lb ,gb);
    ls.writeToBedFile(out);
}
View Full Code Here

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

63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
            String samp= samps[k];
            String tag= samp+"_"+rep;
            File datf= new File(datd+"/norm_wig_chip_"+samp+"_"+rep+"_sacCer1.wig");
            //File datf= new File(datd+"/norm_"+samp+".wig");

            LandscapeScoring ls= new LandscapeScoring(
                new BedGraphReader(datf).getLandscapeBuilder(),
                new BedReader(reg).getGenomicRegionsBuilder()
                );

            File out = new File(outd+"/"+tag+"/set/score-len_"+ FileName.getBaseName(datf)+"_"+FileName.getBaseName(reg)+".bed");
            System.out.println(out);
            ls.writeToBedFile(out);
        }

    }
   
   
View Full Code Here

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

47
48
49
50
51
52
53
54
55
56
57
    od.mkdirs();
    File[] fs= files(org);
    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) /100 / ((gr.getRange()+1)/1000) };
            }
        };
        File out= new File(od+"/"+samp+"_"+FileName.getBaseName(f)+".bed");
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.