Я разработал модель эпидемии и получил результат. Теперь я хочу провести анализ чувствительности, когда я изменяю один параметр (например, переменную бета), а затем исследую его влияние на сгенерированные результаты, сохраняя остальные параметры постоянными (без изменений).
Параметр, который я хотел бы изменить, - это бета. Я хочу отобразить сгенерированные выходные данные моей модели с различными значениями бета: «среднее (бета = 8,4; центр)», «низкое (бета = 5; слева)», «высокое (бета = 15; высокое)», но полосы в моем графике торнадо (справа), который я создал с помощью этих кодов.
Я хотел бы создать сюжет торнадо, подобный этому:
http://protectbenefitrisk.eu/tornadoplot.html
dat <- matrix(c(0.0020, 0.0106,0.0007,0.1196,0.0422,0.2087,0.8730, 0.7015,
0.6460, 0.01691, 0.0081, 0.0246, 0,0,0,0,0,0,0.2776, .7124,0.1204), ncol=7)
dat <- matrix(c(myRes[4:24]), ncol = 7) # you don’t need to run this, if you are using the above data
rownames(dat) <- c('beta=8.4','beta=5', 'beta=15') # Amount of change in beta
colnames(dat) <- c('S', 'LA', 'LB','I', 'T', 'TL', 'Sexp') # Names of output variables
x <- seq(0,1, 0.01) # For plotting the x-axis
#This the output
S LA LB I T TL Sexp
beta=8.4 0.0020776591 0.11966384 0.8730181 0.016916187 0 0 0.2776309
beta=5 0.0106789409 0.04224487 0.7015975 0.008136292 0 0 0.7124848
beta=15 0.0007922731 0.20872385 0.6460100 0.024695246 0 0 0.1204128
# codes for tornado plot
par(mfrow=c(1,1))
barplot(dat[2,], horiz = T, las=1, xlim = c(0,1), xaxt='n', ylab = '', beside=T, col=c('blue')) # low=5
barplot(dat[3,], horiz = T, las=1, xlim = c(0,1), xaxt='n', ylab = '', beside=T, col=c('red'), add = TRUE) # high=15
grid(5, NA, col="lightgray", lwd=2)
axis(1, at=x, labels=x, col.axis="black", las=1)
Я хочу, чтобы сгенерированные выходные данные для бета = 8,4 были центром, а также сгенерированные выходные данные для бета = 5 с левой стороны и сгенерированные выходные данные для бета = 15 с правой стороны графика торнадо. Это подытожит результаты моего анализа чувствительности.
Но я изо всех сил пытаюсь сгенерировать график торнадо для различных значений бета, чтобы показать эффект бета на выходах.
Любое предложение будет с благодарностью!