r - R：GillespieSSA 绘制所有运行图

Question

我正在使用GillespieSSAR 包，并且运行了 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)")

我的问题是我可以绘制一个特定的运行，比如说第三次运行，out[[3]] 给出了数据的第三次运行，但是我希望绘制所有运行以了解每次运行的频率。我不知道如何编写代码。我被困了几个星期，而且我的 R 程序非常薄弱。

Answer 1

如果您只想绘制所有运行的频率，这是您正在寻找的吗？

par(mfrow=c(4,3))

for (i in seq_along(out)){
  hist(out[[i]][,3],xlab="I(t)", main=paste("Run ", i))
}