Markov chain models had proved to be useful tools in many fields, such as
physics, chemistry, information sciences, economics, finances, mathematical
biology, social sciences, and statistics for analyzing data. A discrete time
Markov chain is often used as a statistical model from a random physical
process to fit the observed data. A time-homogeneous Markov chain is a process
that each transition probability from a state to a state does not depend on
time. It is important to test if the assumption of the time-homogeneity of the
chain fits the observed data.
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.
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.
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.
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.
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.
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 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.
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.
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.
An admissible estimator of the eigenvalues of the variance-covariance matrix
is given for multivariate normal distributions with respect to the
scale-invariant squared error loss.
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 $S^n$,
which is used in the directional statistics. The method utilizes the theory and
algorithms of holonomic systems.
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 derive a Markov basis consisting of moves of degree at most three for
two-state toric homogeneous Markov chain model of arbitrary length without
parameters for initial states. Our basis consists of moves of degree three and
degree one, which alter the initial frequencies, in addition to moves of degree
two and degree one for toric homogeneous Markov chain model with parameters for
initial states.
Markov chain models are used in various fields, such behavioral sciences or
econometrics. Although the goodness of fit of the model is usually assessed by
large sample approximation, it is desirable to use conditional tests if the
sample size is not large. We study Markov bases for performing conditional
tests of the toric homogeneous Markov chain model, which is the envelope
exponential family for the usual homogeneous Markov chain model.
We give an expository review of applications of computational algebraic
statistics to design and analysis of fractional factorial experiments based on
our recent works. For the purpose of design, the techniques of Gr\"obner bases
and indicator functions allow us to treat fractional factorial designs without
distinction between regular designs and non-regular designs. For the purpose of
analysis of data from fractional factorial designs, the techniques of Markov
bases allow us to handle discrete observations.
Arrangement theory plays an essential role in the study of the unfolding
model used in many fields. This paper describes how arrangement theory can be
usefully employed in solving the problems of counting (i) the number of
admissible rankings in an unfolding model and (ii) the number of ranking
patterns generated by unfolding models. The paper is mostly expository but also
contains some new results such as simple upper and lower bounds for the number
of ranking patterns in the unidimensional case.
We propose procedures for testing whether stock price processes are
martingales based on limit order type betting strategies. We first show that
the null hypothesis of martingale property of a stock price process can be
tested based on the capital process of a betting strategy. In particular with
high frequency Markov type strategies we find that martingale null hypotheses
are rejected for many stock price processes.
We study Markov bases of decomposable graphical models consisting of
primitive moves (i.e., square-free moves of degree two) by determining the
structure of fibers of sample size two. We show that the number of elements of
fibers of sample size two are powers of two and we characterize primitive moves
in Markov bases in terms of connected components of induced subgraphs of the
independence graph of a hierarchical model. This allows us to derive a complete
description of minimal Markov bases and minimal invariant Markov bases for
decomposable models.
We propose minimum empirical divergence (MED) policy for the multiarmed
bandit problem. We prove asymptotic optimality of the proposed policy for the
case of finite support models. In our setting, Burnetas and Katehakis has
already proposed an asymptotically optimal policy. For choosing an arm our
policy uses a criterion which is dual to the quantity used in Burnetas and
Katehakis. Our criterion is easily computed by a convex optimization technique
and has an advantage in practical implementation.
We give some sufficient conditions of separation of two sets of integer
points by a hyperplane. Our conditions are related to the notion of convexity
of sets of integer points and are weaker than existing notions.
In this paper we propose an investing strategy based on neural network models
combined with ideas from game-theoretic probability of Shafer and Vovk. Our
proposed strategy uses parameter values of a neural network with the best
performance until the previous round (trading day) for deciding the investment
in the current round. We compare performance of our proposed strategy with
various strategies including a strategy based on supervised neural network
models and show that our procedure is competitive with other strategies.
We introduce a new formulation of asset trading games in continuous time in
the framework of the game-theoretic probability established by Shafer and Vovk
(Probability and Finance: It's Only a Game! (2001) Wiley). In our formulation,
the market moves continuously, but an investor trades in discrete times, which
can depend on the past path of the market. We prove that an investor can
essentially force that the asset price path behaves with the variation exponent
exactly equal to two.
We propose a sequential optimizing betting strategy in the multi-dimensional
bounded forecasting game in the framework of game-theoretic probability of
Shafer and Vovk (2001). By studying the asymptotic behavior of its capital
process, we prove a generalization of the strong law of large numbers, where
the convergence rate of the sample mean vector depends on the growth rate of
the quadratic variation process. The growth rate of the quadratic variation
process may be slower than the number of rounds or may even be zero.
For statistical analysis of multiway contingency tables we propose modeling
interaction terms in each maximal compact component of a hierarchical model. By
this approach we can search for parsimonious models with smaller degrees of
freedom than the usual hierarchical model, while preserving conditional
independence structures in the hierarchical model. We discuss estimation and
exacts tests of the proposed model and illustrate the advantage of the proposed
modeling with some data sets.
We discuss connecting tables with zero-one entries by a subset of a Markov
basis. In this paper, as a Markov basis we consider the Graver basis, which
corresponds to the unique minimal Markov basis for the Lawrence lifting of the
original configuration. Since the Graver basis tends to be large, it is of
interest to clarify conditions such that a subset of the Graver basis, in
particular a minimal Markov basis itself, connects tables with zero-one
entries. We give some theoretical results on the connectivity of tables with
zero-one entries.
We discuss connecting tables with zero-one entries by a subset of a Markov
basis. In this paper, as a Markov basis we consider the Graver basis, which
corresponds to the unique minimal Markov basis for the Lawrence lifting of the
original configuration. Since the Graver basis tends to be large, it is of
interest to clarify conditions such that a subset of the Graver basis, in
particular a minimal Markov basis itself, connects tables with zero-one
entries. We give some theoretical results on the connectivity of tables with
zero-one entries.