My favorites | Sign in
Project Home
READ-ONLY: This project has been archived. For more information see this post.
Search
for
Results  
VBA graphical output
Updated Jun 24, 2014 by jean.dau...@gmail.com

This page describes the graphical output of the VBA model inversion.

In particular, we focus on the example of a Q-learning model inversion, given observed choices (see demo_Qlearning.m or this page).

By default, the VBA model inversion terminates by displaying the results in a MATLAB window using tabs, which are described below.

Note: the VBA graphical summary can be reproduced automatically from the posterior and out output structures as follows:

newFig = 1;
[hfp,out] = VBA_ReDisplay(posterior,out,newFig);

The new out structure now contains the inversion diagnostics, in case these had not been derived during the VBA inversion (cf. options.DisplayWin = 0).

Summary

This tab displays general information about the model inversion, including the log-model evidence for Bayesian model comparison.

NB: the "null" log-evidence is given here as well, along with fit accuracy metrics (percentage of variance explained, log-likelihood, etc...)

VB inversion

Sufficient statistics of the posterior densities are displayed here (errorbars depict posterior standard deviations, grey barplots show the difference between the posterior and the prior means). These include the posterior predictive densities (which are plotted against time and against data).

NB: For binary data (as is the case here), the upper-right panel shows the mean and standard deviation of data samples (y-axis) corresponding to six percentiles of the predicted data (x-axis). Ideally, errorbars should align with the red ellipse (along the main diagonal).

diagnostics

Diagnostics include estimated state- and measurement noise time-series, histograms and autocorrelation, as well as the parameters' posterior correlation matrix, etc…

NB: a high posterior correlation between model two model parameters means that they are non-identifiable (they have the same impact on the data). In our example, only the difference between the initial conditions is identifiable (increasing any of these initial conditions has the same effect on data, hence the positive posterior correlation).

kernels

This tab shows the Volterra decomposition of data and system’s states onto the set of the system's input. Would the system be linear, Volterra kernels would be the impulse response function of the system.

NB: Volterra weights (y-axis) are plotted against time lag (x-axis) for each input. Positive weights mean that the system's output increases with the corresponding (lagged) input.

conv

This is useful for checking the convergence of the VB algorithm. One can eyeball the history of free energy optimization over iterations of the VB algorithm.

NB: The "null" log-evidence (+/-3) is also shown, for comparison purposes. This provides a dummy Bayesian p-value test for the model (against chance).

priors

This tab creates a similar display than ‘VB inversion’, but with the sufficient statistics of the prior density (including the prior predictive density, which is also plotted against the data).

NB: this is useful to get an idea of how much was learned during the model inversion (cf. comparison of prior and posterior predictive densities).

Powered by Google Project Hosting