疫学の確率的シミュレーションをシミュレートしています。離散時間でシミュレートするにはどうすればよいですか? 以下のコーディングを使用して、連続時間を取得することができました。
library(GillespieSSA)
parms <- c(beta=0.591,sigma=1/8,gamma=1/7)
x0 <- c(S=50,E=0,I=1,R=0)
a <- c("beta*S*I","sigma*E","gamma*I")
nu <- matrix(c(-1,0,0, 1,-1,0, 0,1,-1, 0,0,1),nrow=4,byrow=TRUE)
set.seed(12345)
out <- lapply(X=1:10,FUN=function(x) ssa(x0,a,nu,parms,tf=50)$data)
out
離散時間を取得するには、コーディングをどのように変更すればよいですか? 事前に感謝します。