Sensitivity R0 comparison

Model 10

syms beta_m beta_h theta_m theta_h mu_m gamma_h mu_h R_0

pars=[beta_m beta_h theta_m theta_h mu_m gamma_h mu_h];

R_0=(beta_m*beta_h*theta_h*theta_m/((theta_m+mu_m)*(gamma_h+mu_h)*(theta_h+mu_h)*(mu_m)))^(1/2)

R_0 =

load Results7.mat

[T,~]=gsua_dpmat(R_0,pars,1,'R_0');

Introduce ranges in the following order:

ans =

common=intersect(T7.Properties.RowNames,T.Properties.RowNames);

T{common,'Range'}=T7{common,'Range'};

T{common,'Nominal'}=T7{common,'Nominal'};

T{'mu_h','Range'}=[4e-4,4e-4];

T{'mu_h','Nominal'}=4e-4

T = 7×2 table 

	Range		Nominal
1 beta_m	0	4.0000	2.0000
2 beta_h	0	4.0000	2.0000
3 theta_m	0.5800	0.8800	0.7300
4 theta_h	0.7000	1.7500	1.2250
5 mu_m	0.1600	0.2000	0.1800
6 gamma_h	0.5000	1.7500	1.1250
7 mu_h	0.0004	0.0004	0.0004

M=gsua_dmatrix(T,100000);
[T2,J,Y]=gsua_sa(M,T,'SensMethod','Xiao');
Progress: 12%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 23%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 34%
Estimated processing time (h:m:s): 0:0:1
Remaining time (h:m:s): 0:0:1
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:39
Number of simulations: 450000
Progress: 45%
Estimated processing time (h:m:s): 0:0:1
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:39
Number of simulations: 450000
Progress: 56%
Estimated processing time (h:m:s): 0:0:1
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 67%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 78%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 89%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000
Progress: 100%
Estimated processing time (h:m:s): 0:0:0
Remaining time (h:m:s): 0:0:0
Elapsed time (h:m:s): 0:0:0
Estimated stop time (h:m:s): 17:23:38
Number of simulations: 450000

clf

gsua_plot('Pie',T2,T2.STi);

gsua_plot('Pie',T2,T2.Si);

gsua_ua(M,T,'xdata',0)

Progress: 100% Estimated processing time (h:m:s): 0:0:0 Remaining time (h:m:s): 0:0:0 Elapsed time (h:m:s): 0:0:0 Estimated stop time (h:m:s): 17:23:55 Number of simulations: 100000

ans = 100000×1

Are there strong interactions?

sum(T2.STi_vec)
ans = 1.0865
sum(T2.Si_vec)
ans = 0.9187

Fiability

sum(abs(T2.Si_vec))
ans = 0.9187
save('Results_r0.mat')