此函数使用简单的连续时间基本细菌模型执行不确定性和敏感性分析。
simulate_usanalysis(
Bmin = 1,
Bmax = 10,
Imin = 1,
Imax = 10,
Bmaxmin = 1e+05,
Bmaxmax = 1e+06,
dBmin = 0.1,
dBmax = 0.1,
kmin = 1e-07,
kmax = 1e-07,
rmin = 0.001,
rmax = 0.001,
dImin = 1,
dImax = 2,
gmean = 0.5,
gvar = 0.1,
samples = 10,
rngseed = 100,
tstart = 0,
tfinal = 200,
dt = 0.1
)
Bmin : :初始细菌数的下限:数字
Bmax : :初始细菌数的上限:数字
Imin : :初始免疫应答的下限:数字
Imax : :初始免疫反应的上限:数字
Bmaxmin : :最大细菌负荷下限:数字
Bmaxmax : :最大细菌负荷上限:数字
dBmin : :细菌死亡率下限:数字
dBmax : :细菌死亡率上限:数字
kmin : :免疫反应致死率下限:数字
kmax : :免疫反应杀死率的上限:数字
rmin : :免疫反应增长率下限:数字
rmax : :免疫反应增长率的上限:数字
dImin : :免疫反应死亡率下限:数字
dImax : :免疫反应死亡率上限:数字
gmean : :细菌生长率平均值:数字
gvar : :细菌生长率的方差:数字
samples : :要运行的LHS样本数:数字
rngseed : :随机数生成器的种子:数字
tstart : :模拟开始时间:数字
tfinal : :最终模拟时间:数字
dt : :返回结果的次数:数字
# To run the simulation with default parameters just call the function:
## Not run: result < - simulate_usanalysis()
# To choose parameter values other than the standard one, specify them, like such:
result < - simulate_usanalysis(dImin = 0.1, dImax = 10, samples = 5, tfinal = 50)
# You should then use the simulation result returned from the function, like this:
plot(result$dat[,"dI"],result$dat[,"Bpeak"],xlab='values for d',ylab='Peak Bacteria',type='l')