How to graph R0 against two parameters on a 3D plot with planes for certain R0 values?

2 次查看（过去 30 天）
Bas123
Bas123 2022-12-20
评论： Bas123 2022-12-20
I have been trying to replicate the following graph from the paper, https://pubmed.ncbi.nlm.nih.gov/32834653/. The authors have plotted R0 against the two parameters gamma_a and beta_s on a 3D plot along with two planes for the values of R0=1 and R0=2. The parameters values and the formula for R0 are given below.
beta_s = 0.274;
alpha_a = 0.4775;
alpha_u = 0.695;
mu = 0.062;
q_i = 0.078;
gamma_a = 0.29;
eta_i = 0.009;
eta_u = 0.05;
R0 = (beta_s*alpha_a)/(gamma_a+mu) + (beta_s*alpha_u*gamma_a*(1-q_i))/((gamma_a+mu)*(eta_u+mu))
I would highly appreciate if someone could help me draw this plot. Thank you so much!

请先登录，再回答此问题。

采纳的回答

Bora Eryilmaz
Bora Eryilmaz 2022-12-20
编辑：Bora Eryilmaz 2022-12-20
Ran in:
[X,Y] = meshgrid(0:0.1:0.4, 0:0.01:0.4);
R0 = 1;
Z = R0 * ones(size(X));
surf(X,Y,Z)
colormap('turbo')
shading interp
xlabel('\beta_s')
ylabel('\gamma_a')
zlabel('R_0')
hold on
R0 = 2;
Z = R0 * ones(size(X));
surf(X,Y,Z)
colormap('turbo')
shading interp
%beta_s = 0.274;
beta_s = X;
alpha_a = 0.4775;
alpha_u = 0.695;
mu = 0.062;
q_i = 0.078;
% gamma_a = 0.29;
gamma_a = Y;
eta_i = 0.009;
eta_u = 0.05;
R0 = (beta_s.*alpha_a)./(gamma_a+mu) + (beta_s.*alpha_u.*gamma_a*(1-q_i))./((gamma_a+mu)*(eta_u+mu));
surf(X,Y,R0)
colormap('turbo')
shading interp

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Discrete Data Plots 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by