随机模型的仿真部分:未感染的目标细胞(U),感染细胞(I),病毒(V)。
simulate_basicvirus_stochastic(
U = 10000,
I = 0,
V = 1,
n = 0,
dU = 0,
b = 1e-04,
dI = 1,
p = 10,
dV = 2,
rngseed = 111,
tfinal = 30
)
U : :目标单元格的初始数目。必须是整数:数字
I : :野生型感染细胞的初始数目。必须是整数。:数字
V : :抗药性病毒的初始数量。必须是整数。:数字
n : :未感染细胞产生率:数字
dU : :未感染细胞清除率:数字
b : :细胞感染水平/感染率:数字
dI : :感染细胞死亡率:数字
p : :病毒产生率:数字
dV : :病毒清除率:数字
rngseed : :允许再现性的随机数生成器种子:数字
tfinal : :最终模拟时间:数字
# To run the simulation with default parameters just call the function:
result < - simulate_basicvirus_stochastic()
# To choose parameter values other than the standard one, specify them, like such:
result < - simulate_basicvirus_stochastic(U = 1e3, dI = 0.1)
# You should then use the simulation result returned from the function, like this:
plot(result$ts[,"time"],result$ts[,"V"],xlab='Time',ylab='Virus',type='l')