1

R用のパッケージを使用しGillespieSSAており、10回実行しています。コードは次のとおりです。

library(GillespieSSA)
init <- c(S=199, I=1)
a <- c("0.03*{S}*{I}","1*{I}")
nu <- matrix(c(-1, +1, +1, -1),nrow=2,byrow=T)
out <- lapply(X=1:10,FUN=function(x) ssa(init, a, nu, tf=10)$data)
plot(out[[3]][,3],type="l",xlab="Time",ylab="I(t)")

私の質問は、特定の実行をプロットできることです。たとえば、3 回目の実行で out[[3]] はデータの 3 回目の実行を示しますが、すべての実行をプロットして各実行の頻度を知りたいと考えています。コードの書き方がわかりません。私は何週間も立ち往生しており、Rプログラムが非常に苦手です。

4

1 に答える 1

0

すべての実行の頻度をプロットしたいだけなら、これはあなたが探しているものですか?

par(mfrow=c(4,3))

for (i in seq_along(out)){
  hist(out[[i]][,3],xlab="I(t)", main=paste("Run ", i))
}
于 2013-04-24T11:37:19.907 に答える