R/rqcCycleAverageQualityCalc.R

In Rqc: Quality Control Tool for High-Throughput Sequencing Data

#' @describeIn rqcCycleAverageQualityPlot calculates necessary statistics
#' @export
rqcCycleAverageQualityCalc <- function(rqcResultSet)
{
    f <- function(x) {
        quality <- mean(Rle(x$score, x$count))
        data.frame(quality)
    }
    
    df <- perCycleQuality(rqcResultSet)
    ddply(df, c("group", "filename", "cycle"), f)
}