r - RのGillespies SSAアルゴリズム

Question

RのGillespie SSAの構文を知っている人はいますか? 私が見つけたコードは次のとおりです。

input n=#iterations,r=#reactions,h=step size (dt)
input x(0) =initial conditions,v=state shift matrix
t←0 
for i= 1 to n
   for j= 1 to r
       p←uniform random number in [0,1]
       If
       p < a(x)h then x(t+h)←x(t)+vj
   end
   t←t+h
end

SSA のこのレイアウトまたは新しいシンプルなレイアウトについてサポートが必要です