/*
* To change this license header, choose License Headers in Project Properties.
* To change this template file, choose Tools | Templates
* and open the template in the editor.
*/
package data.cerevisiae.test.ratio;
import fork.lib.base.file.FileName;
import fork.lib.bio.anno.genomic.BedGraphReader;
import fork.lib.bio.anno.genomic.LandscapeBuilder;
import fork.lib.bio.anno.genomic.LandscapeTransformer;
import fork.lib.bio.anno.genomic.PeakIdentifier;
import java.io.File;
/**
*
* @author man-mqbpjmg4
*/
public class Threshold {
public static void main(String[] args) throws Exception { //debug
File f= new File("r2.wig");
LandscapeBuilder lb= new BedGraphReader(f).getLandscapeBuilder();
LandscapeTransformer lt= new LandscapeTransformer(lb);
lt.subtract(10);
PeakIdentifier pi= new PeakIdentifier(lb, null);
pi.getPeaks().writeToFileBed(new File("peak_"+FileName.getBaseName(f)+".bed"));
}
}