From f737246abc1f1a1acbee4ec125d24745a3c1b29f Mon Sep 17 00:00:00 2001 From: "Cody, Jonathan William" Date: Mon, 27 Jun 2022 19:16:59 -0400 Subject: [PATCH] Rename N803_single.m to ICSolvers/N803_single.m --- N803_single.m => ICSolvers/N803_single.m | 38 ++++++++++++------------ 1 file changed, 19 insertions(+), 19 deletions(-) rename N803_single.m => ICSolvers/N803_single.m (88%) diff --git a/N803_single.m b/ICSolvers/N803_single.m similarity index 88% rename from N803_single.m rename to ICSolvers/N803_single.m index ad44d84..d8b7188 100644 --- a/N803_single.m +++ b/ICSolvers/N803_single.m @@ -8,12 +8,12 @@ % | Pienaar Computational Systems Pharmacology Lab | % \--------------------------------------------------------------/ % -% Nomenclature: V = SIV virions [#/µL] -% T8 = total CD8+ T cells [#/µL] -% E0 = resting SIV-specific CD8+ T cells [#/µL] -% Ea = active SIV-specific CD8+ T cells [#/µL] -% B0 = resting bystander CD8+ T cells [#/µL] -% Ba = active bystander CD8+ T cells [#/µL] +% Nomenclature: V = SIV virions [#/ç„¡] +% T8 = total CD8+ T cells [#/ç„¡] +% E0 = resting SIV-specific CD8+ T cells [#/ç„¡] +% Ea = active SIV-specific CD8+ T cells [#/ç„¡] +% B0 = resting bystander CD8+ T cells [#/ç„¡] +% Ba = active bystander CD8+ T cells [#/ç„¡] % X = N803 at absorption site [pmol/kg] % C = N803 plasma concentration [pM] % R = regulation [] (dimensionless quantity) @@ -52,7 +52,7 @@ % Rename inputed parameters ----------------------------------------------- Vi = AllPars(01) ;% V initial value [log(#/mL)] -EBi = AllPars(02) ;% E+B initial value [#/µL] +EBi = AllPars(02) ;% E+B initial value [#/ç„¡] fE = AllPars(03) ;% initial frequency: E/(E+B) fEA = AllPars(04) ;% initial frequency: Ea/E q = AllPars(05) ;% V growth rate (if E+B were absent) [/d] @@ -61,7 +61,7 @@ V50B = AllPars(08) ;% 50% viral stimulation saturation for B [#/mL] mEn = AllPars(09) ;% normalized Ea reversion rate constant [] mBn = AllPars(10) ;% normalized Ba reversion rate constant [] -EB50 = AllPars(11) ;% 50% E+B proliferation saturation [#/µL] +EB50 = AllPars(11) ;% 50% E+B proliferation saturation [#/ç„¡] pE = AllPars(12) ;% Ea proliferation rate constant [/d] pB = AllPars(13) ;% Ba proliferation rate constant [/d] d = AllPars(14) ;% E0/B0 death rate constant [/d] @@ -85,9 +85,9 @@ %% ------------------------------------------------------------------------ % Calculate some initial conditions & parameters -------------------------- -Vi = 10^(Vi-3) ;% V initial value [#/µL] -V50E = V50E/1000 ;% 50% viral stimulation saturation for E [#/µL] -V50B = V50B/1000 ;% 50% viral stimulation saturation for B [#/µL] +Vi = 10^(Vi-3) ;% V initial value [#/ç„¡] +V50E = V50E/1000 ;% 50% viral stimulation saturation for E [#/ç„¡] +V50B = V50B/1000 ;% 50% viral stimulation saturation for B [#/ç„¡] % restrict mE and mB such that initial activation aE and aB are positive UE = (2*pE/(pE+dA))^7 ; @@ -139,8 +139,8 @@ p0 = pi * (EB50+EBi)/EB50 * (1+p2) ;% E0/B0 proliferation rate con [/d] p1 = pm/p0 ;% E0/B0 proliferation stimulation factor -gE0 = q / ( EA/(1+gE2) + psi*BA/(1+gB2) ) ;% Ea killing rate con [µL/#-d] -gB0 = psi*gE0 ;% Ba killing rate constant [µL/#-d] +gE0 = q / ( EA/(1+gE2) + psi*BA/(1+gB2) ) ;% Ea killing rate con [ç„¡/#-d] +gB0 = psi*gE0 ;% Ba killing rate constant [ç„¡/#-d] sE = dR / ( Vi/(V50E+Vi) + sig*Vi/(V50B+Vi) ) ;% regulation gen by E act sB = sig*sE ;% regulation generation by B activation @@ -149,17 +149,17 @@ % Prepare parameter and initial value vectors and call 'N803_model_2' ----- Pars(01) = q ;% V growth rate (if E+B were absent) [/d] -Pars(02) = gE0 ;% Ea killing rate constant [µL/#-d] -Pars(03) = gB0 ;% Ba killing rate constant [µL/#-d] +Pars(02) = gE0 ;% Ea killing rate constant [ç„¡/#-d] +Pars(03) = gB0 ;% Ba killing rate constant [ç„¡/#-d] -Pars(04) = V50E ;% 50% viral stimulation saturation for E [#/µL] -Pars(05) = V50B ;% 50% viral stimulation saturation for B [#/µL] +Pars(04) = V50E ;% 50% viral stimulation saturation for E [#/ç„¡] +Pars(05) = V50B ;% 50% viral stimulation saturation for B [#/ç„¡] Pars(06) = aE0 ;% E0 activation rate constant [/d] Pars(07) = aB0 ;% B0 activation rate constant [/d] Pars(08) = mE ;% Ea reversion rate constant [/d] Pars(09) = mB ;% Ba reversion rate constant [/d] -Pars(10) = EB50 ;% 50% E+B proliferation saturation [#/µL] +Pars(10) = EB50 ;% 50% E+B proliferation saturation [#/ç„¡] Pars(11) = p0 ;% E0/B0 proliferation rate constant [/d] Pars(12) = pE ;% Ea proliferation rate constant [/d] Pars(13) = pB ;% Ba proliferation rate constant [/d] @@ -193,4 +193,4 @@ Y_OUT = N803_model_2(SoluTimes,DoseTimes,Pars,Yic,varargin) ; -end \ No newline at end of file +end