
We propose an optimal experimental design for a curvilinear regression model
that minimizes the bandwidth of simultaneous confidence bands. Simultaneous
confidence bands for curvilinear regression are constructed by evaluating the
volume of a tube about a curve that is defined as a trajectory of a regression
basis vector (Naiman, 1986). The proposed criterion is constructed based on the
volume of a tube, and the corresponding optimal design that minimizes the
volume of tube is referred to as the tubevolume optimal (TVoptimal) design.
For Fourier and weighted polynomial regressions, the problem is formalized as
one of minimization over the cone of Hankel positive definite matrices, and the
criterion to minimize is expressed as an elliptic integral. We show that the
M\"obius group keeps our problem invariant, and hence, minimization can be
conducted over crosssections of orbits. We demonstrate that for the weighted
polynomial regression and the Fourier regression with three bases, the
tubevolume optimal design forms an orbit of the M\"obius group containing
Doptimal designs as representative elements.

We propose simultaneous confidence bands of the hyperbolictype for the
contrasts between several nonlinear (curvilinear) regression curves. The
critical value of a confidence band is determined from the distribution of the
maximum of a chisquare random process defined on the domain of explanatory
variables. We use the volumeoftube method to derive an upper tail probability
formula of the maximum of a chisquare random process, which is asymptotically
exact and sufficiently accurate in commonly used tail regions. Moreover, we
prove that the formula obtained is equivalent to the expectation of the
EulerPoincare characteristic of the excursion set of the chisquare random
process, and hence conservative. This result is therefore a generalization of
Naiman's inequality for Gaussian random processes. As an illustrative example,
growth curves of consomic mice are analyzed.

We consider one of the most basic multiple testing problems that compares
expectations of multivariate data among several groups. As a test statistic, a
conventional (approximate) $t$statistic is considered, and we determine its
rejection region using a common rejection limit. When there are unknown
correlations among test statistics, the multiplicity adjusted $p$values are
dependent on the unknown correlations. They are usually replaced with their
estimates that are always consistent under any hypothesis. In this paper, we
propose the use of estimates, which are not necessarily consistent and are
referred to as spurious correlations, in order to improve statistical power.
Through simulation studies, we verify that the proposed method asymptotically
controls the familywise error rate and clearly provides higher statistical
power than existing methods. In addition, the proposed and existing methods are
applied to a real multiple testing problem that compares quantitative traits
among groups of mice and the results are compared.

We first review the univariate and bivariate lackofmemory properties
(LMPs). The univariate LMP is a remarkable characterization of the exponential
distribution, while the bivariate LMP is shared by the famous Marshall and
Olkin's, Block and Basu's as well as Freund's bivariate exponential
distributions. We treat all the bivariate lackofmemory (BLM) distributions in
a unified approach and develop some new general properties of the BLM
distributions, including joint moment generating function, product moments and
dependence structure. Necessary and sufficient conditions for the survival
functions of BLM distributions to be totally positive of order two are given.
Some previous results for specific BLM distributions are improved. In
particular, we show that both the MarshallOlkin survival copula and survival
function are totally positive of all orders, regardless of parameters. Besides,
we point out that Slepian's inequality also holds true for the BLM
distributions.

We study zeroforcing detection (ZF) for multipleinput/multipleoutput
(MIMO) spatial multiplexing under transmitcorrelated Rician fading for an N_R
X N_T channel matrix with rank1 lineofsight (LoS) component. By using matrix
transformations and multivariate statistics, our exact analysis yields the
signaltonoise 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
Kfactor, 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 computeralgebra 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.

We show that the distribution of the scalar Schur complement in a noncentral
Wishart matrix is a mixture of central chisquare distributions with different
degrees of freedom. For the case of a rank1 noncentrality matrix, the weights
of the mixture representation arise from a noncentral beta mixture of Poisson
distributions.

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.

Let $V$ be a finite set of indices, and let $B_i$, $i=1,\ldots,m$, be subsets
of $V$ such that $V=\bigcup_{i=1}^{m}B_i$. Let $X_i$, $i\in V$, be independent
random variables, and let $X_{B_i}=(X_j)_{j\in B_i}$. In this paper, we propose
a recursive computation method to calculate the conditional expectation
$E\bigl[\prod_{i=1}^m\chi_i(X_{B_i}) \,\, N\bigr]$ with $N=\sum_{i\in V}X_i$
given, where $\chi_i$ is an arbitrary function. Our method is based on the
recursive summation/integration technique using the Markov property in
statistics. To extract the Markov property, we define an undirected graph whose
cliques are $B_j$, and obtain its chordal extension, from which we present the
expressions of the recursive formula. This methodology works for a class of
distributions including the Poisson distribution (that is, the conditional
distribution is the multinomial). This problem is motivated from the evaluation
of the multiplicityadjusted $p$value of scan statistics in spatial
epidemiology. As an illustration of the approach, we present the real data
analyses to detect temporal and spatial clustering.

For multipleinput multipleoutput (MIMO) spatialmultiplexing transmission,
zeroforcing detection (ZF) is appealing because of its low complexity. Our
recent MIMO ZF performance analysis for RicianRayleigh fading, which is
relevant in heterogeneous networks, has yielded for the ZF outage probability
and ergodic capacity infiniteseries expressions. Because they arose from
expanding the confluent hypergeometric function $ {_1\! F_1} (\cdot, \cdot,
\sigma) $ around 0, they do not converge numerically at realisticallyhigh
Rician $ K $factor values. Therefore, herein, we seek to take advantage of the
fact that $ {_1\! F_1} (\cdot, \cdot, \sigma) $ 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 signaltonoise ratio (SNR) are holonomic. Then, from the
differential equation for $ {_1\! F_1} (\cdot, \cdot, \sigma) $, 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 practicallyrelevant values of $ K $.
Finally, numerical integration of the SNR p.d.f. produced by HGM yields
accurate ZF outage probability and ergodic capacity results.

For multipleinput/multipleoutput (MIMO) spatial multiplexing with
zeroforcing detection (ZF), signaltonoise ratio (SNR) analysis for Rician
fading involves the cumbersome noncentralWishart distribution (NCWD) of the
transmit samplecorrelation (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 SNRdistribution approximation were unknown. Therefore, we have been
attempting to exactly characterize ZF SNR for Rician fading. Our previous
attempts succeeded only for the sole Ricianfading stream under
RicianRayleigh fading, by writing it as scalar Schur complement (SC) in the
Gramian. Herein, we pursue a more general, matrixSCbased analysis to
characterize SNRs when several streams may undergo Rician fading. On one hand,
for fullRician fading, the SC distribution is found to be exactly a CWD if and
only if a channelmeancorrelation \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 SNRdistributions coincide. On the other hand, for RicianRayleigh
fading, the matrixSC distribution is characterized in terms of determinant of
matrix with elementaryfunction entries, which also yields a new
characterization of the ZF SNR. Average error probability results validate our
analysis vs.~simulation.

A method that uses order statistics to construct multivariate distributions
with fixed marginals and which utilizes a representation of the Bernstein
copula in terms of a finite mixture distribution is proposed.
Expectationmaximization (EM) algorithms to estimate the Bernstein copula are
proposed, and a local convergence property is proved. Moreover, asymptotic
properties of the proposed semiparametric estimators are provided. Illustrative
examples are presented using three real data sets and a 3dimensional simulated
data set. These studies show that the Bernstein copula is able to represent
various distributions flexibly and that the proposed EM algorithms work well
for such data.

We analyze the performance of multiple input/multiple output (MIMO)
communications systems employing spatial multiplexing and zeroforcing
detection (ZF). The distribution of the ZF signaltonoise 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 wellknown SNR expression has been
hindered by the noncentrality of the Wishart distribution involved. In
addition, approximation with a centralWishart distribution has not proved
consistently accurate. In contrast, the following exact ZF study proceeds from
a lesserknown SNR expression that separates the intended and interfering
channelgain vectors. By first conditioning on, and then averaging over the
interference, the ZF SNR distribution for RicianRayleigh fading is shown to be
an infinite linear combination of gamma distributions. On the other hand, for
RayleighRician fading, the ZF SNR is shown to be gammadistributed. 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.

Define a chisquare random field on a multidimensional lattice points index
set with a directproduct covariance structure, and consider the distribution
of the maximum of this random field. We provide two approximate formulas for
the upper tail probability of the distribution based on nonlinear renewal
theory and an integralgeometric approach called the volumeoftube method.
This study is motivated by the detection problem of the interactive loci pairs
which play an important role in forming biological species. The joint
distribution of scan statistics for detecting the pairs is regarded as the
chisquare random field above, and hence the multiplicityadjusted $p$value
can be calculated by using the proposed approximate formulas. By using these
formulas, we examine the data of Mizuta, et al. (2010) who reported a new
interactive loci pair of rice intersubspecies.

A polynomial that is nonnegative over a given interval is called a positive
polynomial. The set of such positive polynomials forms a closed convex cone
$K$. In this paper, we consider the likelihood ratio test for the hypothesis of
positivity that the estimand polynomial regression curve is a positive
polynomial. By considering hierarchical hypotheses including the hypothesis of
positivity, we define nested likelihood ratio tests, and derive their null
distributions as mixtures of chisquare distributions by using the
volumeoftubes method. The mixing probabilities are obtained by utilizing the
parameterizations for the cone $K$ and its dual provided in the framework of
Tchebycheff systems for polynomials of degree at most 4. For polynomials of
degree greater than 4, the upper and lower bounds for the null distributions
are provided. Moreover, we propose associated simultaneous confidence bounds
for polynomial regression curves. Regarding computation, we demonstrate that
symmetric cone programming is useful to obtain the test statistics. As an
illustrative example, we conduct data analysis on growth curves of two groups.
We examine the hypothesis that the growth rate (the derivative of growth curve)
of one group is always higher than the other.

Let $K$ be a closed convex polyhedron defined by a finite number of linear
inequalities. In this paper we refine the theory of abstract tubes (Naiman and
Wynn, 1997) associated with $K$ when $K$ is perturbed. In particular, we focus
on the perturbation that is lexicographic and in an outer direction. An
algorithm for constructing the abstract tube by means of linear programming and
its implementation are discussed. Using the abstract tube for perturbed $K$
combined with the recursive integration technique proposed by Miwa, Hayter and
Kuriki (2003), we show that the multidimensional normal probability for a
polyhedral region $K$ can be computed efficiently. In addition, abstract tubes
and the distribution functions of studentized range statistics are exhibited as
numerical examples.

Let $A$ be a real skewsymmetric Gaussian random matrix whose upper
triangular elements are independently distributed according to the standard
normal distribution. We provide the distribution of the largest singular value
$\sigma_1$ of $A$. Moreover, by acknowledging the fact that the largest
singular value can be regarded as the maximum of a Gaussian field, we deduce
the distribution of the standardized largest singular value
$\sigma_1/\sqrt{\mathrm{tr}(A'A)/2}$. These distributional results are utilized
in Scheff\'{e}'s paired comparisons model. We propose tests for the hypothesis
of subtractivity based on the largest singular value of the skewsymmetric
residual matrix. Professional baseball league data are analyzed as an
illustrative example.

We provide formulas for the moments of the real and complex noncentral
Wishart distributions of general degrees. The obtained formulas for the real
and complex cases are described in terms of the undirected and directed graphs,
respectively. By considering degenerate cases, we give explicit formulas for
the moments of bivariate chisquare distributions and $2\times 2$ Wishart
distributions by enumerating the graphs. Noting that the Laguerre polynomials
can be considered to be moments of a noncentral chisquare distributions
formally, we demonstrate a combinatorial interpretation of the coefficients of
the Laguerre polynomials.

The projection pursuit index defined by a sum of squares of the third and the
fourth sample cumulants is known as the moment index proposed by Jones and
Sibson. Limiting distribution of the maximum of the moment index under the null
hypothesis that the population is multivariate normal is shown to be the
maximum of a Gaussian random field with a finite KarhunenLoeve expansion. An
approximate formula for tail probability of the maximum, which corresponds to
the pvalue, is given by virtue of the tube method through determining Weyl's
invariants of all degrees and the critical radius of the index manifold of the
Gaussian random field.

Consider testing normality against a oneparameter family of univariate
distributions containing the normal distribution as the boundary, e.g., the
family of $t$distributions or an infinitely divisible family with finite
variance. We prove that under mild regularity conditions, the sample skewness
is the locally best invariant (LBI) test of normality against a wide class of
asymmetric families and the kurtosis is the LBI test against symmetric
families. We also discuss nonregular cases such as testing normality against
the stable family and some related results in the multivariate cases.

Elliptically contoured distributions can be considered to be the
distributions for which the contours of the density functions are proportional
ellipsoids. We generalize elliptically contoured densities to ``starshaped
distributions'' with concentric starshaped contours and show that many results
in the former case continue to hold in the more general case. We develop a
general theory in the framework of abstract group invariance so that the
results can be applied to other cases as well, especially those involving
random matrices.