如何使用stans integrate_ode x输入？

最近斯坦已经添加了这种integrate_ode方法.不幸的是,我能找到的唯一文件是stan参考手册(第191页).我有一个需要一些驱动信号的模型.据我了解参数x_r,x_i应该用于此.

为了一个具体的例子,我们假设我想从文档中实现以下更改的示例:

 real[] sho(real t,
               real[] y,
               real[] theta,
               real[] x_r,
               int[] x_i) {
      real dydt[2];
      real input_signal; // Change from here!!!
      input_signal <- how_to(t, x_r, x_i); 
      dydt[1] <- y[2] + input_signal; // Change to here!!!
      dydt[2] <- -y[1] - theta[1] * y[2];
      return dydt;
}

输入信号应该是输入的时间序列 - 假设我提交input_signal_vector <- sin(t) + rnorm(T, sd=0.1)(这应该是时间点的信号ts),我打算用于input_signal最接近的时间值.

我能想象的唯一方法是一个可以Concat的ts和input_signal_vector在x_r,然后这个数组中进行搜索.但我无法想象这是这些参数的预期用途.这也是非常低效的.

因此,如果有人能够证明这种情况应该如何解决,我将非常感激.