Rename ICsolvers/N803_shared.m to N803_shared.m

ElsjePienaarGroup · Dec 25, 2022 · 682b4a0 · 682b4a0
1 parent 11c2855
commit 682b4a0
Showing 1 changed file with 20 additions and 20 deletions.
diff --git a/ICsolvers/N803_shared.m → N803_shared.m b/ICsolvers/N803_shared.m → N803_shared.m
@@ -8,12 +8,12 @@
 %     |               Pienaar Computational Systems Pharmacology Lab |
 %     \--------------------------------------------------------------/
 %
-% Nomenclature: V  = SIV virions [#/無]
-%               T8 = total CD8+ T cells [#/無]
-%               S0 = resting SIV-specific CD8+ T cells [#/無]
-%               Sa = active  SIV-specific CD8+ T cells [#/無]
-%               N0 = resting non-SIV-specific CD8+ T cells [#/無]
-%               Na = active  non-SIV-specific CD8+ T cells [#/無]
+% Nomenclature: V  = SIV virions [#/無]
+%               T8 = total CD8+ T cells [#/無]
+%               S0 = resting SIV-specific CD8+ T cells [#/無]
+%               Sa = active  SIV-specific CD8+ T cells [#/無]
+%               N0 = resting non-SIV-specific CD8+ T cells [#/無]
+%               Na = active  non-SIV-specific CD8+ T cells [#/無]
 %               X  = N803 at absorption site [pmol/kg]
 %               C  = N803 plasma concentration [pM]
 %               R  = regulation [] (dimensionless quantity)
@@ -75,8 +75,8 @@
 % Rename inputed parameters -----------------------------------------------
 Vi(1)  = AllPars(01) ;% V initial value [log(#/mL)] (cohort 1)
 Vi(2)  = AllPars(02) ;% V initial value [log(#/mL)] (cohort 2)
-SNi(1) = AllPars(03) ;% S+N initial value [#/無] (cohort 1)
-SNi(2) = AllPars(04) ;% S+N initial value [#/無] (cohort 2)
+SNi(1) = AllPars(03) ;% S+N initial value [#/無] (cohort 1)
+SNi(2) = AllPars(04) ;% S+N initial value [#/無] (cohort 2)
 fS(1)  = AllPars(05) ;% initial frequency: S/(S+N) (cohort 1)
 fSn    = AllPars(06) ;% normalized S/(S+N) (co 2)
 fS(2)  = fS(1)+(0.3-fS(1))*fSn ;% initial frequency: S/(S+N) (cohort 2)
@@ -88,7 +88,7 @@
 mSn  = AllPars(11)  ;% normalized Sa reversion rate constant []
 mNn  = AllPars(12)  ;% normalized Na reversion rate constant []
 
-SN50 = AllPars(13)  ;% 50% S+N proliferation saturation [#/無]
+SN50 = AllPars(13)  ;% 50% S+N proliferation saturation [#/無]
 p0n  = AllPars(14)  ;% nomalized S0/N0 proliferation rate constant [/d]
 pS   = AllPars(15)  ;% Sa proliferation rate constant [/d]
 pN   = AllPars(16)  ;% Na proliferation rate constant [/d]
@@ -112,9 +112,9 @@
 %% ------------------------------------------------------------------------
 % Calculate some initial conditions & parameters --------------------------
 
-Vi = 10.^(Vi - 3) ;% converting V initial value [#/無]
-V50S = V50S/1000  ;% 50% viral stimulation saturation for S [#/無]
-V50N = V50N/1000  ;% 50% viral stimulation saturation for N [#/無]
+Vi = 10.^(Vi - 3) ;% converting V initial value [#/無]
+V50S = V50S/1000  ;% 50% viral stimulation saturation for S [#/無]
+V50N = V50N/1000  ;% 50% viral stimulation saturation for N [#/無]
 YS = Vi ./ (Vi+ V50S) ;% initial V/(V50S+V) [cohort 1,2]
 YN = Vi ./ (Vi+ V50N) ;% initial V/(V50N+V) [cohort 1,2]
 
@@ -166,8 +166,8 @@
 beta = psi * NA ./ (1+gN2*Ri)  ;  z = beta(1) - beta(2)  ;
 a = z*R  ;  b = SA(1)*R - SA(2) + z*(1+R)  ;  c = SA(1) - SA(2) + z  ;
 gS2 = ( -b + sqrt(b^2 - 4*a*c) ) / (2*a)  ;% S killing regulation factor []
-gS0 = q / ( SA(1)/(1+gS2) + beta(1) )  ;% Sa killing rate constant [無/#-d]
-gN0 = psi * gS0                        ;% Na killing rate constant [無/#-d]
+gS0 = q / ( SA(1)/(1+gS2) + beta(1) )  ;% Sa killing rate constant [無/#-d]
+gN0 = psi * gS0                        ;% Na killing rate constant [無/#-d]
 
 %% Do for each cohort (NOT indenting loop) ================================
 for c = 1:2
@@ -187,17 +187,17 @@
 % Prepare parameter and initial value vectors and call 'N803_model_2' -----
 
 Pars(01) = q     ;% V growth rate (if S+N were absent) [/d]
-Pars(02) = gS0   ;% Sa killing rate constant [無/#-d]
-Pars(03) = gN0   ;% Na killing rate constant [無/#-d]
+Pars(02) = gS0   ;% Sa killing rate constant [無/#-d]
+Pars(03) = gN0   ;% Na killing rate constant [無/#-d]
 
-Pars(04) = V50S  ;% 50% viral stimulation saturation for S [#/無]
-Pars(05) = V50N  ;% 50% viral stimulation saturation for N [#/無]
+Pars(04) = V50S  ;% 50% viral stimulation saturation for S [#/無]
+Pars(05) = V50N  ;% 50% viral stimulation saturation for N [#/無]
 Pars(06) = aS0   ;% S0 activation rate constant [/d]
 Pars(07) = aN0   ;% N0 activation rate constant [/d]
 Pars(08) = mS    ;% Sa reversion rate constant [/d]
 Pars(09) = mN    ;% Na reversion rate constant [/d]
 
-Pars(10) = SN50  ;% 50% S+N proliferation saturation [#/無]
+Pars(10) = SN50  ;% 50% S+N proliferation saturation [#/無]
 Pars(11) = p0    ;% S0/N0 proliferation rate constant [/d]
 Pars(12) = pS    ;% Sa proliferation rate constant [/d]
 Pars(13) = pN    ;% Na proliferation rate constant [/d]
@@ -252,4 +252,4 @@
 Y_OUT = [ Y_Cohort{1} , Y_Cohort{2} ] ;
 PARS  = [ P_Cohort{1} ; P_Cohort{2} ] ;
 
-end
+end