随机模型的仿真部分:未感染的目标细胞(U),感染野生型/敏感和未治疗(Is),感染耐药(Ir),野生型病毒(Vs),耐药病毒(Vr)。
simulate_drugresistance_stochastic(
U = 1e+05,
Is = 0,
Ir = 0,
Vs = 10,
Vr = 0,
b = 1e-05,
dI = 1,
e = 0,
m = 1e-04,
p = 100,
c = 4,
f = 0.1,
rngseed = 100,
tfinal = 100
)
U : :目标单元格的初始数目:数字
Is : :野生型感染细胞的初始数目:数字
Ir : :初始耐药感染细胞数:数字
Vs : :野生型病毒的初始数目:数字
Vr : :抗药性病毒的初始数目:数字
b : :细胞感染水平/感染率:数字
dI : :感染细胞死亡率:数字
e : :药物疗效:数字
m : :产生的抗性突变体分数:数字
p : :病毒产生率:数字
c : :病毒清除率:数字
f : :抵抗病毒的适应成本:数字
rngseed : :允许再现性的随机数生成器种子:数字
tfinal : :最大模拟时间:数字
# To run the simulation with default parameters just call the function:
result < - simulate_drugresistance_stochastic()
# To choose parameter values other than the standard one, specify them, like such:
result < - simulate_drugresistance_stochastic(tfinal = 200, e = 0.5)
# You should then use the simulation result returned from the function, like this:
plot(result$ts[,"time"],result$ts[,"Vs"],xlab='Time',ylab='Uninfected cells',type='l')