fitglob 1.3.0 © 2002-2012 by Turku PET Centre
Non-linear fitting of compartmental model to plasma and tissue
time-activity curves (TACs) to estimate K1-k6 and Va using
Iterative Topografical Global Optimization routine (ITGO)
developed by Aimo Torn and Sami Viitanen (1, 2).
Command-line parameters:
1) Parameter file (see below)
2) Plasma file
3) Blood file
4) Tissue TAC file (*.dft)
5) Fit end time (duration)
6) Result file: parameters at global minimum
Options:
-fit=<Filename>
Fitted regional TACs are written in DFT format.
-svg=<Filename>
Fitted and measured TACs are plotted in specified SVG file.
-h or --help
Print this message and exit.
--version or --build
Print software build information and exit.
--silent
Program works silently, printing only error and warning messages.
--verbose
Program prints more information about what it is doing.
Example:
fitglob par.set a919ap.dat a919ab.dat a919.dft 999 a919.res a919fit.dft
The parameter file contents must be exactly as follows:
1) Model identifier:
PR = parallel (receptor) model
with parameters K1, K1/k2, k3, k3/k4, k5, k5/k6, Va
P = parallel (receptor) model
with parameters K1, k2, k3, k4, k5, k6, Va
SR = serial (metabolism) model
with parameters K1, K1/k2, k3, k3/k4, k5, k5/k6, Va
S = serial (metabolism) model
with parameters K1, k2, k3, k4, k5, k6, Va
LR = kLoss model
with parameters K1 K1/k2 k3 kLoss, Va
L = kLoss model
with parameters K1, k2, k3, kLoss, Va
KR = kLoss model2
with parameters K1, K1/k2, k3, k3/k4, k5, k5/k6, kLoss, Va
K = kLoss model2
with parameters K1, k2, k3, k4, k5, k6, kLoss, Va
2) Lower and upper limits for each of the parameters
Parameter can be fixed to a certain value by setting equal limits;
fixing a ratio k3/k4 or k5/k6 to zero sets k4 or k6 to zero.
3) Value of K1/F (0, if blood flow is not considered);
the reduced venous blood activity and its dependence on blood flow can
corrected by specifying the ratio K1/F.
4) Number of points to sample in one iteration (prfmNr)
If you don't know what to put here, set it to 0.
Parameter file example 1:
PR
# K1
0 0.5
# K1/k2
0.001 10.0
# k3
0 0.0
# k3/k4
0 0.0
# k5
0 0.2
# k5/k6
0.001 10.0
# Va
0 0.30
# K1/F
0
#prfmNr
100
Parameter file example 2:
KR
# K1
0 0.5
# K1/k2
0.001 10.0
# k3
0 0.0
# k3/k4
0 0.0
# k5
0 0.2
# k5/k6
0.001 10.0
# kLoss
0 0.2
# Va
0 0.30
# K1/F
0
#prfmNr
100
References:
1. Torn A, Viitanen S. Topografical global optimization. In: C.A. Floudas
and P.M. Pardalos (eds.) Recent advances in Global Optimization,
Princeton University Press, 1992.
2. Sederholm K. Globaali optimointi positroniemissiotomografia-kuvantamiseen
liittyvässä mallintamisessa. Pro Gradu. Turun yliopisto, 2003;
http://www.turkupetcentre.fi/gradu/sederholm.pdf
See also:
fitk2,
fitk3,
fitk4,
fitkloss,
lhsol,
p2t_v3c,
dftweigh
Keywords: global optimization, compartment model, fitting, modelling
This program comes with ABSOLUTELY NO WARRANTY. This is free software, and
you are welcome to redistribute it under GNU General Public License.