That's not a bad idea, and a quick way to do this would be to add pseudo-counts to each fragment. That would make the difference of 1 or 2 reads not nearly as significant.
minor changes
added some basic error checking
Hi, I am using your R script for identification of 4C interactions. It works well...
Hi Joshua, How can I extract the FDR from your pipeline? I want to create a table...
Fixed a problem when there are no reads mapping...
Modified the messages printed when the user inv...
minor typo correction
fixed some bugs with saving output in BED format
added to help with the documentation
Made things easier for "trans" analyses.
improved the functions that generate Bed and Be...
improved documentation for BED and BedGraph files
fixed a problem where files intended for the UC...
We can now read BAM files, and chromosome names...
improved error messages
Fixed a bug in saveSigAsBed - before it could o...
remove .Rhistory
Initial commit
Home