_:n925d078573d042319856b1c7ab13a9b2xb1
<http://schema.org/text>
"I have been trying to replicate the following graph from the paper, https://www.sciencedirect.com/science/article/pii/S2211379722000122 for a class project. \r\n\r\nThe authors draw contour plots of R0 against beta and zeta. The parameter values and the formula for R0 are given below.\r\n\r\n%Basic Reproduction Number\r\n\r\ntheta = 141302; %recruitment rate\r\nmu = 0.001229; %natural death rate\r\ntau = 0.45; %modification factor for A\r\nzeta = 1/14; %influx from Q to S \r\nbeta = 0.88; %transmission coefficient\r\nalpha = 0.75214; %hospitalization rate\r\nq = 0.31167; %influx from Q to I\r\neta_1 = 0.81692; %influx from E to Q\r\neta_2 = 0.02557; %influx from E to A\r\neta_3 = 1/7; %influx from E to I\r\ndelta_1 = 0.16673; %disease death rate for A\r\ndelta_2 = 0.00147; %disease death rate for I\r\ndelta_3 = 0.00038; %disease death rate for J\r\ngamma_1 = 0.00827; %recovery rate for A\r\ngamma_2 = 0.00787; %recovery rate for I\r\ngamma_3 = 0.20186; %recovery rate for J\r\n\r\nK_1 = eta_1 + eta_2 + eta_3 + mu\r\nK_2 = zeta + q + mu\r\nK_3 = gamma_1 + delta_1 + mu\r\nK_4 = alpha + gamma_2 + delta_2 + mu\r\n\r\nR_0 = beta*(tau*eta_2*K_2*K_4 + K_3*(eta_3*K_2 + eta_1*q))/(K_1*K_2*K_3*K_4)\r\n\r\nAny support on drawing the plot would be highly appreciated. Thank you for your time!"
<https://se.mathworks.com/matlabcentral/answers/1881467-plotting-contour-plots-of-r0-against-two-parameters?s_tid=prof_contriblnk>