Compartmental model based spectral analysis of regional PET data. Basis functions (from 1 to n) are formed with a range of k2 from zero to a given maximum, with the provided PTAC as input function. Model Croi(t) = Vb*Cb(t) + Ct[1](t) + Ct[2](t) + ... + Ct[n](t) is fitted to the provided TTACs. ____ _________ | |--K1[i=1]->| Ct[i=1] | | |<-k2[i=1]--|_________| | Cp | ....... | | _________ | |--K1[i=n]->| Ct[i=n] | | |<-k2[i=n]--|_________| |____| Croi = Vb*Cb + K1[i=1]*Ct[i=1] + ... K1[i=n]*Ct[i=n] Usage: tacsa [options] ptacfile btacfile ttacfile maxk2 parfile [safile] Options: -n=<Number of basis functions> Set the number of basis functions; by default 100, minimum 10. -mink2=<Value> Set minimum value for k2 in units 1/min; by default 0.005 min-1. Do not set to zero or very low value, but if you want to add basis function with k2=0 (trapping compartment), enter minimum k2 as negative value. -bf=<filename> Basis function curves are written in specified 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. Parameters Vb and K1 will be written in parameter file as the weights of the basis functions (Vb as the weight of BTAC, and K1 as the sum of weights of other basis functions), and k2 as the weighted average of k2-derived basis functions. In the optional SA file, the weights for each k2-derived basis functions are listed in TAC format. See also: fitdelay, bfmh2o, lhsol, hist2svg Keywords: TAC, modelling, compartmental model, LLSQ

This program comes with ABSOLUTELY NO WARRANTY. This is free software, and you are welcome

to redistribute it under GNU General Public License.

Source codes are available in https://gitlab.utu.fi/vesoik/tpcclib.git.