ETAS 模型的 stan 编程

Question

我是 STAN 的新手。我正在研究时间 ETAS 模型，该模型用于在地震发生时间 t[i] 建模 earthquakes.The 强度，建模为 -

h(t[i]|p,c,mu)=mu+sum((p-1)*(c^(p-1))*(1/((t[i]-t[1:(i-1)]+c)^(p-1))));

其中 t 是时间，p、c、mu 是三个参数。我正在使用 Rstan。我为模型编写了以下 stan 代码：

stan_etas="
data{
  int<lower=0> N;
  real<lower=0> t;
}
parameters{
  real<lower=0> mu;
  real<lower=1.005> p;
  real<lower=0> c;
}

我知道我没有将时间指定为向量。你能帮我把可能性写在模型部分吗？我在写强度时遇到问题。我认为我以前在 R 中写入时间 t[i] 的强度的方式不是在 STAN 中执行此操作的写入方式。

一小部分（只包含20次）数据如下： dat=list(0.0000,310.1907,948.4677,1007.2617,1029.7996,1065.7343,1199.8650, 1234.6809,1298.0234,1316.0350,1381.8400,1413.4311,1546.2059,1591.1326, 1669.5084,1738.9363,1745.5503,1797.9980,1895.6705,1936.3146)

Answer 1

pow函数目前不对向量或数组进行运算，因此您必须循环构建强度。此外，我认为您打算将 t 声明为长度为 N 的实际数组，它看起来像 real<lower=0> t[N];。然后在模型块中，你会有类似的东西：

y[1] <- pow(c, -(p-1));
for (j in 2:N) {
  y[j] <- mu;
  for (i in 1:(j-1))
    y[j] <- y[j] + (p - 1) * c^(p-1) * 
            1 / (t[j]-t[i]+c)^(p-1);
 }

但是，您最终必须使用 increment_log_prob() 函数来注册对数似然。虽然我不熟悉 ETAS 模型，但 ETAS R package 的文档声称它涉及一个积分，目前无法在 Stan 中对其进行数值近似。

ETAS 模型的 stan 编程

stan programming for ETAS model

bayesian

stan

rstan