A novel sparse partial correlation method for ...

5/5/2016

cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

A novel sparse partial correlation method for simultaneous estimation of functional networks in group comparisons Xiaoyun Liang1, David Vaughan2,3, Alan Connelly1,4, and Fernando Calamante1,4 1Imaging

Division, Florey Institute of Neuroscience and Mental Health, Melbourne, Australia, 2Epilepsy Division, Florey Institute of Neuroscience and Mental Health, Melbourne, Australia, 3Department of Neurology, Austin Health, Melbourne, Australia, 4Department of Medicine, University of Melbourne, Melbourne, Australia

Synopsis

We propose a novel approach, Graphical-LAsso with Stability-Selection (GM-GLASS), by employing sparse group penalties for simultaneously estimating networks from healthy control and patient groups. Simulations demonstrate that both GM-GLASS and JGMSS outperform Fisher Z-transform. Our in vivo results further show that GM-GLASS yields highest contrast of network metrics between groups, demonstrating the superiority of GM-GLASS in detecting signi䁀cance group di၀erences over JGMSS and Fisher Z-transform. Overall, by controlling confounding variations between subjects, and therefore enhancing the statistical power, our simulated and in vivo results demonstrate that GMGLASS provides a robust approach for conducting group comparison studies. Purpose

Functional networks at group-level can be used to gain insight into complex brain function, including in group comparison studies. Joint graphical model with stability selection (JGMSS) was shown to greatly control confounding variations when computing group-level connectomic networks [1]. However, JGMSS deals with one single group at a time. Recently, a fused graphical-lasso model (FGM) was proposed to robustly estimate group-level graphs from multiple groups simultaneously [2,3]. We aim to simultaneously estimate individual- and group-level networks from 2 groups using FGM. We propose a novel approach, Group-fused-Multiple-Graphical-LAsso with Stability-Selection (GM-GLASS), by employing sparse group penalties for simultaneously estimating connectivity matrices from within 2 contrast groups, e.g. normal control (NC) and patient (PT) group, while controlling for confounding variations. Methods

GM-GLASS can be described as follows (see Fig. 1 for ᶸowchart): 2K datasets, X

(1)

for 2 groups (ᶠrst K from group 1, last K from group 2), X (k) : matrix, nk : number of time-points, p: number of regions (nodes). The empirical covariance-

,...X

nk × p

(K)

,X

(K+1)

matrix for X (k) is: S (k) [2,3]: max{Θ} (∑k

,...,X

2K

= 1/nk (X

nk (logdetΘ

(k)

,

(k)

)

T

X

(k)

− tr(S

, and the maximum-likelihood estimate of (S (k) )−1 is

(k)

Θ

(k)

,

)) − P ({Θ})) k = 1, . . . , 2K

, with the

constraint that all Θ (k) are positive-deᶠnite; the fused regularization P ({Θ}) [2,3] is reformulated as: P ({Θ}) = λ1 ∑

k

∑

i≠j

|Θ

(k) ij

| + λ2 ∑

k1

∑

i≠j

|Θ

(k1 ) ij

− Θ

(k1 +K) ij

|

with λ1 , λ2 : nonnegative

regularization-parameters, and k = 1, . . . , 2K , and k1 = 1, . . . , K . With stability selection [4], the data are subsampled many times and all positive-correlations that occur in a large fraction of the resulting selection sets are selected. The combination patterns (i.e. the order of subjects) in the second part of P are shuᷘed for each resampling to avoid bias. For a given (λ1 , λ2 ) , λ ,λ S = {(m, n) : 1 ≤ (m, n) ≤ M ; Θ(m, n) < 0} . For a cut-oᷘ Pthr ∈ (0, 1) and a set of regularization parameters, the set of stable positive-correlations are as follows: 1

2

ˆ S stable = {(m, n) : max∀(λ

1 ,λ2 )

λ1 ,λ2

ˆ P m,n

≥ Pthr }

.

Simulation: Simulations were conducted according to [6]. A structural-connectivity matrix [5] was employed to simulate data for NCs (see [6]); by randomly deleting (adding) 25 existing (new) edges, patient data were then simulated, with p = 78 and n = 200 . Ten datasets were generated (i.e. ‘subjects’) for each group. To perform stability selection, n/2 observations were randomly subsampled http://cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

1/4

5/5/2016

cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

without overlapping for each dataset, repeating 100 times to estimate probabilities over the regularization region, with the expected number of falsely-selected connections V bounded by [4], r: number of variables, q: average number of selected connections for a given range of (λ1 , λ2 ) . The per-comparison-error-rate (P C ER = E(V )/r) is employed to control falsely-selected connections [4]. In-vivo data: Resting-state BOLD-fMRI data were acquired from 13 mesial temporal lobe epilepsy (TLE) patients with right hippocampal-sclerosis, as well as 13 age-matched healthy controls, on a 3T Siemens scanner: TR/TE=3000/30 ms, 3mm isotropic voxels, 44 slices, matrix-size: 72 × 72 , 200 time-points. Image analysis: images were preprocessed in a standard way [7]. With individual-level mean time-series within AAL template [8] as inputs, GM-GLASS was employed for estimating individual- and group-level functional networks. Since fMRI signal is not independent and not exchangeable, the time-series were divided into 20 consecutive blocks (200/20 = 10 ). For comparisons, JGMSS [1] and Pearson-cross-correlation followed by Fisher-Z-transformation (Fisher Z-transform) [9] were also employed. For Fisher Z-transform, individual-level networks were averaged across the group, followed by thresholding the group-level networks with multiple-comparison correction. 2

E(V ) = q /(2Pthr ∗ r − r)

Results

Fig. 2 shows estimated networks from simulated normal and patient groups, demonstrating the capability of GM-GLASS in detecting network diᷘerences. The receiver-operating curves (ROCs) in Fig. 3 demonstrate that GM-GLASS and JGMSS estimate true connections from groups more reliably than Fisher Z-transform. Based on the ground truth, simulations have demonstrated that GM-GLASS and JGMSS achieve higher accuracy (75 %) than Fisher Z-transform (50 %) (Fig. 4). In-vivo results show that GM-GLASS has detected 3 signiᶠcantly diᷘerent network metrics between groups, instead of 1 for JGMSS and 0 for Fisher Z-transform (Fig.5), further suggesting that GM-GLASS might have higher sensitivity in detecting group network diᷘerences. Discussion

Simulations demonstrate that GM-GLASS and JGMSS signiᶠcantly outperform Fisher Z-transform (Fig. 3), and they are capable of detecting minor network diᷘerences (Fig. 2); only marginally better performance of GM-GLASS than JGMSS is achieved, likely due to high SNR of simulated data. However, for in vivo data, networks estimated using GM-GLASS are characterized by increased contrast of network metrics between groups compared to the other 2 methods. More speciᶠcally, decreased local-eᶸciency and globaleᶸciency in right TLE are consistent with [10] and decreased network density is in line with [11]; these highlight the capability of GM-GLASS in detecting signiᶠcant diᷘerences between groups that otherwise might not be detected using JGMSS or Fisher Z-transform. Overall, by controlling confounding variations between subjects, and therefore enhancing the statistical power, our simulated and in vivo results demonstrate that GM-GLASS provides a robust approach for conducting group comparison studies. Acknowledgements

We are grateful to the National Health and Medical Research Council (NHMRC) of Australia, the Australian Research Council (ARC), and the Victorian Government's Operational Infrastructure Support Program for their support. References

[1] Liang X, Connelly A, Calamante F. A novel method for robust estimation of group fucntional connectivity based on a joint graphical models approach. Proc. ISMRM p.3059 (2014). [2] Danaher P, Wang P, Witten DM. The joint graphical lasso for inverse covariance estimation across multiple classes. Journal of the Royal Statistical Society Series B-statistical Methodology 76(2): 373-97 (2014).

http://cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

2/4

5/5/2016

cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

[3] Yang S, Lu Z, Shen X, et al. Fused multiple graphical lasso. SIAM Journal on Optimization 25(2): 916-43 (2015). [4] Meinshausen N, Buhlmann P. Stability selection. Journal of the Royal Statistical Society Series BStatistical Methodology 72(4): 417-73 (2010). [5] Gong G, He Y, Concha L et al. Mapping anatomical connectivity patterns of human cerebral cortex using in vivo diᷘusion tensor imaging tractography. Cerebral Cortex 19: 524-36 (2009). [6] Honey CJ, Sporns O, Cammoun L, et al. Predicting human resting-state functional connectivity from structural connectivity. Proc Natl Acad Sci USA 106(6) 2035-40 (2009). [7] Zuo X, Xu T, Jiang L, et al. Toward reliable characterization of functional homogeneity in the human brain: preprocessing, scan duration, imaging resolution and computational space. NeuroImage 65: 37486 (2012). [8] Tzourio-Mazoyer N, Landeau B, Papathanassiou D, et al. Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. NeuroImage 15(1): 273-89 (2002). [9] Fisher RA. On the `probable error' of a coeᶸcient of correlation deduced from a small sample. Metron 1: 3-32 (1921). [10] Liu M, Chen Z, Beaulieu C, et al. Disrupted anatomic white matter network in left mesial temporal lobe epilepsy. Epilepsy 55:674-82 (2014). [11] Besson P, Dinkelacker V, Valabregue R, et al. Structural connectivity diᷘerences in left and right temporal lobe epilepsy. NeuroImage 100: 135-44 (2014). Figures

Figure 1: The 䁀owchart shows the implementation of GM-GLASS in which both sparse individual- and group-level brain functional networks (Sstable) are calculated. Note: the second part of regularization term P was reformulated to incorporate shu䁀ed network di၀erences between groups, achieving the goal of reducing both within- and between-group confounding variability.

Figure 2: Orange and red colors indicate common true edges and truly detected edge di၀erences between groups; blue and yellow colors represent false edges and falsely detected di၀erences between groups (i.e. expected true common edges), with accuracy and sensitivity for NC and PT groups, respectively: 0.9924, 0.9938, 0.9666 and 0.9787.

http://cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

3/4

5/5/2016

cds.ismrm.org/protected/16MPresentations/abstracts/3745.html

Figure 3: The receiver-operating curves (ROC) for NC (left) and PT (right). Estimation of area under the curve (AUC) indicates that GM-GLASS yields the best performance, with accuracy and sensitivity for NC and PT groups, respectively: JGMSS: 0.9911, 0.9928, 0.9605 and 0.9757; Fisher Z-transform: 0.9770, 0.9852, 0.8815 and 0.8298.

Figure 4: Statistical di၀erences of 4 individual-level network metrics estimated from simulated data, i.e. local e䁀ciency (LE), global e䁀ciency (GE), network density and clustering coe䁀cient (CC), between NC and PT groups. Accuracy: 75% for GM-GLASS and JGMSS; 50% for Fisher Z-transform (FisherZ). Note: *p