如何在 rpy2 中将列表传递给 R 并返回结果

2024-04-16

我第一次尝试使用 rpy2。假设我在 python 中有一个列表

l = [1,2,3,4,5,6]

我想用 R 打电话

ks.test(l, pexp)

我怎样才能做到这一点?

我最初的尝试是

#!/usr/bin/python
import rpy2.robjects as robjects

l = [1,2,3,4,5]

f = robjects.r('''
f<-function(l){
    ks.test(l, pexp)
}''')
print f(l)

这显然不是我所得到的正确方法

Error in sort.int(x, na.last = na.last, decreasing = decreasing, ...) : 
  'x' must be atomic
Traceback (most recent call last):
  File "./rpy2-test.py", line 12, in <module>
    print f(l)
  File "/usr/local/lib/python2.7/dist-packages/rpy2/robjects/functions.py", line 166, in __call__
    return super(SignatureTranslatedFunction, self).__call__(*args, **kwargs)
  File "/usr/local/lib/python2.7/dist-packages/rpy2/robjects/functions.py", line 99, in __call__
    res = super(Function, self).__call__(*new_args, **new_kwargs)
rpy2.rinterface.RRuntimeError: Error in sort.int(x, na.last = na.last, decreasing = decreasing, ...) : 
  'x' must be atomic

这样做的正确方法是什么?


我认为你不能直接将 python 列表传递给 R。您可以在使用之前将其转换为 R int 向量。

#!/usr/bin/python
import rpy2.robjects as robjects

l = [1,2,3,4,5]

# get ks.test via execute string as R statement
test = robjects.r('ks.test')
# get a built-in functions variables directly
pexp = robjects.r.pexp

l_vector = robjects.IntVector(l)
result = test(l_vector, pexp)
print result[result.names.index('p.value')]

参考:

  • 使用 R 的 OOP http://rpy.sourceforge.net/rpy2/doc-2.4/html/robjects_oop.html
  • 功能 http://rpy.sourceforge.net/rpy2/doc-2.4/html/robjects_functions.html
  • 向量和数组 http://rpy.sourceforge.net/rpy2/doc-2.2/html/vector.html#assigning-python-style
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

如何在 rpy2 中将列表传递给 R 并返回结果 的相关文章

随机推荐