humanViralDegree: Compare viral and human protein degree distribution, filter...

In vitkl/MItools: Protein interaction data tools: retrieving data from IntAct and most other databases using PSICQUIC service

Description

subset molecular interaction data (cleaned MITAB format in a data.table object) with a list of interactors

Usage

humanViralDegree(data = NULL, directory = "./data_files/",
  Interaction_detection_methods = NULL, Identification_method = NULL,
  PMIDs = NULL, inverse_filter = F, data_name = "")

Arguments

data	a list contatining PPI data as generated by loadHumanViralPPI, if NULL, data will be downloaded internally
directory	where to find / keep PPI data (arg for fullInteractome, interSpeciesInteractome)) and MI ontology (arg for subsetMITABbyMethod)
Interaction_detection_methods	arg for subsetMITABbyMethod
Identification_method	arg for subsetMITABbyMethod
PMIDs	arg for subsetMITABbyPMIDs
inverse_filter	logical, inverse filtering criteria
data_name	to be displayed on plot

Value

data.table to be ggplotted to compare viral and human protein degree distributions

Author(s)

Vitalii Kleshchevnikov

Examples

{
# Load data. This can be done internally by humanViralDegree(),
# however, it takes time (minutes) so it's recommended to load data separately
# if it will be used multiple times
HumanViralPPI = loadHumanViralPPI(directory = "./data_files/")

# prepare all data for the degree distribution analysis
degree_distributions = humanViralDegree(data = HumanViralPPI,
  directory = "./data_files/", data_name = "Full IntAct")

# prepare two-hybrid data for the degree distribution analysis
degree_distributions = humanViralDegree(data = HumanViralPPI,
  directory = "./data_files/", Interaction_detection_methods = "MI:0018",
  data_name = "two-hybrid data")

# prepare AP-MS data for the degree distribution analysis
degree_distributions = humanViralDegree(data = HumanViralPPI,
  directory = "./data_files/", Interaction_detection_methods = "MI:0004",
  Identification_method = "MI:0433", PMIDs = NULL, inverse_filter = F,
  data_name = "AP-MS data")

# to see if viral-interacting human proteins are special in the human network,
# look at the Vidal published and unpublished datasets (only human-human network
# is modified (data for the top 2 plots))
degree_distributions = humanViralDegree(data = HumanViralPPI,
  directory = "./data_files/", PMIDs = c("25416956", "unassigned1304"),
  data_name = "Vidal published and unpublished")
# to see if viral-interacting human proteins are special in the human network,
# do the same for Matthias Mann 2015 paper dataset
degree_distributions = humanViralDegree(data = HumanViralPPI,
  directory = "./data_files/", PMIDs = "26496610",
  data_name = "Matthias Mann 2015 paper")
  }

vitkl/MItools documentation built on May 29, 2019, 2:55 p.m.