Voxel-level mapping of tracer kinetics in pet studies: A statistical approach emphasizing tissue life tables

Finbarr O'Sullivan, Mark Muzi, David A. Mankoff, Janet F. Eary, Alexander M. Spence, Kenneth Krohn

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

Most radiotracers used in dynamic positron emission tomography (PET) scanning act in a linear time-invariant fashion so that the measured timecourse data are a convolution between the time course of the tracer in the arterial supply and the local tissue impulse response, known as the tissue residue function. In statistical terms the residue is a life table for the transit time of injected radiotracer atoms. The residue provides a description of the tracer kinetic information measurable by a dynamic PET scan. Decomposition of the residue function allows separation of rapid vascular kinetics from slower blood-tissue exchanges and tissue retention. For voxel-level analysis, we propose that residues be modeled by mixtures of nonparametrically derived basis residues obtained by segmentation of the full data volume. Spatial and temporal aspects of diagnostics associated with voxel-level model fitting are emphasized. Illustrative examples, some involving cancer imaging studies, are presented. Data from cerebral PET scanning with 18F fluoro-deoxyglucose (FDG) and 15O water (H2O) in normal subjects is used to evaluate the approach. Cross-validation is used to make regional comparisons between residues estimated using adaptive mixture models with more conventional compartmental modeling techniques. Simulations studies are used to theoretically examine mean square error performance and to explore the benefit of voxel-level analysis when the primary interest is a statistical summary of regional kinetics. The work highlights the contribution that multivariate analysis tools and life-table concepts can make in the recovery of local metabolic information from dynamic PET studies, particularly ones in which the assumptions of compartmental-like models, with residues that are sums of exponentials, might not be certain.

Original languageEnglish (US)
Pages (from-to)1065-1094
Number of pages30
JournalAnnals of Applied Statistics
Volume8
Issue number2
DOIs
StatePublished - 2014
Externally publishedYes

Fingerprint

Life Table
Positron emission tomography
Voxel
Kinetics
Tissue
Positron Emission Tomography
Scanning
Impulse response
Convolution
Mean square error
Blood
Decomposition
Imaging techniques
Recovery
Atoms
Life table
Model Fitting
Multivariate Analysis
Impulse Response
Mixture Model

Keywords

  • Kinetic analysis
  • Life-table
  • Mixture modeling
  • PET

ASJC Scopus subject areas

  • Statistics, Probability and Uncertainty
  • Modeling and Simulation
  • Statistics and Probability

Cite this

Voxel-level mapping of tracer kinetics in pet studies : A statistical approach emphasizing tissue life tables. / O'Sullivan, Finbarr; Muzi, Mark; Mankoff, David A.; Eary, Janet F.; Spence, Alexander M.; Krohn, Kenneth.

In: Annals of Applied Statistics, Vol. 8, No. 2, 2014, p. 1065-1094.

Research output: Contribution to journalArticle

O'Sullivan, Finbarr ; Muzi, Mark ; Mankoff, David A. ; Eary, Janet F. ; Spence, Alexander M. ; Krohn, Kenneth. / Voxel-level mapping of tracer kinetics in pet studies : A statistical approach emphasizing tissue life tables. In: Annals of Applied Statistics. 2014 ; Vol. 8, No. 2. pp. 1065-1094.
@article{54575ddee85e4f258b4e22c81d996471,
title = "Voxel-level mapping of tracer kinetics in pet studies: A statistical approach emphasizing tissue life tables",
abstract = "Most radiotracers used in dynamic positron emission tomography (PET) scanning act in a linear time-invariant fashion so that the measured timecourse data are a convolution between the time course of the tracer in the arterial supply and the local tissue impulse response, known as the tissue residue function. In statistical terms the residue is a life table for the transit time of injected radiotracer atoms. The residue provides a description of the tracer kinetic information measurable by a dynamic PET scan. Decomposition of the residue function allows separation of rapid vascular kinetics from slower blood-tissue exchanges and tissue retention. For voxel-level analysis, we propose that residues be modeled by mixtures of nonparametrically derived basis residues obtained by segmentation of the full data volume. Spatial and temporal aspects of diagnostics associated with voxel-level model fitting are emphasized. Illustrative examples, some involving cancer imaging studies, are presented. Data from cerebral PET scanning with 18F fluoro-deoxyglucose (FDG) and 15O water (H2O) in normal subjects is used to evaluate the approach. Cross-validation is used to make regional comparisons between residues estimated using adaptive mixture models with more conventional compartmental modeling techniques. Simulations studies are used to theoretically examine mean square error performance and to explore the benefit of voxel-level analysis when the primary interest is a statistical summary of regional kinetics. The work highlights the contribution that multivariate analysis tools and life-table concepts can make in the recovery of local metabolic information from dynamic PET studies, particularly ones in which the assumptions of compartmental-like models, with residues that are sums of exponentials, might not be certain.",
keywords = "Kinetic analysis, Life-table, Mixture modeling, PET",
author = "Finbarr O'Sullivan and Mark Muzi and Mankoff, {David A.} and Eary, {Janet F.} and Spence, {Alexander M.} and Kenneth Krohn",
year = "2014",
doi = "10.1214/14-AOAS732",
language = "English (US)",
volume = "8",
pages = "1065--1094",
journal = "Annals of Applied Statistics",
issn = "1932-6157",
publisher = "Institute of Mathematical Statistics",
number = "2",

}

TY - JOUR

T1 - Voxel-level mapping of tracer kinetics in pet studies

T2 - A statistical approach emphasizing tissue life tables

AU - O'Sullivan, Finbarr

AU - Muzi, Mark

AU - Mankoff, David A.

AU - Eary, Janet F.

AU - Spence, Alexander M.

AU - Krohn, Kenneth

PY - 2014

Y1 - 2014

N2 - Most radiotracers used in dynamic positron emission tomography (PET) scanning act in a linear time-invariant fashion so that the measured timecourse data are a convolution between the time course of the tracer in the arterial supply and the local tissue impulse response, known as the tissue residue function. In statistical terms the residue is a life table for the transit time of injected radiotracer atoms. The residue provides a description of the tracer kinetic information measurable by a dynamic PET scan. Decomposition of the residue function allows separation of rapid vascular kinetics from slower blood-tissue exchanges and tissue retention. For voxel-level analysis, we propose that residues be modeled by mixtures of nonparametrically derived basis residues obtained by segmentation of the full data volume. Spatial and temporal aspects of diagnostics associated with voxel-level model fitting are emphasized. Illustrative examples, some involving cancer imaging studies, are presented. Data from cerebral PET scanning with 18F fluoro-deoxyglucose (FDG) and 15O water (H2O) in normal subjects is used to evaluate the approach. Cross-validation is used to make regional comparisons between residues estimated using adaptive mixture models with more conventional compartmental modeling techniques. Simulations studies are used to theoretically examine mean square error performance and to explore the benefit of voxel-level analysis when the primary interest is a statistical summary of regional kinetics. The work highlights the contribution that multivariate analysis tools and life-table concepts can make in the recovery of local metabolic information from dynamic PET studies, particularly ones in which the assumptions of compartmental-like models, with residues that are sums of exponentials, might not be certain.

AB - Most radiotracers used in dynamic positron emission tomography (PET) scanning act in a linear time-invariant fashion so that the measured timecourse data are a convolution between the time course of the tracer in the arterial supply and the local tissue impulse response, known as the tissue residue function. In statistical terms the residue is a life table for the transit time of injected radiotracer atoms. The residue provides a description of the tracer kinetic information measurable by a dynamic PET scan. Decomposition of the residue function allows separation of rapid vascular kinetics from slower blood-tissue exchanges and tissue retention. For voxel-level analysis, we propose that residues be modeled by mixtures of nonparametrically derived basis residues obtained by segmentation of the full data volume. Spatial and temporal aspects of diagnostics associated with voxel-level model fitting are emphasized. Illustrative examples, some involving cancer imaging studies, are presented. Data from cerebral PET scanning with 18F fluoro-deoxyglucose (FDG) and 15O water (H2O) in normal subjects is used to evaluate the approach. Cross-validation is used to make regional comparisons between residues estimated using adaptive mixture models with more conventional compartmental modeling techniques. Simulations studies are used to theoretically examine mean square error performance and to explore the benefit of voxel-level analysis when the primary interest is a statistical summary of regional kinetics. The work highlights the contribution that multivariate analysis tools and life-table concepts can make in the recovery of local metabolic information from dynamic PET studies, particularly ones in which the assumptions of compartmental-like models, with residues that are sums of exponentials, might not be certain.

KW - Kinetic analysis

KW - Life-table

KW - Mixture modeling

KW - PET

UR - http://www.scopus.com/inward/record.url?scp=84903793480&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84903793480&partnerID=8YFLogxK

U2 - 10.1214/14-AOAS732

DO - 10.1214/14-AOAS732

M3 - Article

AN - SCOPUS:84903793480

VL - 8

SP - 1065

EP - 1094

JO - Annals of Applied Statistics

JF - Annals of Applied Statistics

SN - 1932-6157

IS - 2

ER -