The R package lglasso is designed to estimate precision matrices (i.e., networks) from high-dimensional longitudinal data under a normal distribution. It can be used to estimate either one-stage models, where a single network underlies all longitudinal measurements, or two-stage models, where the network before treatment differs from the one after treatment (exposure). Both model types can be further classified as homogeneous—treating all individuals as i.i.d. samples—or heterogeneous, which includes individual-level random effects. For details on these definitions and usage, please check the references and links below. The ADMM algorithm is adopted to speed up computation, and the main functions are written in C++. If you have any questions or suggestions, please email chowstat@gmail.com.
You can install the development version of lglasso in R from GitHub with:
First, install R package remotes:
install.packages("remotes")
Then install lglasso :
remotes::install_github("jiezhou-2/lglasso", ref ="main")
Here are some examples for the usage of the package. If you are interested in the underlying statistical model of lglasso, you could check the our manuscript here draft, which is an extension of the paper in the Reference.
Reference
[1] Zhou J et al. Identifying stationary microbial interaction networks based on irregularly spaced longitudinal 16S rRNA gene sequencing data. Front Microbiomes. 2024;3:1366948. doi: 10.3389/frmbi.2024.1366948. Epub 2024 Jun 2. PMID: 40687607; PMCID: PMC12276884.