dataset used this code is almost identical to the searchlight_dsm.py example. Does that one run for you? Did you inspect the dataset that it is somewhat "sane"? What does print blockavged_dataset.summary() say? Michael -- GPG key: 1024D/3144BE0F Michael Hanke http://mih.voxindeserto.de