Model Uncertainty in Panel Vector Autoregressive Models.

We develop methods for Bayesian model averaging (BMA) or selection (BMS) in Panel Vector Autoregressions (PVARs). Our approach allows us to select between or average over all possible combinations of restricted PVARs where the restrictions involve interdependencies between and heterogeneities across cross-sectional units. The resulting BMA framework can find a parsimonious PVAR specification, thus dealing with overparameterization concerns. We use these methods in an application involving the euro area sovereign debt crisis and show that our methods perform better than alternatives. Our findings contradict a simple view of the sovereign debt crisis which divides the euro zone into groups of core and peripheral countries and worries about financial contagion within the latter group.

Large Bayesian VARMAs

Vector Autoregressive Moving Average (VARMA) models have many theoretical properties which should make them popular among empirical macroeconomists. However, they are rarely used in practice due to over-parameterization concerns, difficulties in ensuring identification and computational challenges. With the growing interest in multivariate time series models of high dimension, these problems with VARMAs become even more acute, accounting for the dominance of VARs in this field. In this paper, we develop a Bayesian approach for inference in VARMAs which surmounts these problems. It jointly ensures identification and parsimony in the context of an efficient Markov chain Monte Carlo (MCMC) algorithm. We use this approach in a macroeconomic application involving up to twelve dependent variables. We find our algorithm to work successfully and provide insights beyond those provided by VARs.

Time Varying Dimension Models

Time varying parameter (TVP) models have enjoyed an increasing popularity in empirical macroeconomics. However, TVP models are parameter-rich and risk over-fitting unless the dimension of the model is small. Motivated by this worry, this paper proposes several Time Varying dimension (TVD) models where the dimension of the model can change over time, allowing for the model to automatically choose a more parsimonious TVP representation, or to switch between different parsimonious representations. Our TVD models all fall in the category of dynamic mixture models. We discuss the properties of these models and present methods for Bayesian inference. An application involving US inflation forecasting illustrates and compares the different TVD models. We find our TVD approaches exhibit better forecasting performance than several standard benchmarks and shrink towards parsimonious specifications.

Hierarchical Shrinkage in Time-Varying Parameter Models

In this paper, we forecast EU-area inflation with many predictors using time-varying parameter models. The facts that time-varying parameter models are parameter-rich and the time span of our data is relatively short motivate a desire for shrinkage. In constant coefficient regression models, the Bayesian Lasso is gaining increasing popularity as an effective tool for achieving such shrinkage. In this paper, we develop econometric methods for using the Bayesian Lasso with time-varying parameter models. Our approach allows for the coefficient on each predictor to be: i) time varying, ii) constant over time or iii) shrunk to zero. The econometric methodology decides automatically which category each coefficient belongs in. Our empirical results indicate the benefits of such an approach.

Quantile forecasts of inflation under model uncertainty

Bayesian model averaging (BMA) methods are regularly used to deal with model uncertainty in regression models. This paper shows how to introduce Bayesian model averaging methods in quantile regressions, and allow for different predictors to affect different quantiles of the dependent variable. I show that quantile regression BMA methods can help reduce uncertainty regarding outcomes of future inflation by providing superior predictive densities compared to mean regression models with and without BMA.

A Modulated Closed Form solution for Quantitative Susceptibility Mapping - A thorough evaluation and comparison to iterative methods based on edge prior knowledge.

The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one over-determined multiple orientation method (COSMOS) and anewly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65mm isotropic) brain data acquired at 7T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corrMCF=0.95, r(2)MCF =0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The l-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures.

Quantitative comparison of reconstruction methods for intra-voxel fiber recovery from diffusion MRI.

Validation is arguably the bottleneck in the diffusion magnetic resonance imaging (MRI) community. This paper evaluates and compares 20 algorithms for recovering the local intra-voxel fiber structure from diffusion MRI data and is based on the results of the "HARDI reconstruction challenge" organized in the context of the "ISBI 2012" conference. Evaluated methods encompass a mixture of classical techniques well known in the literature such as diffusion tensor, Q-Ball and diffusion spectrum imaging, algorithms inspired by the recent theory of compressed sensing and also brand new approaches proposed for the first time at this contest. To quantitatively compare the methods under controlled conditions, two datasets with known ground-truth were synthetically generated and two main criteria were used to evaluate the quality of the reconstructions in every voxel: correct assessment of the number of fiber populations and angular accuracy in their orientation. This comparative study investigates the behavior of every algorithm with varying experimental conditions and highlights strengths and weaknesses of each approach. This information can be useful not only for enhancing current algorithms and develop the next generation of reconstruction methods, but also to assist physicians in the choice of the most adequate technique for their studies.

Non-immunologic methods of diagnosis of babesiosis

The diagnosis of tick-borne diseases such as babesiosis still depends on observing the parasite in the infected erythrocyte. Microscopic observation is tedious and often problematic in both early and carrier infections. Better diagnostic methods are needed to prevent clinical disease, especially when susceptible cattle are being moved into disease enzootic areas. This study evaluates two techniques for early diagnosis of Babesia bovis infections in cattle, DNA probes specific for the organism and fluorescent probes specific nucleic acid. The radioisotopically labeled DNA probes are used in slot blot hybridizations whith lysed blood samples, not purified DNA. Thusfar, the probe is specific for B. bovis and can detect as few as 1000 B. bovis parasites in 10µl of blood. The specificity of the fluorescent probe depends on the characteristic morphology of the babesia in whole blood samples, as determined microscopically. The fluorescent probe detects as afew as 10,000 B. bovis parasites in 10*l as blood. The application of each method for alboratory and field use is discussed.

New methods for the diagnosis of Babesia bigemina infection

Accurate diagnosis of Babesia bigemina infection, an economically important tick-transmitted protozoan parasite of cattle, is essential in the management of disease control and in epidemiological studies. The currentlyused methods of diagnosis are blood smear examination and serological tests which include agglutination and immunofluorescence tests. These testes have been used the fild but because they lack sensitivity and specificity, never and improved methods of diagnosis are being developed. The quantitative buffy coat (OBC) method, using microhaematocrit tubes and acridine orange staining allows rapid and quicker diagnosis of B. bigemina and other blood parasites compared to light microscopic examination of stained smears. Parasite specific monoclonal antibodies have been used in antigen/antibody capture enzymelinked immunosorbent assays with grater sensitivity and specificity than previously described serological tests. Similary, DNA probes, derived from a repetitive sequence of the B. bigemina genome, offer a method of detecting very small numbers of parasites which are undetectable by conventional microscopy. An extrachromosomal DNA element, present in all the tick-borne protozoan parasites so far tested, provides an accurate means of diferentiating mixed parasite populations in infected animals. These improved methods will greatly facilitate epidemiological studies.

Endogenous steroid profiling in the athlete biological passport.

The Athlete Biological Passport (ABP) is an individual electronic document that collects data regarding a specific athlete that is useful in differentiating between natural physiologic variations of selected biomarkers and deviations caused by artificial manipulations. A subsidiary of the endocrine module of the ABP, that which here is called Athlete Steroidal Passport (ASP), collects data on markers of an altered metabolism of endogenous steroidal hormones measured in urine samples. The ASP aims to identify not only doping with anabolic-androgenic steroids, but also most indirect steroid doping strategies such as doping with estrogen receptor antagonists and aromatase inhibitors. Development of specific markers of steroid doping, use of the athlete's previous measurements to define individual limits, with the athlete becoming his or her own reference, the inclusion of heterogeneous factors such as the UDPglucuronosyltransferase B17 genotype of the athlete, the knowledge of potentially confounding effects such as heavy alcohol consumption, the development of an external quality control system to control analytical uncertainty, and finally the use of Bayesian inferential methods to evaluate the value of indirect evidence have made the ASP a valuable alternative to deter steroid doping in elite sports. The ASP can be used to target athletes for gas chromatography/combustion/ isotope ratio mass spectrometry (GC/C/IRMS) testing, to withdraw temporarily the athlete from competing when an abnormality has been detected, and ultimately to lead to an antidoping infraction if that abnormality cannot be explained by a medical condition. Although the ASP has been developed primarily to ensure fairness in elite sports, its application in endocrinology for clinical purposes is straightforward in an evidence-based medicine paradigm.

Pharmacokinetic (PK)-based dosage individualization of imatinib: evaluation of efficiency and clinical usefulness

Background: Retrospective analyses suggest that personalized PK-based dosage might be useful for imatinib, as treatment response correlates with trough concentrations (Cmin) in cancer patients. Our objectives were to improve the interpretation of randomly measured concentrations and to confirm its efficiency before evaluating the clinical usefulness of systematic PK-based dosage in chronic myeloid leukemia patients. Methods and Results: A Bayesian method was validated for the prediction of individual Cmin on the basis of a single random observation, and was applied in a prospective multicenter randomized controlled clinical trial. 28 out of 56 patients were enrolled in the systematic dosage individualization arm and had 44 follow-up visits (their clinical follow-up is ongoing). PK-dose-adjustments were proposed in 39% having predicted Cmin significantly away from the target (1000 ng/ml). Recommendations were taken up by physicians in 57%, patients were considered non-compliant in 27%. Median Cmin at study inclusion was 754 ng/ml and differed significantly from the target (p=0.02, Wilcoxon test). On follow-up, Cmin was 984 ng/ml (p=0.82) in the compliant group. CV decreased from 46% to 27% (p=0.02, F-test). Conclusion: PK-based (Bayesian) dosage adjustment is able to bring individual drug exposure closer to a given therapeutic target. Its influence on therapeutic response remains to be evaluated.

Discontinuous Galerkin methods for the one-dimensional Vlasov-Poisson system

We construct a new family of semi-discrete numerical schemes for the approximation of the one-dimensional periodic Vlasov-Poisson system. The methods are based on the coupling of discontinuous Galerkin approximation to the Vlasov equation and several finite element (conforming, non-conforming and mixed) approximations for the Poisson problem. We show optimal error estimates for the all proposed methods in the case of smooth compactly supported initial data. The issue of energy conservation is also analyzed for some of the methods.