我目前用 Rpy2 来拟合一个分布,代码是这样的
%load_ext rpy2.ipython
%R require(movMF)
%%R -i dir_data,n_vM_dir -o theta,alpha
result = movMF(dir_data, n_vM_dir, nruns = 10)
theta = result$theta
alpha = result$alpha
dir_data
,n_vM_dir
theta
,alpha
现在想把这个改写成 Python 函数,主要问题是
dir_data
,n_vM_dir
这两个 Python 中的变量,传到 R 里面去?文档里面的做法
from rpy2.robjects import FloatVector
ctl = FloatVector([4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14])
感觉和我的做法实在是不一样,而且复杂很多。上面用 R magic 一下就解决了 写成函数不止于复杂这么多吧。。。