Dynamic causal modeling


Dynamic causal modeling is a framework for specifying models, fitting them to data and comparing their evidence using Bayesian model comparison. It uses nonlinear state-space models in continuous time, specified using stochastic or ordinary differential equations. DCM was initially developed for testing hypotheses about neural dynamics. In this setting, differential equations describe the interaction of neural populations, which directly or indirectly give rise to functional neuroimaging data e.g., functional magnetic resonance imaging, magnetoencephalography or electroencephalography. Parameters in these models quantify the directed influences or effective connectivity among neuronal populations, which are estimated from the data using Bayesian statistical methods.

Procedure

DCM is typically used to estimate the coupling among brain regions and the changes in coupling due to experimental changes. A model of interacting neural populations is specified, with a level of biological detail dependent on the hypotheses and available data. This is coupled with a forward model describing how neural activity gives rise to measured responses. Estimating the generative model identifies the parameters from the observed data. Bayesian model comparison is used to compare models based on their evidence, which can then be characterised in terms of parameters.
DCM studies typically involve the following stages:
  1. Experimental design. Specific hypotheses are formulated and an experiment is conducted.
  2. Data preparation. The acquired data are pre-processed.
  3. Model specification. One or more forward models are specified for each dataset.
  4. Model estimation. The model are fitted to the data to determine their evidence and parameters.
  5. Model comparison. The evidence for each model is used for Bayesian Model Comparison to select the best model. Bayesian model averaging is used to compute a weighted average of parameter estimates over different models.
The key stages are briefly reviewed below.

Experimental design

Functional neuroimaging experiments are typically either task-based or examine brain activity at rest. In task-based experiments, brain responses are evoked by known deterministic inputs. These experimental variables can change neural activity through direct influences on specific brain regions, such as evoked potentials in the early visual cortex, or via a modulation of coupling among neural populations; for example, the influence of attention. These two types of input - driving and modulatory - are parameterized separately in DCM. To enable efficient estimation of driving and modulatory effects, a 2x2 factorial experimental design is often used - with one factor serving as the driving input and the other as the modulatory input.
Resting state experiments have no experimental manipulations within the period of the neuroimaging recording. Instead, hypotheses are tested about the coupling of endogenous fluctuations in neuronal activity, or in the differences in connectivity between sessions or subjects. The DCM framework includes models and procedures for analysing resting state data, described in the next section.

Model specification

All models in DCM have the following basic form:
The first equality describes the change in neural activity with respect to time, which cannot be directly observed using non-invasive functional imaging modalities. The evolution of neural activity over time is controlled by a neural function with parameters and experimental inputs. The neural activity in turn causes the timeseries , which are generated via an observation function with parameters. Additive observation noise completes the observation model. Usually, the neural parameters are of key interest, which for example represent connection strengths that may change under different experimental conditions.
Specifying a DCM requires selecting a neural model and observation model and setting appropriate priors over the parameters; e.g. selecting which connections should be switched on or off.

Functional MRI

The neural model in DCM for fMRI is a Taylor approximation that captures the gross causal influences between brain regions and their change due to experimental inputs. This is coupled with a detailed biophysical model of the generation of the BOLD response and the MRI signal, based on the Balloon model of Buxton et al., which was supplemented with a model of neurovascular coupling. Additions to the neural model have included interactions between excitatory and inhibitory neural populations and non-linear influences of neural populations on the coupling between other populations.
DCM for resting state studies was first introduced in Stochastic DCM, which estimates both neural fluctuations and connectivity parameters in the time domain, using Generalized Filtering. A more efficient scheme for resting state data was subsequently introduced which operates in the frequency domain, called DCM for Cross-Spectral Density. Both of these can be applied to large-scale brain networks by constraining the connectivity parameters based on the functional connectivity. Another recent development for resting state analysis is Regression DCM implemented in the Tapas software collection. Regression DCM operates in the frequency domain, but linearizes the model under certain simplifications, such as having a fixed haemodynamic response function. The enables rapid estimation of large-scale brain networks.

EEG / MEG

DCM for EEG and MEG data use more biologically detailed neural models than fMRI, due to the higher temporal resolution of these measurement techniques. These can be classed into physiological models, which recapitulate neural circuity, and phenomenological models, which focus on reproducing particular data features. The physiological models can be further subdivided into two classes. derive from the equivalent circuit representation of the cell membrane developed by Hodgkin and Huxley in the 1950s. Convolution models were introduced by Wilson & Cowan and Freeman in the 1970s and involve a convolution of pre-synaptic input by a synaptic kernel function. Some of the specific models used in DCM are as follows:
Model inversion or estimation is implemented in DCM using variational Bayes under the Laplace assumption. This provides two useful quantities: the log marginal likelihood or model evidence is the probability of observing of the data under a given model. Generally, this cannot be calculated explicitly and is approximated by a quantity called the negative variational free energy, referred to in machine learning as the Evidence Lower Bound. Hypotheses are tested by comparing the evidence for different models based on their free energy, a procedure called Bayesian model comparison.
Model estimation also provides estimates of the parameters, for example connection strengths, which maximise the free energy. Where models differ only in their priors, Bayesian Model Reduction can be used to derive the evidence and parameters of nested or reduced models analytically and efficiently.

Model comparison

Neuroimaging studies typically investigate effects that are conserved at the group level, or which differ between subjects. There are two predominant approaches for group-level analysis: random effects Bayesian Model Selection and Parametric Empirical Bayes. Random Effects BMS posits that subjects differ in terms of which model generated their data - e.g. drawing a random subject from the population, there might be a 25% chance that their brain is structured like model 1 and a 75% chance that it is structured like model 2. The analysis pipeline for the BMS approach procedure follows a series of steps:
  1. Specify and estimate multiple DCMs per subject, where each DCM embodies a hypothesis.
  2. Perform Random Effects BMS to estimate the proportion of subjects whose data were generated by each model
  3. Calculate the average connectivity parameters across models using Bayesian Model Averaging. This average is weighted by the posterior probability for each model, meaning that models with greater probability contribute more to the average than models with lower probability.
Alternatively, Parametric Empirical Bayes can be used, which specifies a hierarchical model over parameters. It eschews the notion of different models at the level of individual subjects, and assumes that people differ in the strength of connections. The PEB approach models distinct sources of variability in connection strengths across subjects using fixed effects and between-subject variability. The PEB procedure is as follows:
  1. Specify a single 'full' DCM per subject, which contains all the parameters of interest.
  2. Specify a Bayesian General Linear Model to model the parameters from all subjects at the group level.
  3. Test hypotheses by comparing the full group-level model to reduced group-level models where certain combinations of connections have been switched off.

    Validation

Developments in DCM have been validated using different approaches:
DCM is a hypothesis-driven approach for investigating the interactions among pre-defined regions of interest. It is not ideally suited for exploratory analyses. Although methods have been implemented for automatically searching over reduced models and for modelling large-scale brain networks, these methods require an explicit specification of model space. In neuroimaging, approaches such as psychophysiological interaction analysis may be more appropriate for exploratory use; especially for discovering key nodes for subsequent DCM analysis.
The variational Bayesian methods used for model estimation in DCM are based on the Laplace assumption, which treats the posterior over parameters as Gaussian. This approximation can fail in the context of highly non-linear models, where local minima may preclude the free energy from serving as a tight bound on log model evidence. Sampling approaches provide the gold standard; however, they are time consuming and have typically been used to validate the variational approximations in DCM.

Software implementations

DCM is implemented in the Statistical Parametric Mapping software package, which serves as the canonical or reference implementation. It has been re-implemented and developed in the Tapas software collection and the VBA toolbox.