qsva: A function for computing quality surrogate variables (qSVs)

Description Usage Arguments Value Examples

View source: R/qsva.R

Description

This function computes quality surrogate variables (qSVs) from the library-size- and read-length-normalized degradation matrix for subsequent RNA quality correction

Usage

1
2
3
4
qsva(
  degradationMatrix,
  mod = matrix(1, ncol = 1, nrow = ncol(degradationMatrix))
)

Arguments

degradationMatrix

the normalized degradation matrix, region by sample

mod

(Optional) statistical model used in DE analysis

Value

the qSV adjustment variables

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
## Find files
bwPath <- system.file('extdata', 'bwtool', package = 'sva')

## Read the data
degCovAdj = read.degradation.matrix(
 covFiles = list.files(bwPath,full.names=TRUE),
 sampleNames = list.files(bwPath), readLength = 76, 
 totalMapped = rep(100e6,5),type="bwtool")

## Input data
 head(degCovAdj)

## Results
qsva(degCovAdj)

steveneschrich/msva documentation built on Dec. 23, 2021, 5:33 a.m.