Split-half correlation-based MVPA with group analysis¶
[HGF+01]: first MVPA correlation paper showing evidence for distributed object category representations in ventral temporal cortex.
[PWD06]: illustration of correlation approach to dissociate spatially overlappig but distinct neural populations.
This exercise the most simple and most popular type of MVPA, that was first reported in [HGF+01]. The intuition is to split the data in two halves (e.g. odd and even runs) and estimate the response for each category, voxel and half seperately. If an area represents categories in a distributed manner (over voxels), then correlations (over voxels) of matching categories may be higher than correlations of non-matching categories. If categories are not represented in such a way, one would expect no differences between correlations of matching and non-matching categories.
Briefly, the procedure can be visualized as follows:
In this exercise data we use the handy functionality that cosmo fmri dataset provided us with to compute split-half correlations. Note that this exercise requires you to compute these correlations ‘manually’; another exercise shows how to use a CoSMoMVPA measure (cosmo correlation measure) that computes these values more conveniently.
Part 1 (single subject analysis)¶
Using cosmo fmri dataset, load the t stats for ‘odd’ and ‘even’ runs for s01 (
glm_T_stats_even.nii), while supplying the VT mask. As a reminder, the stimulus labels for each run of the fMRI study were monkey, lemur, mallard, warbler, ladybug, and lunamoth – in that order. These two nifti files contain summary statistics (T statistics from the general linear model
analysis, GLM) for each stimulus for odd and even runs, respectively.
Apply the VT mask to each half, then compute all pairwise correlations between patterns in the first and the second half, resulting in a 6x6 matrix. After applying a Fisher transform, compute the mean difference between values on the diagonal and those off the diagonal. If there is no category information one would expect a difference of zero.
Part 2 (group study)¶
As in the previous exercise, compute the difference between the average on-diagonal versus on the odd-diagonal Fisher-transformed correlations. This time, do this for each subject and for the three ROIs (
brain). For each ROI, gather the correlation difference for each subject and store these, then perform a one-sample t-test against zero to obtain a p-value against the null hypothesis of no correlation information in each region. Use a nested for loop (over ROIs in the outer loop, and subjects in the inner loop) when performing these computations.
Advanced exercise: plot an image of the 6x6 correlation matrix averaged over participants.