我正在使用igraph模拟一个无量纲的SIR模型。
当beta或gamma是0到1之间的分数时,我遇到了内部错误。这是不是在iGraph中允许的?
例如:
g <- erdos.renyi.game(1000, 1/1000)
sir(g,0.4,1)
# Error in .Call("R_igraph_sir", graph, beta, gamma, no.sim, PACKAGE = "igraph") :
# At sir.c:185 : Internal SIR error, Internal error, likely a bug in igraph
对于我的模型,我需要使用以下参数:
b <- 0.0006393279 # beta
gm <- 1/7 # gamma
有没有办法在sir()
上使用上述参数?