peakRefinement: Characterize Euclidean distance peaks using resampling...

In jonathonthill/MMAPPR2: Mutation Mapping Analysis Pipeline for Pooled RNA-Seq

Description

Follows the prePeak step and precedes generateCandidates.

Usage

peakRefinement(mmapprData)

Arguments

mmapprData

The MmapprData object to be analyzed.

Value

A MmapprData object with the peaks slot filled and populated.

Examples

if (requireNamespace('MMAPPR2data', quietly=TRUE)
        & all(Sys.which(c("samtools", "vep")) != "")) {
    mmappr_param <- MmapprParam(refFasta = MMAPPR2data::goldenFasta(),
                               wtFiles = MMAPPR2data::exampleWTbam(),
                               mutFiles = MMAPPR2data::exampleMutBam(),
                               species = "danio_rerio",
                               outputFolder = tempOutputFolder())
}
## Not run: 
md <- new('MmapprData', param = mmappr_param)
postCalcDistMD <- calculateDistance(md)
postLoessMD <- loessFit(postCalcDistMD)
postPrePeakMD <- prePeak(postLoessMD)

postPeakRefMD <- peakRefinement(postPrePeakMD)

## End(Not run)

jonathonthill/MMAPPR2 documentation built on Oct. 2, 2019, 8:16 a.m.