Tube formula approach to testing multivariate normality and testing uniformity on the sphere
|
|
- Dwain Morgan
- 5 years ago
- Views:
Transcription
1 Tube formula approach to testing multivariate normality and testing uniformity on the sphere Akimichi Takemura 1 Satoshi Kuriki 2 1 University of Tokyo 2 Institute of Statistical Mathematics December 11, 2010 December 11, / 24
2 References This talk is based on the following two papers: Kuriki and Takemura (2008). The tube method for the moment index in projection pursuit. Journal of Statistical Planning and Inference, 138, No.9, Kuriki and Takemura (2004). Tail probabilities of the limiting null distributions of the Anderson-Stephens statistics. Journal of Multivariate Analysis, 89, December 11, / 24
3 Contents 1 Tube formula approximation to maximum type test statistics 2 Projection pursuit index and testing multivariate normality 3 Anderson-Stephens statistic for testing uniformity on sphere 4 Summary December 11, / 24
4 Tube formula approximation to maximum type test statistics Tube formula: some historical background Jacob Steiner already had Steiner s formula (1840) for the volume of a tube of a convex set. Minkowski defined mixed volumes. Hotelling (1939) derived the tube formula for a one-dimensional curve and then H.Weyl immediately generalized it to a general dimension. Hotelling s motivation was a nonlinear regression problem. Revival of tube formula in statistics around (Knowles-Siegmund(1989), J.Sun(1991,93) and many other people). December 11, / 24
5 Tube formula approximation to maximum type test statistics Euler characteristic method (independent development) Euler characteristic heuristic was initiated by R.J.Adler for approximating the distribution of the maximum of a random random field (Adler-Hasofer(1976), Adler s book(1981)). This method has been vigorously developed by Adler and Keith Worsley. Some important foundational work was done by Jonathan Taylor (2001 thesis). A standard textbook now is Random Fields and Geometry by Adler and Taylor, 2007, Springer. December 11, / 24
6 Tube formula approximation to maximum type test statistics Two methods are equivalent Around 2000, I and Kuriki were sitting in a talk by Keith Worsley in ISM (Institute of Statistical Mathematics, Tokyo, Japan) and was surprised that he was doing the same computations as us. Takemura and Kuriki (2002) proved the equivalence of these two methods by using Morse theorem (for finite dimensional case). Tube method can be understood as finite dimensional specialization of Euler characteristic method. I should also mention that abstract tube by Naiman and Wynn is a discrete analog of tube formula. December 11, / 24
7 Tube formula approximation to maximum type test statistics Canonical form of tube formula Let z = (z 1,...,z n ) N n (0, I n ). Let M S n 1 be a C 2 -submanifold of dimension d = dimm with piecewise smooth boundaries. Let Z(u) = u z = n u i z i, u = (u 1,...,u n ) M. i=1 Also consider a standardized random field Y (u) = u z/ z, u M, z = z z. December 11, / 24
8 Tube formula approximation to maximum type test statistics Canonical form of tube formula We want to evaluate the distributions of maxima, corresponding to maximum type test statistics: T = max u M Z(u), U = max u M Y (u). The tube method gives an approximation of the tail probabilities P(T x), x, and P(U x), x 1. December 11, / 24
9 Tube formula approximation to maximum type test statistics Spherical tube and its volume Evaluation of the distribution reduces to the evaluation of the volume of a spherical tube around M. M M θ 0 S n-1 Figure: Spherical tube around M December 11, / 24
10 Tube formula approximation to maximum type test statistics Spherical tube and its volume Let M θ = { } v S n 1 min u M cos 1 (u v) θ denote the tube around M with radius θ. Let Vol(M θ ) denote the (n 1)-dimensional spherical volume of M θ. By definition where ( ) P max Y (u) cos θ = Vol(M θ )/Ω n, u M Ω n = Vol(S n 1 ) = 2πn/2 Γ(n/2) and B a,b ( ) denotes the upper probability of beta distribution with parameter (a, b). December 11, / 24
11 Tube formula approximation to maximum type test statistics Tube formula for the volume of a spherical tube Tube formula: For θ smaller than the critical radius, Vol(M θ ) = Ω n {w d+1 B d+1 2, n d w 1 B 1 2, n 1 (cos 2 θ) 2 (cos 2 θ) + w d B d 2, n d (cos 2 θ) 2 }, where w 1,...,w d+1 are geometric invariants of M, which can be evaluated by differential geometric methods. In particular w d+1 = Vol(M)/Ω d+1, w d = Vol( M)/(2Ω d+1 ). We omit explanation of critical radius in this talk. December 11, / 24
12 Tube formula approximation to maximum type test statistics Tail probability for T (non-standardized maximum) For T = max u M Z(u) we need integration of the tube formula in z. By integration on z we have ( ) P max Z(u) x = w d+1 Ḡ d+1 (x 2 ) + w d Ḡ d (x 2 ) + u M +w 1 Ḡ 1 (x 2 ) + O(Ḡ n (x 2 (1 + tan 2 θ c ))), where Ḡ a ( ) is the upper probability of χ 2 distribution with a degrees of freedom and θ c is the critical radius. December 11, / 24
13 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson x t R q, t = 1,...,n: observation vectors S q 1 : the unit sphere in R q. z t = h x t : projection of x t onto the direction h Projection pursuit: looks for the direction h such that the projected data z 1,...,z n do not look like normally distributed. I n (h): projection pursuit index, which measures non-normality of projected data. maximize I n (h) in h S q 1. December 11, / 24
14 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson The null hypothesis: H 0 : x t N q (µ,σ), i.i.d. K k,n (h): the kth sample cumulant of the projected data z 1,...,z n. Skewness and kurtosis: B 1,n (h) = K 3,n (h)/k 2,n (h) 3/2 : the sample skewness B 2,n (h) = K 4,n (h)/k 2,n (h) 2 : the sample kurtosis Jones-Sibson index I n (h) = n 6 B 1,n(h) 2 + n 24 B 2,n(h) 2 Asymptotic distribution of max h S q 1 I n (h)? can be solved by tube formula. December 11, / 24
15 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson Theorem 1 (Asymptotic distribution of the random field) Let ξ 1 R q3, ξ 2 R q4 be random vectors consisting of independent standard normal random variables. For a unit vector h S q 1, let Z 1 (h) = (h h h) ξ 1, Z 2 (h) = (h h h h) ξ 2, where denotes the Kronecker product. Under the null hypothesis H 0 of multivariate normality, as n, max h S q 1 I n (h) converges in distribution to max h S q 1 I(h), where I(h) = Z 1 (h) 2 + Z 2 (h) 2. December 11, / 24
16 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson Theorem 2 (Asymptotic distribution of Jones-Sibson index) As c, ( P max I(h) c2) = h S q 1 where q e=0, e:even κ e ( 3) e/2 (q 1)! κ e = Ω q (q e)! E k = π/2 π/2 Γ( q+1 e 2 ) 2 1+e/2 π (q+1)/2ḡq+1 e(c 2 )(1 + o(1)), e/2 j=0 (3 cos 2 θ + 4 sin 2 θ) k dθ. (q e 2j) (e/2 j)! j! ( 2)j E (q 1 e)/2 j, December 11, / 24
17 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson Note that we are first letting n and then getting the tail probability c n= tube approx tail probability x Figure: Tail probability of limiting distribution (solid line) and its approximation by the tube method (dotted line). December 11, / 24
18 Projection pursuit index and testing multivariate normality Projection pursuit index by Jones and Sibson n= n=3000 n=1000 n= tail probability x Figure: Tail probabilities of finite sample distributions (n = 300,1000,3000, ). December 11, / 24
19 Anderson-Stephens statistic for testing uniformity on sphere Anderson-Stephens statistic H 0 : x t, t = 1,...,n are i.i.d. uniform vectors on S q 1. For h S q 1 let S(h) = 1 n n (h x t ) 2 t=1 Let S max = max h S q 1 S(h), S min = min h S q 1 S(h). Anderson-Stephens test: reject H 0 if S max c or S min c. We also propose to use S range = S max S min December 11, / 24
20 Anderson-Stephens statistic for testing uniformity on sphere Anderson-Stephens statistic S max and S min are the largest and smallest eigenvalues λ 1 (Q) and λ q (Q) of a q q matrix Q = 1 n n z t z t. Let A : q q, symmetric, have the multivariate symmetric normal distribution, i.e., a ii N(0, 1), a ij, i < j, N(0, 1/2), all mutually independent. The limiting null distribution of the eigenvalues of n(q I q /q) is given by the distribution of the eigenvalues of 2 tr(a) (A q(q + 2) q I q) tube formula works again. December 11, / 24 t=1
21 Anderson-Stephens statistic for testing uniformity on sphere Anderson-Stephens statistic Need some more adjustment of constants: Lemma 3 As n, the null distributions of both of n(s max 1/q) and n(s min 1/q) converge to the distribution of 2(q 1)/q 2 (q + 2)T 1, where ( q T 1 = λ 1 (B) with B = A tr(a) ) q 1 q I q. (1) The null distribution of n(s max S min ) converges to the distribution of (2/ q(q + 2))T 2, where T 2 = 1 2 (λ 1 (A) λ q (A)). (2) December 11, / 24
22 Anderson-Stephens statistic for testing uniformity on sphere Anderson-Stephens statistic Theorem 4 When q 3, the asymptotic expansion of the upper tail probability of T 1 = λ 1 (B) is given by where When q = 2, P(T 1 x) = w q e = 1 2 q 1 e=0, e:even w q e Ḡ q e (x 2 )(1 + o(1)), x, (3) ( ) 2q (q 1)/2 ( q + 1 ) e/2 Γ( q+1 2 ) q 1 2q Γ( q e+1 2 )( e (4) 2 )!. P(T 1 x) = Ḡ2(x 2 ), x 0. December 11, / 24
23 Anderson-Stephens statistic for testing uniformity on sphere Anderson-Stephens statistic tail probability limiting distribution n=10,100,1000 approximation by the tube method x Tail probabilities of S max when q = 3. (n = 10, 100, 1000, and approximation by the tube method.) December 11, / 24
24 Summary Summary We gave a brief introduction to tube method. We applied the method to testing multivariate normality based on projection pursuit index and to testing uniformity on the unit sphere. We presented numerical examples to show that tube formula gives a good approximation to the tail probability of test statistics. December 11, / 24
Application of tube formula to distributional problems in multiway layouts
Application of tube formula to distributional problems in multiway layouts Satoshi Kuriki Institute of Statistical Mathematics and Akimichi Takemura Graduate School of Information Science and Technology
More informationON THE EQUIVALENCE OF THE TUBE AND EULER CHARACTERISTIC METHODS FOR THE DISTRIBUTION OF THE MAXIMUM OF GAUSSIAN FIELDS OVER PIECEWISE SMOOTH DOMAINS
The Annals of Applied Probability 2001, Vol. 00, No. 00, 000 000 ON THE EQUIVALENCE OF THE TUBE AND EULER CHARACTERISTIC METHODS FOR THE DISTRIBUTION OF THE MAXIMUM OF GAUSSIAN FIELDS OVER PIECEWISE SMOOTH
More informationTail probability of linear combinations of chi-square variables and its application to influence analysis in QTL detection
Tail probability of linear combinations of chi-square variables and its application to influence analysis in QTL detection Satoshi Kuriki and Xiaoling Dou (Inst. Statist. Math., Tokyo) ISM Cooperative
More informationRANDOM FIELDS AND GEOMETRY. Robert Adler and Jonathan Taylor
RANDOM FIELDS AND GEOMETRY from the book of the same name by Robert Adler and Jonathan Taylor IE&M, Technion, Israel, Statistics, Stanford, US. ie.technion.ac.il/adler.phtml www-stat.stanford.edu/ jtaylor
More informationA test for a conjunction
A test for a conjunction K.J. Worsley a,, K.J. Friston b a Department of Mathematics and Statistics, McGill University, Montreal, Canada H3A 2K6 b Wellcome Department of Cognitive Neurology, Institute
More informationSOME NOTES ON HOTELLING TUBES
SOME NOTES ON HOTELLING TUBES ROGER KOENKER 1. Introduction We begin with the two examples of Johansen and Johnstone. The first is Hotelling s (1939) original example, and the second is closely related
More informationComment. February 1, 2008
Persi Diaconis Stanford University Erich Lehmann UC Berkeley February 1, 2008 Sandy Zabell gives us a good feel for Student and his times. We supplement this with later developments showing that Student
More informationStatistical applications of geometry and random
Statistical applications of geometry and random fields 3 e Cycle romand de statistique et probabilités appliquées Based on joint work with: R. Adler, K. Worsley Jonathan Taylor Stanford University / Université
More informationON THE EQUIVALENCE OF THE TUBE AND EULER CHARACTERISTIC METHODS FOR THE DISTRIBUTION OF THE MAXIMUM OF GAUSSIAN FIELDS OVER PIECEWISE SMOOTH DOMAINS
The Annals of Applied Probability 2002, Vol. 12, No. 2, 768 796 ON THE EQUIVALENCE OF THE TUBE AND EULER CHARACTERISTIC METHODS FOR THE DISTRIBUTION OF THE MAXIMUM OF GAUSSIAN FIELDS OVER PIECEWISE SMOOTH
More informationStatistical Inference
Statistical Inference Classical and Bayesian Methods Revision Class for Midterm Exam AMS-UCSC Th Feb 9, 2012 Winter 2012. Session 1 (Revision Class) AMS-132/206 Th Feb 9, 2012 1 / 23 Topics Topics We will
More informationSpatial Regression. 3. Review - OLS and 2SLS. Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 3. Review - OLS and 2SLS Luc Anselin http://spatial.uchicago.edu OLS estimation (recap) non-spatial regression diagnostics endogeneity - IV and 2SLS OLS Estimation (recap) Linear Regression
More informationAkaike Information Criterion
Akaike Information Criterion Shuhua Hu Center for Research in Scientific Computation North Carolina State University Raleigh, NC February 7, 2012-1- background Background Model statistical model: Y j =
More informationTail probability via the tube formula when the critical radius is zero
Bernoulli 9(3), 2003, 535 558 Tail probability via the tube formula when the critical radius is zero AKIMICHI TAKEMURA 1 and SATOSHI KURIKI 2 1 Department of Mathematical Informatics, Graduate School of
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 informationEvgeny Spodarev WIAS, Berlin. Limit theorems for excursion sets of stationary random fields
Evgeny Spodarev 23.01.2013 WIAS, Berlin Limit theorems for excursion sets of stationary random fields page 2 LT for excursion sets of stationary random fields Overview 23.01.2013 Overview Motivation Excursion
More informationCentral Limit Theorem ( 5.3)
Central Limit Theorem ( 5.3) Let X 1, X 2,... be a sequence of independent random variables, each having n mean µ and variance σ 2. Then the distribution of the partial sum S n = X i i=1 becomes approximately
More informationApplications of the signed distance function to surface geometry. Daniel Mayost
Applications of the signed distance function to surface geometry by Daniel ayost A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Graduate Department of athematics
More information1/24/2008. Review of Statistical Inference. C.1 A Sample of Data. C.2 An Econometric Model. C.4 Estimating the Population Variance and Other Moments
/4/008 Review of Statistical Inference Prepared by Vera Tabakova, East Carolina University C. A Sample of Data C. An Econometric Model C.3 Estimating the Mean of a Population C.4 Estimating the Population
More informationStatistical Inference with Monotone Incomplete Multivariate Normal Data
Statistical Inference with Monotone Incomplete Multivariate Normal Data p. 1/4 Statistical Inference with Monotone Incomplete Multivariate Normal Data This talk is based on joint work with my wonderful
More information11 Survival Analysis and Empirical Likelihood
11 Survival Analysis and Empirical Likelihood The first paper of empirical likelihood is actually about confidence intervals with the Kaplan-Meier estimator (Thomas and Grunkmeier 1979), i.e. deals with
More informationPRE-TEST ESTIMATION OF THE REGRESSION SCALE PARAMETER WITH MULTIVARIATE STUDENT-t ERRORS AND INDEPENDENT SUB-SAMPLES
Sankhyā : The Indian Journal of Statistics 1994, Volume, Series B, Pt.3, pp. 334 343 PRE-TEST ESTIMATION OF THE REGRESSION SCALE PARAMETER WITH MULTIVARIATE STUDENT-t ERRORS AND INDEPENDENT SUB-SAMPLES
More information4.5.1 The use of 2 log Λ when θ is scalar
4.5. ASYMPTOTIC FORM OF THE G.L.R.T. 97 4.5.1 The use of 2 log Λ when θ is scalar Suppose we wish to test the hypothesis NH : θ = θ where θ is a given value against the alternative AH : θ θ on the basis
More informationINFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES
Submitted to the Annals of Statistics INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES By Armin Schwartzman, Walter F. Mascarenhas and Jonathan E. Taylor Harvard School of Public
More informationÜbungen zu RT2 SS (4) Show that (any) contraction of a (p, q) - tensor results in a (p 1, q 1) - tensor.
Übungen zu RT2 SS 2010 (1) Show that the tensor field g µν (x) = η µν is invariant under Poincaré transformations, i.e. x µ x µ = L µ νx ν + c µ, where L µ ν is a constant matrix subject to L µ ρl ν ση
More information1 Appendix A: Matrix Algebra
Appendix A: Matrix Algebra. Definitions Matrix A =[ ]=[A] Symmetric matrix: = for all and Diagonal matrix: 6=0if = but =0if 6= Scalar matrix: the diagonal matrix of = Identity matrix: the scalar matrix
More informationHigh-dimensional asymptotic expansions for the distributions of canonical correlations
Journal of Multivariate Analysis 100 2009) 231 242 Contents lists available at ScienceDirect Journal of Multivariate Analysis journal homepage: www.elsevier.com/locate/jmva High-dimensional asymptotic
More informationStatistical Inference with Monotone Incomplete Multivariate Normal Data
Statistical Inference with Monotone Incomplete Multivariate Normal Data p. 1/4 Statistical Inference with Monotone Incomplete Multivariate Normal Data This talk is based on joint work with my wonderful
More informationSpectra of adjacency matrices of random geometric graphs
Spectra of adjacency matrices of random geometric graphs Paul Blackwell, Mark Edmondson-Jones and Jonathan Jordan University of Sheffield 22nd December 2006 Abstract We investigate the spectral properties
More informationElliptically Contoured Distributions
Elliptically Contoured Distributions Recall: if X N p µ, Σ), then { 1 f X x) = exp 1 } det πσ x µ) Σ 1 x µ) So f X x) depends on x only through x µ) Σ 1 x µ), and is therefore constant on the ellipsoidal
More informationMS-E2112 Multivariate Statistical Analysis (5cr) Lecture 6: Bivariate Correspondence Analysis - part II
MS-E2112 Multivariate Statistical Analysis (5cr) Lecture 6: Bivariate Correspondence Analysis - part II the Contents the the the Independence The independence between variables x and y can be tested using.
More informationInference on distributions and quantiles using a finite-sample Dirichlet process
Dirichlet IDEAL Theory/methods Simulations Inference on distributions and quantiles using a finite-sample Dirichlet process David M. Kaplan University of Missouri Matt Goldman UC San Diego Midwest Econometrics
More informationHypothesis Testing For Multilayer Network Data
Hypothesis Testing For Multilayer Network Data Jun Li Dept of Mathematics and Statistics, Boston University Joint work with Eric Kolaczyk Outline Background and Motivation Geometric structure of multilayer
More informationMinimum distance tests and estimates based on ranks
Minimum distance tests and estimates based on ranks Authors: Radim Navrátil Department of Mathematics and Statistics, Masaryk University Brno, Czech Republic (navratil@math.muni.cz) Abstract: It is well
More informationA significance test for the lasso
1 Gold medal address, SSC 2013 Joint work with Richard Lockhart (SFU), Jonathan Taylor (Stanford), and Ryan Tibshirani (Carnegie-Mellon Univ.) Reaping the benefits of LARS: A special thanks to Brad Efron,
More informationSTAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots. March 8, 2015
STAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots March 8, 2015 The duality between CI and hypothesis testing The duality between CI and hypothesis
More informationStationarity of non-radiating spacetimes
University of Warwick April 4th, 2016 Motivation Theorem Motivation Newtonian gravity: Periodic solutions for two-body system. Einstein gravity: Periodic solutions? At first Post-Newtonian order, Yes!
More informationStatistical Inference On the High-dimensional Gaussian Covarianc
Statistical Inference On the High-dimensional Gaussian Covariance Matrix Department of Mathematical Sciences, Clemson University June 6, 2011 Outline Introduction Problem Setup Statistical Inference High-Dimensional
More informationInstitute of Actuaries of India
Institute of Actuaries of India Subject CT3 Probability & Mathematical Statistics May 2011 Examinations INDICATIVE SOLUTION Introduction The indicative solution has been written by the Examiners with the
More informationLocal Whittle Likelihood Estimators and Tests for non-gaussian Linear Processes
Local Whittle Likelihood Estimators and Tests for non-gaussian Linear Processes By Tomohito NAITO, Kohei ASAI and Masanobu TANIGUCHI Department of Mathematical Sciences, School of Science and Engineering,
More informationStatistical Inference
Statistical Inference Classical and Bayesian Methods Class 7 AMS-UCSC Tue 31, 2012 Winter 2012. Session 1 (Class 7) AMS-132/206 Tue 31, 2012 1 / 13 Topics Topics We will talk about... 1 Hypothesis testing
More informationGEOMETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE
GEOETRICAL TOOLS FOR PDES ON A SURFACE WITH ROTATING SHALLOW WATER EQUATIONS ON A SPHERE BIN CHENG Abstract. This is an excerpt from my paper with A. ahalov [1]. PDE theories with Riemannian geometry are
More informationp(z)
Chapter Statistics. Introduction This lecture is a quick review of basic statistical concepts; probabilities, mean, variance, covariance, correlation, linear regression, probability density functions and
More informationGARCH Models Estimation and Inference
GARCH Models Estimation and Inference Eduardo Rossi University of Pavia December 013 Rossi GARCH Financial Econometrics - 013 1 / 1 Likelihood function The procedure most often used in estimating θ 0 in
More informationOn the conservative multivariate multiple comparison procedure of correlated mean vectors with a control
On the conservative multivariate multiple comparison procedure of correlated mean vectors with a control Takahiro Nishiyama a a Department of Mathematical Information Science, Tokyo University of Science,
More information2.6.3 Generalized likelihood ratio tests
26 HYPOTHESIS TESTING 113 263 Generalized likelihood ratio tests When a UMP test does not exist, we usually use a generalized likelihood ratio test to verify H 0 : θ Θ against H 1 : θ Θ\Θ It can be used
More informationStatistical Inference: Estimation and Confidence Intervals Hypothesis Testing
Statistical Inference: Estimation and Confidence Intervals Hypothesis Testing 1 In most statistics problems, we assume that the data have been generated from some unknown probability distribution. We desire
More informationSummary and discussion of: Exact Post-selection Inference for Forward Stepwise and Least Angle Regression Statistics Journal Club
Summary and discussion of: Exact Post-selection Inference for Forward Stepwise and Least Angle Regression Statistics Journal Club 36-825 1 Introduction Jisu Kim and Veeranjaneyulu Sadhanala In this report
More informationUnified univariate and multivariate random field theory
Unified univariate and multivariate random field theory Keith J. Worsley 1, Jonathan E. Taylor 3, Francesco Tomaiuolo 4, Jason Lerch 1 Department of Mathematics and Statistics, and Montreal Neurological
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More information5 Introduction to the Theory of Order Statistics and Rank Statistics
5 Introduction to the Theory of Order Statistics and Rank Statistics This section will contain a summary of important definitions and theorems that will be useful for understanding the theory of order
More informationRANDOM FIELDS OF MULTIVARIATE TEST STATISTICS, WITH APPLICATIONS TO SHAPE ANALYSIS
Submitted to the Annals of Statistics RANDOM FIELDS OF MULTIVARIATE TEST STATISTICS, WITH APPLICATIONS TO SHAPE ANALYSIS By J.E. Taylor and K.J. Worsley Stanford University and McGill University Our data
More informationGeometric Interpolation by Planar Cubic Polynomials
1 / 20 Geometric Interpolation by Planar Cubic Polynomials Jernej Kozak, Marjeta Krajnc Faculty of Mathematics and Physics University of Ljubljana Institute of Mathematics, Physics and Mechanics Avignon,
More informationAn introduction to General Relativity and the positive mass theorem
An introduction to General Relativity and the positive mass theorem National Center for Theoretical Sciences, Mathematics Division March 2 nd, 2007 Wen-ling Huang Department of Mathematics University of
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationStat 451 Lecture Notes Numerical Integration
Stat 451 Lecture Notes 03 12 Numerical Integration Ryan Martin UIC www.math.uic.edu/~rgmartin 1 Based on Chapter 5 in Givens & Hoeting, and Chapters 4 & 18 of Lange 2 Updated: February 11, 2016 1 / 29
More informationLecture 7. Quaternions
Matthew T. Mason Mechanics of Manipulation Spring 2012 Today s outline Motivation Motivation have nice geometrical interpretation. have advantages in representing rotation. are cool. Even if you don t
More informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationResearch Article A Nonparametric Two-Sample Wald Test of Equality of Variances
Advances in Decision Sciences Volume 211, Article ID 74858, 8 pages doi:1.1155/211/74858 Research Article A Nonparametric Two-Sample Wald Test of Equality of Variances David Allingham 1 andj.c.w.rayner
More informationA CONSTRUCTION OF TRANSVERSE SUBMANIFOLDS
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XLI 2003 A CONSTRUCTION OF TRANSVERSE SUBMANIFOLDS by J. Szenthe Abstract. In case of Riemannian manifolds isometric actions admitting submanifolds
More informationDetecting fmri activation allowing for unknown latency of the hemodynamic response
Detecting fmri activation allowing for unknown latency of the hemodynamic response K.J. Worsley McGill University J.E. Taylor Stanford University January 7, 006 Abstract Several authors have suggested
More informationTesting Statistical Hypotheses
E.L. Lehmann Joseph P. Romano Testing Statistical Hypotheses Third Edition 4y Springer Preface vii I Small-Sample Theory 1 1 The General Decision Problem 3 1.1 Statistical Inference and Statistical Decisions
More informationTesting Restrictions and Comparing Models
Econ. 513, Time Series Econometrics Fall 00 Chris Sims Testing Restrictions and Comparing Models 1. THE PROBLEM We consider here the problem of comparing two parametric models for the data X, defined by
More informationDETECTING SPARSE CONE ALTERNATIVES FOR GAUSSIAN RANDOM FIELDS,
Statistica Sinica 23 2013), 1629-1656 doi:http://dx.doi.org/10.5705/ss.2012-218s DETECTING SPARSE CONE ALTERNATIVES FOR GAUSSIAN RANDOM FIELDS, WITH AN APPLICATION TO fmri Jonathan E. Taylor 1 and Keith
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
More informationThe Canonical Gaussian Measure on R
The Canonical Gaussian Measure on R 1. Introduction The main goal of this course is to study Gaussian measures. The simplest example of a Gaussian measure is the canonical Gaussian measure P on R where
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 informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
More information4sec 2xtan 2x 1ii C3 Differentiation trig
A Assignment beta Cover Sheet Name: Question Done Backpack Topic Comment Drill Consolidation i C3 Differentiation trig 4sec xtan x ii C3 Differentiation trig 6cot 3xcosec 3x iii C3 Differentiation trig
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 informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 9 for Applied Multivariate Analysis Outline Addressing ourliers 1 Addressing ourliers 2 Outliers in Multivariate samples (1) For
More informationLecture 3: Central Limit Theorem
Lecture 3: Central Limit Theorem Scribe: Jacy Bird (Division of Engineering and Applied Sciences, Harvard) February 8, 003 The goal of today s lecture is to investigate the asymptotic behavior of P N (
More informationABOUT PRINCIPAL COMPONENTS UNDER SINGULARITY
ABOUT PRINCIPAL COMPONENTS UNDER SINGULARITY José A. Díaz-García and Raúl Alberto Pérez-Agamez Comunicación Técnica No I-05-11/08-09-005 (PE/CIMAT) About principal components under singularity José A.
More informationAn Introduction to Multivariate Statistical Analysis
An Introduction to Multivariate Statistical Analysis Third Edition T. W. ANDERSON Stanford University Department of Statistics Stanford, CA WILEY- INTERSCIENCE A JOHN WILEY & SONS, INC., PUBLICATION Contents
More informationEngineering Mechanics Statics
Mechanical Systems Engineering -2016 Engineering Mechanics Statics 3. Force Vectors; Position Vector & Dot product Position Vector A position vector is a fixed vector that locates a point in space relative
More informationOn Spatial Involute Gearing
6 th International Conference on Applied Informatics Eger, Hungary, January 27 3, 2004. On Spatial Involute Gearing Hellmuth Stachel Institute of Discrete Mathematics and Geometry, Vienna University of
More information+ Specify 1 tail / 2 tail
Week 2: Null hypothesis Aeroplane seat designer wonders how wide to make the plane seats. He assumes population average hip size μ = 43.2cm Sample size n = 50 Question : Is the assumption μ = 43.2cm reasonable?
More informationCHAPTER 3. Gauss map. In this chapter we will study the Gauss map of surfaces in R 3.
CHAPTER 3 Gauss map In this chapter we will study the Gauss map of surfaces in R 3. 3.1. Surfaces in R 3 Let S R 3 be a submanifold of dimension 2. Let {U i, ϕ i } be a DS on S. For any p U i we have a
More informationGeometric projection of stochastic differential equations
Geometric projection of stochastic differential equations John Armstrong (King s College London) Damiano Brigo (Imperial) August 9, 2018 Idea: Projection Idea: Projection Projection gives a method of systematically
More informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationDaniel M. Oberlin Department of Mathematics, Florida State University. January 2005
PACKING SPHERES AND FRACTAL STRICHARTZ ESTIMATES IN R d FOR d 3 Daniel M. Oberlin Department of Mathematics, Florida State University January 005 Fix a dimension d and for x R d and r > 0, let Sx, r) stand
More informationHigher order moments of the estimated tangency portfolio weights
WORKING PAPER 0/07 Higher order moments of the estimated tangency portfolio weights Farrukh Javed, Stepan Mazur, Edward Ngailo Statistics ISSN 403-0586 http://www.oru.se/institutioner/handelshogskolan-vid-orebro-universitet/forskning/publikationer/working-papers/
More informationLikelihood Ratio Tests. that Certain Variance Components Are Zero. Ciprian M. Crainiceanu. Department of Statistical Science
1 Likelihood Ratio Tests that Certain Variance Components Are Zero Ciprian M. Crainiceanu Department of Statistical Science www.people.cornell.edu/pages/cmc59 Work done jointly with David Ruppert, School
More informationA nonparametric two-sample wald test of equality of variances
University of Wollongong Research Online Faculty of Informatics - Papers (Archive) Faculty of Engineering and Information Sciences 211 A nonparametric two-sample wald test of equality of variances David
More informationCanonical Correlation Analysis of Longitudinal Data
Biometrics Section JSM 2008 Canonical Correlation Analysis of Longitudinal Data Jayesh Srivastava Dayanand N Naik Abstract Studying the relationship between two sets of variables is an important multivariate
More informationHypothesis Testing One Sample Tests
STATISTICS Lecture no. 13 Department of Econometrics FEM UO Brno office 69a, tel. 973 442029 email:jiri.neubauer@unob.cz 12. 1. 2010 Tests on Mean of a Normal distribution Tests on Variance of a Normal
More informationLarge sample covariance matrices and the T 2 statistic
Large sample covariance matrices and the T 2 statistic EURANDOM, the Netherlands Joint work with W. Zhou Outline 1 2 Basic setting Let {X ij }, i, j =, be i.i.d. r.v. Write n s j = (X 1j,, X pj ) T and
More informationLARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS. S. G. Bobkov and F. L. Nazarov. September 25, 2011
LARGE DEVIATIONS OF TYPICAL LINEAR FUNCTIONALS ON A CONVEX BODY WITH UNCONDITIONAL BASIS S. G. Bobkov and F. L. Nazarov September 25, 20 Abstract We study large deviations of linear functionals on an isotropic
More informationBayesian Econometrics
Bayesian Econometrics Christopher A. Sims Princeton University sims@princeton.edu September 20, 2016 Outline I. The difference between Bayesian and non-bayesian inference. II. Confidence sets and confidence
More informationNeural Network Training
Neural Network Training Sargur Srihari Topics in Network Training 0. Neural network parameters Probabilistic problem formulation Specifying the activation and error functions for Regression Binary classification
More informationEigenvalues and eigenfunctions of the Laplacian. Andrew Hassell
Eigenvalues and eigenfunctions of the Laplacian Andrew Hassell 1 2 The setting In this talk I will consider the Laplace operator,, on various geometric spaces M. Here, M will be either a bounded Euclidean
More informationGaussian Measure of Sections of convex bodies
Gaussian Measure of Sections of convex bodies A. Zvavitch Department of Mathematics, University of Missouri, Columbia, MO 652, USA Abstract In this paper we study properties of sections of convex bodies
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationOn Spatial Involute Gearing
6 th International Conference on Applied Informatics Eger, Hungary, January 27 31, 2004. On Spatial Involute Gearing Hellmuth Stachel Institute of Discrete Mathematics and Geometry, Vienna University of
More informationSTAT 730 Chapter 5: Hypothesis Testing
STAT 730 Chapter 5: Hypothesis Testing Timothy Hanson Department of Statistics, University of South Carolina Stat 730: Multivariate Analysis 1 / 28 Likelihood ratio test def n: Data X depend on θ. The
More informationTesting equality of two mean vectors with unequal sample sizes for populations with correlation
Testing equality of two mean vectors with unequal sample sizes for populations with correlation Aya Shinozaki Naoya Okamoto 2 and Takashi Seo Department of Mathematical Information Science Tokyo University
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 informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
More information