r - 如何在python中使用rpy包调用ltm函数

Question

我正在尝试以下代码：

from rpy import *
r.library("ltm")

dat= #some data frame or matrix
r.ltm(r('dat~z1'))

即将出现的错误是--- RPy_RException：eval（expr，envir，enclos）中的错误：找不到对象'dat'

请告诉我使用 rpy 库调用ltm函数的正确方法

Answer 1

我会尝试使用 rpy2 和类似以下内容的通用方法：

from rpy2.robjects import *
r("library('ltm')")
r.assign('r_var_name',py_var_name)
r("r_var_name<-as.desired.data.type(r_var_name)")

然后在进一步的 r("blah") 语句中使用 'ltm' 包函数对 'r_var_name' 执行任何命令。

例如，获取 ltm 包示例之一的系数：

In [30]: py_obj = r("coef(ltm(Abortion ~ z1, control = list(GHk = 20, iter.em = 20)))")
In [32]: py_obj
Out[32]: 
<Matrix - Python:0x4db0290 / R:0x52f04f0>
[0.188998, -0.256378, -0.367623, ..., 4.542567, 5.840821, 3.243826]