Invertible Network for Classification and Biomarker Selection for ASD

07/23/2019 ∙ by Juntang Zhuang, et al. ∙ 0

Determining biomarkers for autism spectrum disorder (ASD) is crucial to understanding its mechanisms. Recently deep learning methods have achieved success in the classification task of ASD using fMRI data. However, due to the black-box nature of most deep learning models, it's hard to perform biomarker selection and interpret model decisions. The recently proposed invertible networks can accurately reconstruct the input from its output, and have the potential to unravel the black-box representation. Therefore, we propose a novel method to classify ASD and identify biomarkers for ASD using the connectivity matrix calculated from fMRI as the input. Specifically, with invertible networks, we explicitly determine the decision boundary and the projection of data points onto the boundary. Like linear classifiers, the difference between a point and its projection onto the decision boundary can be viewed as the explanation. We then define the importance as the explanation weighted by the gradient of prediction w.r.t the input, and identify biomarkers based on this importance measure. We perform a regression task to further validate our biomarker selection: compared to using all edges in the connectivity matrix, using the top 10% important edges we generate a lower regression error on 6 different severity scores. Our experiments show that the invertible network is both effective at ASD classification and interpretable, allowing for discovery of reliable biomarkers.



There are no comments yet.


page 1

page 2

page 3

page 4

This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

1 Introduction

Autism spectrum disorder (ASD) is a neurodevelopmental disorder that affects social interaction and communication, yet the causes for ASD are still unknown [16]. Functional MRI (fMRI) can measure the contrast dependent on blood oxygenation [15] and reflect brain activities, and therefore has the potential to help in understanding ASD.

Recent research efforts have applied machine learning and deep learning methods on fMRI data to classify ASD versus control groups

[1, 11] and predict treatment outcomes [19, 18]. However, deep learning models are typically hard to interpret, and thus are difficult to use for identifying biomarkers for ASD.

Various methods have been proposed to interpret a deep neural network. Bach et al. proposed to assign the decision of a neural network to its input with layer-wise relevance propagation

[2]; Mahendran et al. proposed to approximately invert the network to explain its decision [13]; Sundararajan et al. proposed integrated gradient to explain a model’s decision [17]. However, all these methods only generate an approximation to the inversion of neural networks, and can not unravel the black-box representation of deep learning models.

The recently proposed invertible network can accurately reconstruct the input from its output [12, 3]. Based on this property, we propose a novel method to interpret ASD classification on fMRI data. As shown in Fig. 2, an invertible network first maps data from the input domain (e.g. connectivity matrix calculated from fMRI data) to the feature domain, then applies a fully-connected layer to classify ASD from control group. Since a fully-connected layer is equivalent to a linear classifier in the feature domain, we can determine the decision boundary as a high-dimensional plane, and calculate projection of a point onto the boundary. Since our network is invertible, we can invert the decision boundary in the feature domain to the input domain. As shown in Fig. 3, the difference between the input and its projection on the decision boundary can be viewed as the explanation for the model’s decision. We applied the proposed method on ASD classification of the ABIDE dataset achieving 71% accuracy, and then identified biomarkers for ASD based on the proposed interpretation method. We further validated the selected biomarkers (edges in the connectivity matrix as in Fig. 6 and ROIs as in Fig. 7) in a regression task: compared with using all edges, the selected edges generate more accurate predictions for ASD severity (Table 2).

Our contributions can be summarized as:

  • Based on the invertible network, we proposed a novel method to intepret model decision and rank feature importance.

  • We applied the proposed method on an ASD classification task, achieved a high classification accuracy (71% on the whole ABIDE dataset), and identified biomarkers.

  • We demonstrated effectiveness of selected biomarkers in a regression task.

2 Methods

The classification pipeline is summarized in Fig. 2. In this section, we first discuss input data and pre-processing of fMRI, then introduce the structure of invertible networks as in Fig. 2. Next, we propose a novel method to determine the decision boundary and identify biomarkers for ASD.

2.1 Dataset and Inputs

The ABIDE dataset [9] consists of fMRI data for 530 subjects with ASD and 505 control subjects. We use the pre-processed data by the Connectome Computation System (CCS) [20] pipeline, which includes registration, 4D global mean-based intensity normalization, nuisance correction and band-pass filtering. We extract the mean time series for each region of interest (ROI) defined by the CC200 atlas [8]

, consisting of 200 ROIs. We compute the Pearson Correlation between the time series of every pair of ROIs, and reshape the connectivity into a vector of length

. This vector is the input to the invertible network.

Figure 1: Classification pipeline for fMRI and structure of invertible network.
Figure 2: Structure of invertible blocks.

2.2 Invertible Networks

The structure of the invertible network is shown in Fig. 2. An invertible network is composed of a stack of invertible blocks and a final fully-connected (FC) layer to perform classification. Note that inversion does not update parameters, thus is different from the backward propagation method; invertible networks can be trained with the backward propagation method as normal networks.

2.2.1 Invertible Blocks

An invertible block is shown in Fig. 2, where the input is split by channel into two parts and , and the outputs are denoted as and . Feature maps have the same shape. and are non-linear functions with parameters to learn: for 1D input, and

can be a sequence of FC layers, 1D batch normalization layers and activation layers; for 2D input,

and can be a sequence of convolutional layers, 2D batch normalization layers and activation layers. and are required to generate outputs of the same shape as input. The invertible block can accurately recover the input from its output, where the forward pass and inverse of an invertible block are:


2.2.2 Notations

The invertible network classifier can be viewed as a 2-stage model:


where is the invertible transform and is the final FC layer which is a classifier. We denote as the input domain (e.g. connectivity matrix from fMRI reshaped to 1D vector), and as the feature domain. A data point is mapped from the input domain to the feature domain by an invertible transform , then a linear classifer is applied.

2.3 Model Decision Interpretation and Biomarker Selection

In this section we propose a novel method to interpret decisions of an invertible network classifier. We begin with a linear classifier, then generalize to non-linear classifier such as neural networks. An example is shown in Fig. 3.

2.3.1 Interpret Decision of a Linear Classifier

A linear classifier is calculated as:


where is the weight vector, and is the bias. The decision boundary is a high dimensional plane, and can be denoted as .

For a data point , we calculate its projection onto the decision boundary as


and define the and as:


where is element-wise product, the difference can be viewed as the explanation for the linear classifier, and the of each input dimension is defined as the absolute value of the weighted by . An example is shown in Fig. 3(a).

2.3.2 Determine Projection onto Decision Boundary of Invertible Networks

As in Equation 2, a point in the input domain is mapped to the feature domain, denoted as . Since the classifier is linear in the feature domain, we can calculate the projection of onto the boundary as in Equation 4, denoted as . Note that since is invertible, we can map it to the input domain, denoted as where is the inverse operation of as in Equation 1. Furthermore, we can invert the decision boundary from the feature domain to the input domain.

An example is shown in Fig. 3. Fig. 3(a) shows a data point and its projection in the feature domain, Fig. 3(b) inverts the projection back to the input domain. Fig. 3(c) and (d) show the decision boundary in the feature and input domain, respectively.

2.3.3 Feature Importance and Biomarker Selection

We expand the network around current point using Taylor expansion:


Similar to Equation 5, we approximate locally as a linear classifier and define the and as:


Equation 7 defines individualized importance for each data point, and we calculate the mean importance across the entire dataset. We select edges with top importance values as the biomarkers for ASD.

Note that two classes may be close in the input but separated in the explanation. We show two examples in Fig.  4. In the top row, the decision boundary is a horizontal line, thus is useful to distinguish the two clusters while is not. In this case, the distribution of the two clusters are overlapped in the axis but separated in the axis for both input and the explanation. In the bottom row, two clusters are separated by line : neither nor axis can distinguish two clusters using the input (Fig. 4(b),(d)), but both axes have a large separation margin using the explanation(Fig. 4(c),(e)).

Figure 3: Illustration on a simulation dataset. Fig. (a): data point and its projection onto the decision boundary. In the feature domain the model is a linear classifier. Fig. (b) Corresponding points and in the input domain, calculated as inversion of points from (a). Fig. (c) and (d): Points from two classes are sampled around two interleaving half circles. Decision boundary is a line in the feature domain (c), and is a curve in the input domain (d).
Figure 4: Explanation method on simulation datasets. Columns (a) to (e) represent data distribution, input value in axis, explanation in axis, input value in axis and explanation in axis. In the top row, both input and explanation fail to distinguish two classes in axis, but succeed in axis. In the bottom row, both and values are useful; two clusters are overlapped in distribution of () input values, but separated in the explanation for () axis.

3 Experiments

3.0.1 Classification Accuracy

We perform a 10-fold cross validation on the entire ABIDE dataset to classify ASD from control group. We compare the invertible network classifier with other methods including SVM, random forest (RF) with 5,000 trees, a 1-layer multi-layer perceptron (MLP) and a 2-layer MLP. Our invertible network (InvNet) has 2 invertible blocks, and the


in an invertible block is FC-ReLU-FC-ReLU, where the first FC layer maps a 19,900 (edges in a connectivity matrix using CC200 atlas with 200 ROIs


) dimensional vector to a 512 dimensional vector, and the second FC maps a vector from 512 dimension to 19,900. For each of the 1035 subjects, we augment data 50 times by bootstrapping voxels within each ROI, then calculating the connectome from the bootstrapped mean time series. This results in 51570 examples. All deep learning models are trained with SGD optimizer for 50 epochs, with a learning rate of 1e-5 and a cross-entropy loss.

Numerical results are summarized in Table 1. Compared with other methods including a Deep Neural Network (DNN) model in [11], the proposed InvNet generates the highest accuracy (0.71), recall (0.71) and F1 score (0.71).

Table 1: Classification results
Accuracy 0.67 0.66 0.66 0.68 0.70 0.71
Precision 0.68 0.65 0.66 0.68 0.74 0.72
Recall 0.68 0.65 0.66 0.67 0.63 0.71
F1 0.68 0.68 0.67 0.69 - 0.71
Table 2: Regression results for ASD severity scores
# Edges Total Awareness Cognition Comm Motivation Mannerism
MSE 100% 44.2 3.86 6.90 11.95 7.59 6.37
RF(top 10%) 71.8 8.2 12.6 22.1 12.7 11.8
ours 42.8 3.72 6.84 11.70 7.17 6.19
Cor 100% 0.14 0.33 0.38 0.20 -0.03 0.22
RF(top 10%) 0.09 -0.19 0.01 -0.01 -0.12 -0.08
ours 0.18 0.42 0.40 0.28 0.20 0.32
Figure 5: Top row: distribution of connectivity edge for two groups (ASD v.s. Control). Bottom row: distribution of explanation for two groups. From left to right: results on edges with top 1, top 2 and top 3 importance defined in Equation 7.

3.0.2 Explanation and Biomarker Selection

We plot the histogram of explanation for edges with top 3 importance in Fig. 5. The histogram of edge values (top row) can not distinguish ASD from control group, while the distribution of explanations (bottom row) for the two groups are separated. The proposed explanation method can be viewed as a naive embedding method, mapping data to a low-dimensional space where two classes are separated.

For each edge, we calculate the importance as defined in Equation 7 and select the top 20 connections, as shown in Fig. 6, and plot the corresponding ROIs in Fig. 7. The proposed method found many regions that are shown to be closely related to autism in the literature, including: superior temporal gyrus [4], frontal cortex [7, 5], precentral gyrus [14], insular cortex [10] and other regions. The selected edges and ROIs can be viewed as biomarkers for ASD. Detailed results are in the appendix.

Figure 6: Top 20 connections selected by the proposed method. Figure 7: Top ROIs selected by the proposed method.

3.0.3 Validate Biomarker Selection in Regression

We validate our biomarker selection in a regression task on the ABIDE dataset. The input is the connectivity matrix reshaped to a vector, and the output is different subscores of the social responsiveness scale (SRS) [6]. SRS provides a continuous measure of different aspects of social ability, including a total score and subscores for awareness, cognition, communication, motivation and mannerism. We use the top 10% of edges selected by the proposed method, and compared its performance with using 100% of edges. We perform a 10-fold cross validation with linear support vector regression (SVR) using penalty, and within each fold the penalty parameter is chosen by nested cross-validation (choices for penalty strength are ).

We validate our model by comparison with other feature selection methods. We first selected top 10% important features with RF based on out-of-bag prediction importance with default parameters in MATLAB, then refitted a RF using selected features. Results are marked as RF(top 10%) in Table 2.

We calculate the mean squared error (MSE) and cross correlation (Cor) between predictions and measurements. Results are summarized in Table 2. Compared with other methods, our selected biomarkers consistently generates a smaller MSE and a higher Cor, validating our biomarker selection method.

3.0.4 Generalization to Convolutional Invertible Networks

We generalize the proposed model decision interpretation method to convolutional networks, and validate it on an image classification task with the MNIST dataset. Results are summarized in the appendix. The proposed method generates more intuitive explanation results on 2D images.

4 Conclusions

We introduced a novel decision interpretation and feature importance ranking method based on invertible networks. We then applied the proposed method on a classification task to classify ASD from control group based on fMRI scans. We selected important connections and ROIs as biomarkers for ASD, and validated these biomarkers in a regression task on the ABIDE dataset. Our invertible network generates a high classification accuracy, and our biomarkers consistently generate a smaller MSE and a higher Cor compared with using all edges for regression tasks on different severity measures. The proposed interpretation method is generic, and has the potential in other aspects of interpretable deep learning such as 2D image classification.


This research was funded by the National Institutes of Health (NINDS-R01NS035193).


  • [1] Anderson, J.S., et al.: Functional connectivity magnetic resonance imaging classification of autism. Brain (2011)
  • [2] Bach, S., et al.: On pixel-wise explanations for non-linear classifier decisions by layer-wise relevance propagation. PloS one (2015)
  • [3] Behrmann, J., Duvenaud, D., Jacobsen, J.H.: Invertible residual networks. arXiv preprint arXiv:1811.00995 (2018)
  • [4] Bigler, E.D., et al.: Superior temporal gyrus, language function, and autism. Developmental neuropsychology (2007)
  • [5] Carper, R.A., Courchesne, E.: Localized enlargement of the frontal cortex in early autism. Biological psychiatry 57(2), 126–133 (2005)
  • [6] Constantino, J.N.: Social responsiveness scale. Springer (2013)
  • [7] Courchesne, E., Pierce, K.: Why the frontal cortex in autism might be talking only to itself: local over-connectivity but long-distance disconnection. Current opinion in neurobiology 15(2), 225–230 (2005)
  • [8]

    Craddock, R.C., et al.: A whole brain fmri atlas generated via spatially constrained spectral clustering. Human brain mapping (2012)

  • [9] Di Martino, A., et al.: The autism brain imaging data exchange: towards a large-scale evaluation of the intrinsic brain architecture in autism. Mol. Psychiatry (2014)
  • [10]

    Gogolla, N., et al.: Sensory integration in mouse insular cortex reflects gaba circuit maturation. Neuron (2014)

  • [11] Heinsfeld, A.S., Franco, A.R., Craddock, R.C., Buchweitz, A., Meneguzzi, F.: Identification of autism spectrum disorder using deep learning and the abide dataset. NeuroImage: Clinical 17, 16–23 (2018)
  • [12] Jacobsen, J.H., Smeulders, A., Oyallon, E.: i-revnet: Deep invertible networks. arXiv preprint arXiv:1802.07088 (2018)
  • [13] Mahendran, A., Vedaldi, A.: Understanding deep image representations by inverting them. In: CVPR (2015)
  • [14] Nebel, M.B., et al.: Precentral gyrus functional connectivity signatures of autism. Frontiers in systems neuroscience (2014)
  • [15] Ogawa, S., et al.: Brain magnetic resonance imaging with contrast dependent on blood oxygenation. Proceedings of the National Academy of Sciences (1990)
  • [16] Speaks, A.: What is autism. Retrieved on November 17,  2011 (2011)
  • [17] Sundararajan, M., Taly, A., Yan, Q.: Axiomatic attribution for deep networks. arXiv preprint arXiv:1703.01365 (2017)
  • [18] Zhuang, J., et al.: Prediction of pivotal response treatment outcome with task fmri using random forest and variable selection. In: ISBI 2018 (2018)
  • [19] Zhuang, J., et al.: Prediction of severity and treatment outcome for asd from fmri. In: International Workshop on PRedictive Intelligence In MEdicine (2018)
  • [20] Zuo, X.N., et al.: Toward reliable characterization of functional homogeneity in the human brain: preprocessing, scan duration, imaging resolution and computational space. Neuroimage 65 (2013)