BZPEERSource author record
Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.
Catalog footprint
Published work
We show that the cohomology intersection number of a twisted Gauss-Manin connection with regularization condition is a rational function. As an application, we obtain a new quadratic relation associated to period integrals of a certain family of K3 surfaces.
We give an approximate formula for the distribution of the largest eigenvalue of real Wishart matrices by the expected Euler characteristic method for the general dimension. The formula is expressed in terms of a definite integral with parameters. We derive a differential equation satisfied by the integral for the 2×2 matrix case and perform a numerical analysis of it.
We give a bijection between a quotient space of the parameters and the space of moments for any A-hypergeometric distribution. An algorithmic method to compute the inverse image of the map is proposed utilizing the holonomic gradient method and an asymptotic equivalence of the map and the iterative proportional scaling. The algorithm gives a method to solve a conditional maximum likelihood estimation problem in statistics. Our interplay between the theory of hypergeometric functions and statistics gives some new formulas of A-hypergeometric polynomials.
We give two efficient methods to derive Pfaffian systems for A-hypergeometric systems for the application to the holonomic gradient method for statistics. We utilize the Hilbert driven Buchberger algorithm and Macaulay type matrices in the two methods.
This is the third revision. We study bases of Pfaffian systems for A-hypergeometric system. Gröbner deformations give bases. These bases also give those for twisted cohomology groups. For hypergeometric system associated to a class of order polytopes, these bases have a combinatorial description. The size of the bases associated to a subclass of the order polytopes have the growth rate of the polynomial order. Bases associated to two chain posets and bouquets are studied.
We propose an accelerated version of the holonomic gradient descent and apply it to calculating the maximum likelihood estimate (MLE) of the Fisher-Bingham distribution on a d-dimensional sphere. We derive a Pfaffian system (an integrable connection) and a series expansion associated with the normalizing constant with an error estimation. These enable us to solve some MLE problems up to dimension d=7 with a specified accuracy.
We apply the holonomic gradient method introduced by Nakayama et al.(2011) to the evaluation of the exact distribution function of the largest root of a Wishart matrix, which involves a hypergeometric function 1F1 of a matrix argument. Numerical evaluation of the hypergeometric function has been one of the longstanding problems in multivariate distribution theory. The holonomic gradient method offers a totally new approach, which is complementary to the infinite series expansion around the origin in terms of zonal polynomials. It allows us to move away from the origin by the use of partial differential equations satisfied by the hypergeometric function. From numerical viewpoint we show that the method works well up to dimension 10. From theoretical viewpoint the method offers many challenging problems both to statistics and D-module theory.
A modified A-hypergeometric system is a system of differential equations for the function f(tw⋅x) where f(y) is a solution of an A-hypergeometric system in n variables and w is an n dimensional integer vector, which is called the weight vector. We study the irregularity of modified systems by adapting to this case the notion of umbrella introduced by M. Schulze and U. Walther. Especially, we study slopes and Gevrey series solutions. We develop some applications of this study. Under some conditions we give Laplace integral representations of divergent series solutions of the modified system and we show that certain Gevrey series solutions of the original A-hypergeometric system along coordinate varieties are Gevrey asymptotic expansions of holomorphic solutions of the A-hypergeometric system.
We study properties of Fisher distribution (von Mises-Fisher distribution, matrix Langevin distribution) on the rotation group SO(3). In particular we apply the holonomic gradient descent, introduced by Nakayama et al. (2011), and a method of series expansion for evaluating the normalizing constant of the distribution and for computing the maximum likelihood estimate. The rotation group can be identified with the Stiefel manifold of two orthonormal vectors. Therefore from the viewpoint of statistical modeling, it is of interest to compare Fisher distributions on these manifolds. We illustrate the difference with an example of near-earth objects data.
The Fisher-Bingham system is a system of linear partial differential equations satisfied by the Fisher-Bingham integral for the n-dimensional sphere Sn. The system is given in [Nakayama et al. (2011), Theorem 2] and it is shown that it is a holonomic system [Koyama]. We show that the holonomic rank of the system is equal to 2n+2.
The incomplete beta function is an important special function in statistics. In modern theory of hypergeometric functions, we regard hypergeometric functions as pairings of twisted cycles and twisted cocycles. However, the incomplete beta function cannot be understood in this scheme; in other words, the domain of the integration is not cycle (incomplete). We will generalize the theory of A-hypergeometric systems for incomplete functions. We give a general study as well as a detailed study on an incomplete Gauss hypergeometric function.
We determine the singular locus of the Lauricella function FC by utilizing the theory of D-modules and Gröbner basis. The A-hypergeometric system associated to FC is also discussed.
We give an approximate algorithm of computing holonomic systems of linear differential equations for definite integrals with parameters. We show that this algorithm gives a correct answer in finite steps, but we have no general stopping condition. We apply the approximate method to find differential equations for integrals associated to smooth Fano polytopes. They are interested in the study of K3 surfaces and the toric mirror symmetry. In this class of integrals, we can apply Stienstra's rank formula to our algorithm, which gives a stopping condition of the approximate algorithm.
We give a new algorithm to find local maximum and minimum of a holonomic function and apply it for the Fisher-Bingham integral on the sphere Sn, which is used in the directional statistics. The method utilizes the theory and algorithms of holonomic systems.
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.