用于循环和在R中绘图

该程序是一个简单的流行病模拟器，除了一件事外，大部分都是完成的。我想得到一个图上每个时间步的感染人数（它可以是任何类型的图，最好是线形图）

以下是整个脚本的关键部分：

 for (i in 1:simlength) {
  discordant.links <- which(infected[links[,1]]+infected[links[,2]]==1)
  transmit <- rbinom(length(discordant.links),1,p.t)
  transmitter.links <- discordant.links[transmit==1]
  nodes.of.transmitter.links <-   unique(as.vector(links[transmitter.links,1:2]))
  infected[nodes.of.transmitter.links] <- 1
  infected[infected==1] <- sample(c(1,2),length(infected[infected==1]),replace=TRUE,prob=c(1-p.r,p.r))
  if (plot.spread) {
    node.colour[infected==1] <- "red"
    node.colour[infected==2] <- "green"
    plot(network.i,layout=fixlayout, main=paste("Time =", i), vertex.color=node.colour)
 }
}

for（1中的i:simlength）{
discordant.links如果您的示例受到影响肯定会有帮助。
n这可能会受到影响，但是我需要得到一个图。我应该给出一个更好的示例…我想得到一个图，显示每个时间步有多少人受到感染。例如，如果simlength==5，那么我可以得到一个显示受感染pe总数的图在每一个时间步都有人感染，在这个例子中是：1，2，3，4和5。就像在时间步1有10人感染，在时间步2有25人感染等等。对不起，我是一个不专业的人，但我在这方面还是比较新的。