[pymvpa] regression analysis to predict subject-specific score

Zhen Zonglei zonglei.fsl at gmail.com
Fri Aug 12 02:27:46 UTC 2011


Thanks very much. I will try them in the next days.


Zonglei Zhen

On Fri, Aug 12, 2011 at 12:00 AM, Yaroslav Halchenko
<debian at onerussian.com>wrote:

> indeed we have no good tutorial/example for regressions yet besides one
> for GPR
>
> doc/examples/gpr.py
>
> also we interface to SVR regressions (from libsvm and shogun) and started
> to add interfaces to scikits.learn.  Some samples of them are available
> from
> regrswh, so  on my system:
>
> *In [5]: print '\n'.join([str(x) for x in regrswh[:]])
> <libsvm epsilon-SVR>
> <libsvm nu-SVR>
> <sg.LinSVMR()/libsvr>
> <skl.PLSRegression_1d()>
> <skl.LARS()>
> <skl.LassoLARS()>
> <GPR(kernel='linear')>
> <GPR(kernel='sqexp')>
>
> As for "howto" -- just the same way you use classifiers -- then
> ConfusionMatrix
> in .stats would be replaced with RegressionStatistics. e.g.
>
> In [5]: from mvpa.suite import *
>
> In [6]: from mvpa.testing.datasets import datasets as testing_datasets
>
> In [7]: cve = CrossValidation(regrswh[:][0], NFoldPartitioner(),
> postproc=mean_sample(), errorfx=corr_error, enable_ca=['training_stats',
> 'stats'])
>
> In [8]: print cve(testing_datasets['chirp_linear'])
> <Dataset: 1x1 at float64, <sa: cvfolds>>
>
> In [9]: print cve.ca.stats
> Statistics  Mean  Std   Min       Max
> ---------- ----- ----- -----     -----
> Data:
>   RMP_t   0.668 0.015 0.639     0.681
>   STD_t   0.661 0.015 0.631     0.675
>   RMP_p   0.644 0.043 0.593     0.731
>   STD_p   0.637 0.042 0.583     0.721
> Results:
>    CCe     0.06 0.016 0.036     0.084
>   RMSE    0.232 0.027 0.184     0.266
> RMSE/RMP_t 0.348 0.043  0.27      0.4
> Summary:
>    CCe     0.06         p=  3.65268e-137
>   RMSE     0.23
> RMSE/RMP_t  0.35
>  # of sets   6
>
>
> On Mon, 08 Aug 2011, Zhen Zonglei wrote:
>
> >    Hi Guys:
>
>
>
> >    I am trying to do multivariate regression analysis to predict
> >    subject-specific score with pymvpa 0.6. But, I did not find some
> >    examples about this in the manual. What regressions are implemented in
> >    the toolbox? Could you please show me how to do regression analysis in
> >    the toolbox?
>
>
> >    Best
>
>
> >    Zonglei Zhen
>
> > _______________________________________________
> > 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
>
> _______________________________________________
> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.alioth.debian.org/pipermail/pkg-exppsy-pymvpa/attachments/20110812/a53af349/attachment.html>


More information about the Pkg-ExpPsy-PyMVPA mailing list