Description Usage Arguments Details Note References Examples
This function provide filterByExpr with two customized options. See 'Arguments'. Different strategies can be useful either if you are trying to compare different approaches of normalization or if you want to analyze a particular biotype where a different variation of expression is expected under certain conditions. A BioInsight data frame will return with your new count matrix where you can proceed with your Differential Expression Analysis.
1 2 3 4 5 6 7 8 | describeRNA(
counts,
biotypes,
groups,
report = FALSE,
verbose = FALSE,
filter = 1
)
|
counts |
data.frame where you have gene counts |
biotypes |
data.frame where you have a gene_biotype column |
groups |
factor with groups and samples (see Examples) |
report |
if TRUE will generate a .pdf file in /tmp/ folder ('wordcloud' and 'RColorBrewer' are necessary) |
verbose |
if TRUE will print your result in your console ('knitr' is necessary) |
filter |
Your threshold (see edgeR::filterByExpr) |
For filter you may want to use "1" if you want the default option from filterByExpr. If 2 "Slightly" above the default will be applied. If "3" A more restrictive option will be applied.
You can use trace(describeRNA, edit=T) to set different values as threshold for "filter" option.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139-140. doi:10.1093/bioinformatics/btp616
1 2 3 4 5 6 7 8 9 10 11 | counts = system.file("extdata", "count_matrix.tsv", package="BioInsight")
counts = read.table(counts, row.names=1, header=TRUE)
biotypes = system.file("extdata", "Rattus_Norvegicus_biomart.tsv", package="BioInsight")
biotypes = read.table(biotypes, row.names=1, header=TRUE)
groups = rep(as.factor(c("1","2")), each=5)
describeRNA(counts=counts,
biotypes=biotypes,
groups=groups,
filter=2)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.