Compute log2(counts per mil reads) and library size for each sample
log2CPM(qcounts, lib.size = NULL)
quantile normalized counts
default is colsums(qcounts)
list containing log2(quantile counts per mil reads) and library sizes
example_data_dir <- system.file("example/data", package = "PathoStat")
pathoreport_file_suffix <- "-sam-report.tsv"
datlist <- readPathoscopeData(example_data_dir, pathoreport_file_suffix)
countdat <- datlist$countdata
lcpm <- log2CPM(countdat)
Questions? Problems? Suggestions? Tweet to @rdrrHQian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.