interactome: Create an interactome from BioGRID data sets

In HTSanalyzeR: Gene set over-representation, enrichment and network analyses for high-throughput screens

Description

This is a generic function.

When implemented as the S4 method of class NWA, this function creates an interactome before conducting network analysis.

To use this function for objects of class NWA:

interactome(object, interactionMatrix, species, link, reportDir = "HTSanalyzerReport", genetic=FALSE, verbose=TRUE)

Usage

interactome(object, ...)

Arguments

object	an object. When this function is implemented as the S4 method of class NWA, this argument is an object of class 'NWA'
...	other arguments (see below for the arguments supported by the method of class NWA)

interactionMatrix:

an interaction matrix including columns 'InteractionType', 'InteractorA' and 'InteractorB'. If this matrix is available, the interactome can be directly built based on it.

species:

a single character value specifying the species for which the data should be read. The current version supports one of the following species: "Dm" ("Drosophila_melanogaster"), "Hs" ("Homo_sapiens"), "Rn" ("Rattus_norvegicus"), "Mm" ("Mus_musculus"), "Ce" ("Caenorhabditis_elegans").

link:

the link (url) where the data should be downloaded (in tab2 format). The default link is version 3.1.71 (valid on Dec. 5 2010).

reportDir:

a single character value specifying the directory to store reports. The BioGRID data set will be downloaded and stored in a subdirectory called 'Data' in 'reportDir'.

genetic:

a single logical value. If TRUE, genetic interactions will be kept; otherwise, they will be removed from the data set.

verbose:

a single logical value indicating to display detailed messages (when verbose=TRUE) or not (when verbose=FALSE)

Details

This function provides two options to create an interactome for network analysis. The user can either input an interaction matrix including columns 'InteractionType', 'InteractionA' and 'InteractionB', or set 'species', 'link' and 'genetic' to download data set from BioGRID and extract corresponding interactions to build the interactome.

Another way to set up the interactome is to input a graphNEL object when the NWA object is created (i.e. nwa=new("NWA", pvalues, phenotypes, interactome)).

Value

In the end, this function will return an updated object with slot 'interactome' as an object of class graphNEL.

Author(s)

Xin Wang [e-mail]

See Also

biogridDataDownload, NWA

Examples

## Not run: 
##example 1, input interactome when initializing an 'NWA' object
##load p-values and phenotypes
data("KcViab_PVals","KcViab_Data4Enrich")
##load BioGRID interactome for Drosophila Melanogaster
data("Biogrid_DM_Interactome")
##create a NWA (NetWork Analysis) object
nwa <- new("NWA",pvalues=KcViab_PVals, phenotypes=KcViab_Data4Enrich, 
interactome=Biogrid_DM_Interactome)
##print nwa
nwa

library(BioNet)
##example 2, build an interactome from Biogrid data base
##create a NWA (NetWork Analysis) object
nwa <- new("NWA", pvalues=KcViab_PVals, phenotypes=KcViab_Data4Enrich)
##print nwa
nwa
##download data from BioGRID and build the interactome
nwa <- interactome(nwa, species="Dm", reportDir="NWATest")
##print nwa again
nwa

## End(Not run)

HTSanalyzeR documentation built on Oct. 31, 2019, 7:10 a.m.