R/rqcCycleAverageQualityPcaCalc.R

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

#' @describeIn rqcCycleAverageQualityPcaPlot calculates necessary statistics
#' @export
rqcCycleAverageQualityPcaCalc <- function(rqcResultSet)
{
  df <- rqcCycleAverageQualityCalc(rqcResultSet)
  df.wide <- dcast(df, cycle ~ filename, value.var = "quality")
  fit <- prcomp(~ ., data=df.wide[,-1, drop=FALSE])
  return(fit)
}