我正在嘗試創建一個內部帶有 expression() 的 function，但出現錯誤

Question

我想創建這個 function

library(sde)

myfunction <- function(r,a,K,vi){
  set.seed(123)
  d <- expression(r*x*(K-x))
  s <- expression(a*x*(K-x))
  sde.sim(X0=1,drift=d, sigma=s,M=3) -> X
  plot(X, main="Multiple")
}

但是當我運行 function 時：

myfunction(r=0.5,a=0.35,K=100,vi=1)

我得到

Error in eval(drift) : object 'r' not found

為什么R找不到expression()里面的對象？

我在 function 之外得到了想要的結果並分配了值； 但我想要一個 function 來做到這一點。 問題是 R 在 function 的 expression() 中找不到值。

r <- 0.5  
a <- 0.35  
K <- 20  
vi <- 1  
set.seed(123)  
d <- expression(r*x*(K-x))
s <- expression(ax(K-x))
sde.sim(X0=1,drift=d, sigma=s,M=3) -> X
plot(X, main="Multiple")`

Answer 1

一種解決方法是將bquote與.()一起使用。 bquote評估包含在.()中的表達式：

library(sde)
myfunction <- function(r,a,K,vi){
  set.seed(123)
  d <- as.expression(bquote(.(r)*x*(.(K)-x)))
  s <- as.expression(bquote(.(a)*x*(.(K)-x)))
  sde.sim(X0=1,drift=d, sigma=s,M=3) -> X
  plot(X, main="Multiple")
}

myfunction(r=0.5,a=0.35,K=100,vi=1)