[pymvpa] RPy and ENET
michael.hanke at gmail.com
Mon Jan 11 14:33:04 UTC 2010
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:
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?
GPG key: 1024D/3144BE0F Michael Hanke
More information about the Pkg-ExpPsy-PyMVPA