## Get some sequence files
sq.files = list.files(system.file('extdata', package='qrqc'),
pattern="test.*fastq", full.names=TRUE)
names(sq.files) <- gsub("(.*)\\.fastq", "\1", basename(sq.files))
sq <- lapply(sq.files, readSeqFile)
## Take the FASTQSummary objects, extract quality data from each of
## the, and combine.
d <- list2df(sq, getQual)
## Look at difference in average quality
aggregate(d$mean, list(sample=d$sample), mean)
## Look at difference in variance - this is where we really see a
## change.
aggregate(d$mean, list(sample=d$sample), var)
Run the code above in your browser using DataLab