Nothing
`mean.outbreaks.series` <-
function (network, num.reps=10, num.points=10, make.plot=0){
t.crit <- get.t.crit(network=network)
if(make.plot!=0)
curve(get.mean.outbreak(network=network, T=x), from=0, to=c(t.crit-t.crit/10))
df <- data.frame()
for(point in 1:num.points){
print(paste("point",point,"of",num.points))
trans <- (point/num.points)*(9*t.crit/10)
sizes <- list()
for(rep in 1:num.reps){
sim.res <- simulate.epidemic(network=network, initial.infectious=c(1), T=trans)
sizes[[rep]] <- sum(sim.res %v% "status" %in% "recovered")
}
df <- rbind(df, c(trans, mean(as.numeric(sizes))))
if(make.plot!=0)
points(trans, mean(as.numeric(sizes)))
}
df
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.