Skip to content

Commit

Permalink
cleanup some comments from vignette
Browse files Browse the repository at this point in the history
  • Loading branch information
gpertea committed Sep 22, 2024
1 parent 63abfee commit 4294381
Showing 1 changed file with 1 addition and 4 deletions.
5 changes: 1 addition & 4 deletions vignettes/Intro_qsvaR.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -222,11 +222,8 @@ dim(qsvs_wrapper)
Next we can use a standard limma package approach to do differential expression on the data. The key here is that we add our qsvs to the `model.matrix()`. Here we input our `RangedSummarizedExperiment` object and our `model.matrix()` with qSVs. Note here that the `RangedSummarizedExperiment` object is the original object loaded with the full list of transcripts, not the the one we subsetted for qSVs. This is because while PCs can be generated from a subset of genes, differential expression is best done on the full dataset. The expected output is a `sigTx` object that shows the results of differential expression.

```{r "perform DE"}
### should be done by cbinding mod to pcs
## subset to an expression cutoff
#rse_tx <- rse_tx[rowData(rse_tx)$passExprsCut, ]
# create a model.matrix with demographic infor and qsvs
# create a model.matrix with demographic info and qsvs
mod_qSVA <- cbind(mod, qsvs)
# log tranform transcript expression
Expand Down

0 comments on commit 4294381

Please sign in to comment.