[pymvpa] Regression estimates?

Per B. Sederberg psederberg at gmail.com
Thu Mar 31 02:58:47 UTC 2011


Hi folks:

Long time no chat  :)

I've been struggling for hours to do something extremely easy (with
latest github master), so I'm breaking down to ask.  I'm running a
cross-validated regression and sensitivity analysis:

clfr = GLMNET_R(alpha=.5, model_type='naive', enable_ca=['estimates'])
sclf = SplitClassifier(clfr,
                       NFoldPartitioner(),
                       enable_ca=['stats','estimates'])

# Compute sensitivity, which internally trains the classifier
print "Running CV"
analyzer = sclf.get_sensitivity_analyzer()
senses = analyzer(dat_all)
print sclf.ca.stats

I'm having trouble understanding the performance in the stats:

-------> print(sclf.ca.stats)
Statistics  Mean  Std   Min      Max
---------- ----- ----- -----    -----
Data:
   RMP_t   2.285 0.696  1.3       4
   STD_t     0     0     0        0
   RMP_p   2.357 0.627 1.494    4.502
   STD_p     0     0     0        0
Results:
    CCe      1     0     1        1
   RMSE    0.392 0.267 0.019    1.117
RMSE/RMP_t  0.18 0.125 0.009    0.477
Summary:
    CCe     0.25         p=  0.000125809
   RMSE     0.47
RMSE/RMP_t  0.2
 # of sets   20


Also, I really just want to get the predictions/estimates of the
classifier from each fold, but I don't know how to get that out.
Every place I look for estimates just gives me:

UnknownStateError: Unknown yet value of estimates


Any pointers would be greatly appreciated as I have about one more
hour to get this working.

Thanks,
Per



More information about the Pkg-ExpPsy-PyMVPA mailing list