Processing Overview

<< Click to Display Table of Contents >>

Navigation:  PKIN Data Processing >

Processing Overview

Data processing of studies with PET or SPECT tracers typically consists of the following parts:

1. Data Loading

Create a new workspace (tab) in the PKIN tool.

Load the whole-blood curve (total activity concentration of tracer in  the blood samples).

Load the plasma activity curve. Depending on the preparation of the blood data the plasma activity curve can represent (1) the unchanged tracer (parent) in plasma, or (2) total tracer in the plasma of the blood samples including the metabolites. In the latter case a metabolite correction has to be performed in PKIN.

An alternative to loading the plasma activity curve is loading the plasma fraction curve. In this case the total plasma activity is calculated by multiplying whole blood activity with the plasma fraction.

Load the parent fraction curve (fraction of unchanged tracer in plasma). This is only necessary if the loaded plasma curve represents total tracer activity in plasma including metabolites.

Load the average time-activity curves of one or multiple tissue regions.

2. Input Curve Configuration

Define interpolation functions (also called models) for the different types of blood data and fit them to the data. The default is linear interpolation between measurements, and a parent fraction identical to 1.

The input curve is obtained as the multiplication of the plasma model with the parent fraction model.

3. Kinetic Model Configuration

Select a region with a "typical" time-activity curve (TAC) from the Region list.

Select a simple model from the Model list.

Enable the parameters to be fitted by checking their boxes; unchecked parameters remain at the initial values entered.

Select a weighting scheme of the residuals; variable or constant (default) weighting is available.

4. Kinetic Model Fitting

Fit the model parameters by activating the Fit model, current region button.

Consult the residuals to check whether the model is adequate; there should ideally be no bias in the residuals, just random noise.

If the model fit is fine, configure Copy the model to all regions to Model & Par and activate the button to establish the same initial model configuration for all TACs. Then activate Fit model, all regions to adjust the parameters to the regional measurements.

Check the fit result for all the TACs by stepping through the regions. Note that mouse scrolling is supported when placing the cursor on top of the Region list.

If the model is not yet fine, test more complex models.

5. Kinetic Model Comparison

Switch between compartment models of different complexity and fit. The parameters are either maintained for each model type, or converted, according to the Parameter conversion setting in the Extras Panel panel.

The results of the different fits are stored in the Model History which can be used for inspecting the results per region.

Check the residuals for judging model adequacy.

Check the different criteria on the Details tab (Schwartz Criterion SC, Akaike Information Criterion AIC, Model Selection Criterion MSC) to decide whether a more complex model is supported by the data.

Check for parameter identifiability. As an indicator of the parameter identifiability standard errors (%SE) are returned from the fit. They should remain limited for all relevant parameters. Additionally, Monte Carlo simulations can be performed to obtain distribution statistics of the parameter estimates.

If justified by physiology, try to improve the stability of parameter estimation by enforcing common parameters among regions in a Coupled fitting procedure.

Compare the outcome of compartment models with that of other models, such as reference models or graphical plots.

6. Saving

Save all model information together with the data in a composite text .km file which allows restoring a session including the fitting history.

Save a summary of all regional parameters in an EXCEL-ready text file .kinPar.

7. Batch Processing

Lengthy calculations such as coupled fits with many regional TACs or Monte Carlo simulations can be run in a batch job. The results are saved both as .km files and in a .kinPar text file.