Statistical Analysis of Functions on Surfaces, With an Application to Medical Imaging

In functional data analysis, data are commonly assumed to be smooth functions on a fixed interval of the real line. In this work, we introduce a comprehensive framework for the analysis of functional data, whose domain is a two-dimensional manifold and the domain itself is subject to variability from sample to sample. We formulate a statistical model for such data, here called functions on surfaces, which enables a joint representation of the geometric and functional aspects, and propose an associated estimation framework. We assess the validity of the framework by performing a simulation study and we finally apply it to the analysis of neuroimaging data of cortical thickness, acquired from the brains of different subjects, and thus lying on domains with different geometries. Supplementary materials for this article are available online

Representation and reconstruction of covariance operators in linear inverse problems

Abstract We introduce a framework for the reconstruction and representation of functions in a setting where these objects cannot be directly observed, but only indirect and noisy measurements are available, namely an inverse problem setting. The proposed methodology can be applied either to the analysis of indirectly observed functional images or to the associated covariance operators, representing second-order information, and thus lying on a non-Euclidean space. To deal with the ill-posedness of the inverse problem, we exploit the spatial structure of the sample data by introducing a flexible regularizing term embedded in the model. Thanks to its efficiency, the proposed model is applied to MEG data, leading to a novel approach to the investigation of functional connectivity.</jats:p

A Functional Approach to Deconvolve Dynamic Neuroimaging Data.

Positron emission tomography (PET) is an imaging technique which can be used to investigate chemical changes in human biological processes such as cancer development or neurochemical reactions. Most dynamic PET scans are currently analyzed based on the assumption that linear first-order kinetics can be used to adequately describe the system under observation. However, there has recently been strong evidence that this is not the case. To provide an analysis of PET data which is free from this compartmental assumption, we propose a nonparametric deconvolution and analysis model for dynamic PET data based on functional principal component analysis. This yields flexibility in the possible deconvolved functions while still performing well when a linear compartmental model setup is the true data generating mechanism. As the deconvolution needs to be performed on only a relative small number of basis functions rather than voxel by voxel in the entire three-dimensional volume, the methodology is both robust to typical brain imaging noise levels while also being computationally efficient. The new methodology is investigated through simulations in both one-dimensional functions and 2D images and also applied to a neuroimaging study whose goal is the quantification of opioid receptor concentration in the brain.The research of Ci-Ren Jiang is supported in part by NSC 101-2118-M-001-013-MY2 (Taiwan); the research of Jane-Ling Wang is supported by NSF grants, DMS-09-06813 and DMS-12-28369. JA is supported by EPSRC grant EP/K021672/2. The authors would like to thank SAMSI and the NDA programme where some of this research was carried out.This is the final version of the article. It first appeared from Taylor & Francis via http://dx.doi.org/10.1080/01621459.2015.106024

An Introduction to Applications of Wavelet Benchmarking with Seasonal Adjustment

Summary Before adjustment, low and high frequency data sets from national accounts are frequently inconsistent. Benchmarking is the procedure used by economic agencies to make such data sets consistent. It typically involves adjusting the high frequency time series (e.g. quarterly data) so that they become consistent with the lower frequency version (e.g. annual data). Various methods have been developed to approach this problem of inconsistency between data sets. The paper introduces a new statistical procedure, namely wavelet benchmarking. Wavelet properties allow high and low frequency processes to be jointly analysed and we show that benchmarking can be formulated and approached succinctly in the wavelet domain. Furthermore the time and frequency localization properties of wavelets are ideal for handling more complicated benchmarking problems. The versatility of the procedure is demonstrated by using simulation studies where we provide evidence showing that it substantially outperforms currently used methods. Finally, we apply this novel method of wavelet benchmarking to official data from the UK's Office for National Statistics.Engineering and Physical Sciences Research CouncilThis is the final version of the article. It first appeared from Wiley via https://doi.org/10.1111/rssa.1224

Smooth Principal Component Analysis over two-dimensional manifolds with an application to neuroimaging

Motivated by the analysis of high-dimensional neuroimaging signals located over the cortical surface, we introduce a novel Principal Component Analysis technique that can handle functional data located over a two-dimensional manifold. For this purpose a regularization approach is adopted, introducing a smoothing penalty coherent with the geodesic distance over the manifold. The model introduced can be applied to any manifold topology, can naturally handle missing data and functional samples evaluated in different grids of points. We approach the discretization task by means of finite element analysis and propose an efficient iterative algorithm for its resolution. We compare the performances of the proposed algorithm with other approaches classically adopted in literature. We finally apply the proposed method to resting state functional magnetic resonance imaging data from the Human Connectome Project, where the method shows substantial differential variations between brain regions that were not apparent with other approaches.Engineering and Physical Sciences Research Council (Grants EP/K021672/2, EP/N014588/1)This is the author accepted manuscript. The final version is available from the Institute of Mathematical Statistics via https://doi.org/10.1214/16-AOAS97

Stress-Induced Reinstatement of Drug Seeking: 20 Years of Progress

In human addicts, drug relapse and craving are often provoked by stress. Since 1995, this clinical scenario has been studied using a rat model of stress-induced reinstatement of drug seeking. Here, we first discuss the generality of stress-induced reinstatement to different drugs of abuse, different stressors, and different behavioral procedures. We also discuss neuropharmacological mechanisms, and brain areas and circuits controlling stress-induced reinstatement of drug seeking. We conclude by discussing results from translational human laboratory studies and clinical trials that were inspired by results from rat studies on stress-induced reinstatement. Our main conclusions are (1) The phenomenon of stress-induced reinstatement, first shown with an intermittent footshock stressor in rats trained to self-administer heroin, generalizes to other abused drugs, including cocaine, methamphetamine, nicotine, and alcohol, and is also observed in the conditioned place preference model in rats and mice. This phenomenon, however, is stressor specific and not all stressors induce reinstatement of drug seeking. (2) Neuropharmacological studies indicate the involvement of corticotropin-releasing factor (CRF), noradrenaline, dopamine, glutamate, kappa/dynorphin, and several other peptide and neurotransmitter systems in stress-induced reinstatement. Neuropharmacology and circuitry studies indicate the involvement of CRF and noradrenaline transmission in bed nucleus of stria terminalis and central amygdala, and dopamine, CRF, kappa/dynorphin, and glutamate transmission in other components of the mesocorticolimbic dopamine system (ventral tegmental area, medial prefrontal cortex, orbitofrontal cortex, and nucleus accumbens). (3) Translational human laboratory studies and a recent clinical trial study show the efficacy of alpha-2 adrenoceptor agonists in decreasing stress-induced drug craving and stress-induced initial heroin lapse

Search for rare quark-annihilation decays, B --> Ds(*) Phi

We report on searches for B- --> Ds- Phi and B- --> Ds*- Phi. In the context of the Standard Model, these decays are expected to be highly suppressed since they proceed through annihilation of the b and u-bar quarks in the B- meson. Our results are based on 234 million Upsilon(4S) --> B Bbar decays collected with the BABAR detector at SLAC. We find no evidence for these decays, and we set Bayesian 90% confidence level upper limits on the branching fractions BF(B- --> Ds- Phi) Ds*- Phi)<1.2x10^(-5). These results are consistent with Standard Model expectations.Comment: 8 pages, 3 postscript figues, submitted to Phys. Rev. D (Rapid Communications

Measurement of the branching fraction for B−→D0K∗−B^- \to D^0 K^{*-}

We present a measurement of the branching fraction for the decay B- --> D0 K*- using a sample of approximately 86 million BBbar pairs collected by the BaBar detector from e+e- collisions near the Y(4S) resonance. The D0 is detected through its decays to K- pi+, K- pi+ pi0 and K- pi+ pi- pi+, and the K*- through its decay to K0S pi-. We measure the branching fraction to be B.F.(B- --> D0 K*-)= (6.3 +/- 0.7(stat.) +/- 0.5(syst.)) x 10^{-4}

Observation of a significant excess of π0π0\pi^{0}\pi^{0} events in B meson decays

We present an observation of the decay $B^{0} \to \pi^{0} \pi^{0}$ based on a sample of 124 million $B\bar{B}$ pairs recorded by the BABAR detector at the PEP-II asymmetric-energy $B$ Factory at SLAC. We observe $46 \pm 13 \pm 3$ events, where the first error is statistical and the second is systematic, corresponding to a significance of 4.2 standard deviations including systematic uncertainties. We measure the branching fraction \BR(B^{0} \to \pi^{0} \pi^{0}) = (2.1 \pm 0.6 \pm 0.3) \times 10^{-6}, averaged over $B^{0}$ and $\bar{B}^{0}$ decays