Description Usage Arguments Details Value References Examples
The R implementation of the FASJEM method, which is introduced in the paper "A Fast and Scalable Joint Estimator for Learning Multiple Related Sparse Gaussian Graphical Models". Please run demo(fasjem) to learn the basic functions provided by this package. For more details, please see <http://proceedings.mlr.press/v54/wang17e/wang17e.pdf>.
1 |
X |
A List of input matrices. They can be either data matrices or covariance matrices. If every matrix in the X is a symmetric matrix, the input matrices are assumed to be the covariance matrices from the multiple related tasks. |
method |
By using two different regularization functions as the second norm in the objective, this package provides two different options for regularizing the sparsity pattern shared among multiple graphs. This parameter decides which function to use for the second regularization norm. When When |
lambda |
A positive number. This hyperparameter controls the sparsity level of the matrices. The λ_n in the following section: Details. |
epsilon |
A positive number. This hyperparameter represents the ratio between the l1 norm and the second group norm. The ε in the following section: Details. |
gamma |
A positive number. This hyperparameter is used in calculating each proximity during optimization. Please check the Algorithm 1 in our paper for more details. |
rho |
A positive number. This hyperparameter controls the learning rate of the proximal gradient method. Please check the Algorithm 1 in our paper for more details. |
iterMax |
An integer. The max number of iterations in the optimization of fasjem. |
The FASJEM algorithm is a fast and scalable method to estimate multiple related sparse Gaussian Graphical models. It solves the following equation:
\min\limits_{Ω_{tot}} ||Ω_{tot}||_1 + ε \mathcal{R}'(Ω_{tot})
Subject to :
||Ω_{tot} - inv(T_v(\hat{Σ}_{tot}))||_{∞} ≤ λ_n
\mathcal{R}'^*(Ω_{tot} - inv(T_v(\hat{Σ}_{tot}))) ≤ ελ_n
More details are provided in the equation (3.1) of our original paper.
The λ_n in the above equation represents the hyperparameter lambda
who controls the sparsity level of the target precision matrices.
The ελ_n in the above equation represents the regularization parameter of the second norm who controls how multiple graphs share a certain pattern. Here ε represents the input parameter epsilon
whose default value is 0.1.
Other parameters in the fasjem function are described in details by the Algorithm 1 in our paper.
When method = "fasjem-g"
, \mathcal{R}'(\cdot) = ||\cdot||_{\mathcal{G},2}.
When method = "fasjem-i"
, \mathcal{R}'(\cdot) = ||\cdot||_{\mathcal{G},∞}.
Please run demo(fasjem)
to learn the basic functions provided by this package. For more details, please see <http://proceedings.mlr.press/v54/wang17e/wang17e.pdf>.
Graphs |
A list of the estimated inverse covariance matrices. |
Beilun Wang, Ji Gao, Yanjun Qi (2017). A Fast and Scalable Joint Estimator for Learning Multiple Related Sparse Gaussian Graphical Models. <http://proceedings.mlr.press/v54/wang17e/wang17e.pdf>
1 2 3 | data(exampleData)
fasjem(X = exampleData, method = "fasjem-g", 0.1, 0.1, 0.1, 0.05, 10)
fasjem(X = exampleData, method = "fasjem-i", 0.1, 0.1, 0.1, 0.05, 10)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.