Description Usage Arguments Details Value References See Also Examples
Infer networks from genomics data using Markov Random Fields specified by node-conditional univariate exponential family distributions.
1 2 3 |
X |
a pxn data matrix. |
method |
specification of the type of MRF model, default to " |
stability |
specification of the stability method, default to " |
N |
number of iterations for stability selection, default to 100. |
beta |
threshold value on sparsity of the network, default to 0.01. |
lmin |
ratio of minimum lambda value from the maximum lambda value, default to 0.01. |
nlams |
number of lambda for regularization, default to 20. |
lambda.path |
vector lambda used for regularization, default ot |
parallel |
logical value to indicate if the process should be run parallelly in multiple threads, default to |
nCpus |
number of (maximum) cores to use for parallel execution, default to 4. |
sym |
logical value to indicate if symmetry is enforced on the inferred edges, default to |
th |
threshold value for the estimated edge coefficient, default to 0.005. |
sth |
an inferred edge is retained only if its stability score is greater than |
R |
truncation level for classes " |
R0 |
lower-bound truncation level for " |
This is the main function of the package that fits exponential family Markov Networks to genomics data. To estimate the network structures using native distribution of the genomics data,
specify the MRF family types in the "method
" parameter. For genomic networks based on next-generation sequencing data, we recommend using the LPGM
family.
The table at the beginning of the document lists the family type recommended for each of the genomic data platforms.
An object of class GMS
will be returned representing the inferred Markov networks over the regularization path. See GMS
for details.
Allen, G.I., and Liu, Z. (2012). A Log-Linear graphical model for inferring genetic networks from high-throughput sequencing data. The IEEE International Conference on Bioinformatics and Biomedicine (BIBM 2012).
Allen, G. I., and Liu, Z. (2013). A Local Poisson Graphical Model for Inferring Genetic Networks from Next Generation Sequencing Data. IEEE Transactions on NanoBioscience, 12(3), pp.1-10
Liu, H., Roeder, K., and Wasserman, L. (2010). Stability approach to regularization selection (stars) for high dimensional graphical models. NIPS 23, pp.1432?1440.
Meinshausen, N. and Buhlmann, P. (2006). High-dimensional graphs and variable selection with the lasso. The Annals of Statistics, 34(3), pp.1436?1462.
Meinshausen, N. and Buhlmann, P. (2010). Stability selection. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 72(4), pp.417?473.
Ravikumar, P., Wainwright, M., and Lafferty, J. (2010). High-dimensional ising model selection using l1-regularized logistic regression. The Annals of Statistics, 38(3), pp.1287?1319.
Yang, E., Ravikumar, P.K., Allen, G.I., and Liu, Z. (2012). Graphical models via generalized linear models. NIPS, 25, pp.1367–1375.
Yang, E., Ravikumar, P.K., Allen, G.I., and Liu, Z. (2013a). On graphical models via univariate exponential family distributions. arXiv preprint arXiv:1301.4183.
Yang, E., Ravikumar, P.K., Allen, G.I., and Liu, Z. (2013b). On Poisson graphical models. NIPS, pp.1718-1726.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | # Example for LPGM
# Refer to the package's introduction for identical example
## Not run: n = 100
## Not run: p = 20
## Not run: sim <- XMRF.Sim(n=n, p=p, model="LPGM", graph.type="scale-free")
## Not run: simDat <- sim$X
## Not run: # Compute the optimal lambda
## Not run: lmax = lambdaMax(t(simDat))
## Not run: lambda = 0.01* sqrt(log(p)/n) * lmax
## Not run: # Run LPGM
## Not run: lpgm.fit <- XMRF(simDat, method="LPGM", N=10, lambda.path=lambda)
## Not run: ml = plotNet(sim$B, fn="simDat.netPlot.pdf")
## Not run: ml = plot(lpgm.fit, fn="lpgm.netPlot_1.pdf", i=1, mylayout=ml)
## Not run: plot(lpgm.fit, fn="lpgm.fit.net.pdf")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.