In order to determine the statistical significance of the mutual information values between genes we test for each pair of genes the following null hypothesis.
H_0^I: The mutual information between gene i and j is zero.
Because we are using a nonparametric test we need to obtain the corresponding null distribution for H_0^I from a randomization of the data.
'makenull' performs a randomization with the formulated null hypothesis by permuting the sample and gene labels for all genes of the entire expression matrix at once. The vector of the mutual information null distribution is obtained from repeated randomizations for a given number of iterations.
1 |
dataset |
gene expression data matrix |
nullit |
nullit defines the size of the generated null distribution vector used for hypothesis testing of significant edges inferred by c3net. The null distribution of mutual information is generated from sample and gene label randomization. default number of iterations: nullit=ceiling(10^5/(((genes*genes)/2)-genes)) genes: number of genes |
estimator |
estimators for continuous variables "pearson" (default), "spearman", "kendall", "spearman" estimators for discrete variables "emp", "mm","sg","shrink" |
disc |
only required for discrete estimators, method for discretize function (infotheo package) "equalfreq" (default), "equalwidth", "globalequalwith" |
The function 'makenull' generates a numeric vector of the mutual information null distribution estimated from a given number of randomized datasets.
Ricardo de Matos Simoes <r.dematossimoes@qub.ac.uk> Frank Emmert-Streib <f.emmert-streib@qub.ac.uk>
de Matos Simoes R, Emmert-Streib F., Bagging statistical network inference from large-scale gene expression data., PLoS One. 2012;7(3):e33624. Epub 2012 Mar 30.
c3mtc bc3net
1 2 |
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.