nem.calcSignificance: Statistical significance of network hypotheses

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/nem.calcSignificance.R

Description

Assess statistical significance of a network hypothesis by comparing it to a null hypothesis.

Usage

1
	nem.calcSignificance(D, x, N=1000, which.test=c("permutation", "rand.net", "rand.modify"), seed=1, mc.cores=8)

Arguments

D

data matrix with experiments in the columns (binary or continious)

x

nem object

N

number of random networks to sample

which.test

type of significance test to conduct (see details)

seed

random seed

mc.cores

number of cores to be used on a multicore processor

Details

Given data, N random network hypotheses from a null distribution are drawn as follows: For each S-gene S_{k} we randomly choose a number o of outgoing edges between 0 and 3. We then select o S-genes having at most 1 ingoing edge, connected S_{k} to them and transitively closed the graph. For all random network hypotheses it is counted, how often their likelihood is bigger than that of the given network. This yields an exact p-value.

Another way of assessing the statistical significance of the network hypothesis is to draw random permutations of the node labels. Note that in this case the node degree distribution is the same as in the given network. Again, we can obtain an exact p-value by counting, how often the likelihood of the permuted network is bigger than that of the given network.

Finally, comparison to randomly perturbed networks (insertion, deletion or reversion of 1 edge) yields an exact p-value describing the stability of the network. For dynoNEMs network modification operations are edge weight increase, decrease and swap.

Value

For "random" network: p-value of the network according to the null hypothesis that it is random. For permuted nework: p-value of the network according to the null hypothesis that a network with permuted node labels is at least as good. For perturbed / modified network: p-value of the network according to the null hypothesis a randomly peturbed network is at least as good.

Author(s)

Holger Froehlich

See Also

nem.consensus, nem.jackknife, nem.bootstrap, nem

Examples

1
2
3
4
5
6
7
## Not run: 
   data("BoutrosRNAi2002")
   D <- BoutrosRNAiDiscrete[,9:16]   
   res = nem(D, control=set.default.parameters(unique(colnames(D)), para=c(0.13,0.05))) # get best network
   nem.calcSignificance(D,res) # assess its statistical significance

## End(Not run)

cbg-ethz/pcNEM documentation built on Sept. 27, 2019, 8:58 a.m.