JMGM: Joint Mixed Graphical Models

In equSA: Learning High-Dimensional Graphical Models

Description

Infer network structures from multiple datasets with mixed types of variables and edge restrictions option available.

Usage

JMGM(data,ALPHA1=0.05,ALPHA2=0.01,restrict=FALSE,parallel=FALSE,nCPUs)

Arguments

data	a list of nxp data matrices. n can be different for each dataset but p should be the same.
ALPHA1	The significance level of correlation screening. In general, a high significance level of correlation screening will lead to a slightly large separator set S_{ij}, which reduces the risk of missing some important variables in the conditioning set. Including a few false variables in the conditioning set will not hurt much the accuracy of the ψ-partial correlation coefficient.
ALPHA2	The significance level of ψ screening.
restrict	Should edge restriction applied? (logical). If TRUE, we assume that there should be no edge among binary variables. The default is FALSE.
parallel	Should parallelization be used? (logical), default is FALSE.
nCPUs	Number of cores used for parallelization. Recommend to be equal to the number of datasets.

Details

This is the function that can jointly estimate multiple graphical models with mixed types of data and also consider the edge restriction scenarios. The method has three novelties: First, the proposed method resolves the conditional independence information using a p-learning algorithm and therefore can be applied to the mixed types of random variables. Second, the proposed method can construct networks with restricted edges determined by some preliminary knowledges. Third, the proposed method involves a Fast Bayesian joint estimation method which works on edge-wise scores and can achieve both fast and accurate integration performance for constructing multiple networks. See Jia and Liang (2018).

Value

A list of three elements:

A	An array of multiple adjacency matrices of networks which is a Mxpxp array. M is the number of dataset groups, p is the dimension of variables in each group.
score.sep	Separately estimated ψ scores matrix for all pairs in multiple datasets. The first two columns denote the pair indices of variables i and j and the rest columns denote the estimated ψ scores for this pair in different groups.
score.joint	Estimated integrative ψ scores matrix for all pairs in multiple datasets. The first two columns denote the pair indices of variables i and j and the rest columns denote the estimated integrative ψ scores for this pair in different groups.

Author(s)

Bochao [e-mail] and Faming Liang

References

Jia, B., and Liang, F. (2018) Joint Estimation of Restricted Mixed Graphical Models. manuscript.

Examples

library(equSA)
data1 <- DAGsim(n=200, p=100, type="AR(2)")$data
data2 <- DAGsim(n=200, p=100, type="AR(2)")$data
data_all <- vector("list",2)
data_all[[1]] <- data1
data_all[[2]] <- data2
JMGM(data_all,ALPHA1=0.1,ALPHA2=0.05,parallel=TRUE,nCPUs=2)




         

equSA documentation built on May 6, 2019, 1:06 a.m.