重复使用其输出作为自己的输入的 R 函数

Question

给定 R 函数 auto（如下），我想知道是否有可能从第二次运行 for 循环输出 ps 被用作 pr？

例如，如果循环要运行 3 次（即，length(n) = 3），那么在循环的第一次运行中，pr 将按原样使用，但从第二次运行开始，ps 的结果（即，ps(x)) 从之前的运行中取代pr 直到length(n) 的数量。

因此，在for 循环的第一次运行之后，每次，上一次运行中的ps(x) 都会在下一次运行中扮演pr 的角色。我的最终目标是curve最终以这种方式获得的ps。

auto <- function(n, dat){

for(i in 1:length(n)){           
 pr = function(x) dbeta(x, 1, 1)
 lk = function(x) dbinom(dat[i], n[i], x) # So, here first `n = 100` and 
 ps = function(x) pr(x)*lk(x)             # `dat = 55` will go thru first 
  }                                       # round of the loop and produce
 curve(ps)                                # a `ps`. But in the second run of
}                                         # the loop, the `ps` just
# Example of use:                         # produced will be used as `pr`
auto(n = c(100, 50), dat = c(55, 60) )    # for `n = 50` and `dat = 60`
                                          # to produce a new `ps`.

Answer 1

让我试着给出一个答案。

所以，首先，我认为您需要有限状态自动机函数。这种方法需要使用环境使用，您需要在其中存储函数调用的状态和值。 但是，您只需要根据您的周期指数值进行选择。 这要容易得多，因为您只需要 switch 运算符，仅此而已。

auto <- function(n, dat, n.loop){
x <- 0.1;
for(i in 1:n.loop){
val <- switch(i, dbeta(x, 1, 1), dbinom(dat, n, x), dbeta(x, 1, 1)*dbinom(dat, n, x))
print(val);
  }
}

auto(n = 100, dat = 55, n.loop = 3)

我有 2 个问题应该很好地重新组织代码。 1) 您在汽车中使用的x 值是什么？这不是您的代码中的参数，因此，您应该直接确定auto 中x 的值。也许这是一个论点？ 2) 你的函数需要返回什么值？

因此，添加了x <- 0.1 和print(val) 以检查代码执行情况。

Answer 2

我想用不同的方式。虽然，我还没有完成它，但仍然想与他人分享改进的想法。我正在寻找非常基本和粗略的方法来解决它。

auto <- function(n, dat){
  pr = function(x) dbeta(x, 1, 1)
  lk = function(x) dbinom(dat[1], n[1], x)
  ps = function(x) pr(x)*lk(x) 
  #keep the funtion till now in psold 
  psold = ps
  #loop through 2nd to 2nd last element
  for(i in 2:length(n)-1){
    lk = function(x) dbinom(dat[i], n[i], x)
    #Use psold to evaluate new function ps
    ps = function(x) psold(x)*lk(x)
    psold = ps
  }
  lk = function(x) dbinom(dat[length(n)], n[length(n)], x)
  #Finaly function to be returned.
  ps = function(x) psold(x)*lk(x)  
}
# Example of use:
auto(n = c(100, 50), dat = c(55, 60) )

Answer 3

不确定我是否完全理解这个问题。我在 for 循环中添加了一个 if 语句来检查它是否是第一次迭代。如果是，它将通过您的原始语句定义pr。如果不是，则通过 pr(x) * 上一次迭代的 lk(x) 函数定义它，该函数由i-1 子集。我确定我在这里遗漏了一些东西..

我也很困惑 x 输入的来源。

auto <- function(n, dat){

  for(i in 1:length(n)){
    if(i == 1){pr = function(x) dbeta(x, 1, 1)} else{pr = function(x) dbeta(x, 1, 1) * function(x) dbinom(dat[i-1], n[i-1], x)}
    lk = function(x) dbinom(dat[i], n[i], x)
    ps = function(x) pr(x)*lk(x)
  }
  curve(ps)
}
# Example of use:
auto(n = c(100, 50), dat = c(55, 60) )