Holonomic gradient method for distribution function of a weighted sum of noncentral chi-square random variables
|
|
- Edwina Foster
- 6 years ago
- Views:
Transcription
1 Holonomic gradient method for distribution function of a weighted sum of noncentral chi-square random variables T.Koyama & A.Takemura March 3,
2 The Ball Probability (1). Let X be a d-dimensional normal random vector with the mean vector µ and the covariance matrix Σ: X N(µ, Σ) Then the probability P (X X2 d r2 ) of the ball with radius r is written as 1 x x2 d r2 (2π) d/2 exp Σ 1/2 ( 1 ) 2 (x µ) Σ 1 (x µ) dx. 2
3 The Ball Probability (2). By rotation, we can assume that Σ is a diagonal matrix without loss of generality. Hence X 2 = X X2 d is the weighted sum of independent noncentral chi-square random variables. The integral for the infinitesimal interval is the Fisher-Bingham integral. r < x x2 d < r + dr In statistical interpretation, the conditional distribution of X given X = r is the Fisher-Bingham distribution. 3
4 Fisher Bingham Integral. The Fisher Bingham Integral is an integral on the sphere S d 1 (r) = { t R d t t2 d = r2} with parameters λ 1,..., λ d, τ 1,..., τ d defined by f(λ, τ, r) = d exp S d 1 (r) i=1 λ i t 2 i + Here, dt is the volume element on S d 1 (r) with d i=1 τ i t i dt. Sd 1(r) dt = rd 1 S d 1, S d 1 = Vol(S d 1 (1)) = 2πd/2 Γ(d/2). 4
5 Related works with the Fisher Bingham Integral [KW] K. Kume, S. G. Walker, On the Fisher-Bingham distribution, Statistics and Computing 19, , (2009) H. Nakayama, K. Nishiyama, M. Noro, K. Ohara, T. Sei, N. Takayama, A. Takemura, Holonomic Gradient Descent and its Application to the Fisher-Bingham Integral, Advances in Applied Mathematics 47, , (2011) [KNNT] T. Koyama, H. Nakayama, K. Nishiyama, and N. Takayama, Holonomic gradient descent for the Fisher-Bingham distribution on the d-dimensional sphere, Computational Statistics, 29, , (2014) T.Sei, A.Kume, Calculating the Normalising Constant of the Bingham Distribution on the Sphere using the Holonomic Gradient Method, Statistics and Computing, (2013) 5
6 The Ball Probability and the Fisher Bingham Integral. Let Σ = diag(σ 2 1,..., σ2 d ), µ = (µ 1,..., µ d ). Put new parameters as λ i := 1 2σi 2, τ i := µ i σi 2. Then the ball probability is written by the Fisher Bingham integral: P (X X2 d r2 ) di=1 λi = exp 1 4 π d/2 d i=1 τi 2 r λ i 0 f(λ, τ, s)ds. 6
7 The Ordinary Differential Equation w.r.t r. Let f(λ, τ, r) be the Fisher Bingham integral. Put F = ( f,..., f, f,..., f ). τ 1 τ d λ 1 λ d Then, the vector F satisfies the ODE[KNNT]: r F = P r F. The matrix P r = (p ij ) is explicitly given by rp ij = (2λ i r 2 + 1)δ ij + d k=1 τ i δ j(k+d) (1 i d), rp (i+d)j = τ i r 2 δ ij + (2λ i r 2 + 2)δ j(i+d) + k i δ j(k+d) (1 i d) for 1 j 2d. 7
8 P r = 1 r 2r 2 x O y 1 y O 2r 2 x d + 1 y d y d r 2 y 1 O 2r 2 x O r 2 y d 1 2r 2 x d + 2 where O denotes off-diagonal block of 0 s and 1 denotes off-diagonal block of 1 s., 8
9 Remark. Since the Fisher Bingham integral f(λ, τ, r) satisfies the equation f(λ, τ, r) = S d 1 (r) 1 r 2(t t2 d ) exp = 1 r 2 ( f λ f λ d ), d i=1 λ i t 2 i + d i=1 τ i t i dt we can obtain the value of f(λ, τ, r) from the value of the vector F. 9
10 Series Expansion of the Fisher Bingham Integral. The Fisher Bingham integral has the following series expansion [KW]: f(λ, τ, r) = r d 1 S d 1 Here, α,β N d 0 S d 1 = r 2 α+β (d 2)!! d i=1 (2α i + 2β i 1)!! (d α + 2 β )!!α!(2β)! λα τ 2β. Sd 1(1) dt = 2πd/2 Γ(d/2) and N 0 = {0, 1, 2,... }. For a multi-index α N d 0, we define α! = d i=1 α i!, α!! = d i=1 α i!!, and α = d i=1 α i. 10
11 Laplace Approximation (1). Suppose σ 1 > σ 2 σ d, i.e., λ 1 > λ 2 λ d, then we have the Laplace approximation for the Fisher Bingham integral as follows: f(λ, τ, r) f(λ, τ, r). Here, we put f(λ, τ, r) = (e rτ 1 + e rτ 1) exp r 2 λ 1 d i=2 π (d 1)/2 1 di=2 (λ 1 λ i ) 1/2. τ 2 i 4(λ i λ 1 ) 11
12 Laplace Approximation (2). Under the same assumption, we have the Laplace approximation for derivatives of the Fisher Bingham integral: f λ 1 r 2 f, f τ 1 r erτ 1 e rτ 1 e rτ 1 + e rτ 1 where j = 2,..., d. f, f λ j f τ j ( τ j 2(λ j λ 1 ) τ j 2(λ 1 λ j ) f, ) (λ 1 λ j ) f, 12
13 Computation of Initial Values. By the power series expansion of the Fisher Bingham integral, we obtain an approximation for 0 < r 1: f = r d+1 τ i + O(r d+2 ) τ i (i = 1,..., d), f = r d+1 + O(r d+2 ) λ i (i = 1,..., d) In our implementation of the HGM, we take r = for computing the initial value. 13
14 HGM for 0 < r 1. When 0 < r 1, we utilize the ordinary differential equation F r = P rf (1) to evaluate the Fisher Bingham integral by HGM. Here, we put ( f F =,..., f, f,..., f ), τ 1 τ d λ 1 λ d and P r is an explicitly given 2d 2d matrix with rational function elements, as shown above. 14
15 Procedure. Compute the value of r (d+1) F for r = by the approximate expression. The reason of computing r (d+1) F is that the values of elements of F is too small for floating point number. Solve the differential equation (1) numerically. In our implementation, we utilize Explicit embedded Runge-Kutta Prince-Dormand (8, 9) method, and we put the accuracy In order to prevent the elements of F becoming too large, we multiply a constant to F several times when we applying Runge- Kutta method. We use different set of standard monomials for small r and large r. 15
16 HGM for 1 r. When parameter r is large, we utilize the following ordinary differential equation: ( ) Q D r = r D 1 (2rλ 1 + τ 1 )I 2d + DP r D 1 Q (2) where ( 1 Q = exp( r 2 f λ 1 r τ 1 ), f,..., f, 1 f r τ 1 τ 2 τ d r 2, f,..., f λ 1 λ 2 λ d I 2d is the identity matrix with size 2d, and we put ), D = diag( 1 r, 1,..., 1, 1 r2, 1,..., 1). By the Laplace approximation, we expect that Q converges to a vector when r. 16
17 Procedure. Compute the value of F for r = 1.0. Compute the value of Q from F. Solve the differential equation (2) numerically. Compute the value of F from Q. Tamio s code will be made open. 17
18 Cumulative Distribution Function. Let d = 3 and σ 1 = 3.0, σ 2 = 2.0, σ 3 = 1.0, µ 1 = 0.01, µ 2 = 0.02, µ 3 = 0.03, (3) i.e., λ 1 = , λ 2 = 0.125, λ 3 = 0.5, τ 1 = , τ 2 = 0.005, τ 3 = By the HGM, we can compute the probability P ( di=1 X 2 i < r 2) for each r. 18
19 Ball Probability Figure 1: A graph of r and P ( di=1 X 2 i < r 2) r 19
20 Ratio of FB Integral and It s Laplace Approximation. When r 1, the Laplace approximation for the Fisher Bingham integral is f(λ, τ, r) (e rτ 1 + e rτ 1) exp r 2 λ 1 π (d 1)/2 1 di=2 (λ 1 λ i ) 1/2. d i=2 τ 2 i 4(λ i λ 1 ) If HGM is numerically accurate enough, the ratio of both sides of the above expression goes to 1 when r goes to infinity. 20
21 The following numerical calculations confirm this conjecture experimentally. For f/ τ i and f/ λ i, we also obtained the Laplace approximations. We check by numerical calculations that the ratio of the differential of f and it s Laplace approximation goes to 1 for r. 21
22 For the parameter (3), the graph of the ratio are the following: (Fisher Bingham Integral)/(Laprace Approximation) r Figure 2: Ratio of FB Integral and It s Laplace Approximation 22
23 The graph of the ratio of differentials are the following: (Differantial of FB)/(Laprace Approximation) df/d\tau_1 df/d\tau_2 df/d\tau_3 df/d\lambda_1 df/d\lambda_2 df/d\lambda_ r Figure 3: Ratio of differentials of FB and It s LA 23
24 Let d be dimension, σ (1) and σ (2) be diagonal matrices. Suppose the diagonal elements of σ (1) and σ (2) are σ (1) d + 1 i = (1 k d) (4) k(k + 1) and σ (2) i = 2(d + 2)(d + 3) k(k + 1)(k + 2)(k + 3) respectively (based on cumulative χ 2 statistic ). Let µ (1) = 0, µ (2) = ( (d 1) ). (1 k d) (5) 24
25 For each d, we computed the probability P ( d i=1 X 2 i < 40.0 ) and measured the computational time. We considered the following four patterns of parameters: (Σ (1), µ (1) ), (Σ (1), µ (2) ), (Σ (2), µ (1) ), (Σ (2), µ (2) ). 25
26 For Σ (1), we obtained dim µ = µ (1) = 0 µ = µ (2) 0 p time(s) p time(s)
27 For Σ (2), we obtained dim µ = µ (1) = 0 µ = µ (2) 0 p time(s) p time(s)
28 In the case where the dimension is large. Let d be the dimension, and parameter σ k = 1 k(k + 1), µ k = 0 (1 k d). This corresponds to Anderson-Darling statistic, where in fact d =. We computed the probability and measured it s computational time. Here, we fixed the radius as r =
29 dim p time(s) Computational time r 29
30 Summary and Discussion. HGM is again amazingly accurate. It is important to re-scale the elements of the standard monomial vector as needed, so that they do not become too small or too large. Asymptotic evaluation is useful for re-scaling. Question: can we start from infinity and integrate in the reverse direction? In this problem it was not practical. 30
Use of Asymptotics for Holonomic Gradient Method
Use of Asymptotics for Holonomic Gradient Method Akimichi Takemura Shiga University July 25, 2016 A.Takemura (Shiga U) HGM and asymptotics 2015/3/3 1 / 30 Some disclaimers In this talk, I consider numerical
More informationHolonomic Gradient Method for Multivariate Normal Distribution Theory. Akimichi Takemura, Univ. of Tokyo August 2, 2013
Holonomic Gradient Method for Multivariate Normal Distribution Theory Akimichi Takemura, Univ. of Tokyo August 2, 2013 Items 1. First example (non-central t-distribution, by T.Sei) 2. Principles of HGM
More informationHolonomic gradient method for hypergeometric functions of a matrix argument Akimichi Takemura, Shiga University
Holonomic gradient method for hypergeometric functions of a matrix argument Akimichi Takemura, Shiga University (joint with H.Hashiguchi and N.Takayama) June 20, 2016, RIMS, Kyoto Items 1. Summary of multivariate
More informationAnderson-Darling Type Goodness-of-fit Statistic Based on a Multifold Integrated Empirical Distribution Function
Anderson-Darling Type Goodness-of-fit Statistic Based on a Multifold Integrated Empirical Distribution Function S. Kuriki (Inst. Stat. Math., Tokyo) and H.-K. Hwang (Academia Sinica) Bernoulli Society
More informationEffective Sample Size in Spatial Modeling
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS029) p.4526 Effective Sample Size in Spatial Modeling Vallejos, Ronny Universidad Técnica Federico Santa María, Department
More information44 CHAPTER 2. BAYESIAN DECISION THEORY
44 CHAPTER 2. BAYESIAN DECISION THEORY Problems Section 2.1 1. In the two-category case, under the Bayes decision rule the conditional error is given by Eq. 7. Even if the posterior densities are continuous,
More informationLecture 5 : Projections
Lecture 5 : Projections EE227C. Lecturer: Professor Martin Wainwright. Scribe: Alvin Wan Up until now, we have seen convergence rates of unconstrained gradient descent. Now, we consider a constrained minimization
More informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Tuesday, 5 June, 2012 9:00 am to 12:00 pm PAPER 1 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question
More informationAs usual, these notes are intended for use by class participants only, and are not for circulation. Week 6: Lectures 11, 12
As usual, these notes are intended for use by class participants only, and are not for circulation Week 6: Lectures, The Dirac equation and algebra March 5, 0 The Lagrange density for the Dirac equation
More informationDifferential equations
Differential equations Math 7 Spring Practice problems for April Exam Problem Use the method of elimination to find the x-component of the general solution of x y = 6x 9x + y = x 6y 9y Soln: The system
More information3. (a) (8 points) There is more than one way to correctly express the null hypothesis in matrix form. One way to state the null hypothesis is
Stat 501 Solutions and Comments on Exam 1 Spring 005-4 0-4 1. (a) (5 points) Y ~ N, -1-4 34 (b) (5 points) X (X,X ) = (5,8) ~ N ( 11.5, 0.9375 ) 3 1 (c) (10 points, for each part) (i), (ii), and (v) are
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationEigenvectors. Prop-Defn
Eigenvectors Aim lecture: The simplest T -invariant subspaces are 1-dim & these give rise to the theory of eigenvectors. To compute these we introduce the similarity invariant, the characteristic polynomial.
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationStat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2
Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationMATH 3795 Lecture 10. Regularized Linear Least Squares.
MATH 3795 Lecture 10. Regularized Linear Least Squares. Dmitriy Leykekhman Fall 2008 Goals Understanding the regularization. D. Leykekhman - MATH 3795 Introduction to Computational Mathematics Linear Least
More informationAsymptotic Statistics-VI. Changliang Zou
Asymptotic Statistics-VI Changliang Zou Kolmogorov-Smirnov distance Example (Kolmogorov-Smirnov confidence intervals) We know given α (0, 1), there is a well-defined d = d α,n such that, for any continuous
More informationDifferential Equations for Dyson Processes
Differential Equations for Dyson Processes Joint work with Harold Widom I. Overview We call Dyson process any invariant process on ensembles of matrices in which the entries undergo diffusion. Dyson Brownian
More informationMARN 5898 Regularized Linear Least Squares.
MARN 5898 Regularized Linear Least Squares. Dmitriy Leykekhman Spring 2010 Goals Understanding the regularization. D. Leykekhman - MARN 5898 Parameter estimation in marine sciences Linear Least Squares
More informationDepartment of Mathematics, University of Wisconsin-Madison Math 114 Worksheet Sections (4.1),
Department of Mathematics, University of Wisconsin-Madison Math 114 Worksheet Sections (4.1), 4.-4.6 1. Find the polynomial function with zeros: -1 (multiplicity ) and 1 (multiplicity ) whose graph passes
More information1.1. Fields Partial derivatives
1.1. Fields A field associates a physical quantity with a position A field can be also time dependent, for example. The simplest case is a scalar field, where given physical quantity can be described by
More informationDifferentiation of functions of covariance
Differentiation of log X May 5, 2005 1 Differentiation of functions of covariance matrices or: Why you can forget you ever read this Richard Turner Covariance matrices are symmetric, but we often conveniently
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
More informationMATH 162. Midterm 2 ANSWERS November 18, 2005
MATH 62 Midterm 2 ANSWERS November 8, 2005. (0 points) Does the following integral converge or diverge? To get full credit, you must justify your answer. 3x 2 x 3 + 4x 2 + 2x + 4 dx You may not be able
More informationChapter Parabolic Partial Differential Equations
hapter. Parabolic Partial Differential Equations After reading this chapter, you should be able to:. Use numerical methods to solve parabolic partial differential equations by explicit, implicit, and rank-nicolson
More informationOn the Fisher Bingham Distribution
On the Fisher Bingham Distribution BY A. Kume and S.G Walker Institute of Mathematics, Statistics and Actuarial Science, University of Kent Canterbury, CT2 7NF,UK A.Kume@kent.ac.uk and S.G.Walker@kent.ac.uk
More informationMATRIX INTEGRALS AND MAP ENUMERATION 2
MATRIX ITEGRALS AD MAP EUMERATIO 2 IVA CORWI Abstract. We prove the generating function formula for one face maps and for plane diagrams using techniques from Random Matrix Theory and orthogonal polynomials.
More informationExam TFY4230 Statistical Physics kl Wednesday 01. June 2016
TFY423 1. June 216 Side 1 av 5 Exam TFY423 Statistical Physics l 9. - 13. Wednesday 1. June 216 Problem 1. Ising ring (Points: 1+1+1 = 3) A system of Ising spins σ i = ±1 on a ring with periodic boundary
More informationGlobal Optimization of Ordinary Differential Equations Models
Global Optimization of Ordinary Differential Equations Models Angelo Lucia*, Meghan L. Bellows and Leah M. Octavio Department of Chemical Engineering, University of Rhode Island, Kingston, RI 02881 Abstract
More informationSYSTEMTEORI - ÖVNING 1. In this exercise, we will learn how to solve the following linear differential equation:
SYSTEMTEORI - ÖVNING 1 GIANANTONIO BORTOLIN AND RYOZO NAGAMUNE In this exercise, we will learn how to solve the following linear differential equation: 01 ẋt Atxt, xt 0 x 0, xt R n, At R n n The equation
More informationOn the Diagonal Approximation of Full Matrices
On the Diagonal Approximation of Full Matrices Walter M. Lioen CWI P.O. Box 94079, 090 GB Amsterdam, The Netherlands ABSTRACT In this paper the construction of diagonal matrices, in some
More informationEcon 508B: Lecture 5
Econ 508B: Lecture 5 Expectation, MGF and CGF Hongyi Liu Washington University in St. Louis July 31, 2017 Hongyi Liu (Washington University in St. Louis) Math Camp 2017 Stats July 31, 2017 1 / 23 Outline
More information8-5. A rational inequality is an inequality that contains one or more rational expressions. x x 6. 3 by using a graph and a table.
A rational inequality is an inequality that contains one or more rational expressions. x x 3 by using a graph and a table. Use a graph. On a graphing calculator, Y1 = x and Y = 3. x The graph of Y1 is
More informationSerco Assurance. Resonance Theory and Transport Theory in WIMSD J L Hutton
Serco Assurance Resonance Theory and Transport Theory in WIMSD J L Hutton 2 March 2004 Outline of Talk Resonance Treatment Outline of problem - pin cell geometry U 238 cross section Simple non-mathematical
More informationNon white sample covariance matrices.
Non white sample covariance matrices. S. Péché, Université Grenoble 1, joint work with O. Ledoit, Uni. Zurich 17-21/05/2010, Université Marne la Vallée Workshop Probability and Geometry in High Dimensions
More informationLecture 20: Linear model, the LSE, and UMVUE
Lecture 20: Linear model, the LSE, and UMVUE Linear Models One of the most useful statistical models is X i = β τ Z i + ε i, i = 1,...,n, where X i is the ith observation and is often called the ith response;
More informationChapter 7, continued: MANOVA
Chapter 7, continued: MANOVA The Multivariate Analysis of Variance (MANOVA) technique extends Hotelling T 2 test that compares two mean vectors to the setting in which there are m 2 groups. We wish to
More informationConvex Optimization & Parsimony of L p-balls representation
Convex Optimization & Parsimony of L p -balls representation LAAS-CNRS and Institute of Mathematics, Toulouse, France IMA, January 2016 Motivation Unit balls associated with nonnegative homogeneous polynomials
More informationIYGB Mathematical Methods 1
IYGB Mathematical Methods Practice Paper A Time: 3 hours Candidates may use any non programmable, non graphical calculator which does not have the capability of storing data or manipulating algebraic expressions
More informationStat 5101 Notes: Brand Name Distributions
Stat 5101 Notes: Brand Name Distributions Charles J. Geyer September 5, 2012 Contents 1 Discrete Uniform Distribution 2 2 General Discrete Uniform Distribution 2 3 Uniform Distribution 3 4 General Uniform
More informationLecture 10: Dimension Reduction Techniques
Lecture 10: Dimension Reduction Techniques Radu Balan Department of Mathematics, AMSC, CSCAMM and NWC University of Maryland, College Park, MD April 17, 2018 Input Data It is assumed that there is a set
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
More informationThe purpose of this section is to derive the asymptotic distribution of the Pearson chi-square statistic. k (n j np j ) 2. np j.
Chapter 9 Pearson s chi-square test 9. Null hypothesis asymptotics Let X, X 2, be independent from a multinomial(, p) distribution, where p is a k-vector with nonnegative entries that sum to one. That
More informationMass, quasi-local mass, and the flow of static metrics
Mass, quasi-local mass, and the flow of static metrics Eric Woolgar Dept of Mathematical and Statistical Sciences University of Alberta ewoolgar@math.ualberta.ca http://www.math.ualberta.ca/~ewoolgar Nov
More information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationMATH 5640: Functions of Diagonalizable Matrices
MATH 5640: Functions of Diagonalizable Matrices Hung Phan, UMass Lowell November 27, 208 Spectral theorem for diagonalizable matrices Definition Let V = X Y Every v V is uniquely decomposed as u = x +
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationEQ: What are limits, and how do we find them? Finite limits as x ± Horizontal Asymptote. Example Horizontal Asymptote
Finite limits as x ± The symbol for infinity ( ) does not represent a real number. We use to describe the behavior of a function when the values in its domain or range outgrow all finite bounds. For example,
More informationProgress in the method of Ghosts and Shadows for Beta Ensembles
Progress in the method of Ghosts and Shadows for Beta Ensembles Alan Edelman (MIT) Alex Dubbs (MIT) and Plamen Koev (SJS) Oct 8, 2012 1/47 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationForces of Constraint & Lagrange Multipliers
Lectures 30 April 21, 2006 Written or last updated: April 21, 2006 P442 Analytical Mechanics - II Forces of Constraint & Lagrange Multipliers c Alex R. Dzierba Generalized Coordinates Revisited Consider
More informationChapter 5 continued. Chapter 5 sections
Chapter 5 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationMath 115 Spring 11 Written Homework 10 Solutions
Math 5 Spring Written Homework 0 Solutions. For following its, state what indeterminate form the its are in and evaluate the its. (a) 3x 4x 4 x x 8 Solution: This is in indeterminate form 0. Algebraically,
More informationContents. Preface. Notation
Contents Preface Notation xi xv 1 The fractional Laplacian in one dimension 1 1.1 Random walkers with constant steps.............. 1 1.1.1 Particle number density distribution.......... 2 1.1.2 Numerical
More informationSolutions Ph 236b Week 1
Solutions Ph 236b Week 1 Page 1 of 7 Solutions Ph 236b Week 1 Kevin Barkett and Mark Scheel January 19, 216 Contents Problem 1................................... 2 Part (a...................................
More informationHYPERBOLIC SETS WITH NONEMPTY INTERIOR
HYPERBOLIC SETS WITH NONEMPTY INTERIOR TODD FISHER, UNIVERSITY OF MARYLAND Abstract. In this paper we study hyperbolic sets with nonempty interior. We prove the folklore theorem that every transitive hyperbolic
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More informationThe Laplacian and Mean and Extreme Values
Mathematical Assoc. of America American Mathematical Monthly 2: September 29, 205 2:9 p.m. LaplaceMeanValueFinal.tex page The Laplacian and Mean and Extreme Values Jeffrey S. Ovall Abstract. The Laplace
More informationVariations on a Theme: Fields of Definition, Fields of Moduli, Automorphisms, and Twists
Variations on a Theme: Fields of Definition, Fields of Moduli, Automorphisms, and Twists Michelle Manes (mmanes@math.hawaii.edu) ICERM Workshop Moduli Spaces Associated to Dynamical Systems 17 April, 2012
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 15 Special Relativity (Chapter 7) What We Did Last Time Defined Lorentz transformation Linear transformation of 4-vectors that conserve the length in Minkowski space Derived
More informationBIOS 2083 Linear Models Abdus S. Wahed. Chapter 2 84
Chapter 2 84 Chapter 3 Random Vectors and Multivariate Normal Distributions 3.1 Random vectors Definition 3.1.1. Random vector. Random vectors are vectors of random variables. For instance, X = X 1 X 2.
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationLECTURE III Asymmetric Simple Exclusion Process: Bethe Ansatz Approach to the Kolmogorov Equations. Craig A. Tracy UC Davis
LECTURE III Asymmetric Simple Exclusion Process: Bethe Ansatz Approach to the Kolmogorov Equations Craig A. Tracy UC Davis Bielefeld, August 2013 ASEP on Integer Lattice q p suppressed both suppressed
More informationErrata, Mahler Study Aids for Exam S, 2017 HCM, 10/24/17 Page 1. 1, page 2, in Table of Contents: Section 33 is Specific Examples of Markov Processes.
Errata, Mahler Study Aids for Exam S, 2017 HCM, 10/24/17 Page 1 1, page 2, in Table of Contents: Section 33 is Specific Examples of Markov Processes. 1, page 24, solution ot the exercise: (60%)(100) +
More informationRandom Vectors and Multivariate Normal Distributions
Chapter 3 Random Vectors and Multivariate Normal Distributions 3.1 Random vectors Definition 3.1.1. Random vector. Random vectors are vectors of random 75 variables. For instance, X = X 1 X 2., where each
More informationKinetic Theory 1 / Probabilities
Kinetic Theory 1 / Probabilities 1. Motivations: statistical mechanics and fluctuations 2. Probabilities 3. Central limit theorem 1 Reading check Main concept introduced in first half of this chapter A)Temperature
More informationWeek 6: Differential geometry I
Week 6: Differential geometry I Tensor algebra Covariant and contravariant tensors Consider two n dimensional coordinate systems x and x and assume that we can express the x i as functions of the x i,
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More information7.3. Determinants. Introduction. Prerequisites. Learning Outcomes
Determinants 7.3 Introduction Among other uses, determinants allow us to determine whether a system of linear equations has a unique solution or not. The evaluation of a determinant is a key skill in engineering
More informationProposed Proof of Riemann Hypothesis
Proposed Proof of Riemann Hypothesis Aron Palmer Abstract Proposed proof of the Riemann hypothesis showing that positive decreasing continuous function which tends to zero as t goes to infinity can t have
More informationMATHEMATICS 217 NOTES
MATHEMATICS 27 NOTES PART I THE JORDAN CANONICAL FORM The characteristic polynomial of an n n matrix A is the polynomial χ A (λ) = det(λi A), a monic polynomial of degree n; a monic polynomial in the variable
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationML estimation: Random-intercepts logistic model. and z
ML estimation: Random-intercepts logistic model log p ij 1 p = x ijβ + υ i with υ i N(0, συ) 2 ij Standardizing the random effect, θ i = υ i /σ υ, yields log p ij 1 p = x ij β + σ υθ i with θ i N(0, 1)
More information1 Adeles over Q. 1.1 Absolute values
1 Adeles over Q 1.1 Absolute values Definition 1.1.1 (Absolute value) An absolute value on a field F is a nonnegative real valued function on F which satisfies the conditions: (i) x = 0 if and only if
More informationSymplectic integration. Yichao Jing
Yichao Jing Hamiltonian & symplecticness Outline Numerical integrator and symplectic integration Application to accelerator beam dynamics Accuracy and integration order Hamiltonian dynamics In accelerator,
More informationthe renormalization group (RG) idea
the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed
More informationSemidefinite programming lifts and sparse sums-of-squares
1/15 Semidefinite programming lifts and sparse sums-of-squares Hamza Fawzi (MIT, LIDS) Joint work with James Saunderson (UW) and Pablo Parrilo (MIT) Cornell ORIE, October 2015 Linear optimization 2/15
More informationLimits at Infinity. Horizontal Asymptotes. Definition (Limits at Infinity) Horizontal Asymptotes
Limits at Infinity If a function f has a domain that is unbounded, that is, one of the endpoints of its domain is ±, we can determine the long term behavior of the function using a it at infinity. Definition
More informationGradient Methods Using Momentum and Memory
Chapter 3 Gradient Methods Using Momentum and Memory The steepest descent method described in Chapter always steps in the negative gradient direction, which is orthogonal to the boundary of the level set
More informationEXERCISES IN MODULAR FORMS I (MATH 726) (2) Prove that a lattice L is integral if and only if its Gram matrix has integer coefficients.
EXERCISES IN MODULAR FORMS I (MATH 726) EYAL GOREN, MCGILL UNIVERSITY, FALL 2007 (1) We define a (full) lattice L in R n to be a discrete subgroup of R n that contains a basis for R n. Prove that L is
More informationHyperuniformity on the Sphere
Hyperuniformity on the Sphere Peter Grabner Institute for Analysis and Number Theory Graz University of Technology Problem Session, ICERM, February 12, 2018 Two point distributions Quantify evenness For
More informationDepartment of Mathematics, University of California, Berkeley. GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016
Department of Mathematics, University of California, Berkeley YOUR 1 OR 2 DIGIT EXAM NUMBER GRADUATE PRELIMINARY EXAMINATION, Part A Fall Semester 2016 1. Please write your 1- or 2-digit exam number on
More informationPhysics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =...
Physics 606, Quantum Mechanics, Final Exam NAME Please show all your work. (You are graded on your work, with partial credit where it is deserved.) All problems are, of course, nonrelativistic. 1. Consider
More informationConvolutions and fluctuations: free, finite, quantized.
Convolutions and fluctuations: free, finite, quantized. Vadim Gorin MIT (Cambridge) and IITP (Moscow) February 2018 Hermitian matrix operations a 1 0 b 1 0 0 a 2 0 0 b 2 0 A =. 0.. B = 0. 0.. 0 0 a N 0
More informationab = c a If the coefficients a,b and c are real then either α and β are real or α and β are complex conjugates
Further Pure Summary Notes. Roots of Quadratic Equations For a quadratic equation ax + bx + c = 0 with roots α and β Sum of the roots Product of roots a + b = b a ab = c a If the coefficients a,b and c
More informationUPPER DEVIATIONS FOR SPLIT TIMES OF BRANCHING PROCESSES
Applied Probability Trust 7 May 22 UPPER DEVIATIONS FOR SPLIT TIMES OF BRANCHING PROCESSES HAMED AMINI, AND MARC LELARGE, ENS-INRIA Abstract Upper deviation results are obtained for the split time of a
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationCommutative Matrix Multiplication
Commutative Matrix Multiplication Daniel Beer and Henning Thielemann 3th August 009 Introduction Matrix interpolation for fractals geometric interpolation, which one is the right: A p B p or B p A p? Definition
More informationMatrices of Dirac Characters within an irrep
Matrices of Dirac Characters within an irrep irrep E : 1 0 c s 2 c s D( E) D( C ) D( C ) 3 3 0 1 s c s c 1 0 c s c s D( ) D( ) D( ) a c b 0 1 s c s c 2 1 2 3 c cos( ), s sin( ) 3 2 3 2 E C C 2 3 3 2 3
More informationCS 257: Numerical Methods
CS 57: Numerical Methods Final Exam Study Guide Version 1.00 Created by Charles Feng http://www.fenguin.net CS 57: Numerical Methods Final Exam Study Guide 1 Contents 1 Introductory Matter 3 1.1 Calculus
More informationThe Yau-Tian-Donaldson Conjectuture for general polarizations
The Yau-Tian-Donaldson Conjectuture for general polarizations Toshiki Mabuchi, Osaka University 2015 Taipei Conference on Complex Geometry December 22, 2015 1. Introduction 2. Background materials Table
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Approximations Simo Särkkä Aalto University, Finland November 18, 2014 Simo Särkkä (Aalto) Lecture 4: Numerical Solution of SDEs November
More informationChecking the Radioactive Decay Euler Algorithm
Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present
More informationQuantum Chaos and Nonunitary Dynamics
Quantum Chaos and Nonunitary Dynamics Karol Życzkowski in collaboration with W. Bruzda, V. Cappellini, H.-J. Sommers, M. Smaczyński Phys. Lett. A 373, 320 (2009) Institute of Physics, Jagiellonian University,
More informationLecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations
Lecture 4: Numerical Solution of SDEs, Itô Taylor Series, Gaussian Process Approximations Simo Särkkä Aalto University Tampere University of Technology Lappeenranta University of Technology Finland November
More informationJoint work with Nottingham colleagues Simon Preston and Michail Tsagris.
/pgf/stepx/.initial=1cm, /pgf/stepy/.initial=1cm, /pgf/step/.code=1/pgf/stepx/.expanded=- 10.95415pt,/pgf/stepy/.expanded=- 10.95415pt, /pgf/step/.value required /pgf/images/width/.estore in= /pgf/images/height/.estore
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationWorkshop on Discrete Harmonic Analysis Newton Institute, March 2011
High frequency criteria for Boolean functions (with an application to percolation) Christophe Garban ENS Lyon and CNRS Workshop on Discrete Harmonic Analysis Newton Institute, March 2011 C. Garban (ENS
More information