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

조회 수: 3 (최근 30일)
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 2022년 12월 20일
편집: Bora Eryilmaz 님. 2022년 12월 20일
[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')
xlabel('\beta_s')
ylabel('\gamma_a')
zlabel('R_0')
hold on
R0 = 2;
Z = R0 * ones(size(X));
surf(X,Y,Z)
colormap('turbo')
%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') ##### 댓글 수: 1없음 표시없음 숨기기
Bas123 2022년 12월 20일
It worked great. I do sincerely appreciate your support. Thanks a lot!

댓글을 달려면 로그인하십시오.

### 카테고리

Help CenterFile Exchange에서 Lighting, Transparency, and Shading에 대해 자세히 알아보기

### Community Treasure Hunt

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

Start Hunting!