[pymvpa] normalization: zscore by example?

Yaroslav Halchenko debian at onerussian.com
Mon Nov 7 22:29:40 UTC 2011


nothing is impossible... but why would you like to do that? quick
situation/answers pairs

* I care only about classification performance and don't do any feature
  (voxel) diagnosticity analysis such as searchlight

  well -- then zscoring is legit... but I don't think we have builtin
  mapper/function for it, since if you don't have any groups of voxels
  to care about zscoring independently, then it is as simple to in
  pure python as

  ds.samples -= np.mean(ds, axis=1)[:, None]
  ds.samples /= np.std(ds, axis=1)[:, None]

* I do care about per-voxel estimates of diagnosticity

  DO NOT ZSCORE ACROSS FEATURES (VOXELS)! ;)

On Mon, 07 Nov 2011, Mike E. Klein wrote:

>    Hi all,
>    I realize now (after a boneheaded mistake), that PyMVPA's zscore()
>    function is normalizing for all the values across a chunk (looking at each
>    voxel individually). I'm wondering if it's possible to, instead, find the
>    zscores for all the voxels in a single volume (against each voxel's
>    resting state), considering each 3d volume separately.�
>    Thanks,
>    Mike

> _______________________________________________
> Pkg-ExpPsy-PyMVPA mailing list
> Pkg-ExpPsy-PyMVPA at lists.alioth.debian.org
> http://lists.alioth.debian.org/cgi-bin/mailman/listinfo/pkg-exppsy-pymvpa


-- 
=------------------------------------------------------------------=
Keep in touch                                     www.onerussian.com
Yaroslav Halchenko                 www.ohloh.net/accounts/yarikoptic



More information about the Pkg-ExpPsy-PyMVPA mailing list