Showing 1–50 of 112 results for author: Takemura, A

MOD-CL: Multi-label Object Detection with Constrained Loss

Authors: Sota Moriyama, Koji Watanabe, Katsumi Inoue, Akihiro Takemura

Abstract: We introduce MOD-CL, a multi-label object detection framework that utilizes constrained loss in the training process to produce outputs that better satisfy the given requirements. In this paper, we use $\mathrm{MOD_{YOLO}}$, a multi-label object detection model built upon the state-of-the-art object detection model YOLOv8, which has been published in recent years. In Task 1, we introduce the Corre… ▽ More We introduce MOD-CL, a multi-label object detection framework that utilizes constrained loss in the training process to produce outputs that better satisfy the given requirements. In this paper, we use $\mathrm{MOD_{YOLO}}$, a multi-label object detection model built upon the state-of-the-art object detection model YOLOv8, which has been published in recent years. In Task 1, we introduce the Corrector Model and Blender Model, two new models that follow after the object detection process, aiming to generate a more constrained output. For Task 2, constrained losses have been incorporated into the $\mathrm{MOD_{YOLO}}$ architecture using Product T-Norm. The results show that these implementations are instrumental to improving the scores for both Task 1 and Task 2. △ Less

Submitted 31 January, 2024; originally announced March 2024.

Towards end-to-end ASP computation

Authors: Taisuke Sato, Akihiro Takemura, Katsumi Inoue

Abstract: We propose an end-to-end approach for answer set programming (ASP) and linear algebraically compute stable models satisfying given constraints. The idea is to implement Lin-Zhao's theorem \cite{Lin04} together with constraints directly in vector spaces as numerical minimization of a cost function constructed from a matricized normal logic program, loop formulas in Lin-Zhao's theorem and constraint… ▽ More We propose an end-to-end approach for answer set programming (ASP) and linear algebraically compute stable models satisfying given constraints. The idea is to implement Lin-Zhao's theorem \cite{Lin04} together with constraints directly in vector spaces as numerical minimization of a cost function constructed from a matricized normal logic program, loop formulas in Lin-Zhao's theorem and constraints, thereby no use of symbolic ASP or SAT solvers involved in our approach. We also propose precomputation that shrinks the program size and heuristics for loop formulas to reduce computational difficulty. We empirically test our approach with programming examples including the 3-coloring and Hamiltonian cycle problems. As our approach is purely numerical and only contains vector/matrix operations, acceleration by parallel technologies such as many-cores and GPUs is expected. △ Less

Submitted 13 June, 2023; v1 submitted 11 June, 2023; originally announced June 2023.

Comments: 29 pages, 9 figures

ACM Class: I.2.4

Non-minimaxity of debiased shrinkage estimators

Authors: Yuzo Maruyama, Akimichi Takemura

Abstract: We consider the estimation of the $p$-variate normal mean of $X\sim N_p(θしーた,I)$ under the quadratic loss function. We investigate the decision theoretic properties of debiased shrinkage estimator, the estimator which shrinks towards the origin for smaller $\|x\|^2$ and which is exactly equal to the unbiased estimator $X$ for larger $\|x\|^2$. Such debiased shrinkage estimator seems superior to the u… ▽ More We consider the estimation of the $p$-variate normal mean of $X\sim N_p(θしーた,I)$ under the quadratic loss function. We investigate the decision theoretic properties of debiased shrinkage estimator, the estimator which shrinks towards the origin for smaller $\|x\|^2$ and which is exactly equal to the unbiased estimator $X$ for larger $\|x\|^2$. Such debiased shrinkage estimator seems superior to the unbiased estimator $X$, which implies minimaxity. However we show that it is not minimax under mild conditions. △ Less

Submitted 7 June, 2023; originally announced June 2023.

MSC Class: 62C20

Generating Explainable Rule Sets from Tree-Ensemble Learning Methods by Answer Set Programming

Authors: Akihiro Takemura, Katsumi Inoue

Abstract: We propose a method for generating explainable rule sets from tree-ensemble learners using Answer Set Programming (ASP). To this end, we adopt a decompositional approach where the split structures of the base decision trees are exploited in the construction of rules, which in turn are assessed using pattern mining methods encoded in ASP to extract interesting rules. We show how user-defined constr… ▽ More We propose a method for generating explainable rule sets from tree-ensemble learners using Answer Set Programming (ASP). To this end, we adopt a decompositional approach where the split structures of the base decision trees are exploited in the construction of rules, which in turn are assessed using pattern mining methods encoded in ASP to extract interesting rules. We show how user-defined constraints and preferences can be represented declaratively in ASP to allow for transparent and flexible rule set generation, and how rules can be used as explanations to help the user better understand the models. Experimental evaluation with real-world datasets and popular tree-ensemble algorithms demonstrates that our approach is applicable to a wide range of classification tasks. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: In Proceedings ICLP 2021, arXiv:2109.07914

Journal ref: EPTCS 345, 2021, pp. 127-140

The volume-of-tube method for Gaussian random fields with inhomogeneous variance

Authors: Satoshi Kuriki, Akimichi Takemura, Jonathan E. Taylor

Abstract: The tube method or the volume-of-tube method approximates the tail probability of the maximum of a smooth Gaussian random field with zero mean and unit variance. This method evaluates the volume of a spherical tube about the index set, and then transforms it to the tail probability. In this study, we generalize the tube method to a case in which the variance is not constant. We provide the volume… ▽ More The tube method or the volume-of-tube method approximates the tail probability of the maximum of a smooth Gaussian random field with zero mean and unit variance. This method evaluates the volume of a spherical tube about the index set, and then transforms it to the tail probability. In this study, we generalize the tube method to a case in which the variance is not constant. We provide the volume formula for a spherical tube with a non-constant radius in terms of curvature tensors, and the tail probability formula of the maximum of a Gaussian random field with inhomogeneous variance, as well as its Laplace approximation. In particular, the critical radius of the tube is generalized for evaluation of the asymptotic approximation error. As an example, we discuss the approximation of the largest eigenvalue distribution of the Wishart matrix with a non-identity matrix parameter. The Bonferroni method is the tube method when the index set is a finite set. We provide the formula for the asymptotic approximation error for the Bonferroni method when the variance is not constant. △ Less

Submitted 9 September, 2021; v1 submitted 4 August, 2021; originally announced August 2021.

Comments: 30 pages, 3 figures

MSC Class: 62H10 (Primary); 60G60 (Secondary)

Game-theoretic derivation of upper hedging prices of multivariate contingent claims and submodularity

Authors: Takeru Matsuda, Akimichi Takemura

Abstract: We investigate upper and lower hedging prices of multivariate contingent claims from the viewpoint of game-theoretic probability and submodularity. By considering a game between "Market" and "Investor" in discrete time, the pricing problem is reduced to a backward induction of an optimization over simplexes. For European options with payoff functions satisfying a combinatorial property called subm… ▽ More We investigate upper and lower hedging prices of multivariate contingent claims from the viewpoint of game-theoretic probability and submodularity. By considering a game between "Market" and "Investor" in discrete time, the pricing problem is reduced to a backward induction of an optimization over simplexes. For European options with payoff functions satisfying a combinatorial property called submodularity or supermodularity, this optimization is solved in closed form by using the Lovász extension and the upper and lower hedging prices can be calculated efficiently. This class includes the options on the maximum or the minimum of several assets. We also study the asymptotic behavior as the number of game rounds goes to infinity. The upper and lower hedging prices of European options converge to the solutions of the Black-Scholes-Barenblatt equations. For European options with submodular or supermodular payoff functions, the Black-Scholes-Barenblatt equation is reduced to the linear Black-Scholes equation and it is solved in closed form. Numerical results show the validity of the theoretical results. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Journal ref: Japan Journal of Industrial and Applied Mathematics, 37, 213--248, 2020

Distribution of ratio of two Wishart matrices and evaluation of cumulative probability by holonomic gradient method

Authors: Hiroki Hashiguchi, Nobuki Takayama, Akimichi Takemura

Abstract: We study the distribution of the ratio of two central Wishart matrices with different covariance matrices. We first derive the density function of a particular matrix form of the ratio and show that its cumulative distribution function can be expressed in terms of the hypergeometric function 2F1 of a matrix argument. Then we apply the holonomic gradient method for numerical evaluation of the hyper… ▽ More We study the distribution of the ratio of two central Wishart matrices with different covariance matrices. We first derive the density function of a particular matrix form of the ratio and show that its cumulative distribution function can be expressed in terms of the hypergeometric function 2F1 of a matrix argument. Then we apply the holonomic gradient method for numerical evaluation of the hypergeometric function. This approach enables us to compute the power function of Roy's maximum root test for testing the equality of two covariance matrices. △ Less

Submitted 28 October, 2016; originally announced October 2016.

MSC Class: 62H10; 13N10

Journal ref: Journal of Multivariate Analysis 165 (2018) 270-278

Differential relations for the largest root distribution of complex non-central Wishart matrices

Authors: Raimundas Vidunas, Akimichi Takemura

Abstract: A holonomic system for the probability density function of the largest eigenvalue of a non-central complex Wishart distribution with identity covariance matrix is derived. Furthermore a new determinantal formula for the probability density function is derived (for m=2,3) or conjectured. A holonomic system for the probability density function of the largest eigenvalue of a non-central complex Wishart distribution with identity covariance matrix is derived. Furthermore a new determinantal formula for the probability density function is derived (for m=2,3) or conjectured. △ Less

Submitted 6 September, 2016; originally announced September 2016.

Comments: 23 pages

MSC Class: 60E05; 13N10; 62H10

Relation between the rate of convergence of strong law of large numbers and the rate of concentration of Bayesian prior in game-theoretic probability

Authors: Ryosuke Sato, Kenshi Miyabe, Akimichi Takemura

Abstract: We study the behavior of the capital process of a continuous Bayesian mixture of fixed proportion betting strategies in the one-sided unbounded forecasting game in game-theoretic probability. We establish the relation between the rate of convergence of the strong law of large numbers in the self-normalized form and the rate of divergence to infinity of the prior density around the origin. In parti… ▽ More We study the behavior of the capital process of a continuous Bayesian mixture of fixed proportion betting strategies in the one-sided unbounded forecasting game in game-theoretic probability. We establish the relation between the rate of convergence of the strong law of large numbers in the self-normalized form and the rate of divergence to infinity of the prior density around the origin. In particular we present prior densities ensuring the validity of Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm. △ Less

Submitted 26 April, 2016; originally announced April 2016.

MSC Class: 60G42

Journal ref: Stochastic Processes and their Applications 128 (2018) 1466-1484

Chi-Square Mixture Representations for the Distribution of the Scalar Schur Complement in a Noncentral Wishart Matrix

Authors: Constantin Siriteanu, Satoshi Kuriki, Donald Richards, Akimichi Takemura

Abstract: We show that the distribution of the scalar Schur complement in a noncentral Wishart matrix is a mixture of central chi-square distributions with different degrees of freedom. For the case of a rank-1 noncentrality matrix, the weights of the mixture representation arise from a noncentral beta mixture of Poisson distributions. We show that the distribution of the scalar Schur complement in a noncentral Wishart matrix is a mixture of central chi-square distributions with different degrees of freedom. For the case of a rank-1 noncentrality matrix, the weights of the mixture representation arise from a noncentral beta mixture of Poisson distributions. △ Less

Submitted 26 December, 2015; originally announced December 2015.

Journal ref: Statistics and Probability Letters 115 (2016) 79-87

Exponential decay rate of partial autocorrelation coefficients of ARMA and short-memory processes

Authors: Akimichi Takemura

Abstract: We present a short proof of the fact that the exponential decay rate of partial autocorrelation coefficients of a short-memory process, in particular an ARMA process, is equal to the exponential decay rate of the coefficients of its infinite autoregressive representation. We present a short proof of the fact that the exponential decay rate of partial autocorrelation coefficients of a short-memory process, in particular an ARMA process, is equal to the exponential decay rate of the coefficients of its infinite autoregressive representation. △ Less

Submitted 4 January, 2016; v1 submitted 22 November, 2015; originally announced November 2015.

MSC Class: 60G12

Journal ref: Statistics and Probability Letters 110 (2016) 207-210

$A$-Hypergeometric Distributions and Newton Polytopes

Authors: Nobuki Takayama, Satoshi Kuriki, Akimichi Takemura

Abstract: 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 estim… ▽ More 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. △ Less

Submitted 11 November, 2015; v1 submitted 8 October, 2015; originally announced October 2015.

MSC Class: 33C65; 33C90; 62H12; 62H17

Exact ZF Analysis and Computer-Algebra-Aided Evaluation in Rank-1 LoS Rician Fading

Authors: Constantin Siriteanu, Akimichi Takemura, Christoph Koutschan, Satoshi Kuriki, Donald St. P. Richards, Hyundong Shin

Abstract: We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distr… ▽ More We study zero-forcing detection (ZF) for multiple-input/multiple-output (MIMO) spatial multiplexing under transmit-correlated Rician fading for an N_R X N_T channel matrix with rank-1 line-of-sight (LoS) component. By using matrix transformations and multivariate statistics, our exact analysis yields the signal-to-noise ratio moment generating function (m.g.f.) as an infinite series of gamma distribution m.g.f.'s and analogous series for ZF performance measures, e.g., outage probability and ergodic capacity. However, their numerical convergence is inherently problematic with increasing Rician K-factor, N_R , and N_T. We circumvent this limitation as follows. First, we derive differential equations satisfied by the performance measures with a novel automated approach employing a computer-algebra tool which implements Groebner basis computation and creative telescoping. These differential equations are then solved with the holonomic gradient method (HGM) from initial conditions computed with the infinite series. We demonstrate that HGM yields more reliable performance evaluation than by infinite series alone and more expeditious than by simulation, for realistic values of K , and even for N_R and N_T relevant to large MIMO systems. We envision extending the proposed approaches for exact analysis and reliable evaluation to more general Rician fading and other transceiver methods. △ Less

Submitted 19 May, 2016; v1 submitted 24 July, 2015; originally announced July 2015.

Comments: Accepted for publication by the IEEE Transactions on Wireless Communications, on April 7th, 2016; this is the final revision before publication

Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for self-normalized martingales: a game-theoretic approach

Authors: Takeyuki Sasai, Kenshi Miyabe, Akimichi Takemura

Abstract: We prove an Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for self-normalized martingales. Our proof is given in the framework of the game-theoretic probability of Shafer and Vovk. As many other game-theoretic proofs, our proof is self-contained and explicit. We prove an Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for self-normalized martingales. Our proof is given in the framework of the game-theoretic probability of Shafer and Vovk. As many other game-theoretic proofs, our proof is self-contained and explicit. △ Less

Submitted 24 April, 2015; originally announced April 2015.

MSC Class: 60G42

Journal ref: The Annals of Probability 47 (2019) 1136-1161

Holonomic gradient method for distribution function of a weighted sum of noncentral chi-square random variables

Authors: Tamio Koyama, Akimichi Takemura

Abstract: We apply the holonomic gradient method to compute the distribution function of a weighted sum of independent noncentral chi-square random variables. It is the distribution function of the squared length of a multivariate normal random vector. We treat this distribution as an integral of the normalizing constant of the Fisher-Bingham distribution on the unit sphere and make use of the partial diffe… ▽ More We apply the holonomic gradient method to compute the distribution function of a weighted sum of independent noncentral chi-square random variables. It is the distribution function of the squared length of a multivariate normal random vector. We treat this distribution as an integral of the normalizing constant of the Fisher-Bingham distribution on the unit sphere and make use of the partial differential equations for the Fisher-Bingham distribution. △ Less

Submitted 29 August, 2015; v1 submitted 1 March, 2015; originally announced March 2015.

MSC Class: 32C38; 16S32; 62H10

A power-law decay model with autocorrelation for posting data to social networking services

Authors: Toshifumi Fujiyama, Chihiro Matsui, Akimichi Takemura

Abstract: We propose a power-law decay model with autocorrelation for posting data to social networking services concerning particular events such as national holidays or major sport events. In these kinds of events we observe people's interest both before and after the events. In our model the number of postings has a Poisson distribution whose expected value decays as a power law. Our model also incorpora… ▽ More We propose a power-law decay model with autocorrelation for posting data to social networking services concerning particular events such as national holidays or major sport events. In these kinds of events we observe people's interest both before and after the events. In our model the number of postings has a Poisson distribution whose expected value decays as a power law. Our model also incorporates autocorrelations by autoregressive specification of the expected value. We show that our proposed model well fits the data from social networking services. △ Less

Submitted 9 November, 2014; originally announced November 2014.

MSC Class: 62M10

Journal ref: PLoS ONE 11(8): e0160592. 2016

An objective look at obtaining the plotting positions for QQ-plots

Authors: Reza Hosseini, Akimichi Takemura

Abstract: Choosing the plotting positions for the QQ-plot has been a subject of much debate in the statistical and engineering literature. This paper looks at this problem objectively by considering three frameworks: distribution-theoretic; decision-theoretic; game-theoretic. In each framework, we derive the plotting positions and show that there are more than one legitimate solution depending on the practi… ▽ More Choosing the plotting positions for the QQ-plot has been a subject of much debate in the statistical and engineering literature. This paper looks at this problem objectively by considering three frameworks: distribution-theoretic; decision-theoretic; game-theoretic. In each framework, we derive the plotting positions and show that there are more than one legitimate solution depending on the practitioner's objective. This work clarifies the choice of the plotting positions by allowing one to easily find the mathematical equivalent of their view and choose the corresponding solution. This work also discusses approximations to the plotting positions when no closed form is available. △ Less

Submitted 24 September, 2014; originally announced September 2014.

A game-theoretic proof of Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for fair-coin tossing

Authors: Takeyuki Sasai, Kenshi Miyabe, Akimichi Takemura

Abstract: We give a game-theoretic proof of the celebrated Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for fair coin tossing. Our proof, based on Bayesian strategy, is explicit as many other game-theoretic proofs of the laws in probability theory. We give a game-theoretic proof of the celebrated Erdos-Feller-Kolmogorov-Petrowsky law of the iterated logarithm for fair coin tossing. Our proof, based on Bayesian strategy, is explicit as many other game-theoretic proofs of the laws in probability theory. △ Less

Submitted 9 November, 2014; v1 submitted 8 August, 2014; originally announced August 2014.

MSC Class: 60G42

A framework for fitting sparse data

Authors: Reza Hosseini, Akimichi Takemura, Kiros Berhane

Abstract: This paper develops a framework for fitting functions with domains in the Euclidean space, when data are sparse but a slow variation allows for a useful fit. We measure the variation by Lipschitz Bound (LB). Functions which admit smaller LB are considered to vary more slowly. Since most functions in practice are wiggly and do not admit a small LB, we extend this framework by approximating a wiggly… ▽ More This paper develops a framework for fitting functions with domains in the Euclidean space, when data are sparse but a slow variation allows for a useful fit. We measure the variation by Lipschitz Bound (LB). Functions which admit smaller LB are considered to vary more slowly. Since most functions in practice are wiggly and do not admit a small LB, we extend this framework by approximating a wiggly function, f, by ones which admit a smaller LB and do not deviate from f by more than a specified Bound Deviation (BD). In fact for any positive LB, one can find such a BD, thus defining a trade-off function (LB-BD function) between the variation measure (LB) and the deviation measure (BD). We show that the LB-BD function satisfies nice properties: it is non-increasing and convex. We also present a method to obtain it using convex optimization. For a function with given LB and BD, we find the optimal fit and present deterministic bounds for the prediction error of various methods. Given the LB-BD function, we discuss picking an appropriate LB-BD pair for fitting and calculating the prediction errors. The developed methods can naturally accommodate an extra assumption of periodicity to obtain better prediction errors. Finally we present the application of this framework to air pollution data with sparse observations over time. △ Less

Submitted 4 July, 2014; originally announced July 2014.

Conformal Geometry of Sequential Test in Multidimensional Curved Exponential Family

Authors: Masayuki Kumon, Akimichi Takemura, Kei Takeuchi

Abstract: This article presents a differential geometrical method for analyzing sequential test procedures. It is based on the primal result on the conformal geometry of statistical manifold developed in Kumon, Takemura and Takeuchi (2011). By introducing curvature-type random variables, the condition is first clarified for a statistical manifold to be an exponential family under an appropriate sequential t… ▽ More This article presents a differential geometrical method for analyzing sequential test procedures. It is based on the primal result on the conformal geometry of statistical manifold developed in Kumon, Takemura and Takeuchi (2011). By introducing curvature-type random variables, the condition is first clarified for a statistical manifold to be an exponential family under an appropriate sequential test procedure. This result is further elaborated for investigating the efficient sequential test in a multidimensional curved exponential family. The theoretical results are numerically examined by using von Mises-Fisher and hyperboloid models. △ Less

Submitted 30 June, 2014; originally announced June 2014.

Journal ref: Sequential Analysis 35 (2016) 30-68

Properties of powers of functions satisfying second-order linear differential equations with applications to statistics

Authors: Naoki Marumo, Toshinori Oaku, Akimichi Takemura

Abstract: We derive properties of powers of a function satisfying a second-order linear differential equation. In particular we prove that the n-th power of the function satisfies an (n+1)-th order differential equation and give a simple method for obtaining the differential equation. Also we determine the exponents of the differential equation and derive a bound for the degree of the polynomials, which are… ▽ More We derive properties of powers of a function satisfying a second-order linear differential equation. In particular we prove that the n-th power of the function satisfies an (n+1)-th order differential equation and give a simple method for obtaining the differential equation. Also we determine the exponents of the differential equation and derive a bound for the degree of the polynomials, which are coefficients in the differential equation. The bound corresponds to the order of differential equation satisfied by the n-fold convolution of the Fourier transform of the function. These results are applied to some probability density functions used in statistics. △ Less

Submitted 17 May, 2014; originally announced May 2014.

Journal ref: Japan Journal of Industrial and Applied Mathematics 32 (2015) 553-572

Markov degree of configurations defined by fibers of a configuration

Authors: Takayuki Koyama, Mitsunori Ogawa, Akimichi Takemura

Abstract: We consider a series of configurations defined by fibers of a given base configuration. We prove that Markov degree of the configurations is bounded from above by the Markov complexity of the base configuration. As important examples of base configurations we consider incidence matrices of graphs and study the maximum Markov degree of configurations defined by fibers of the incidence matrices. In… ▽ More We consider a series of configurations defined by fibers of a given base configuration. We prove that Markov degree of the configurations is bounded from above by the Markov complexity of the base configuration. As important examples of base configurations we consider incidence matrices of graphs and study the maximum Markov degree of configurations defined by fibers of the incidence matrices. In particular we give a proof that the Markov degree for two-way transportation polytopes is three. △ Less

Submitted 25 July, 2014; v1 submitted 12 May, 2014; originally announced May 2014.

Comments: 28 pages

MSC Class: 05B20

Estimation of exponential-polynomial distribution by holonomic gradient descent

Authors: Jumpei Hayakawa, Akimichi Takemura

Abstract: We study holonomic gradient decent for maximum likelihood estimation of exponential-polynomial distribution, whose density is the exponential function of a polynomial in the random variable. We first consider the case that the support of the distribution is the set of positive reals. We show that the maximum likelihood estimate (MLE) can be easily computed by the holonomic gradient descent, even t… ▽ More We study holonomic gradient decent for maximum likelihood estimation of exponential-polynomial distribution, whose density is the exponential function of a polynomial in the random variable. We first consider the case that the support of the distribution is the set of positive reals. We show that the maximum likelihood estimate (MLE) can be easily computed by the holonomic gradient descent, even though the normalizing constant of this family does not have a closed-form expression and discuss determination of the degree of the polynomial based on the score test statistic. Then we present extensions to the whole real line and to the bivariate distribution on the positive orthant. △ Less

Submitted 15 September, 2014; v1 submitted 30 March, 2014; originally announced March 2014.

MIMO Zero-Forcing Performance Evaluation Using the Holonomic Gradient Method

Authors: Constantin Siriteanu, Akimichi Takemura, Satoshi Kuriki, Hyundong Shin, Christoph Koutschan

Abstract: For multiple-input multiple-output (MIMO) spatial-multiplexing transmission, zero-forcing detection (ZF) is appealing because of its low complexity. Our recent MIMO ZF performance analysis for Rician--Rayleigh fading, which is relevant in heterogeneous networks, has yielded for the ZF outage probability and ergodic capacity infinite-series expressions. Because they arose from expanding the conflue… ▽ More For multiple-input multiple-output (MIMO) spatial-multiplexing transmission, zero-forcing detection (ZF) is appealing because of its low complexity. Our recent MIMO ZF performance analysis for Rician--Rayleigh fading, which is relevant in heterogeneous networks, has yielded for the ZF outage probability and ergodic capacity infinite-series expressions. Because they arose from expanding the confluent hypergeometric function $ {_1\! F_1} (\cdot, \cdot, σしぐま) $ around 0, they do not converge numerically at realistically-high Rician $ K $-factor values. Therefore, herein, we seek to take advantage of the fact that $ {_1\! F_1} (\cdot, \cdot, σしぐま) $ satisfies a differential equation, i.e., it is a \textit{holonomic} function. Holonomic functions can be computed by the \textit{holonomic gradient method} (HGM), i.e., by numerically solving the satisfied differential equation. Thus, we first reveal that the moment generating function (m.g.f.) and probability density function (p.d.f.) of the ZF signal-to-noise ratio (SNR) are holonomic. Then, from the differential equation for $ {_1\! F_1} (\cdot, \cdot, σしぐま) $, we deduce those satisfied by the SNR m.g.f. and p.d.f., and demonstrate that the HGM helps compute the p.d.f. accurately at practically-relevant values of $ K $. Finally, numerical integration of the SNR p.d.f. produced by HGM yields accurate ZF outage probability and ergodic capacity results. △ Less

Submitted 15 April, 2015; v1 submitted 15 March, 2014; originally announced March 2014.

Comments: This manuscript was accepted in December 2014

Journal ref: IEEE Transactions on Wireless Communications, vol. 14, no. 4, April 2015, pp. 2322-2335

Derandomization in Game-theoretic Probability

Authors: Kenshi Miyabe, Akimichi Takemura

Abstract: We give a general method for constructing a deterministic strategy of Reality from a randomized strategy in game-theoretic probability. The construction can be seen as derandomization in game-theoretic probability. We give a general method for constructing a deterministic strategy of Reality from a randomized strategy in game-theoretic probability. The construction can be seen as derandomization in game-theoretic probability. △ Less

Submitted 4 August, 2014; v1 submitted 22 February, 2014; originally announced February 2014.

Comments: 19 pages

MSC Class: 60G42

Schur Complement Based Analysis of MIMO Zero-Forcing for Rician Fading

Authors: Constantin Siriteanu, Akimichi Takemura, Satoshi Kuriki, Donald St. P. Richards, Hyundong Shin

Abstract: For multiple-input/multiple-output (MIMO) spatial multiplexing with zero-forcing detection (ZF), signal-to-noise ratio (SNR) analysis for Rician fading involves the cumbersome noncentral-Wishart distribution (NCWD) of the transmit sample-correlation (Gramian) matrix. An \textsl{approximation} with a \textsl{virtual} CWD previously yielded for the ZF SNR an approximate (virtual) Gamma distribution.… ▽ More For multiple-input/multiple-output (MIMO) spatial multiplexing with zero-forcing detection (ZF), signal-to-noise ratio (SNR) analysis for Rician fading involves the cumbersome noncentral-Wishart distribution (NCWD) of the transmit sample-correlation (Gramian) matrix. An \textsl{approximation} with a \textsl{virtual} CWD previously yielded for the ZF SNR an approximate (virtual) Gamma distribution. However, analytical conditions qualifying the accuracy of the SNR-distribution approximation were unknown. Therefore, we have been attempting to exactly characterize ZF SNR for Rician fading. Our previous attempts succeeded only for the sole Rician-fading stream under Rician--Rayleigh fading, by writing it as scalar Schur complement (SC) in the Gramian. Herein, we pursue a more general, matrix-SC-based analysis to characterize SNRs when several streams may undergo Rician fading. On one hand, for full-Rician fading, the SC distribution is found to be exactly a CWD if and only if a channel-mean--correlation \textsl{condition} holds. Interestingly, this CWD then coincides with the \textsl{virtual} CWD ensuing from the \textsl{approximation}. Thus, under the \textsl{condition}, the actual and virtual SNR-distributions coincide. On the other hand, for Rician--Rayleigh fading, the matrix-SC distribution is characterized in terms of determinant of matrix with elementary-function entries, which also yields a new characterization of the ZF SNR. Average error probability results validate our analysis vs.~simulation. △ Less

Submitted 26 September, 2014; v1 submitted 2 January, 2014; originally announced January 2014.

Comments: 32 pages, 4 figures, 1 table

Optimality of Thompson Sampling for Gaussian Bandits Depends on Priors

Authors: Junya Honda, Akimichi Takemura

Abstract: In stochastic bandit problems, a Bayesian policy called Thompson sampling (TS) has recently attracted much attention for its excellent empirical performance. However, the theoretical analysis of this policy is difficult and its asymptotic optimality is only proved for one-parameter models. In this paper we discuss the optimality of TS for the model of normal distributions with unknown means and va… ▽ More In stochastic bandit problems, a Bayesian policy called Thompson sampling (TS) has recently attracted much attention for its excellent empirical performance. However, the theoretical analysis of this policy is difficult and its asymptotic optimality is only proved for one-parameter models. In this paper we discuss the optimality of TS for the model of normal distributions with unknown means and variances as one of the most fundamental example of multiparameter models. First we prove that the expected regret of TS with the uniform prior achieves the theoretical bound, which is the first result to show that the asymptotic bound is achievable for the normal distribution model. Next we prove that TS with Jeffreys prior and reference prior cannot achieve the theoretical bound. Therefore the choice of priors is important for TS and non-informative priors are sometimes risky in cases of multiparameter models. △ Less

Submitted 8 November, 2013; originally announced November 2013.

Exact MIMO Zero-Forcing Detection Analysis for Transmit-Correlated Rician Fading

Authors: Constantin Siriteanu, Steven Blostein, Akimichi Takemura, Hyundong Shin, Shahram Yousefi, Satoshi Kuriki

Abstract: We analyze the performance of multiple input/multiple output (MIMO) communications systems employing spatial multiplexing and zero-forcing detection (ZF). The distribution of the ZF signal-to-noise ratio (SNR) is characterized when either the intended stream or interfering streams experience Rician fading, and when the fading may be correlated on the transmit side. Previously, exact ZF analysis ba… ▽ More We analyze the performance of multiple input/multiple output (MIMO) communications systems employing spatial multiplexing and zero-forcing detection (ZF). The distribution of the ZF signal-to-noise ratio (SNR) is characterized when either the intended stream or interfering streams experience Rician fading, and when the fading may be correlated on the transmit side. Previously, exact ZF analysis based on a well-known SNR expression has been hindered by the noncentrality of the Wishart distribution involved. In addition, approximation with a central-Wishart distribution has not proved consistently accurate. In contrast, the following exact ZF study proceeds from a lesser-known SNR expression that separates the intended and interfering channel-gain vectors. By first conditioning on, and then averaging over the interference, the ZF SNR distribution for Rician-Rayleigh fading is shown to be an infinite linear combination of gamma distributions. On the other hand, for Rayleigh-Rician fading, the ZF SNR is shown to be gamma-distributed. Based on the SNR distribution, we derive new series expressions for the ZF average error probability, outage probability, and ergodic capacity. Numerical results confirm the accuracy of our new expressions, and reveal effects of interference and channel statistics on performance. △ Less

Submitted 2 January, 2014; v1 submitted 10 July, 2013; originally announced July 2013.

Comments: 14 pages, two-colum, 1 table, 10 figures

Report number: METR 2013-07

Markov degree of the Birkhoff model

Authors: Takashi Yamaguchi, Mitsunori Ogawa, Akimichi Takemura

Abstract: We prove the conjecture by Diaconis and Eriksson (2006) that the Markov degree of the Birkhoff model is three. In fact, we prove the conjecture in a generalization of the Birkhoff model, where each voter is asked to rank a fixed number, say r, of candidates among all candidates. We also give an exhaustive characterization of Markov bases for small r. We prove the conjecture by Diaconis and Eriksson (2006) that the Markov degree of the Birkhoff model is three. In fact, we prove the conjecture in a generalization of the Birkhoff model, where each voter is asked to rank a fixed number, say r, of candidates among all candidates. We also give an exhaustive characterization of Markov bases for small r. △ Less

Submitted 18 November, 2013; v1 submitted 3 April, 2013; originally announced April 2013.

MSC Class: 62H17

Journal ref: Journal of Algebraic Combinatorics, 40 (2014) pp. 293-311

Calculation of orthant probabilities by the holonomic gradient method

Authors: Tamio Koyama, Akimichi Takemura

Abstract: We apply the holonomic gradient method (HGM) introduced by [9] to the calculation of orthant probabilities of multivariate normal distribution. The holonomic gradient method applied to orthant probabilities is found to be a variant of Plackett's recurrence relation ([14]). However an implementation of the method yields recurrence relations more suitable for numerical computation than Plackett's re… ▽ More We apply the holonomic gradient method (HGM) introduced by [9] to the calculation of orthant probabilities of multivariate normal distribution. The holonomic gradient method applied to orthant probabilities is found to be a variant of Plackett's recurrence relation ([14]). However an implementation of the method yields recurrence relations more suitable for numerical computation than Plackett's recurrence relation. We derive some theoretical results on the holonomic system for the orthant probabilities. These results show that multivariate normal orthant probabilities possess some remarkable properties from the viewpoint of holonomic systems. Finally we show that numerical performance of our method is comparable or superior compared to existing methods. △ Less

Submitted 18 January, 2013; v1 submitted 29 November, 2012; originally announced November 2012.

Comments: 17 pages

MSC Class: 32C38; 16S32; 62H10

The law of the iterated logarithm in game-theoretic probability with quadratic and stronger hedges

Authors: Kenshi Miyabe, Akimichi Takemura

Abstract: We prove both the validity and the sharpness of the law of the iterated logarithm in game-theoretic probability with quadratic and stronger hedges. We prove both the validity and the sharpness of the law of the iterated logarithm in game-theoretic probability with quadratic and stronger hedges. △ Less

Submitted 13 April, 2013; v1 submitted 24 August, 2012; originally announced August 2012.

Journal ref: Stochastic Processes and their Applications 123 (2013) 3132-3152

Bayesian logistic betting strategy against probability forecasting

Authors: Masayuki Kumon, Jing Li, Akimichi Takemura, Kei Takeuchi

Abstract: We propose a betting strategy based on Bayesian logistic regression modeling for the probability forecasting game in the framework of game-theoretic probability by Shafer and Vovk (2001). We prove some results concerning the strong law of large numbers in the probability forecasting game with side information based on our strategy. We also apply our strategy for assessing the quality of probabilit… ▽ More We propose a betting strategy based on Bayesian logistic regression modeling for the probability forecasting game in the framework of game-theoretic probability by Shafer and Vovk (2001). We prove some results concerning the strong law of large numbers in the probability forecasting game with side information based on our strategy. We also apply our strategy for assessing the quality of probability forecasting by the Japan Meteorological Agency. We find that our strategy beats the agency by exploiting its tendency of avoiding clear-cut forecasts. △ Less

Submitted 16 April, 2012; originally announced April 2012.

MSC Class: 60G42

Journal ref: Stochastic Analysis and Applications 31 (2013) 214-234

Markov degree of the three-state toric homogeneous Markov chain model

Authors: David Haws, Abraham Martín del Campo, Akimichi Takemura, Ruriko Yoshida

Abstract: We consider the three-state toric homogeneous Markov chain model (THMC) without loops and initial parameters. At time $T$, the size of the design matrix is $6 \times 3\cdot 2^{T-1}$ and the convex hull of its columns is the model polytope. We study the behavior of this polytope for $T\geq 3$ and we show that it is defined by 24 facets for all $T\ge 5$. Moreover, we give a complete description of t… ▽ More We consider the three-state toric homogeneous Markov chain model (THMC) without loops and initial parameters. At time $T$, the size of the design matrix is $6 \times 3\cdot 2^{T-1}$ and the convex hull of its columns is the model polytope. We study the behavior of this polytope for $T\geq 3$ and we show that it is defined by 24 facets for all $T\ge 5$. Moreover, we give a complete description of these facets. From this, we deduce that the toric ideal associated with the design matrix is generated by binomials of degree at most 6. Our proof is based on a result due to Sturmfels, who gave a bound on the degree of the generators of a toric ideal, provided the normality of the corresponding toric variety. In our setting, we established the normality of the toric variety associated to the THMC model by studying the geometric properties of the model polytope. △ Less

Submitted 17 September, 2013; v1 submitted 13 April, 2012; originally announced April 2012.

Comments: 26 pages, 1 figure

Finite-time Regret Bound of a Bandit Algorithm for the Semi-bounded Support Model

Authors: Junya Honda, Akimichi Takemura

Abstract: In this paper we consider stochastic multiarmed bandit problems. Recently a policy, DMED, is proposed and proved to achieve the asymptotic bound for the model that each reward distribution is supported in a known bounded interval, e.g. [0,1]. However, the derived regret bound is described in an asymptotic form and the performance in finite time has been unknown. We inspect this policy and derive a… ▽ More In this paper we consider stochastic multiarmed bandit problems. Recently a policy, DMED, is proposed and proved to achieve the asymptotic bound for the model that each reward distribution is supported in a known bounded interval, e.g. [0,1]. However, the derived regret bound is described in an asymptotic form and the performance in finite time has been unknown. We inspect this policy and derive a finite-time regret bound by refining large deviation probabilities to a simple finite form. Further, this observation reveals that the assumption on the lower-boundedness of the support is not essential and can be replaced with a weaker one, the existence of the moment generating function. △ Less

Submitted 17 February, 2012; v1 submitted 10 February, 2012; originally announced February 2012.

MSC Class: 60F10; 93E35

Holonomic gradient method for the distribution function of the largest root of a Wishart matrix

Authors: Hiroki Hashiguchi, Yasuhide Numata, Nobuki Takayama, Akimichi Takemura

Abstract: 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 off… ▽ More 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. △ Less

Submitted 20 January, 2013; v1 submitted 2 January, 2012; originally announced January 2012.

MSC Class: 62H10; 13N10

Journal ref: Journal of Multivariate Analysis 117 (2031) 296-312

Markov Bases for Typical Block Effect Models of Two-way Contingency Tables

Authors: Mitsunori Ogawa, Akimichi Takemura

Abstract: Markov basis for statistical model of contingency tables gives a useful tool for performing the conditional test of the model via Markov chain Monte Carlo method. In this paper we derive explicit forms of Markov bases for change point models and block diagonal effect models, which are typical block-wise effect models of two-way contingency tables, and perform conditional tests with some real data… ▽ More Markov basis for statistical model of contingency tables gives a useful tool for performing the conditional test of the model via Markov chain Monte Carlo method. In this paper we derive explicit forms of Markov bases for change point models and block diagonal effect models, which are typical block-wise effect models of two-way contingency tables, and perform conditional tests with some real data sets. △ Less

Submitted 19 December, 2011; originally announced December 2011.

Comments: 16 pages

MSC Class: 62H17

Journal ref: Journal of Multivariate Analysis 112 (2012) 219--229

Properties and applications of Fisher distribution on the rotation group

Authors: Tomonari Sei, Hiroki Shibata, Akimichi Takemura, Katsuyoshi Ohara, Nobuki Takayama

Abstract: 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 ide… ▽ More 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. △ Less

Submitted 3 February, 2013; v1 submitted 4 October, 2011; originally announced October 2011.

Comments: 25 pages, 1 figures

On computation of the characteristic polynomials of the discriminantal arrangements and the arrangements generated by generic points

Authors: Yasuhide Numata, Akimichi Takemura

Abstract: In this article we give a computational study of combinatorics of the discriminantal arrangements. The discriminantal arrangements are parametrized by two positive integers n and k such that n>k. The intersection lattice of the discriminantal arrangement with the parameter (n,k) is isomorphic to the intersection lattice of the hyperplane arrangement generated by n generic points in the d-dimension… ▽ More In this article we give a computational study of combinatorics of the discriminantal arrangements. The discriminantal arrangements are parametrized by two positive integers n and k such that n>k. The intersection lattice of the discriminantal arrangement with the parameter (n,k) is isomorphic to the intersection lattice of the hyperplane arrangement generated by n generic points in the d-dimensional vector space where d=n-k-1. The combinatorics of the discriminantal arrangements is very hard, except for the special cases of the Boolean arrangements (k=0) and the braid arrangements (k=1). We review some results on the intersection lattices of the arrangements generated by generic points and use them to obtain some computational results on the characteristic polynomials of the discriminantal arrangements. △ Less

Submitted 2 October, 2011; originally announced October 2011.

Journal ref: Harmony of Grobner Bases and the Modern Industrial Society, (Takayuki Hibi, editor), World Scientific, pp.228-252. 2012

Cones of elementary imsets and supermodular functions: a review and some new results

Authors: Takuya Kashimura, Tomonari Sei, Akimichi Takemura, Kentaro Tanaka

Abstract: In this paper we give a review of the method of imsets introduced by Studeny (2005) from a geometric point of view. Elementary imsets span a polyhedral cone and its dual cone is the cone of supermodular functions. We review basic facts on the structure of these cones. Then we derive some new results on the following topics: i) extreme rays of the cone of standardized supermodular functions, ii) fa… ▽ More In this paper we give a review of the method of imsets introduced by Studeny (2005) from a geometric point of view. Elementary imsets span a polyhedral cone and its dual cone is the cone of supermodular functions. We review basic facts on the structure of these cones. Then we derive some new results on the following topics: i) extreme rays of the cone of standardized supermodular functions, ii) faces of the cones, iii) small relations among elementary imsets, and iv) some computational results on Markov basis for the toric ideal defined by elementary imsets. △ Less

Submitted 12 September, 2011; originally announced September 2011.

MSC Class: 62H17

Journal ref: Harmony of Grobner Bases and the Modern Industrial Society, World Scientific, pp.117-152, 2012

Running Markov chain without Markov basis

Authors: Hisayuki Hara, Satoshi Aoki, Akimichi Takemura

Abstract: The methodology of Markov basis initiated by Diaconis and Sturmfels(1998) stimulated active research on Markov bases for more than ten years. It also motivated improvements of algorithms for Grobner basis computation for toric ideals, such as those implemented in 4ti2. However at present explicit forms of Markov bases are known only for some relatively simple models, such as the decomposable model… ▽ More The methodology of Markov basis initiated by Diaconis and Sturmfels(1998) stimulated active research on Markov bases for more than ten years. It also motivated improvements of algorithms for Grobner basis computation for toric ideals, such as those implemented in 4ti2. However at present explicit forms of Markov bases are known only for some relatively simple models, such as the decomposable models of contingency tables. Furthermore general algorithms for Markov bases computation often fail to produce Markov bases even for moderate-sized models in a practical amount of time. Hence so far we could not perform exact tests based on Markov basis methodology for many important practical problems. In this article we propose to use lattice bases for performing exact tests, in the case where Markov bases are not known. Computation of lattice bases is much easier than that of Markov bases. With many examples we show that the approach with lattice bases is practical. We also check that its performance is comparable to Markov bases for the problems where Markov bases are known. △ Less

Submitted 31 August, 2011; originally announced September 2011.

Comments: 18 pages

Properties of semi-elementary imsets as sums of elementary imsets

Authors: Takuya Kashimura, Tomonari Sei, Akimichi Takemura, Kentaro Tanaka

Abstract: We study properties of semi-elementary imsets and elementary imsets introduced by Studeny (2005). The rules of the semi-graphoid axiom (decomposition, weak union and contraction) for conditional independence statements can be translated into a simple identity among three semi-elementary imsets. By recursively applying the identity, any semi-elementary imset can be written as a sum of elementary im… ▽ More We study properties of semi-elementary imsets and elementary imsets introduced by Studeny (2005). The rules of the semi-graphoid axiom (decomposition, weak union and contraction) for conditional independence statements can be translated into a simple identity among three semi-elementary imsets. By recursively applying the identity, any semi-elementary imset can be written as a sum of elementary imsets, which we call a representation of the semi-elementary imset. A semi-elementary imset has many representations. We study properties of the set of possible representations of a semi-elementary imset and prove that all representations are connected by relations among four elementary imsets. △ Less

Submitted 30 May, 2011; originally announced May 2011.

MSC Class: 62H17

Journal ref: Journal of Algebraic Statistics 2 (2011) 14-35

Stochastic Bandit Based on Empirical Moments

Authors: Junya Honda, Akimichi Takemura

Abstract: In the multiarmed bandit problem a gambler chooses an arm of a slot machine to pull considering a tradeoff between exploration and exploitation. We study the stochastic bandit problem where each arm has a reward distribution supported in a known bounded interval, e.g. [0,1]. For this model, policies which take into account the empirical variances (i.e. second moments) of the arms are known to perf… ▽ More In the multiarmed bandit problem a gambler chooses an arm of a slot machine to pull considering a tradeoff between exploration and exploitation. We study the stochastic bandit problem where each arm has a reward distribution supported in a known bounded interval, e.g. [0,1]. For this model, policies which take into account the empirical variances (i.e. second moments) of the arms are known to perform effectively. In this paper, we generalize this idea and we propose a policy which exploits the first d empirical moments for arbitrary d fixed in advance. The asymptotic upper bound of the regret of the policy approaches the theoretical bound by Burnetas and Katehakis as d increases. By choosing appropriate d, the proposed policy realizes a tradeoff between the computational complexity and the expected regret. △ Less

Submitted 14 May, 2011; originally announced May 2011.

MSC Class: 62L05; 60G40

Journal ref: JMLR Workshop and Conference Proceedings, Volume 22: AISTATS 2012, 529-537. 2012

Arrangements stable under the Coxeter groups

Authors: Hidehiko Kamiya, Akimichi Takemura, Hiroaki Terao

Abstract: Let B be a real hyperplane arrangement which is stable under the action of a Coxeter group W. Then B acts naturally on the set of chambers of B. We assume that B is disjoint from the Coxeter arrangement A=A(W) of W. In this paper, we show that the W-orbits of the set of chambers of B are in one-to-one correspondence with the chambers of C=A\cup B which are contained in an arbitrarily fixed chamber… ▽ More Let B be a real hyperplane arrangement which is stable under the action of a Coxeter group W. Then B acts naturally on the set of chambers of B. We assume that B is disjoint from the Coxeter arrangement A=A(W) of W. In this paper, we show that the W-orbits of the set of chambers of B are in one-to-one correspondence with the chambers of C=A\cup B which are contained in an arbitrarily fixed chamber of A. From this fact, we find that the number of W-orbits of the set of chambers of B is given by the number of chambers of C divided by the order of W. We will also study the set of chambers of C which are contained in a chamber b of B. We prove that the cardinality of this set is equal to the order of the isotropy subgroup W_b of b. We illustrate these results with some examples, and solve an open problem in Kamiya, Takemura and Terao [Ranking patterns of unfolding models of codimension one, Adv. in Appl. Math. (2010)] by using our results. △ Less

Submitted 17 October, 2011; v1 submitted 26 March, 2011; originally announced March 2011.

MSC Class: 20F55; 32S22; 52C35

Journal ref: Configuration Spaces: Geometry, Combinatorics and Topology, Scuola Normale Superiore Pisa, pp.327-354, 2012

Convergence of random series and the rate of convergence of the strong law of large numbers in game-theoretic probability

Authors: Kenshi Miyabe, Akimichi Takemura

Abstract: We give a unified treatment of the convergence of random series and the rate of convergence of strong law of large numbers in the framework of game-theoretic probability of Shafer and Vovk (2001). We consider games with the quadratic hedge as well as more general weaker hedges. The latter corresponds to existence of an absolute moment of order smaller than two in the measure-theoretic framework. W… ▽ More We give a unified treatment of the convergence of random series and the rate of convergence of strong law of large numbers in the framework of game-theoretic probability of Shafer and Vovk (2001). We consider games with the quadratic hedge as well as more general weaker hedges. The latter corresponds to existence of an absolute moment of order smaller than two in the measure-theoretic framework. We prove some precise relations between the convergence of centered random series and the convergence of the series of prices of the hedges. When interpreted in measure-theoretic framework, these results characterize convergence of a martingale in terms of convergence of the series of conditional absolute moments. In order to prove these results we derive some fundamental results on deterministic strategies of Reality, who is a player in a protocol of game-theoretic probability. It is of particular interest, since Reality's strategies do not have any counterparts in measure-theoretic framework, ant yet they can be used to prove results, which can be interpreted in measure-theoretic framework. △ Less

Submitted 1 November, 2011; v1 submitted 7 March, 2011; originally announced March 2011.

Journal ref: Stochastic Processes and their Applications 122 (2012) 1-30

A lower bound for the Graver complexity of the incidence matrix of a complete bipartite graph

Authors: Taisei Kudo, Akimichi Takemura

Abstract: We give an exponential lower bound for the Graver complexity of the incidence matrix of a complete bipartite graph of arbitrary size. Our result is a generalization of the result by Berstein and Onn (2009) for 3xr complete bipartite graphs, r \ge 3. We give an exponential lower bound for the Graver complexity of the incidence matrix of a complete bipartite graph of arbitrary size. Our result is a generalization of the result by Berstein and Onn (2009) for 3xr complete bipartite graphs, r \ge 3. △ Less

Submitted 23 February, 2011; originally announced February 2011.

MSC Class: 62H17

Journal ref: Journal of Combinatorics 3 (2012) 695-708

Conformal geometry of statistical manifold with application to sequential estimation

Authors: Masayuki Kumon, Akimichi Takemura, Kei Takeuchi

Abstract: We present a geometrical method for analyzing sequential estimating procedures. It is based on the design principle of the second-order efficient sequential estimation provided in Okamoto, Amari and Takeuchi (1991). By introducing a dual conformal curvature quantity, we clarify the conditions for the covariance minimization of sequential estimators. These conditions are further elabolated for the… ▽ More We present a geometrical method for analyzing sequential estimating procedures. It is based on the design principle of the second-order efficient sequential estimation provided in Okamoto, Amari and Takeuchi (1991). By introducing a dual conformal curvature quantity, we clarify the conditions for the covariance minimization of sequential estimators. These conditions are further elabolated for the multidimensional curved exponential family. The theoretical results are then numerically examined by using typical statistical models, von Mises-Fisher and hyperboloid models. △ Less

Submitted 16 February, 2011; originally announced February 2011.

MSC Class: 62L12

Journal ref: Sequential Analysis 30 (2011) 308-337

Standard imsets for undirected and chain graphical models

Authors: Takuya Kashimura, Akimichi Takemura

Abstract: We derive standard imsets for undirected graphical models and chain graphical models. Standard imsets for undirected graphical models are described in terms of minimal triangulations for maximal prime subgraphs of the undirected graphs. For describing standard imsets for chain graphical models, we first define a triangulation of a chain graph. We then use the triangulation to generalize our result… ▽ More We derive standard imsets for undirected graphical models and chain graphical models. Standard imsets for undirected graphical models are described in terms of minimal triangulations for maximal prime subgraphs of the undirected graphs. For describing standard imsets for chain graphical models, we first define a triangulation of a chain graph. We then use the triangulation to generalize our results for the undirected graphs to chain graphs. △ Less

Submitted 27 July, 2015; v1 submitted 14 February, 2011; originally announced February 2011.

Comments: Published at http://dx.doi.org/10.3150/14-BEJ611 in the Bernoulli (http://isi.cbs.nl/bernoulli/) by the International Statistical Institute/Bernoulli Society (http://isi.cbs.nl/BS/bshome.htm)

Report number: IMS-BEJ-BEJ611

Journal ref: Bernoulli 2015, Vol. 21, No. 3, 1467-1493

Graver basis for an undirected graph and its application to testing the beta model of random graphs

Authors: Mitsunori Ogawa, Hisayuki Hara, Akimichi Takemura

Abstract: In this paper we give an explicit and algorithmic description of Graver basis for the toric ideal associated with a simple undirected graph and apply the basis for testing the beta model of random graphs by Markov chain Monte Carlo method. In this paper we give an explicit and algorithmic description of Graver basis for the toric ideal associated with a simple undirected graph and apply the basis for testing the beta model of random graphs by Markov chain Monte Carlo method. △ Less

Submitted 30 November, 2011; v1 submitted 13 February, 2011; originally announced February 2011.

MSC Class: 05C80; 62H17

Journal ref: Annals of the Institute of Statistical Mathematics 65 (2013) 191-212

On intersection lattices of hyperplane arrangements generated by generic points

Authors: Hiroshi Koizumi, Yasuhide Numata, Akimichi Takemura

Abstract: We consider hyperplane arrangements generated by generic points and study their intersection lattices. These arrangements are known to be equivalent to discriminantal arrangements. We show a fundamental structure of the intersection lattices by decomposing the poset ideals as direct products of smaller lattices corresponding to smaller dimensions. Based on this decomposition we compute the Möbius… ▽ More We consider hyperplane arrangements generated by generic points and study their intersection lattices. These arrangements are known to be equivalent to discriminantal arrangements. We show a fundamental structure of the intersection lattices by decomposing the poset ideals as direct products of smaller lattices corresponding to smaller dimensions. Based on this decomposition we compute the Möbius functions of the lattices and the characteristic polynomials of the arrangements up to dimension six. △ Less

Submitted 19 September, 2010; originally announced September 2010.

Comments: 17 pages, 4 figures

MSC Class: 52C35; 05A99

Journal ref: Annals of Combinatorics 16 (2012) 789-813

Approximations and asymptotics of upper hedging prices in multinomial models

Authors: Ryuichi Nakajima, Masayuki Kumon, Akimichi Takemura, Kei Takeuchi

Abstract: We give an exposition and numerical studies of upper hedging prices in multinomial models from the viewpoint of linear programming and the game-theoretic probability of Shafer and Vovk. We also show that, as the number of rounds goes to infinity, the upper hedging price of a European option converges to the solution of the Black-Scholes-Barenblatt equation. We give an exposition and numerical studies of upper hedging prices in multinomial models from the viewpoint of linear programming and the game-theoretic probability of Shafer and Vovk. We also show that, as the number of rounds goes to infinity, the upper hedging price of a European option converges to the solution of the Black-Scholes-Barenblatt equation. △ Less

Submitted 29 June, 2011; v1 submitted 25 July, 2010; originally announced July 2010.

MSC Class: 91G20

Journal ref: Japan Journal of Industrial and Applied Mathematics 25 (2012) 1-21