Hi Michael,<br><br>I tried running the script on other datasets and they seem to work fine. I've run into this sort of thing several times so far with PyMVPA (i.e., cases where a classifier will run for one subject but crash for others). <br>
<br><blockquote style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;" class="gmail_quote">Did you inspect the dataset that it is somewhat "sane"? What does<br>
print blockavged_dataset.summary()
say?<br></blockquote><br>Yes, the data looks fine. In fact, I've run the same data with several other classifiers and they worked fine. Further, the data will run with the searchlight for several spheres, and only crash when it hits an especially small one (4 voxels). I'm not sure if it's an issue with the number of voxels or not. (I imagine it must've hit similarly small feature sets for the other subjects it worked for.)<br>
<br>Any idea what could be causing this?<br><br>Michael<br><br><br><br><br><br><div class="gmail_quote">On Wed, Dec 16, 2009 at 4:15 PM, Michael Hanke <span dir="ltr"><<a href="mailto:michael.hanke@gmail.com" target="_blank">michael.hanke@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hey,<br>
<br>
[ sorry for the delay ]<br>
<div><br>
On Wed, Dec 16, 2009 at 12:04:28PM +0100, Valentin Haenel wrote:<br>
> Would anyone of the senior developers like to comment on this, please.<br>
<br>
</div>You're senior too -- no? ;-)<br>
<div><div></div><div><br>
> * Michael W. Cole <<a href="mailto:mwcole@gmail.com" target="_blank">mwcole@gmail.com</a>> [091214]:<br>
> > Hello Valentin,<br>
> ><br>
> > Thanks for the reply.<br>
> ><br>
> > I'm including the requested info:<br>
> > PyMVPA 0.4.3 on Fedora 10 x86_64.<br>
> ><br>
> > Here is the relevant part of the script:<br>
> ><br>
> > blockavged_dataset = NiftiDataset(samples=niftiFile,<br>
> > > labels=attrs.labels, chunks=attrs.chunks, labels_map={'cat1': 0,<br>
> > > 'cat2': 1, 'cat3': 2, 'cat4': 3, 'cat5': 4},<br>
> > > mask=maskFile)<br>
> > ><br>
> > > # remove first category from dataset<br>
> > > blockavged_dataset = blockavged_dataset.selectSamples(N.array([l != 0<br>
> > > for l in blockavged_dataset.labels],<br>
> > > dtype='bool'))<br>
> > ><br>
> > > """Create a dissimilarity matrix based on the labels of the data points<br>
> > > in our test dataset. This will allow us to see if there is a<br>
> > > correlation<br>
> > > between any given searchlight sphere and the experimental<br>
> > > conditions."""<br>
> > ><br>
> > > # create dissimilarity matrix using the 'confusion' distance<br>
> > > # metric<br>
> > > dsm = DSMatrix(blockavged_dataset.labels, 'confusion')<br>
> ><br>
> ><br>
> > #<br>
> > > # Run Searchlight<br>
> > > #<br>
> > ><br>
> > > # choose classifier<br>
> > > dsmetric = DSMDatasetMeasure(dsm, 'pearson', 'pearson')<br>
> > ><br>
> > > for radius in [6.4,9.6]:<br>
> > ><br>
> > > # tell which one we are doing<br>
> > > print "Running searchlight with radius: %.2fmm..." % (radius)<br>
> > ><br>
> > > # setup Searchlight with a custom radius<br>
> > > # radius has to be in the same unit as the nifti file's pixdim<br>
> > > # property.<br>
> > > sl = Searchlight(dsmetric, radius=radius)<br>
> > ><br>
> > > # run searchlight on example dataset and retrieve error map<br>
> > > sl_map = sl(blockavged_dataset)<br>
<br>
</div></div>From just looking at it I have no clue what is wrong. Except for the<br>
dataset used this code is almost identical to the searchlight_dsm.py<br>
example. Does that one run for you?<br>
<br>
Did you inspect the dataset that it is somewhat "sane"? What does<br>
<br>
print blockavged_dataset.summary()<br>
<br>
say?<br>
<br>
<br>
Michael<br>
<font color="#888888"><br>
<br>
--<br>
GPG key: 1024D/3144BE0F Michael Hanke<br>
<a href="http://mih.voxindeserto.de" target="_blank">http://mih.voxindeserto.de</a><br>
</font><div><div></div><div><br>
_______________________________________________<br>
Pkg-ExpPsy-PyMVPA mailing list<br>
<a href="mailto:Pkg-ExpPsy-PyMVPA@lists.alioth.debian.org" target="_blank">Pkg-ExpPsy-PyMVPA@lists.alioth.debian.org</a><br>
<a href="http://lists.alioth.debian.org/mailman/listinfo/pkg-exppsy-pymvpa" target="_blank">http://lists.alioth.debian.org/mailman/listinfo/pkg-exppsy-pymvpa</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Michael W. Cole, Ph.D.<br>Post-doctoral researcher<br>Cognitive Control & Psychopathology Lab<br>Washington University in St. Louis<br>