Skip to content

Commit

Permalink
Rename HIV_model_5.m to Legacy Code/HIV_model_5.m
Browse files Browse the repository at this point in the history
  • Loading branch information
cody2 authored Jun 30, 2019
1 parent 8f3ff16 commit 5a66639
Showing 1 changed file with 17 additions and 17 deletions.
34 changes: 17 additions & 17 deletions HIV_model_5.m → Legacy Code/HIV_model_5.m
Original file line number Diff line number Diff line change
Expand Up @@ -36,8 +36,8 @@

% Initial Conditions
% P(01) = virion i.c. [log(/mL)]
% P(02) = CD8+ T cell (R+A+E) i.c. [/�L]
% P(03) = NK cell i.c. [/�L]
% P(02) = CD8+ T cell (R+A+E) i.c. [/]
% P(03) = NK cell i.c. [/]

% Pharmacokinetic Parameters
% P(04) = N803 injection site i.c. [pmol/kg]
Expand All @@ -58,7 +58,7 @@
% Parameters governing R,K proliferation and death
% P(16) = R death rate constant [/d]
% P(17) = K "
% P(18) = R max proliferating concentration [/�L]
% P(18) = R max proliferating concentration [/]
% P(19) = K "
% P(20) = R drug-induced proliferation factor [/d]
% P(21) = K "
Expand All @@ -70,12 +70,12 @@
% P(25) = E initial killing rate (a_E * E_0) [/d]
% P(26) = K "
% P(27) = R drug-induced activation factor [/d]
% P(28) = E drug-induced killing factor [�L/d]
% P(28) = E drug-induced killing factor [/d]
% P(29) = K "
% P(30) = R activation suppression factor []
% P(31) = E killing suppression factor []
% P(32) = K "
% P(33) = R half-activation virus concentration [/�L]
% P(33) = R half-activation virus concentration [/]

% Parameters governing CD8+ T cell clonal expansion and retraction
% P(34) = A death rate constant [/d]
Expand All @@ -85,10 +85,10 @@
% P(38) = E deactivation rate constant [/d]

% Parameters governing CD4+ T cells
% P(39) = CD4+ T cell i.c. [/�L]
% P(39) = CD4+ T cell i.c. [/]
% P(40) = T death rate constant [/d]
% P(41) = T max proliferating concentration [/�L]
% P(42) = T infection rate constant [/�L-d]
% P(41) = T max proliferating concentration [/]
% P(42) = T infection rate constant [/-d]
% P(43) = infected T decay rate constant [*(P(40)]

%% ========================================================================
Expand All @@ -97,9 +97,9 @@

% Y_MAIN = model output at each time point in 'SolTimes'
% Y_MAIN(:,1) = virion [log(#/mL)]
% Y_MAIN(:,2) = CD8+ T cell (R+A+E) [/�L]
% Y_MAIN(:,3) = NK cell [/�L]
% Y_MAIN(:,4) = CD4+ T cell [/�L]
% Y_MAIN(:,2) = CD8+ T cell (R+A+E) [/]
% Y_MAIN(:,3) = NK cell [/]
% Y_MAIN(:,4) = CD4+ T cell [/]

% Y_ALL = cell array of extra outputs (see 'output' section)

Expand Down Expand Up @@ -133,7 +133,7 @@
% -------------------------------------------------------------------------
% Declare initial conditions & calculated parameters ----------------------

V0 = 10^(P(01)-3) ;% convert V i.c. to [/�L]
V0 = 10^(P(01)-3) ;% convert V i.c. to [/]

% calculate i.c. for CD8+ T cell subsets (R,A,E)
CD8TCOEFS = [ ones(1,P(36)+2) ; zeros(P(36)+1,P(36)+2) ] ;
Expand All @@ -156,7 +156,7 @@
% redefine P(25:26) as E,K per-capita killing rate constants
P(25:26) = P(25:26)./ Yic(6:7) ;

% (infection rate constant) * (infected cell to virus conversion) [/�L-d]
% (infection rate constant) * (infected cell to virus conversion) [/-d]
beta = ( P(40)*P(43) + sum(P(25:26).*Yic(6:7)) ) / Yic(8) ;

% calculate [R;K;T] proliferation rate [/d]
Expand Down Expand Up @@ -305,9 +305,9 @@
Y = abs(Y) ;% removing negatives resulting from numerical error

Y_MAIN(:,1) = log10(Y(:,9))+3 ;% virion [log(#/mL)]
Y_MAIN(:,2) = sum(Y(:,[5,6,10:in(1)]),2) ;% CD8+ T cell (R+A+E) [/�L]
Y_MAIN(:,3) = Y(:,7) ;% NK cell [/�L]
Y_MAIN(:,4) = Y(:,8) ;% CD4+ T cell [/�L]
Y_MAIN(:,2) = sum(Y(:,[5,6,10:in(1)]),2) ;% CD8+ T cell (R+A+E) [/]
Y_MAIN(:,3) = Y(:,7) ;% NK cell [/]
Y_MAIN(:,4) = Y(:,8) ;% CD4+ T cell [/]

% -------------------------------------------------------------------------
% Prepare extra outputs 'Y_ALL' -----------------------------------------
Expand Down Expand Up @@ -467,4 +467,4 @@
Y_ALL{20}{3,1}(n) = Yn(4) ;
end

end
end

0 comments on commit 5a66639

Please sign in to comment.