The three models share the same R0
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=sqrt(beta_m*beta_h*theta_h*theta_m/((theta_m+mu_m)*(gamma_h+mu_h)*(theta_h+mu_h)*(mu_m)))
R_0 =
[T,~]=gsua_dpmat(R_0,pars,1,
'R_0'
)
Introduce ranges in the following order:
ans =
T =
7×2 table
Range
Nominal
1 beta_m
0
0
0
2 beta_h
0
0
0
3 theta_m
0
0
0
4 theta_h
0
0
0
5 mu_m
0
0
0
6 gamma_h
0
0
0
7 mu_h
0
0
0
load
Results7.mat
common=intersect(T7.Properties.RowNames,T.Properties.RowNames)
common =
6×1 cell array
{'beta_h' } {'beta_m' } {'gamma_h'} {'mu_m' } {'theta_h'} {'theta_m'}
T{
'mu_h'
,
'Range'
}=[4e-4,4e-4];
T{
'mu_h'
,
'Nominal'
}=4e-4;
T{common,
'Range'
}=T7{common,
'Range'
};
%T.Nominal=[];
T{common,
'T7'
}=T7{common,
'Estlsqc'
}(:,1:476);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
load
Results8_complete.mat
T{common,
'T8'
}=T8{common,
'Estlsqc'
}(:,1:158);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
load
Results10_complete.mat
T{common,
'T10'
}=T10{common,
'Estlsqc'
}(:,1:136);
Warning: The new variables being added to the table have fewer rows than the table. They have been extended with rows containing default values.
T.T7(end,:)=repmat(4e-4,1,476);
T.T8(end,:)=repmat(4e-4,1,158);
T.T10(end,:)=repmat(4e-4,1,136);
clf
gsua_eval(T.T7,T);
clf
gsua_eval(T.T8,T);
clf
gsua_eval(T.T10,T);