ADAnalysis: Group of Functionally Associated Genes (GFAG) partial...

In jrybarczyk/ADAM: ADAM: Activity and Diversity Analysis Module

Group of Functionally Associated Genes (GFAG) partial analysis

Description

Analysis of functionally associated gene groups, based on gene diversity and activity, for different species according to existing annotation packages or user's personal annotations. ADAnalysis function allows to run a partial analysis, where is calculated just gene diversity and activity of each GFAG with no signicance by bootrstrap, Wilcoxon or Fisher.

Usage

ADAnalysis(ComparisonID, ExpressionData, MinGene, MaxGene, 
DBSpecies, AnalysisDomain, GeneIdentifier)

Arguments

ComparisonID	Sample comparisons identification. It must be a vector in which each element corresponds to 2 sample columns from the expression data. The data's sample columns in each element from the vector are comma separated. The first one refers to the control sample, while the second refers to the experiment. This argument must be informed by the user. There is no default value for it.
ExpressionData	Gene expression data (microarray or RNA-seq, for example). It must be a SummarizedExperiment object, a data frame or a path for a text file tab separated containing at least 3 columns. First column mandatory corresponds to the gene identification, according to GeneIdentifier argument. Second, third and the other columns correspond to the gene expression values realated to the genes in the first column and each of these columns correspond to a different sample (control versus experiment). This argument must be informed by the user. There is no default value for it.
MinGene	Minimum number of genes per GFAG. It must be a positive integer value different from zero and lower than MaxGene argument. Default is 3.
MaxGene	Maximum number of genes per GFAG. It must be a positive integer value different from zero and greater than MinGene argument. Default is 2000.
DBSpecies	A string corresponding to the name of an OrgDb species gene annotation package: org.Ag.eg.db (Anopheles gambiae), org.At.tair.db (Arabdopsis thaliana), org.Bt.eg.db (Bos taurus), org.Ce.eg.db (Caenorhabditis elegans), org.Cf.eg.db (Canis familiaris), org.Dm.eg.db (Drosophila melanogaster), org.Dr.eg.db (Danio rerio), org.EcK12.eg.db (Escherichia coli K12), org.EcSakai.eg.db (Escherichia coli Sakai), org.Gg.eg.db (Gallus gallus), org.Hs.eg.db (Homo sapiens), org.Mm.eg.db (Mus musculus), org.Mmu.eg.db (Macaca mulatta), org.Pf.plasmo.db (Plasmodium falciparum), org.Pt.eg.db (Pan troglodytes), org.Rn.eg.db (Rattus norvegicus), org.Sc.sgd.db (Saccharomyces cerevisiae), org.Ss.eg.db (Sus scrofa) and org.Xl.eg.db (Xenopus laevis). If there is no package, it's possible for the user to create a personal gene annotation file, tab separated, containing 3 columns: gene, term annotation code and description of the term annotation. So, istead of a string with an OrgDb name, inform a data frame or a path for the file. This argument must be informed by the user. There is no default value for it.
AnalysisDomain	Analysis domain to be considered for building GFAGs, according: gobp (Gene Ontology - Biological Processes), gocc (Gene Ontology - Celular Components), gomf (Gene Ontology - Molecular Functions), kegg (KEGG Pathways) or own (if there is no annotation package - the annotations were defined in a file by the user). This argument must be informed by the user. There is no default value for it.
GeneIdentifier	Gene nomenclature to be used: entrez or tairID for Arabdopsis thaliana, entrez or orfID for Saccharomyces cerevisiae, symbol or orfID for Plasmodium falciparum and symbol or entrez for the other species. If there is no annotation package, just put the gene nomenclature present in the user's personal annotations. This argument must be informed by the user. There is no default value for it.

Details

The genes present in the expression data are grouped by their respective functions according to the domains described by AnalysisDomain argument. The relationship between genes and functions are made based on the species annotation package. If there is no annotation package, a three column file (gene, function and function description) must be provided. For each GFAG, gene diversity and activity in each sample are calculated. As the package always compare two samples (control versus experiment), relative gene diversity and activity for each GFAG are calculated.

Value

Return a list with two elements. The first one refers to a data frame with the GFAGs and their respective genes. The second one is a a list where each position is a data frame presenting the result of the analysis, according to ComparisonID argument.

Author(s)

André Luiz Molan ([e-mail])

References

CASTRO, M. A., RYBARCZYK-FILHO, J. L., et al. Viacomplex: software for landscape analysis of gene expression networks in genomic context. Bioinformatics, Oxford Univ Press, v. 25, n. 11, p. 1468–1469, 2009.

Examples

##
#Partial Analysis with Aedes aegypti through ADAnalysis function
##
data(ExpressionAedes)
data(KeggPathwaysAedes)
ResultAnalysis <- ADAnalysis(ComparisonID = c("control1,experiment1"),
ExpressionData = ExpressionAedes, MinGene = 3L, MaxGene = 20L, 
DBSpecies = KeggPathwaysAedes, AnalysisDomain = "own",
GeneIdentifier = "geneStableID")
## Not run: 
head(ResultAnalysis[[1]]) #Relation between genes and functions
head(ResultAnalysis[[2]][1]) #Result comparison 1

## End(Not run)

jrybarczyk/ADAM documentation built on March 3, 2024, 1:08 a.m.