0

我需要将我的数据拟合到 Beta 分布中并检索alpha参数。

我一直在用 Python 编写代码，但 SciPy 中似乎没有任何 beta 拟合函数。要么我在我不太熟悉的 Matlab 中做所有事情，要么在带有 R 及其fitdistr功能的 Python 中做所有事情。所以我选择了后者。

from rpy2.robjects.packages import importr

MASS = importr('MASS')

然后我在 [0,1) 范围内获取我的浮点数向量，并将其传递给fitdistr：

myVector = myVector.tolist()
MASS.fitdistr(myVector,"beta")

太糟糕了，它需要某种其他向量。rpy 和 rpy2 不是应该为我做所有的转换吗？

Error in function (x, densfun, start, ...)  : 
  'x' must be a non-empty numeric vector
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "/usr/lib/python2.7/dist-packages/rpy2/robjects/functions.py", line 82, in __call__
    return super(SignatureTranslatedFunction, self).__call__(*args, **kwargs)
  File "/usr/lib/python2.7/dist-packages/rpy2/robjects/functions.py", line 34, in __call__
    res = super(Function, self).__call__(*new_args, **new_kwargs)
rpy2.rinterface.RRuntimeError: Error in function (x, densfun, start, ...)  : 
  'x' must be a non-empty numeric vector

我需要在这里做什么？

4

1 回答 1

1

答案就在这里： Converting python objects for rpy2

import rpy2.robjects.numpy2ri
rpy2.robjects.numpy2ri.activate()

请参阅http://rpy.sourceforge.net/rpy2/doc/html/numpy.html：

于 2012-07-18T14:05:58.113 回答