nbGLM: Identify Differential Expression Modules Based on the...
In SeqMADE: Network Module-Based Model in the Differential Expression Analysis for RNA-Seq

Description Usage Arguments Details Value Author(s) See Also Examples

The algorithm identify differential expression modules using Generalized Linear Model (GLM) for differential expression analysis in RNA-Seq data, and in the model two indicator variables Group and Module are adopted to fit the GLM.

1	nbGLM(factors, N, networkModule, modulematrix, distribution = c("poisson", "NB")[1])

`factors`	Factors with three variables including Count, Group, Direction.
`N`	The total sample sizes.
`networkModule`	NetworkModule is the gene sets or modules in the biological network or metabolic pathway, with the 1th column as the module names and the 2th columnn as the gene symbol constituting the module.
`modulematrix`	Modulematrix is a matrix, in which the indicator variables 1 or 0 represent whether a gene belong to a given module or not.
`distribution`	a character string indicating the distribution of RNA-Seq count value, default is 'NB'.

The GLM method was determined by the distribution of RNA-Seq count value including Poisson and Negative Binomial distribution, and there are two indicator variables Group and Module, Module=1 when a gene belongs to the module and Module= 0 otherwise; Group=1 for case values and Group=0 for control values. Group * Module represents the interaction effects between Group and Module, and the significance of a module is decided by the interaction and adjusted p-values are calculated to correct for multiple testing.

The nominal pvalue and FDR for the significance of each gene set or module.

Mingli Lei, Jia Xu, Li-Ching Huang, Lily Wang, Jing Li

glm()

data(exprs)
data(networkModule)
case <- c("A1","A2","A3","A4","A5","A6","A7")
control <- c("B1","B2","B3","B4","B5","B6","B7")
factors <- Factor(exprs, case, control) 
modulematrix <- moduleMatrix(exprs,networkModule)
Result <- nbGLM(factors, 14, networkModule, modulematrix, distribution = "NB")