考慮使用rpy2的SignatureTranslatedAnonymousPackage (STAP)在Python環境中導入任意R代碼作爲可用包。爲了演示,寫入中的R到Python的Rdtq github使用rpy2
以下轉換:
ř
# Loading required package: Rdtq
require(Rdtq)
# Assigning drift and diff functions
mydrift = function(x) { -x }
mydiff = function(x) { rep(1,length(x)) }
# Running rdtq()
test = rdtq(h=0.1, k=0.01, bigm=250, init=0, fT=1,
drift=mydrift, diffusion=mydiff, method="sparse")
# Plotting output
plot(test$xvec, test$pdf, type='l')
的Python
from rpy2 import robjects
from rpy2.robjects.packages import STAP
from rpy2.robjects.packages import importr
# Loading required package: Rdtq
Rdtq = importr('Rdtq')
fct_string = """
my_drift <- function(x) { -x }
my_diff <- function(x) { rep(1,length(x)) }
"""
# Creating package with above drift and diff methods
my_fcts = STAP(fct_string, "my_fcts")
# Running rdtq() --notice per Python's model: all methods are period qualified
test = Rdtq.rdtq(h=0.1, k=0.01, bigm=250, init=0, fT=1,
drift=my_fcts.my_drift(), diffusion=my_fcts.my_diff(), method="sparse")
# Load plot function
plot = robjects.r.plot
# Plotting by name index
plot(test[test.names.index('xvec')], test[test.names.index('pdf')], type='l')
這是一個[XY問題](https://開頭meta.stackexchange.com/questions/66377/what-is-the-xy-problem)。您只告訴我們您的Y解決方案,但不解釋X問題。請給出真實的,完整的方案,並提供具體或可重複的例子。雖然你可能會這樣想,但你可能需要'lambda'。 – Parfait
夠公平的,我已經調整了這個問題。 – user56643