fitk2 - tpcclib 0.8.0 © 2023 by Turku PET Centre

Non-linear fitting of one-tissue compartment model to plasma input, blood,
and tissue time-activity curves (PTAC, BTAC, and TTAC) to estimate
parameters K1, k2, and Vb. Sample times must be in minutes.
To fit only parameters K1 and Vb, assuming k2=0, create limit file (see options),
and set the lower and upper limit for K1/k2 to zero.
    ______        ______    
   |      |  K1  |      |   
   |  Ca  | ---> |  C1  |   
   |______|      |______|   
                    | k2    
Usage: fitk2 [Options] ptacfile btacfile ttacfile endtime resultfile
     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.
     Without filename the default values are printed on screen.
     Standard deviations are calculated and saved in results (Y, default),
     or not calculated (n).
     Program runs a lot faster if SD and CL are not calculated.
     95% Confidence limits are calculated and saved in results (y), or
     not calculated (N, default).
     Enter a fixed Vb; fitted by default.
 -<BPnd|BPp|DVR>=<Reference region id or filename>
     Optional reference region is used to calculate BPnd, BPp, or DVR;
     BPnd=DVroi/DVref-1, BPp=DVroi-DVref, and DVR=DVroi/DVref
     Fitted regional TACs are written in DFT format.
     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.
Example 1: estimate K1, K1/k2 and Vb, with estimation of parameter SD
     fitk2 -sd=y a919ap.bld a919ab.bld a919.tac 60 a919k2.res
Example 2: estimate K1 and DV (=K1/k2); Vb is constrained to 0%;
DVRs are calculated by dividing DVs by the DV of region 'cer'
     fitk2 -Vb=0 -R=cer p25apc.bld none p25.tac 60 p25_k2.res
See also: logan, lhsoldv, fitk2di, fitk4, fit_h2o, p2t_v3c, tacweigh, taccbv
Keywords: TAC, modelling, distribution volume, reversible uptake, 1TCM