p2t_v3c - tpcclib 0.8.0 © 2023 by Turku PET Centre

Deprecated, use sim_3tcm instead.
 
Simulation of PET tissue time-radioactivity concentration curve (TAC)
from arterial plasma (Ca) and blood (Cb) TACs, based on three-tissue
compartmental model, where the compartments are in series (default):
 
  ____    K1   ____    k3   ____    k5   ____     
 | Ca | ----> | C1 | ----> | C2 | ----> | C3 |    
 |____| <---- |____| <---- |____| <---- |____|    
          k2           k4           k6            
 
  dC1(t)/dt = K1*Ca(T) - (k2+k3)*C1(T) + k4*C2(T) 
  dC2(t)/dt = k3*C1(T) - (k4+k5)*C2(T) + k6*C3(T) 
  dC3(t)/dt = k5*C2(T) - k6*C3(T)                 
  Ct(T) = C1(T) + C2(T) + C3(T)                   
  Cvb(T) = Cab(T) - dCt(t)/dt / f                 
  Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T) 
 
, or, optionally, three-tissue compartmental model, where the 2nd and 3rd
tissue compartments are parallel, often used to represent specific and
non-specific binding:
 
  ____    K1   ____    k3   ____   
 | Ca | ----> | C1 | ----> | C2 |  
 |____| <---- |____| <---- |____|  
          k2           k4          
                | ^                
             k5 | | k6             
                v |                
               ____                
              | C3 |               
              |____|               
 
  dC1(t)/dt = K1*Ca(T) - (k2+k3+k5)*C1(T) + k4*C2(T) + k6*C3(T) 
  dC2(t)/dt = k3*C1(T) - k4*C2(T)  
  dC3(t)/dt = k5*C1(T) - k6*C3(T)  
  Ct(T) = C1(T) + C2(T) + C3(T)    
  Cvb(T) = Cab(T) - dCt(t)/dt / f                 
  Cpet(T)= Vb*fA*Cab(T) + Vb*(1-fA)*Cvb(T) + (1-Vb)*Ct(T) 
 
Usage: p2t_v3c [options] plasmafile bloodfile K1 k2 k3 k4 k5 k6 f Vb fA simfile
 
Options:
 -paral[lel]
     Model with parallel compartments C2 and C3 is applied.
 -ser[ies]
     Model with compartments C1, C2, and C3 in series is applied (default).
 -sub | -nosub
     TACs of sub-compartments (C1, C2 and C3) are written (-sub, default)
     or not written (-nosub) into the simfile.
 -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.
 
Enter blood flow (f) in units (ml/(min*ml)); set f=0 to assume that f>>K1
and Cvb=Cab. Enter the vascular volume fraction (Vb) and its arterial
fraction (fA) as percentages.
If the times in input files are in seconds, the units of rate constants
(k's) and blood flow (f) must also be specified as 1/sec.
For accurate results, plasma TAC should have very short sampling intervals.
To reduce the model, k5 or k3 can be set to 0.
 
Example:
     p2t_v3c -nosub plasma.dat blood.dat 0.3 0.2 0.1 0.2 0 0 0.6 5 30 sim.dat
 
Simulated TACs are written in ASCII format with columns:
  1) Sample time
  2) Total tissue activity concentration (Cpet)
  3) Activity concentration in 1st tissue compartment, (1-Vb)*C1
  4) Activity concentration in 2nd tissue compartment, (1-Vb)*C2
  5) Activity concentration in 3rd tissue compartment, (1-Vb)*C3
  6) Arterial contribution to tissue activity, Vb*fA*Cab
  7) Venous contribution to tissue activity, Vb*(1-fA)*Cvb
 
References:
1. TPCMOD0001.
 
See also: sim_3tcm, p2t_3c, p2t_di, tacadd, tacren, simframe, dft2img
 
Keywords: TAC, simulation, modelling, compartmental model