R/rqcCycleGCCalc.R

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

#' @describeIn rqcCycleGCPlot calculates necessary statistics
#' @export
rqcCycleGCCalc <- function(rqcResultSet)
{
    f <- function(x) {
        gc <- sum(x$count[c(2,3)]) / sum(x$count) * 100
        data.frame(gc)
    }
    
    df <- perCycleBasecall(rqcResultSet)
    ddply(df, c("filename", "cycle"), f)
}