RNAseq_normalization: Normalize RNAseq data according to DESeq2 method.

View source: R/trscr_analysis.R

RNAseq_normalizationR Documentation

Normalize RNAseq data according to DESeq2 method.

Description

Perform normalization of transcriptomic data (RNAseq) using DESeq2 R package.

Usage

RNAseq_normalization(data_trscr, exp_grp, gene_list,
  filter_indiv = "no_filter", alpha = 0.05, contrast = c("tissue_status",
  "patho", "normal"), fitType = "parametric")

Arguments

data_trscr

A data matrix that contains transcriptome information (RNAseq counts from HTseq). Columns correspond to indivuals, row correspond to genes.

exp_grp

A exp_grp dataframe that contains metadatas on data_trscr individuals.

gene_list

A gene_list bedfile containing the genes to screen for differential expression.

filter_indiv

A vector of individual names to be screened for differential expression. Optionnal (set on "no_filter" by default).

alpha

A parameter to indicate the significance cutoff used by DESeq2::results funtion for optimizing the independent filtering (by default 0.05). If the adjusted p-value cutoff (FDR) will be a value other than 0.1, alpha should be set to that value.

contrast

A vector containing the constrast to be used to estimate the logarithmic fols change. By default: c("tissue_status","patho","normal")

fitType

A DESeq fitting paramater, by default set to "parametric"

Value

A list composed of a data_ntrscr matrix of normalized counts and a size_factor vector of normalization factors.


magrichard/dmprocr documentation built on July 21, 2023, 11:01 p.m.