Practical Priors for Bayesian Inference of Latent Biomarkers.

Latent biomarkers are quantities that strongly relate to patient's disease diagnosis and prognosis, but are difficult to measure or even not directly observable. The objective of this study was to develop, analyze and validate new priors for Bayesian inference of such biomarkers. Theoretical analysis revealed a relationship between the estimates inferred from the model and the true values of measured quantities, and the impact of the priors. This led to a new prior encoding scheme that incorporates objectively measurable domain knowledge, i.e. by performing two measurements with a reference method, which imply scale of the prior distribution. Second, priors on parameters of systematic error are non-informative, which enables biomarker estimation from a set of different quantities. Analysis showed that the volume of nucleus basalis of Meynert, which is reduced in early stages of Alzheimer's dementia and Parkinson's disease, is inter-related and could be inferred from compartmental brain volume measurements performed on routine clinical MR scans. Another experiment showed that total lesion load, associated to future disability progression in multiple sclerosis patients, could be inferred from lesion volume measurements based on multiple automated MR scan segmentations. Besides, figures of merit derived from the estimates could, without comparing against reference gold standard segmentations, identify the best performing lesion segmentation method. The proposed new priors substantially simplify the application of Bayesian inference for latent biomarkers and thus open an avenue for clinical implementation of new biomarkers, which may ultimately advance the evidence-based medicine.

Reference-free error estimation for multiple measurement methods.

We present a computational framework to select the most accurate and precise method of measurement of a certain quantity, when there is no access to the true value of the measurand. A typical use case is when several image analysis methods are applied to measure the value of a particular quantitative imaging biomarker from the same images. The accuracy of each measurement method is characterized by systematic error (bias), which is modeled as a polynomial in true values of measurand, and the precision as random error modeled with a Gaussian random variable. In contrast to previous works, the random errors are modeled jointly across all methods, thereby enabling the framework to analyze measurement methods based on similar principles, which may have correlated random errors. Furthermore, the posterior distribution of the error model parameters is estimated from samples obtained by Markov chain Monte-Carlo and analyzed to estimate the parameter values and the unknown true values of the measurand. The framework was validated on six synthetic and one clinical dataset containing measurements of total lesion load, a biomarker of neurodegenerative diseases, which was obtained with four automatic methods by analyzing brain magnetic resonance images. The estimates of bias and random error were in a good agreement with the corresponding least squares regression estimates against a reference.

A framework for automatic creation of gold-standard rigid 3D-2D registration datasets.

PURPOSE: Advanced image-guided medical procedures incorporate 2D intra-interventional information into pre-interventional 3D image and plan of the procedure through 3D/2D image registration (32R). To enter clinical use, and even for publication purposes, novel and existing 32R methods have to be rigorously validated. The performance of a 32R method can be estimated by comparing it to an accurate reference or gold standard method (usually based on fiducial markers) on the same set of images (gold standard dataset). Objective validation and comparison of methods are possible only if evaluation methodology is standardized, and the gold standard  dataset is made publicly available. Currently, very few such datasets exist and only one contains images of multiple patients acquired during a procedure. To encourage the creation of gold standard 32R datasets, we propose an automatic framework. METHODS: The framework is based on rigid registration of fiducial markers. The main novelty is spatial grouping of fiducial markers on the carrier device, which enables automatic marker localization and identification across the 3D and 2D images. RESULTS: The proposed framework was demonstrated on clinical angiograms of 20 patients. Rigid 32R computed by the framework was more accurate than that obtained manually, with the respective target registration error below 0.027 mm compared to 0.040 mm. CONCLUSION: The framework is applicable for gold standard setup on any rigid anatomy, provided that the acquired images contain spatially grouped fiducial markers. The gold standard datasets and software will be made publicly available.