[pymvpa] Using ridge regression in a searchlight analysis

VINCE TASCHEREAU-DUMOUCHEL vincenttd at ucla.edu
Tue Jan 17 22:36:26 UTC 2017


Hi everyone,

I am trying to run a searchlight with a linear outcome ranging from 0 to 6
(I recoded the Targets of my dataset accordingly). I have 1 average pattern
for each of these target values for each of my 6 chunks (I should probably
include more trials instead of using averages, but this is what I am
playing with so far to have an even dataset) .

I thought of using ridge regression and here is the code I have so far:

    clfsd = RidgeReg(lm=None, enable_ca=['stats'])

    cv = CrossValidation(clfsd,
                     NFoldPartitioner(),
                     errorfx=lambda x,y: pearsonr(x,y)[0],
                     enable_ca=['training_stats','stats'])

    sl = sphere_searchlight(cv, radius=3, postproc=mean_sample())
    sl_results = sl(ds_alll)

I am having a hard time figuring out what the pearson r means in this
context. Is it the correlation between Y  and Y-Hat? I end up with some
negative values. Any chance it could be interpreted has a negative relation
between the voxels within the searchlight and the continuous target values?

Also, since I am fairly new to pyMVPA, anyone has any suggestion on how to
figure that out by myself in the code? I tried running in debug mode but I
seem to run into an error that I don't have if I run the code normally.

Thanks!

V
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.alioth.debian.org/pipermail/pkg-exppsy-pymvpa/attachments/20170117/2aeba6c8/attachment.html>


More information about the Pkg-ExpPsy-PyMVPA mailing list