Руководство по использованию выражения rpy2 eval - PullRequest
Новая
       6

Руководство по использованию выражения rpy2 eval

3 голосов
/ 12 августа 2011

Я новичок как в R, так и в rpy2. Я пытаюсь портировать следующий пример 

library(MASS)  # for eqscplot  
data(topo, package="MASS")  
topo.kr <- surf.ls(2, topo)  
trsurf <- trmat(topo.kr, 0, 6.5, 0, 6.5, 50)

до rpy2.

Пока у меня есть 

import  rpy2.robjects as robjects
robjects.r('library(spatial)')  
f1 = robjects.r['surf.ls']  
x = robjects.IntVector([1,2,3])  
y = robjects.IntVector([1,2,3])  
z = robjects.IntVector([1,30,3])  
res = f1(2, x,y,z)

Я предполагаю, что результат должен быть res. Однако, когда я печатаю res, используя print(res.r_repr()), я получаю выражение, которое не могу оценить. Буду признателен за любую помощь в этом.

Ответы [ 2 ]

2 голосов
/ 16 августа 2011

Тогда возникает вопрос: как мне преобразовать список R в словарь Python (мало общего с оценкой выражений).С rpy2-2.2.x (и 2.3-dev): 

from rpy2.robjects.vectors import ListVector
# make an R list
l = ListVector({'a': 1, 'b': 'b'})

# get a Python dict out of an R list
dict(l.iteritems())
1 голос
/ 15 августа 2011

Ваш код работает нормально.Я думаю, что у вас просто проблемы с доступом к результатам.Ваш результирующий объект res по сути является списком R.Я бы преобразовал его в соответствующий словарь Python. 

rListObj = {}
for key,val in zip(robjects.r.names(res),res):
  rListObj[key] = [i for i in val] #R Vector to List

Результаты: 

{'f': [1.0, 1.0, 1.0, -1.0, 0.0, 1.0, 1.0, 0.0, 1.0, -1.0, 0.0, 1.0, 1.0, 0.0, 1.0, 1.0, 0.0, 1.0], 'rx': [1, 3], 'ry': [1, 3], 'np': [2], 'beta': [0.0, 0.0, 0.0, 0.0, 0.0, 0.0], 'r': [-1.7320508075688772, -1.6729239521451333e-16, -1.4142135623730951, -1.1547005383792512, -5.187907395343139e-17, -0.8164965809277259, -1.6729239521451333e-16, -1.4142135623730951, 3.415236843329339e-17, nan, -1.1547005383792512, -5.187907395343139e-17, -0.8164965809277259, nan, 0.0, -1.1547005383792512, -5.187907395343139e-17, -0.8164965809277259, nan, 0.0, 0.0], 'call': [<SignatureTranslatedFunction - Python:0xb7539dec / R:0xa686cec>, <IntVector - Python:0xb7534cac / R:0xa69e788>, <IntVector - Python:0xb7534d2c / R:0xa5f72f8>, <IntVector - Python:0xb7534c2c / R:0xa5f7320>, <IntVector - Python:0xb7534bac / R:0xa5f7348>], 'y': [1, 2, 3], 'x': [1, 2, 3], 'z': [1, 30, 3], 'wz': [0.0, 0.0, 0.0]}

Я проверил это на несколько старой версии rpy2 (2.1.9), вероятноУдивительные способы сделать это с более свежими версиями.

...