DOCSLIB.ORG

Arxiv:2006.03086V2 [Quant-Ph] 12 Nov 2020

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Arxiv:2006.03086V2 [Quant-Ph] 12 Nov 2020 Augmented fidelities for single qubit gates Filip Wudarski,1, 2, ∗ Jeffrey Marshall,1, 2, y Andre Petukhov,1, 3 and Eleanor Rieffel1 1QuAIL, NASA Ames Research Center, Moffett Field, California 94035, USA 2USRA Research Institute for Advanced Computer Science, Mountain View, California 94043, USA 3Google Inc., Santa Barbara, California 93117, USAz An average gate fidelity is a standard performance metric to quantify deviation between an ideal unitary gate transformation and its realistic experimental implementation. The average is taken with respect to states uniformly distributed over the full Hilbert space. We analytically (single-qubit) and numerically (two-qubit) show how this average changes if the uniform distribution condition is relaxed, replaced by parametrized distributions – polar cap and von Mises-Fisher distributions – and how the resulting fidelities can differentiate certain noise models. In particular, we demonstrate that Pauli channels with different noise rates along the three axes can be faithfully distinguished using these augmented fidelities. I. INTRODUCTION less, they are all valuable benchmarking tools that allow researchers to capture and quantify some of the most rel- Impressive progress in quantum technologies has taken evant aspects of the behavior of quantum devices and quantum computing from a theoretical framework to an their building blocks - qubits and gates. experimental playground, where basic proof-of-principle Special attention should be given to the average fidelity concepts can be tested and verified. The most promi- - currently the figure of merit of performance metric. It nent recent example of the latter is the demonstration describes an error between an ideal and experimental re- of quantum advantage [1]: that even the imperfect cur- alization of a gate, and is averaged over all possible states rently available quantum hardware can perform tasks in- uniformly distributed in the Hilbert space (according to tractable for the most powerful supercomputers. Further the Haar measure). Unlike the diamond norm and the advances toward more capable and robust quantum hard- min fidelity, the average gate fidelity can be efficiently ware depend on gaining a better understanding of un- estimated through protocols like randomized benchmark- derlying physical effects, including the characterization ing (RB) [9–12], cross-entropy benchmark [1, 13] or di- of noise in actual quantum hardware. rect fidelity estimation [14–16]. However, being only a In theory, quantum process tomography (QPT) [2] can single parameter the metric cannot distinguish various be used to exhaustively benchmark a quantum device, noise models; it reports only a single element of the pro- identifying all of its imperfections. QPT reconstructs cess matrix, the χ element, which is associated with n n 00 the full process matrix χ (of size 2 × 2 , where n is the depolarizing rate. That is, it effectively identifies all the number of qubits), a matrix that encodes complete channels as depolarizing channels. information about underlying quantum transformation (including unwanted effects caused by noise). Unfortu- In this article, we propose to relax the uniform distri- nately, QPT scales exponentially with system size, be- bution condition and introduce augmented fidelity met- coming impractical for systems larger than a few qubits rics via parametrized distributions. In particular, we an- [3]. Intuitively it seems plausible that well-controlled sys- alytically investigate what information about noise pro- tems will have only few dominating error sources - i.e. the cesses can be extracted from average fidelity with respect χ matrix will be sparse up to some accuracy. Therefore, to a von Mises-Fisher distribution (a normal distribu- lower parameter approximations and associated metrics tion in directional statistics) and a polar cap distribu- and protocols that could assess the performance of quan- tion, i.e. a uniform distribution over a subset of states tum devices, and identify the crucial elements of χ, are parametrized by polar angle (colatitude). This approach arXiv:2006.03086v2 [quant-ph] 12 Nov 2020 promising approaches to noise characterization and quan- augments standard uniform-average fidelity by adding tum hardware benchmarking. Currently, the most com- extra tunable parameters to the metric, that depend on mon figures of merit considered are: diamond norm [4], distribution properties. This work provides a partial so- minimum fidelity [5] or average fidelity [6–8]. All three lution to the problem posed by Nielsen in [7] regarding techniques yield a single value that characterizes devia- gate fidelities over non-uniform distributions. Addition- tions from the ideal transformation. A parameter count ally, we derive the maximal spread in fidelity the value shows that all three methods provide only limited infor- (the difference between minimum and maximum attain- mation about the process matrix for a device. Neverthe- able values) and provide error bars (based on the stan- dard deviation derived from the considered distributions) for processes that share the same depolarizing rate. In particular, we show how to identify noise biases [17, 18] ∗ fi[email protected] in Pauli channels. Our analysis focuses mainly on sin- y jeff[email protected] gle qubit gates, where analytical formulas are derived, z Permanent address but we also open the discussion for similar approaches 2 for multi-qubit gates, in particular we numerically show A. Polar Cap Distribution local distribution for two-qubit gates and how they differ from the uniform-average fidelity. First let us define a single state gate fidelity as y Fj ih j(U; EU ) = Tr Uj ih jU EU (j ih j) ; (4) II. AVERAGE GATE FIDELITY where U is the investigated gate, EU its CPTP (imper- Consider the fidelity between a state transformed ac- fect) realization and j i is a state upon which the gate cording to a given unitary (gate) transformation (ideal acts. Since measuring and computing Eq. (4) for all pos- action) U and the same state transformed with noisy re- sible states is infeasible, one usually reports the average alization of U, which is a completely positive and trace fidelity value, which is taken over j i distributed uni- preserving (CPTP) map EU . The average fidelity of the formly in the entire Hilbert space. noisy realization is the average of this fidelity over all Following derivation from [6] we define the restricted initial (pure) states distributed uniformly in the entire ¯ ¯ Hilbert space. Bowdrey et al. [6] introduced a simple for- average gate fidelity FΘ(U; EU ) ≡ FΘ as mula for calculating average gate fidelity for single qubit Z ¯ gate, which was later generalized to multi-qubit gates and FΘ = Fj ih j(U; EU )dΩ = connected with entanglement fidelity [7,8]. The average gate fidelity (from now on referred as “uniform-average Z Θ Z 2π 3 1 X σj fidelity”) for n-qubit gates is therefore given by = Tr U c (θ; φ) U y × S(Θ) j 2 θ=0 φ=0 j=0 P22n−1 y y 2n Tr UV U EU (Vk) + 2 3 ¯ ¯ y k=0 k X σk F (U; EU ) = F (U ◦ EU ) = ; × E c (θ; φ) sin θdφdθ; (5) 22n(2n + 1) U k 2 (1) k=0 where Vk are traceless unitary matrices forming an or- where S(Θ) := 2π(1 − cos Θ) is the solid angle for thonormal basis with respect to Hilbert-Schmidt inner normalization of the distribution, and cj(θ; φ) are pure y n product (Tr(VkVj ) = 2 δkj, Tr(Vk) = 0 for k = state’s Bloch vector coefficients: c0(θ; φ) = 1; c1(θ; φ) = 2n 1;:::; 2 − 1 and V0 = 1). By writing the composed sin θ cos φ, c2(θ; φ) = sin θ sin φ and c3(θ; φ) = cos θ. map These correspond respectively to the identity matrix σ0, and the three x; y; z Pauli matrices σ1;2;3. Note that 22n−1 Eq. (4) assumes coordinate system where θ = 0 corre- y X y E(ρ) ≡ U ◦ EU (ρ) = χklVkρVl ; (2) sponds to j0i state, and the polar cap distribution is cen- k;l=0 tered around it. However, transformation to an arbitrary ~ central state is straightforward via rotation j i = UR j i. where the χ matrix is called the process matrix for E, it Now performing the integration leaves us with can be demonstrated that the average fidelity only de- 2 Θ 2 pends on the χ00 element corresponding to a “depolariz- 1 (2 + cos Θ) sin X F¯ = + 2 Tr Uσ U yE σ + ing” rate, a unitary invariant element, i.e. Θ 2 12 k U k k=1 2nχ + 1 1 + cos Θ + cos2 Θ F¯(U y ◦ E ) = 00 : (3) + Tr Uσ U yE σ + U 2n + 1 12 3 U 3 1 + cos Θ + Tr Uσ U yE σ : (6) The above formula demonstrates inability to distinguish 8 3 U 0 different noise processes that differ in other χ parameters, and this limitation stems from the averaging procedure It is transparent that for Θ = π one recovers result and properties of the Haar measure. In order to have a Eq. (1) for uniform distribution over the entire Hilbert more sensitive metric, we propose to use several differ- space. Expressing the composed gate-noisy-gate map y ent initial state distributions. In addition to averaging E = U ◦ EU in the form Eq. (2), one can show that over all pure states distributed uniformly, we explore two 2 models: i) uniform distribution parametrized by a po- 1 X Tr σkE σk = χ0;0 − χ3;3; (7) lar angle Θ 2 [0; π], which we call polar cap distribution 4 k=1 (e.g. for Θ = π=2 we have a distribution over the north- 1 ern hemisphere, while for Θ = π we recover the entire Tr σ Eσ = χ − χ − χ + χ ; (8) 2 3 3 0;0 1;1 2;2 3;3 space distribution), and ii) von Mises-Fisher distribution 1 around a state j i (without loss of generality we can fix Tr σ3E σ0 = Re(χ0;3); (9) it to j0i) parametrized by “variance” parameter κ.
Load more

Recommended publications
  • Contributions to Directional Statistics Based Clustering Methods
    University of California Santa Barbara Contributions to Directional Statistics Based Clustering Methods A dissertation submitted in partial satisfaction of the requirements for the degree Doctor of Philosophy in Statistics & Applied Probability by Brian Dru Wainwright Committee in charge: Professor S. Rao Jammalamadaka, Chair Professor Alexander Petersen Professor John Hsu September 2019 The Dissertation of Brian Dru Wainwright is approved. Professor Alexander Petersen Professor John Hsu Professor S. Rao Jammalamadaka, Committee Chair June 2019 Contributions to Directional Statistics Based Clustering Methods Copyright © 2019 by Brian Dru Wainwright iii Dedicated to my loving wife, Carmen Rhodes, without whom none of this would have been possible, and to my sons, Max, Gus, and Judah, who have yet to know a time when their dad didn’t need to work from early till late. And finally, to my mother, Judith Moyer, without your tireless love and support from the beginning, I quite literally wouldn’t be here today. iv Acknowledgements I would like to offer my humble and grateful acknowledgement to all of the wonderful col- laborators I have had the fortune to work with during my graduate education. Much of the impetus for the ideas presented in this dissertation were derived from our work together. In particular, I would like to thank Professor György Terdik, University of Debrecen, Faculty of Informatics, Department of Information Technology. I would also like to thank Professor Saumyadipta Pyne, School of Public Health, University of Pittsburgh, and Mr. Hasnat Ali, L.V. Prasad Eye Institute, Hyderabad, India. I would like to extend a special thank you to Dr Alexander Petersen, who has held the dual role of serving on my Doctoral supervisory commit- tee as well as wearing the hat of collaborator.
    [Show full text]
  • A Guide on Probability Distributions
    powered project A guide on probability distributions R-forge distributions Core Team University Year 2008-2009 LATEXpowered Mac OS' TeXShop edited Contents Introduction 4 I Discrete distributions 6 1 Classic discrete distribution 7 2 Not so-common discrete distribution 27 II Continuous distributions 34 3 Finite support distribution 35 4 The Gaussian family 47 5 Exponential distribution and its extensions 56 6 Chi-squared's ditribution and related extensions 75 7 Student and related distributions 84 8 Pareto family 88 9 Logistic ditribution and related extensions 108 10 Extrem Value Theory distributions 111 3 4 CONTENTS III Multivariate and generalized distributions 116 11 Generalization of common distributions 117 12 Multivariate distributions 132 13 Misc 134 Conclusion 135 Bibliography 135 A Mathematical tools 138 Introduction This guide is intended to provide a quite exhaustive (at least as I can) view on probability distri- butions. It is constructed in chapters of distribution family with a section for each distribution. Each section focuses on the tryptic: definition - estimation - application. Ultimate bibles for probability distributions are Wimmer & Altmann (1999) which lists 750 univariate discrete distributions and Johnson et al. (1994) which details continuous distributions. In the appendix, we recall the basics of probability distributions as well as \common" mathe- matical functions, cf. section A.2. And for all distribution, we use the following notations • X a random variable following a given distribution, • x a realization of this random variable, • f the density function (if it exists), • F the (cumulative) distribution function, • P (X = k) the mass probability function in k, • M the moment generating function (if it exists), • G the probability generating function (if it exists), • φ the characteristic function (if it exists), Finally all graphics are done the open source statistical software R and its numerous packages available on the Comprehensive R Archive Network (CRAN∗).
    [Show full text]
  • Recent Advances in Directional Statistics
    Recent advances in directional statistics Arthur Pewsey1;3 and Eduardo García-Portugués2 Abstract Mainstream statistical methodology is generally applicable to data observed in Euclidean space. There are, however, numerous contexts of considerable scientific interest in which the natural supports for the data under consideration are Riemannian manifolds like the unit circle, torus, sphere and their extensions. Typically, such data can be represented using one or more directions, and directional statistics is the branch of statistics that deals with their analysis. In this paper we provide a review of the many recent developments in the field since the publication of Mardia and Jupp (1999), still the most comprehensive text on directional statistics. Many of those developments have been stimulated by interesting applications in fields as diverse as astronomy, medicine, genetics, neurology, aeronautics, acoustics, image analysis, text mining, environmetrics, and machine learning. We begin by considering developments for the exploratory analysis of directional data before progressing to distributional models, general approaches to inference, hypothesis testing, regression, nonparametric curve estimation, methods for dimension reduction, classification and clustering, and the modelling of time series, spatial and spatio- temporal data. An overview of currently available software for analysing directional data is also provided, and potential future developments discussed. Keywords: Classification; Clustering; Dimension reduction; Distributional
    [Show full text]
  • A General Approach for Obtaining Wrapped Circular Distributions Via Mixtures
    UC Santa Barbara UC Santa Barbara Previously Published Works Title A General Approach for obtaining wrapped circular distributions via mixtures Permalink https://escholarship.org/uc/item/5r7521z5 Journal Sankhya, 79 Author Jammalamadaka, Sreenivasa Rao Publication Date 2017 Peer reviewed eScholarship.org Powered by the California Digital Library University of California 1 23 Your article is protected by copyright and all rights are held exclusively by Indian Statistical Institute. This e-offprint is for personal use only and shall not be self- archived in electronic repositories. If you wish to self-archive your article, please use the accepted manuscript version for posting on your own website. You may further deposit the accepted manuscript version in any repository, provided it is only made publicly available 12 months after official publication or later and provided acknowledgement is given to the original source of publication and a link is inserted to the published article on Springer's website. The link must be accompanied by the following text: "The final publication is available at link.springer.com”. 1 23 Author's personal copy Sankhy¯a:TheIndianJournalofStatistics 2017, Volume 79-A, Part 1, pp. 133-157 c 2017, Indian Statistical Institute ! AGeneralApproachforObtainingWrappedCircular Distributions via Mixtures S. Rao Jammalamadaka University of California, Santa Barbara, USA Tomasz J. Kozubowski University of Nevada, Reno, USA Abstract We show that the operations of mixing and wrapping linear distributions around a unit circle commute, and can produce a wide variety of circular models. In particular, we show that many wrapped circular models studied in the literature can be obtained as scale mixtures of just the wrapped Gaussian and the wrapped exponential distributions, and inherit many properties from these two basic models.
    [Show full text]
  • Multivariate Statistical Functions in R
    Multivariate statistical functions in R Michail T. Tsagris [email protected] College of engineering and technology, American university of the middle east, Egaila, Kuwait Version 6.1 Athens, Nottingham and Abu Halifa (Kuwait) 31 October 2014 Contents 1 Mean vectors 1 1.1 Hotelling’s one-sample T2 test ............................. 1 1.2 Hotelling’s two-sample T2 test ............................ 2 1.3 Two two-sample tests without assuming equality of the covariance matrices . 4 1.4 MANOVA without assuming equality of the covariance matrices . 6 2 Covariance matrices 9 2.1 One sample covariance test .............................. 9 2.2 Multi-sample covariance matrices .......................... 10 2.2.1 Log-likelihood ratio test ............................ 10 2.2.2 Box’s M test ................................... 11 3 Regression, correlation and discriminant analysis 13 3.1 Correlation ........................................ 13 3.1.1 Correlation coefficient confidence intervals and hypothesis testing us- ing Fisher’s transformation .......................... 13 3.1.2 Non-parametric bootstrap hypothesis testing for a zero correlation co- efficient ..................................... 14 3.1.3 Hypothesis testing for two correlation coefficients . 15 3.2 Regression ........................................ 15 3.2.1 Classical multivariate regression ....................... 15 3.2.2 k-NN regression ................................ 17 3.2.3 Kernel regression ................................ 20 3.2.4 Choosing the bandwidth in kernel regression in a very simple way . 23 3.2.5 Principal components regression ....................... 24 3.2.6 Choosing the number of components in principal component regression 26 3.2.7 The spatial median and spatial median regression . 27 3.2.8 Multivariate ridge regression ......................... 29 3.3 Discriminant analysis .................................. 31 3.3.1 Fisher’s linear discriminant function ....................
    [Show full text]
  • Handbook on Probability Distributions
    R powered R-forge project Handbook on probability distributions R-forge distributions Core Team University Year 2009-2010 LATEXpowered Mac OS' TeXShop edited Contents Introduction 4 I Discrete distributions 6 1 Classic discrete distribution 7 2 Not so-common discrete distribution 27 II Continuous distributions 34 3 Finite support distribution 35 4 The Gaussian family 47 5 Exponential distribution and its extensions 56 6 Chi-squared's ditribution and related extensions 75 7 Student and related distributions 84 8 Pareto family 88 9 Logistic distribution and related extensions 108 10 Extrem Value Theory distributions 111 3 4 CONTENTS III Multivariate and generalized distributions 116 11 Generalization of common distributions 117 12 Multivariate distributions 133 13 Misc 135 Conclusion 137 Bibliography 137 A Mathematical tools 141 Introduction This guide is intended to provide a quite exhaustive (at least as I can) view on probability distri- butions. It is constructed in chapters of distribution family with a section for each distribution. Each section focuses on the tryptic: definition - estimation - application. Ultimate bibles for probability distributions are Wimmer & Altmann (1999) which lists 750 univariate discrete distributions and Johnson et al. (1994) which details continuous distributions. In the appendix, we recall the basics of probability distributions as well as \common" mathe- matical functions, cf. section A.2. And for all distribution, we use the following notations • X a random variable following a given distribution, • x a realization of this random variable, • f the density function (if it exists), • F the (cumulative) distribution function, • P (X = k) the mass probability function in k, • M the moment generating function (if it exists), • G the probability generating function (if it exists), • φ the characteristic function (if it exists), Finally all graphics are done the open source statistical software R and its numerous packages available on the Comprehensive R Archive Network (CRAN∗).
    [Show full text]
  • Bivariate Angular Estimation Under Consideration of Dependencies Using Directional Statistics
    Bivariate Angular Estimation Under Consideration of Dependencies Using Directional Statistics Gerhard Kurz1, Igor Gilitschenski1, Maxim Dolgov1 and Uwe D. Hanebeck1 Abstract— Estimation of angular quantities is a widespread issue, but standard approaches neglect the true topology of the problem and approximate directional with linear uncertainties. In recent years, novel approaches based on directional statistics have been proposed. However, these approaches have been unable to consider arbitrary circular correlations between multiple angles so far. For this reason, we propose a novel recursive filtering scheme that is capable of estimating multiple angles even if they are dependent, while correctly describing their circular correlation. The proposed approach is based on toroidal probability distributions and a circular correlation coefficient. We demonstrate the superiority to a standard approach based on the Kalman filter in simulations. Fig. 1: A bivariate wrapped normal probability distribution Index Terms— recursive filtering, wrapped normal, circular on the torus shown as a heat map. correlation coefficient, moment matching. I. INTRODUCTION our knowledge, this is the first work on recursive estimation There are many applications that require estimation of based on toroidal probability distributions. angular quantities. These applications include, but are not It should be noted that there are some directional filters that, limited to, robotics, augmented reality, and aviation, as well in a sense, take correlation between angles into account. We as biology, geology, and medicine. In many cases, not just have performed research on a filter based on the hyperspheri- one, but several angles have to be estimated. Furthermore, cal Bingham distribution [8], [9], which captures correlations correlations may exist between those angles and have to be when estimating rotations represented as quaternions.
    [Show full text]
  • A New Unified Approach for the Simulation of a Wide Class of Directional Distributions
    This is a repository copy of A New Unified Approach for the Simulation of a Wide Class of Directional Distributions. White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/123206/ Version: Accepted Version Article: Kent, JT orcid.org/0000-0002-1861-8349, Ganeiber, AM and Mardia, KV orcid.org/0000-0003-0090-6235 (2018) A New Unified Approach for the Simulation of a Wide Class of Directional Distributions. Journal of Computational and Graphical Statistics, 27 (2). pp. 291-301. ISSN 1061-8600 https://doi.org/10.1080/10618600.2017.1390468 © 2018 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America. This is an author produced version of a paper published in Journal of Computational and Graphical Statistics. Uploaded in accordance with the publisher's self-archiving policy. Reuse Items deposited in White Rose Research Online are protected by copyright, with all rights reserved unless indicated otherwise. They may be downloaded and/or printed for private study, or other acts as permitted by national copyright laws. The publisher or other rights holders may allow further reproduction and re-use of the full text version. This is indicated by the licence information on the White Rose Research Online record for the item. Takedown If you consider content in White Rose Research Online to be in breach of UK law, please notify us by emailing [email protected] including the URL of the record and the reason for the withdrawal request. [email protected] https://eprints.whiterose.ac.uk/ A new unified approach for the simulation of a wide class of directional distributions John T.
    [Show full text]
  • A Logistic Normal Mixture Model for Compositions with Essential Zeros
    A Logistic Normal Mixture Model for Compositions with Essential Zeros Item Type text; Electronic Dissertation Authors Bear, John Stanley Publisher The University of Arizona. Rights Copyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction or presentation (such as public display or performance) of protected items is prohibited except with permission of the author. Download date 28/09/2021 15:07:32 Link to Item http://hdl.handle.net/10150/621758 A Logistic Normal Mixture Model for Compositions with Essential Zeros by John S. Bear A Dissertation Submitted to the Faculty of the Graduate Interdisciplinary Program in Statistics In Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy In the Graduate College The University of Arizona 2 0 1 6 2 THE UNIVERSITY OF ARIZONA GRADUATE COLLEGE As members of the Dissertation Committee, we certify that we have read the disserta- tion prepared by John Bear, titled A Logistic Normal Mixture Model for Compositions with Essential Zeros and recommend that it be accepted as fulfilling the dissertation requirement for the Degree of Doctor of Philosophy. Date: August 18, 2016 Dean Billheimer Date: August 18, 2016 Joe Watkins Date: August 18, 2016 Bonnie LaFleur Date: August 18, 2016 Keisuke Hirano Final approval and acceptance of this dissertation is contingent upon the candidate's submission of the final copies of the dissertation to the Graduate College. I hereby certify that I have read this dissertation prepared under my direction and recommend that it be accepted as fulfilling the dissertation requirement.
    [Show full text]
  • High-Precision Computation of Uniform Asymptotic Expansions for Special Functions Guillermo Navas-Palencia
    UNIVERSITAT POLITÈCNICA DE CATALUNYA Department of Computer Science High-precision Computation of Uniform Asymptotic Expansions for Special Functions Guillermo Navas-Palencia Supervisor: Argimiro Arratia Quesada A dissertation submitted in fulfillment of the requirements for the degree of Doctor of Philosophy in Computing May 7, 2019 i Abstract In this dissertation, we investigate new methods to obtain uniform asymptotic ex- pansions for the numerical evaluation of special functions to high-precision. We shall first present the theoretical and computational fundamental aspects required for the development and ultimately implementation of such methods. Applying some of these methods, we obtain efficient new convergent and uniform expan- sions for numerically evaluating the confluent hypergeometric functions 1F1(a; b; z) and U(a, b, z), and the Lerch transcendent F(z, s, a) at high-precision. In addition, we also investigate a new scheme of computation for the generalized exponential integral En(z), obtaining one of the fastest and most robust implementations in double-precision arithmetic. In this work, we aim to combine new developments in asymptotic analysis with fast and effective open-source implementations. These implementations are com- parable and often faster than current open-source and commercial state-of-the-art software for the evaluation of special functions. ii Acknowledgements First, I would like to express my gratitude to my supervisor Argimiro Arratia for his support, encourage and guidance throughout this work, and for letting me choose my research path with full freedom. I also thank his assistance with admin- istrative matters, especially in periods abroad. I am grateful to Javier Segura and Amparo Gil from Universidad de Cantabria for inviting me for a research stay and for their inspirational work in special func- tions, and the Ministerio de Economía, Industria y Competitividad for the financial support, project APCOM (TIN2014-57226-P), during the stay.
    [Show full text]
  • Package 'Directional'
    Package ‘Directional’ May 26, 2021 Type Package Title A Collection of R Functions for Directional Data Analysis Version 5.0 URL Date 2021-05-26 Author Michail Tsagris, Giorgos Athineou, Anamul Sajib, Eli Amson, Micah J. Waldstein Maintainer Michail Tsagris <[email protected]> Description A collection of functions for directional data (including massive data, with millions of ob- servations) analysis. Hypothesis testing, discriminant and regression analysis, MLE of distribu- tions and more are included. The standard textbook for such data is the ``Directional Statis- tics'' by Mardia, K. V. and Jupp, P. E. (2000). Other references include a) Phillip J. Paine, Si- mon P. Preston Michail Tsagris and Andrew T. A. Wood (2018). An elliptically symmetric angu- lar Gaussian distribution. Statistics and Computing 28(3): 689-697. <doi:10.1007/s11222-017- 9756-4>. b) Tsagris M. and Alenazi A. (2019). Comparison of discriminant analysis meth- ods on the sphere. Communications in Statistics: Case Studies, Data Analysis and Applica- tions 5(4):467--491. <doi:10.1080/23737484.2019.1684854>. c) P. J. Paine, S. P. Pre- ston, M. Tsagris and Andrew T. A. Wood (2020). Spherical regression models with general co- variates and anisotropic errors. Statistics and Computing 30(1): 153--165. <doi:10.1007/s11222- 019-09872-2>. License GPL-2 Imports bigstatsr, parallel, doParallel, foreach, RANN, Rfast, Rfast2, rgl RoxygenNote 6.1.1 NeedsCompilation no Repository CRAN Date/Publication 2021-05-26 15:20:02 UTC R topics documented: Directional-package . .4 A test for testing the equality of the concentration parameters for ciruclar data . .5 1 2 R topics documented: Angular central Gaussian random values simulation .
    [Show full text]
  • Spherical-Homoscedastic Distributions: the Equivalency of Spherical and Normal Distributions in Classification
    Journal of Machine Learning Research 8 (2007) 1583-1623 Submitted 4/06; Revised 1/07; Published 7/07 Spherical-Homoscedastic Distributions: The Equivalency of Spherical and Normal Distributions in Classification Onur C. Hamsici [email protected] Aleix M. Martinez [email protected] Department of Electrical and Computer Engineering The Ohio State University Columbus, OH 43210, USA Editor: Greg Ridgeway Abstract Many feature representations, as in genomics, describe directional data where all feature vectors share a common norm. In other cases, as in computer vision, a norm or variance normalization step, where all feature vectors are normalized to a common length, is generally used. These repre- sentations and pre-processing step map the original data from Rp to the surface of a hypersphere p 1 S − . Such representations should then be modeled using spherical distributions. However, the difficulty associated with such spherical representations has prompted researchers to model their spherical data using Gaussian distributions instead—as if the data were represented in Rp rather p 1 than S − . This opens the question to whether the classification results calculated with the Gaus- sian approximation are the same as those obtained when using the original spherical distributions. In this paper, we show that in some particular cases (which we named spherical-homoscedastic) the answer to this question is positive. In the more general case however, the answer is negative. For this reason, we further investigate the additional error added by the Gaussian modeling. We conclude that the more the data deviates from spherical-homoscedastic, the less advisable it is to employ the Gaussian approximation.
    [Show full text]