# Compartmental model ODEs

Compartmental models for PET data can be described in terms of a set of linear, first-order, constant-coefficient, ordinary differential equations (ODEs). Change of tracer concentration in one of the tissue compartments is a linear function of the concentrations in all other compartments:

From the linearity follows that kinetic measurements are the *convolution* of the tracer
input function and the response function of the system.

Convolution:

For example, for two-tissue compartment model with
tissue compartments C_{1} and C_{2}, input function C_{0}(t), and
rate constants
*K _{1}*,

*k*,

_{2}*k*, and

_{3}*k*, can be described mathematically with a system of two ODEs:

_{4}The simulated tissue TACs are used to

- test and validate analysis methods and software
- estimate compartmental model rate constants (model fitting).

## Analytic solutions using Laplace transformation

Analytic solutions for linear first-order differential
equations can be derived using *Laplace transformation*.
For this 2TCM example the solution is
(Phelps et al., 1979):

, where ⊗ denotes the operation of convolution. For one-tissue compartmental model the solution is a lot simpler:

For simplified reference tissue model the solution is:

## Direct solutions for discrete-time data

Direct numerical approaches based on simple (Eulerian) integration are conceptually simpler and
computationally more effective than analytic solutions requiring convolution operations
(Graham, 1985).
Adams-Bashforth Predictor requires present and past values of *x* (*x ^{m}*
and

*x*) to predict

^{m-1}*x*:

^{m+1}Δt is the sample time difference.
Fourth-order Adams Predictor requires four previous values to predict *x ^{m+1}*:

As a fourth-order algorithm its errors are similar to the fourth-order Runge-Kutta algorithm (DiStefano, 2013). The implicit Second-order Adams-Moulton Corrector requires only one value of the past:

The Second-order Adams-Moulton method has been applied to PET data modelling (Kuwabara et al., 1993). Solutions for the most common compartmental models using this method are shown in TPC modelling reports.

As an example, the solution of the one-tissue compartmental model ODE,

starts by integrating the equation, assuming that at time zero all concentrations are zero:

Integrals are calculated using trapezoidal method.
Concentrations in tissue compartments are calculated starting at time zero when all concentrations
and integrals are assumed to be zero, and then proceeding one input function sample at a time, using
the results from previous sample to calculate the next ones.
Definite integral of *n* th compartment is implicitly estimated
(Kuwabara et al., 1993) based on second-order
Adams-Moulton method with Euler integration:

In this example we have just one compartment (*n*=1).
After substitution and rearrangement:

For simplified reference tissue model the solution is:

Using this method the TACs from compartmental models can be easily simulated even in spreadsheet software. For an example, you can download an Excel file where PET tissue data can be simulated using two-tissue compartmental model from measured plasma data and user definable rate constants. Excel files for simulating tissue data with three-tissue compartmental models are also available (compartments in parallel and in series).

## Linearization

Some compartment model ODEs can be linearized into multilinear equations from which the model parameters can be solved using standard linear optimization methods. These methods may produce somewhat biased results with noisy data, but the computations are very fast, enabling pixel-by-pixel calculation of parametric images.

## Comparison of solution methods

The convolution method provides precise results only if the sampling time interval is very short
(relative to the rate constants). The alternative method is computationally faster, and produces
correct results (as far as it is possible from discrete data) without any interpolation;
the convolution method only approaches the correct result with shorter interpolated sampling
intervals (that is, smaller *Δt* is required in the convolution method).
As an example, convolution and the alternative Adams-Moulton method are
compared in a simulation of one-tissue
compartmental model.

With data sampled at discrete times some errors in simulation will always be present. The longer the sampling time frame durations are, the more important it is to appropriately account for the frame durations in calculating AUC (integral).

## See also:

- Fitting compartmental models
- Compartmental model analysis of regional TACs
- Reference tissue input compartmental models
- PET data
- Input function
- Dispersion
- ODEs for CM with dual plasma input: compartments in parallel and in series

## Literature

Budinger TF, Huesman RH, Knittel B, Friedland RP, Derenzo SE (1985):
Physiological modeling of dynamic measurements of metabolism using positron emission tomography.
*In:* The Metabolism of the Human Brain Studied with Positron Emission Tomography.
(Eds: Greitz T et al.) Raven Press, New York, 165-183.

DiStefano III J. *Dynamic Systems Biology Modeling and Simulation*. Academic Press, 2013.
ISBN: 9780124104112.

Graham MM. Model simplification: complexity versus reduction. *Circulation* 1985;
72(5 Pt 2):IV63-IV68. PMID: 4053329.

Gunn RN, Gunn RS, Cunningham VJ. Positron emission tomography compartmental models.
*J Cereb Blood Flow Metab.* 2001; 21: 635-652.
doi: 00004647-200106000-00002.

Huang SC, Phelps ME (1986): Principles of tracer kinetic modeling in positron emission tomography
and autoradiography. *In:* Positron Emission Tomography and Autoradiography: Principles and
Applications for the Brain and Heart.
(Eds: Phelps,M; Mazziotta,J; Schelbert,H) Raven Press, New York, 287-346.

Kuwabara H, Cumming P, Reith J, Léger G, Diksic M, Evans AC, Gjedde A. Human striatal L-DOPA
decarboxylase activity estimated in vivo using 6-[^{18}F]fluoro-DOPA and positron emission
tomography: error analysis and application to normal subjects.
*J Cereb Blood Flow Metab.* 1993; 13: 43-56.
doi: 10.1038/jcbfm.1993.7.

Phelps ME, Huang S-C, Hoffman EJ, Selin C, Sokoloff L, Kuhl DE.
Tomographic measurement of local cerebral glucose metabolic rate in humans with
[F-18]2-fluoro-2-deoxy-D-glucose: validation of method. *Ann Neurol.* 1979; 6: 371-388.
doi: 10.1002/ana.410060502.

Tags: Modeling, Compartmental model, Simulation, ODE, Convolution

Updated at: 2023-06-14

Created at: 2004-06-03

Written by: Vesa Oikonen