********** MODEL NAME model_3cpt_linear_tc2 ********** MODEL NOTES PK model for simulation of drug concentration in central compartment with following characteristics: Compartments: 3 Elimination : linear Absorption : 3 transit compartments Unit convention Dose: mg Concentration: ug/mL Time: hours The annotation of the parameter units is consistent with the given unit convention. Units of the inputs (dose) and outputs (concentration) in the dataset for parameter estimation need to match the unit convention. ********** MODEL STATES d/dt(Ad) = -kabs*Ad + Fabs1*INPUT1 d/dt(TC1) = +kabs*Ad - ktr*TC1 d/dt(TC2) = + ktr*TC1 - ktr*TC2 d/dt(TC3) = + ktr*TC2 - ktr*TC3 d/dt(Ac) = + ktr*TC3 - Q1/Vc*Ac + Q1/Vp1*Ap1 - Q2/Vc*Ac + Q2/Vp2*Ap2 - CL/Vc*Ac d/dt(Ap1) = + Q1/Vc*Ac - Q1/Vp1*Ap1 d/dt(Ap2) = + Q2/Vc*Ac - Q2/Vp2*Ap2 Ad(0) = 0 TC1(0) = 0 TC2(0) = 0 TC3(0) = 0 Ac(0) = 0 Ap1(0) = 0 Ap2(0) = 0 ********** MODEL PARAMETERS Fabs1 = 1 # Relative bioavailability (-) kabs = 2 # Absorption rate parameter (1/hour) CL = 3 # Apparent clearance (L/hour) Vc = 32 # Apparent central volume (L) Q1 = 1 # Apparent intercompartmental clearance to first peripheral compartment (L/hour) Vp1 = 10 # Apparent first peripheral volume (L) Q2 = 1 # Apparent intercompartmental clearance to second peripheral compartment (L/hour) Vp2 = 10 # Apparent second peripheral volume (L) ktr = 0.5 # Transit compartment rate constant (1/hours) ********** MODEL VARIABLES % Calculation of concentration in central compartment Cc = Ac/Vc % Defining an output (only needed when interfacing with NLME % parameter estimation tools such as NONMEM and MONOLIX) OUTPUT1 = Cc # Compound concentration (ug/mL) ********** MODEL REACTIONS ********** MODEL FUNCTIONS ********** MODEL EVENTS