R 未知变量的积分

我想在R上编写如下代码：

 \int_{t_{r-1}}^y  \left\{ \int_{t_{r-1}}^{\tau} w(z+x) q(\tau-x\mid z+x)dx\right\} \times  g_i(y-\tau \mid z+\tau) d\tau

所以我试着这样做：

library(cubature)
f <- function(x){ W(Z+x[1])*h(x[2]-x[1],Z+x[1])*g_i(y-x[2],Z+x[2])}
adaptIntegrate(f, lowerLimit=c(t_r-1,t_r-1), upperLimit=c(x[2],y), y=y, Z=Z)

我怎么办

---

谢谢

你的例子是不可复制的。什么是t_r*？我认为您需要将x定义为一个值，并且不允许上限失效。极限中的x不能与函数定义中的任意积分x相同。我不确定，但我担心你试图对τtau进行解析积分，然后是数值积分，这是这个软件包无法做到的。我不能定义x，因为我不知道它是什么。我只想在x上积分，然后在tau上积分，但我很惊讶，在R上不可能这样做，这不是t_R，这是t_R-1，你忽略了一点：你不能对一个伪变量tau进行数值积分。这在任何语言中都是正确的，不仅仅是R。你需要一个计算机代数系统来计算第一个积分，然后你可以做一个数值积分，假设你有y和tr-1的值

object 'x' not found.