Robustness analysis table: effect of individual parameter variation on dynamics of the whole system. Here CRD and non-CRD indicate control response dynamics and non-control response dynamics respectively. By varying each parameter value individually, we observed and recorded the ranges for CRD model dynamics and non-CRD model dynamics

ParameterRangeDynamics
α0 < α < 0.08stable CRD
β0 < β < 10.40stable CRD
γ0 < γ < 0.0002stable CRD
δ0 < δ < 0.0028stable CRD
ξ0 < ξ < 1.33stable non-CRD
1.33 < ξ < 9.8stable CRD
9.8 < ξ < 14stable non-CRD
p10 < p1 < 0.009stable non-CRD
0.009 < p1 < 0.031stable CRD
0.031 < p1 < 0.035stable non-CRD
q10 < q1 < 0.1036stable CRD
p20 < p2 < 7.2336stable CRD
q20 < q2 < 0.2stable CRD
p30 < p3 < 14.68stable CRD
q30 < q3 < 0.6stable non-CRD
0.6 < q3 < 3stable CRD
ces0 < ces < 2000.3624stable CRD
ks0 < ks < 100stable CRD
τs0 < τs < 60.3504stable CRD