qSVA | R Documentation |
A wrapper function used to perform qSVA in one step.
qSVA(
rse_tx,
type = c("cell_component", "standard", "top1500"),
sig_transcripts = select_transcripts(type),
mod,
assayname
)
rse_tx |
A RangedSummarizedExperiment-class object containing the transcript data desired to be studied. |
type |
a character string specifying which model you would like to use when selecting a degradation matrix. |
sig_transcripts |
A list of transcripts that are associated with
degradation signal. Use |
mod |
Model Matrix with necessary variables the you would model for in differential expression |
assayname |
character string specifying the name of the assay desired in rse_tx |
matrix with k principal components for each sample
## First we need to define a statistical model. We'll use the example
## covComb_tx_deg data. Note that the model you'll use in your own data
## might look different from this model.
mod <- model.matrix(~ mitoRate + Region + rRNA_rate + totalAssignedGene + RIN,
data = colData(covComb_tx_deg)
)
## To ensure that the results are reproducible, you will need to set a
## random seed with the set.seed() function. Internally, we are using
## sva::num.sv() which needs a random seed to ensure reproducibility of the
## results.
set.seed(20230621)
qSVA(rse_tx = covComb_tx_deg, type = "cell_component", mod = mod, assayname = "tpm")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.