Integrator: AbHAC Internal Enrichment Calculator

Description Usage Arguments Value Author(s)

Description

Internal Function that performs abhac analysis inside other functions

Usage

1
2
3
4
Integrator(ppi.database = NULL, list.categories = NULL, fac = NULL,
  id.conversion.set = NULL, fisher.fdr = "Permutation.FDR",
  fisher.fdr.cutoff = 0.05, num.permuted.ppi = 10,
  method.permuted.ppi = "AsPaper", bins.permuted.ppi = 4, num.cores = 6)

Arguments

ppi.database

2 column whole protein interaction network. Either loaded by data(ppi.database)(filtering is recommended based on types of interactions) or bu used.

list.categories

(Internal) is list of proteins for each enrichment category with accurate names as of enrichment.categories

fac

is all the proteins that exist in protein interaction network. If not using data(ppi.database), it is necessary to specify.

id.conversion.set

A dataframe for ID conversions provided as global variable id.conversion.set. Columns represent Entrez gene ID, Uniprot Accession, Gene Symbol, Ensembl gene ID and refseq protein ID (all human)

fisher.fdr

Can be either "Permutation.FDR", "Permutation.FWER" or any of the methods parsed into p.adjust. Type ?p.adjust for more details.

fisher.fdr.cutoff

Cutoff to be used for fisher's exact test false discovery rates.

num.permuted.ppi

If you have selected any of the two permutation based methods, the number of permuted networks to be used for multiple testing correction must be specified.

method.permuted.ppi

If you have selected any of the two permutation based methods, the method for binning proteins by their edge degree for creting permuted networks for multiple testing correction must be specified. It should be one of ("AsPaper", "equal", "ByDegree").

bins.permuted.ppi

If you have selected any of the two permutation based methods, specified the number of bins for proteins to be grouped into. If you have selected "AsPaper", you would better leave this as 4 (default). For the two other methods, we advise a number between 10-20.

num.cores

Note that a parallel for loop using foreach package calculates the p-values for all the different permuted networks. The number of processors to be used for this foreach has to be set by you using the "registerDoMC(cores=4)". However this parameters determined the number of processors to be used for calculation of pvalues by an mclapply feature. So if you are using registerDoMC(cores=4) and you want to limit the analysis to 12 processors, you must specify num.cores = 3.

Value

A dataframe with Benjamini-Hochberg FDR bases on fisher's one tail exact test for all enrichment categories for all proteins that interact with at least one of the input proteins

Author(s)

Mehran Karimzadeh mehran.karimzadehreghbati at mail dot mcgill dot ca


mehrankr/AbHAC documentation built on May 22, 2019, 6:49 p.m.