Description Usage Arguments Value Note Author(s) See Also Examples
RaMWAS functions for estimation and plotting of the fragment size distribution.
1 2 3 4 5 6 | estimateFragmentSizeDistribution(frdata, seqLength)
plotFragmentSizeDistributionEstimate(
frdata,
estimate,
col1 = "blue",
col2 = "red")
|
frdata |
Distribution of distances from the starts of isolated reads to the respective CpGs. |
seqLength |
The length of sequenced part of the fragments. |
estimate |
Fragment size distribution estimate. |
col1 |
Color of |
col2 |
Color of |
The function estimateFragmentSizeDistribution
returns the estimate of the fragment size distribution.
If the length of frdata
is equal to seqLength
,
the fragments are assumed to all be of length seqLength
.
Andrey A Shabalin andrey.shabalin@gmail.com
See vignettes: browseVignettes("ramwas")
.
1 2 3 4 5 6 7 8 9 10 | # Simulate data
x = 0:250
truemean = 1 - pnorm(x, mean = 150, sd = 50)
frdata = rpois(n = length(x), lambda = truemean*300)
# Estimate fragment size distribution
estimate = estimateFragmentSizeDistribution(frdata, seqLength = 50)
# Plot fragment size distribution estimate
plotFragmentSizeDistributionEstimate(frdata, estimate)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.