RandomWalkRestart: Random Walk Restart (RWR)

Description Usage Arguments Value Examples

View source: R/PhenoGeneRankerFunctions.R

Description

This method runs the random walk with restart on the provided walk matrix. It returns a data frame including ranked genes and phenotypes, and RWR scores of the genes and phenotypes. If generatePvalue is TRUE then it generates p-values along with the ranks.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
RandomWalkRestart(
  walkMatrix,
  geneSeeds,
  phenoSeeds,
  generatePValue = TRUE,
  numCores = 1,
  r = 0.7,
  eta = 0.5,
  tau = NULL,
  phi = NULL,
  S = 1000
)

Arguments

walkMatrix

This is the walk matrix generated by the function CreateWalkMatrix.

geneSeeds

This is a vector for storing the ids of the genes that RWR starts its walk. The final ranks show the proximity of the genes/phenotypes to the seed genes.

phenoSeeds

This is a vector for storing the ids of the phenotypes that RWR starts its walk. The final ranks show the proximity of the genes/phenotypes to the seed phenotypes.

generatePValue

If this is TRUE, it will generate the probability values for each of the gene/phenotype rankings. If it is FALSE then the function will only return the ranks of genes/phenotype.

numCores

This is the number of cores used for parallel processing.

r

This parameter controls the global restart probability of RWR, and it has a default value of 0.7.

eta

This parameter controls restarting of RWR either to a gene seed or phenotype seeds, higher eta means utilizing gene seeds more than phenotype seeds, and it has a default value of 0.5.

tau

This is a vector that stores weights for each of the 'gene' layer in the complex gene and phenotype network. Each value of the vector corresponds to the order of the network files in the input file of CreateWalkMatrix function. The weights must sum up to the same number of gene layers. Default value gives equal weight to gene layers.

phi

This is a vector that stores weights for each of the 'phenotype' layer in the complex gene and phenotype network. Each value of the vector corresponds to the order of the network files in the input file of CreateWalkMatrix function. The weights must sum up to the same number of phenotype layers. Default value gives equal weight to phenotype layers.

S

This is the number of random samples to be used for p-value calculation It is highly recommended to use S=1000.

Value

If the parameter generatePValue is TRUE, then this function returns a data frame of ranked genes/phenotypes with p-values with three columns; Gene/Phenotype ID, score, p-value. If generatePValue is FALSE, then it returns a data frame of ranked genes/phenotypes with two columns; Gene/Phenotype ID, score.

Examples

1
2
3
4
5
wm <- CreateWalkMatrix('input_file.txt')
ranksWithoutPVal <- RandomWalkRestart(wm, c('g1', 'g2'), c('p1'), FALSE)
ranksWithPVal <- RandomWalkRestart(wm, c('g1', 'g2'), c(), TRUE, S=10)
ranksWithoutPval <- RandomWalkRestart(wm, c('g1'), c(),
       FALSE, 1, 0.5, 0.6, tau=c(1.5,0.5), phi=c(0.5,1.5))

bozdaglab/PhenoGeneRanker documentation built on March 11, 2021, 5:37 a.m.