diff --git a/ICsolvers/N803_shared.m b/N803_shared.m similarity index 90% rename from ICsolvers/N803_shared.m rename to N803_shared.m index 0510cb1..de3be4f 100644 --- a/ICsolvers/N803_shared.m +++ b/N803_shared.m @@ -8,12 +8,12 @@ % | Pienaar Computational Systems Pharmacology Lab | % \--------------------------------------------------------------/ % -% Nomenclature: V = SIV virions [#/µL] -% T8 = total CD8+ T cells [#/µL] -% S0 = resting SIV-specific CD8+ T cells [#/µL] -% Sa = active SIV-specific CD8+ T cells [#/µL] -% N0 = resting non-SIV-specific CD8+ T cells [#/µL] -% Na = active non-SIV-specific CD8+ T cells [#/µL] +% 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 [#/µL] (cohort 1) -SNi(2) = AllPars(04) ;% S+N initial value [#/µL] (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 [#/µL] +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 [#/µL] -V50S = V50S/1000 ;% 50% viral stimulation saturation for S [#/µL] -V50N = V50N/1000 ;% 50% viral stimulation saturation for N [#/µL] +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 [µL/#-d] -gN0 = psi * gS0 ;% Na killing rate constant [µL/#-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 [µL/#-d] -Pars(03) = gN0 ;% Na killing rate constant [µL/#-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 [#/µL] -Pars(05) = V50N ;% 50% viral stimulation saturation for N [#/µL] +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 [#/µL] +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 \ No newline at end of file +end