fit_trtm - tpcclib 0.7.8 © 2022 by Turku PET Centre

NLLSQ estimation of R1 (=K1/K1'), k2', and k3 applying the transport-limited
reference tissue compartment model, TRTM (1). This model is based on
the reference tissue compartment model (2), but here it is assumed that
metabolism of tracer is irreversible (k4=0) during the PET scanning, and that
in (positive) reference tissue k3'>>k2', and thus the uptake in it is limited
only by transport into tissue (3,4).
 
Usage: fit_trtm [Options] ttacfile reference endtime resultfile
 
TTAC file can be in DFT or PMOD format. Sample times must be in minutes.
If TTAC file contains weights, those are used in the NLLSQ fitting.
Reference region TAC can be given separate TAC file or as the name or number
of the reference region in TTAC file.
 
Options:
 -rk2=<<value>|mean|median>
     Constrain k2' to specified <value>, or to mean or median of regional
     k2' values excluding reference region(s)
 -lim=<filename>
     Specify the constraints for model parameters;
     This file with default values can be created by giving this option
     as the only command-line argument to this program.
 -SD[=<y|N>]
     Standard deviations are calculated and saved in results (y), or
     not calculated (n).
 -CL[=<y|N>]
     95% Confidence limits are calculated and saved in results (y), or
     not calculated (n).
 -w1
     All weights are set to 1.0 (no weighting); by default, weights in
     TTAC file are used, if available.
 -wf
     Weight by sampling interval.
 -fit=<Filename>
     Fitted regional TACs are written in file.
 -svg=<Filename>
     Fitted and measured TACs are plotted in specified SVG file.
 -h, --help
     Display usage information on standard output and exit.
 -v, --version
     Display version and compile information on standard output and exit.
 -d[n], --debug[=n], --verbose[=n]
     Set the level (n) of debugging messages and listings.
 -q, --quiet
     Suppress displaying normal results on standard output.
 -s, --silent
     Suppress displaying anything except errors.
 
 
Values of R1, k2, and k3 are written in the specified result file.
Fitted curves are written in DFT format, if file name is given.
 
Example 1: file a789.tac contains regions-of-interest and positive reference
region, with name 'putam'. The whole time range is used in the fit.
Fitted TACs are plotted in SVG format in file a789trtm.svg.
     fit_trtm -svg=a789trtm.svg a789.tac putam 999 a789.res
 
Example 2: Reference region TAC is in a separate file, a789ref.tac.
Standard deviations and confidence limits are also estimated.
TAC data from injection to 60 min is used in the fitting.
     fit_trtm -SD=y -CL=y a789.tac a789ref.tac 60 a789.res
 
Example 3a: Create a file containing default parameter limits:
     fit_trtm -lim=trtm.lim
 
Example 3b: Apply user-defined parameter constraints specified in trtm.lim:
     fit_trtm -lim=trtm.lim a789.tac 'putam mean' 999 a789.tac
 
References:
1. Oikonen V. Model equations for reference tissue compartmental models.
   http://www.turkupetcentre.net/reports/tpcmod0002.pdf
2. Cunningham VJ, Hume SP, Price GR, Ahier RG, Cremer JE, Jones AKP.
   Compartmental analysis of diprenorphine binding to opiate receptors
   in the rat in vivo and its comparison with equilibrium data in vitro.
   J Cereb Blood Flow Metab 1991;11:1-9.
3. Herholz K, Lercher M, Wienhard K, Bauer B, Lenz O, Heiss W-D.
   PET measurement of cerebral acetylcholine esterase activity without
   blood sampling. Eur J Nucl Med 2001;28:472-477.
4. Nagatsuka S, Fukushi K, Shinotoh H, Namba H, Iyo M, Tanaka N, Aotsuka A,
   Ota T, Tanada S, Irie T. Kinetic analysis of [11C]MP4A using a high-
   radioactivity brain region that represents an integrated input function
   for measurement of cerebral acetylcholinesterase activity without
   arterial blood sampling. J Cereb Blood Flow Metab 2001; 21: 1354-1366.
 
See also: fitk3, tacweigh, fit_rrtm, lhtrtm, sim_rtcm, rescoll
 
Keywords: TAC, modelling, irreversible uptake, RTCM, reference input