Description Usage Arguments Details Value Author(s) See Also Examples
In preparation for a differential expression analysis, you will have to choose how to scale the raw counts provided by the recount project. Note that the raw counts are the sum of the base level coverage so you have to take into account the read length or simply the total coverage for the given sample (default option). You might want to do some further scaling to take into account the gene or exon lengths. If you prefer to calculate read counts without scaling check the function read_counts.
1 2 3 4 5 6 7 8 | scale_counts(
rse,
by = "auc",
targetSize = 4e+07,
L = 100,
factor_only = FALSE,
round = TRUE
)
|
rse |
A RangedSummarizedExperiment-class object as downloaded with download_study. |
by |
Either |
targetSize |
The target library size in number of single end reads. |
L |
The target read length. Only used when |
factor_only |
Whether to only return the numeric scaling factor or
to return a RangedSummarizedExperiment-class
object with the counts scaled. If set to |
round |
Whether to round the counts to integers or not. |
Rail-RNA http://rail.bio uses soft clipping when aligning
which is why we recommed using by = 'auc'
.
If the reads are from a paired-end library, then the avg_read_length
is the average fragment length. This is taken into account when using
by = 'mapped_reads'
.
If factor_only = TRUE
it returns a numeric vector with the
scaling factor for each sample. If factor_only = FALSE
it returns a
RangedSummarizedExperiment-class object with
the counts already scaled.
Leonardo Collado-Torres
download_study, read_counts
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | ## Load an example rse_gene object
rse_gene <- rse_gene_SRP009615
## Scale counts
rse <- scale_counts(rse_gene)
## Find the project used as an example
project_info <- abstract_search("GSE32465")
## See some summary information for this project
project_info
## Use the following code to re-download this file
## Not run:
## Download
download_study(project_info$project)
## Load file
load(file.path(project_info$project, "rse_gene.Rdata"))
identical(rse_gene, rse_gene_SRP009615)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.