PLNet is a package for calculating the precision matrix of poission log-normal model using moment method and dtrace loss. We first calculate the covariance matrix by moment method and then calculate the precision matrix by dtrace loss using R package EQUAL.
PLNet is available on Github.
- Before installing PLNet package, please install cescwang85-EQUAL-da9cbbb.tar.gz on Github.
- For the development version, use the github install
remotes::install_github("https://github.com/XiDsLab/PLNet")
The package comes with a simulation data and a single-cell RNA sequence data about Interferon β-Stimulated PBMCs to present the functionality of main function.
library(PLNet)
## simulation data for testing
data("obs_mat")
PLNet_res_sim<-PLNet::PLNet(obs_mat = PLNet::obs_mat,
Sd_est = "GMPR",lambda_vec = NULL , n_lambda = 50 , lambda_valuemax0 = 10,penalize.diagonal = TRUE)
## real data analysis
data("realdatamatrix")
PLNet_res_real<-PLNet::PLNet(obs_mat = PLNet::realdatamatrix,
Sd_est = "GMPR",lambda_vec = NULL , n_lambda = 50 , lambda_valuemax0 = 10,penalize.diagonal = TRUE)
library(MASS)
data_generator_res<-PLNet::data_generator(n = 100, p = 50, densy_degree = 0.05,
sd_ls = 0.1, mean_ls = log(10),value_nondiagonal = 0.3)
Please cite our work using the following references: