GraphSPME: Markov Precision Matrix Estimation and Asymptotic Stein-Type Shrinkage
| Autor: | Lunde, Berent Ånund Strømnes, Curic, Feda, Sortland, Sondre |
|---|---|
| Rok vydání: | 2022 |
| Předmět: | |
| Druh dokumentu: | Working Paper |
| Popis: | GraphSPME is an open source Python, R and C++ header-only package implement-ing non-parametric sparse precision matrix estimation along with asymptotic Stein-type shrinkage estimation of the covariance matrix. The user defines a potential neighbourhood structure and provides data that potentially are p >> n. This paper introduces a novel approach for finding the optimal order (that data allows to estimate) of a potential Markov property. The algorithm is implemented in the package, alleviating the problem of users making Markov assumptions and implementing corresponding complex higher-order neighbourhood structures. Estimation is made accurate and stable by simultaneously utilising both Markov properties and Stein-type shrinkage. Asymptotic results on Stein-type shrinkage ensure that non-singular well conditioned matrices are obtained in an automatic manner. Final symmetry conversion creates symmetric positive definite estimates. Furthermore, the estimation routine is made efficient and scalable to very high-dimensional problems (~10^7) by utilising the sparse nature of the precision matrix under Markov assumptions. Implementation wise, the sparsity is exploited by employing the sparsity possibilities made available by the Eigen C++ linear-algebra library. The package and examples are available at https://github.com/equinor/GraphSPME Comment: 17 pages, 6 figures |
| Databáze: | arXiv |
| Externí odkaz: |
|