postMerge {RJMCMCNucleosomes} | R Documentation |
rjmcmc
functionA internal helper function which merges closely positioned nucleosomes to rectify the over splitting and provide a more conservative approach. Beware that each chromosome must be treated separatly.
The function uses the Bioconductor package
consensusSeeker
to
group closely positioned nucleosomes.
postMerge(reads, resultRJMCMC, extendingSize, chrLength, minReads = 5, seqName = NULL)
reads |
a |
resultRJMCMC |
an object of class 'rjmcmcNucleosomes' or 'rjmcmcNucleosomesMerge' containing informations about nucleosomes. |
extendingSize |
a positive |
chrLength |
a positive |
minReads |
a positive |
a array
of numeric
, the updated values of the
nucleosome positions. When no nucleosome is present, NULL
is
returned.
Pascal Belleau, Astrid Deschenes
## Loading dataset data(RJMCMC_result) data(reads_demo_02) ## Results before post-treatment RJMCMC_result$mu ## Post-treatment function which merged closely positioned nucleosomes postResult <- RJMCMCNucleosomes:::postMerge(reads = reads_demo_02, resultRJMCMC = RJMCMC_result, extendingSize = 80, chrLength = 73500) ## Results after post-treatment postResult