1

我正在使用 rpy2-2.2.2 和新的免费 Enthought python 发行版,其中包括 numpy 1.6.0 和 python 2.7.2。我 easy_installed rpy2 导致 v. 2.2.2 被安装并且所有测试都成功。

我遇到的问题是我编写的代码在 rpy2 2.1.8 和 python 2.6 上运行良好。问题在于将数组从 numpy 转换为 R。

以下是相关代码的片段:

import rpy2                    
import rpy2.rinterface as rinterface
import rpy2.robjects as rob    
import rpy2.rlike.container as rlc
import numpy as np            
import rpy2.robjects.numpy2ri  
r = rob.r
...
    HGr     = rob.conversion.py2ri(HG_reg)
    RHSr    = rob.conversion.py2ri(RHS)
    #
    CalData = rlc.TaggedList([HGr,RHSr],tags=('hg','rhs'))
    CalData = rob.DataFrame(CalData)

    r('''library(pls)''')
    #rob.globalEnv["HGr"] = HGr
    #rob.globalEnv["RHSr"] = RHSr
    rob.globalenv["CalData"] = CalData


    # perform the PLS regression
    if wetlflag:
        HGresults = r.plsr(r("hg ~ rhs.1 + rhs.2 + rhs.3 + rhs.4"),data=CalData,validation="LOO")

我很乐意承认这不是最优雅的做事方式,但它在我需要提供结果之前和现在都有效(!)。我得到的错误如下:

   Traceback (most recent call last):
  File "Mercury_PLS_WL_DF.py", line 224, in <module>
    HGr     = rob.conversion.py2ri(HG_reg)
  File "/Library/Frameworks/Python.framework/Versions/7.1/lib/python2.7/site-packages/rpy2-2.2.2dev_20110726-py2.7-macosx-10.5-i386.egg/rpy2/robjects/__init__.py", line 134, in default_py2ri
    raise(ValueError("Nothing can be done for the type %s at the moment." %(type(o))))
ValueError: Nothing can be done for the type <type 'numpy.ndarray'> at the moment.

我在这里找到了讨论并得到了 numpy 数组现在自动转换为 R 数组的印象,但是注释掉rob.conversion.py2ri(HG_reg)语句并直接使用 numpy 数组似乎也失败了。我错过了一些明显的东西吗?为什么会在 2.1.8 和 2.2.2 之间中断?

4

1 回答 1

3

来自http://rpy.sourceforge.net/rpy2/doc-2.2/html/numpy.html#from-numpy-to-rpy2

警告

在 rpy2 的早期版本中,只需要导入即可进行转换。导入模块时的副作用可能会导致问题,现在有一个额外的步骤来激活转换:调用函数 rpy2.robjects.activate()。

所以放在rpy2.robjects.activate()导入之后,你应该没问题。

于 2011-07-26T23:26:56.080 回答