Description Usage Arguments Value
This function can be used to prepare artificial replicate maps (with no differential interactions) from any 2 Hi-C contact maps. To do so one must first prepare HiCglm object from any pair of contact maps (can be technical or biological replicates) and then feed it into this function. Resulting maps can be used to benchmark DI detection methods (for example by adding some artificial DI and then testing it). This approach has an advantage of producing artificial maps resembling real Hi-C data, because it keeps decay, compartment and TADs unaffected (at least to some extent).
1 2 | ## S3 method for class 'HiCglm'
simulate_map(hic.glm, dpools = "all", remove.na.contacts = TRUE)
|
hic.glm |
HiCglm object |
dpools |
character or vector, if character "all" then use all pooled diagonals in background model, otherwise use specified subset |
remove.na.contacts |
logical if true remove NAs |
list of list of data frames with expected interactions, first list enumerates different maps (usually chromosomes), second list contains 2 elements refering to Y (simulate Y condition from X) and X (simulate X condition from Y) contact maps in sparse format.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.