PythonRPy“x$terms:$operator中的错误对于原子向量无效”

Question

我是RPy的新人，如果我的问题很小，请原谅。我试图用Python编写这个主题中的顶级解决方案：Screening (multi)collinearity in a regression model，但我得到以下错误：

rpy.RPy_RException: Error in x$terms : $ operator is invalid for atomic vectors

我写的代码：

from rpy import *
r.set_seed(42)
a=r.rnorm(100)
b=r.rnorm(100)
m=r.model_matrix('~a+b')

我做错什么了？在

编辑： 使用agstudy撰写的回复（感谢您的帮助！）我为rpy2准备了解决方案

^{pr2}$

Answer 1

这应该行得通

fmla = r.Formula('~a+b')
env = fmla.environment
env['a'] = a
env['b'] = b
r.model_matrix(fmla)

在R中，可以重现错误

^{pr2}$