bayesian - ETAS モデルの stan プログラミング

Question

私はスタンの新人です。私は、地震のモデル化に使用されるモデルである時間 ETAS モデルに取り組んでいます。地震発生時刻 t[i] での震度は、

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 は 3 つのパラメーターです。Rstanを使用しています。モデル用に次のスタンコードを書きました。

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 でこれを行うための書き込み方法ではないと思います。

A small part (containing 20 times only) of the data is as follows: 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)