[pymvpa] Loading multiple ROIs with fmri_dataset

Yaroslav O Halchenko yoh at onerussian.com
Wed Aug 14 16:25:19 BST 2019

On Wed, 14 Aug 2019, Lyam Bailey wrote:

>    Dear Yaroslav,

>    Thanks so much for the quick reply. I tried implementing this with:

>    fa_fname = {a: "haxby2001/sub001/masks/25mm/%s.nii.gz" % a for a in ["vt,
>    hoc"]}
>    temp_ds = fmri_dataset(bold_fname, mask=None, targets = targets,
>    add_fa=fa_fname, chunks = int(i)-1)

>    The code ran without errors, however when I continued through the example
>    RSA script (http://www.pymvpa.org/examples/rsa_fmri.html) the
>    dissimilarity matrices appeared very different to those generated when
>    using load_tutorial_data. I experimented with
>    fa_fname = {a: "haxby2001/sub001/masks/25mm/%s.nii.gz" % a for a in
>    ["white"]}

>    ...and this turned out to generate exactly the same results. I also get
>    the same results with add_fa = None. So I get the feeling that add_fa is
>    not actually being implemented. Can you see a problem in the way I've
>    called fmri_dataset?

add_fa just adds those feature attributes.  You would still need to
subselect those features if you decide to work with specific ROIs. 

for se3lection could do smth like

    roi_ds = temp_ds[:, temp_ds.fa.vt != 0]

or am I missing what you are missing? ;)

Yaroslav O. Halchenko
Center for Open Neuroscience     http://centerforopenneuroscience.org
Dartmouth College, 419 Moore Hall, Hinman Box 6207, Hanover, NH 03755
Phone: +1 (603) 646-9834                       Fax: +1 (603) 646-1419
WWW:   http://www.linkedin.com/in/yarik        

More information about the Pkg-ExpPsy-PyMVPA mailing list