[pymvpa] RPy and ENET

[Michael Hanke]

On Mon, Jan 11, 2010 at 09:31:45AM -0500, Michael Hanke wrote:
> On Mon, Jan 11, 2010 at 08:27:15AM -0500, Per B. Sederberg wrote:
> > exists('rpy2', raiseException=True); import rpy2.robjects;
> > rpy2.robjects.r.library('glmnet')
> > 
> > You just ran the last two and it works.  What happens when you run the
> > first?  Perhaps it is something in the exists method...  I assume if
> > you can run the last two, then you can also run:
> > 
> > import rpy2 as __
> > 
> > which is the check for rpy2...
> > 
> > What happens if you change the glmnet check to be just:
> > 
> > import rpy2.robjects; rpy2.robjects.r.library('glmnet')
> 
> The problem was RPy1 -- they didn't play nicely together. Once I removed
> RPy1, things were fine.
> 
> That raises the question, what we should do with the RPy1 dependent
> code? Per, could you port ithe rest to RPy2 as well (possibly removing ENET, if
> it is fully replaced by GLMNET now)?
> 
> 
> Sidenote: If refactored the glmnet tests a bit and included them into
> the main battery. That revealed a few problems that I am addressing
> right now -- will push what I have now nevertheless and push the fixes
> later on.

One more: I get lots of this:

Warning message:
from glmnet fortran code - Convergence for 50th lambda value not reached after maxit=100 iterations; solutions for larger lambdas returned 

Is there anything we can do about that?


michael

-- 
GPG key:  1024D/3144BE0F Michael Hanke
http://mih.voxindeserto.de