plotFragmentSizeDistributionEstimate {ramwas} | R Documentation |
RaMWAS functions for estimation and plotting of the fragment size distribution.
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")
.
# 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)