A new approach to distribution free tests in contingency tables

We suggest an extremely wide class of asymptotically distribution free goodness of fit tests for testing independence in two-way contingency tables, or equivalently, independence of two discrete random variables. The nature of these tests is that the test statistics can be viewed as definite functions of the transformation of \(\widehat{T}_n = (\widehat{T}_{ij})=\Big (\frac{\nu _{ij}- n\hat{a}_i\hat{b}_j}{\sqrt{n\hat{a}_i\hat{b}_j}}\Big )\) where \(\nu _{ij}\) are frequencies and \(\hat{a}_i, \hat{b}_j\) are estimated marginal distributions. Our method is also applicable for testing independence of two discrete random vectors. We make some comparisons on statistical powers of the new tests with the conventional chi-square test and suggest some cases in which this class is significantly more powerful.     

Notes on consistency of some minimum distance estimators with simulation results

We focus on the minimum distance density estimators \({\widehat{f}}_n\) of the true probability density \(f_0\) on the real line. The consistency of the order of \(n^{-1/2}\) in the (expected) L\(_1\)-norm of Kolmogorov estimator (MKE) is known if the degree of variations of the nonparametric family \(\mathcal {D}\) is finite. Using this result for MKE we prove that minimum Lévy and minimum discrepancy distance estimators are consistent of the order of \(n^{-1/2}\) in the (expected) L\(_1\)-norm under the same assumptions. Computer simulation for these minimum distance estimators, accompanied by Cramér estimator, is performed and the function \(s(n)=a_0+a_1\sqrt{n}\) is fitted to the L\(_1\)-errors of \({\widehat{f}}_n\) leading to the proportionality constant \(a_1\) determination. Further, (expected) L\(_1\)-consistency rate of Kolmogorov estimator under generalized assumptions based on asymptotic domination relation is studied. No usual continuity or differentiability conditions are needed.     

Ordering properties of the smallest order statistics from generalized Birnbaum–Saunders models with associated random shocks

Let \(X_{1},\ldots , X_{n}\) be lifetimes of components with independent non-negative generalized Birnbaum–Saunders random variables with shape parameters \(\alpha _{i}\) and scale parameters \(\beta _{i},~ i=1,\ldots ,n\), and \(I_{p_{1}},\ldots , I_{p_{n}}\) be independent Bernoulli random variables, independent of \(X_{i}\)’s, with \(E(I_{p_{i}})=p_{i},~i=1,\ldots ,n\). These are associated with random shocks on \(X_{i}\)’s. Then, \(Y_{i}=I_{p_{i}}X_{i}, ~i=1,\ldots ,n,\) correspond to the lifetimes when the random shock does not impact the components and zero when it does. In this paper, we discuss stochastic comparisons of the smallest order statistic arising from such random variables \(Y_{i},~i=1,\ldots ,n\). When the matrix of parameters \((h({\varvec{p}}), {\varvec{\beta }}^{\frac{1}{\nu }})\) or \((h({\varvec{p}}), {\varvec{\frac{1}{\alpha }}})\) changes to another matrix of parameters in a certain mathematical sense, we study the usual stochastic order of the smallest order statistic in such a setup. Finally, we apply the established results to two special cases: classical Birnbaum–Saunders and logistic Birnbaum–Saunders distributions.     

Unbiased estimates for a lognormal regression problem and a nonparametric alternative

Given a normal sample with means \({{\bf x}_{1}^{\prime} {\bf \varphi}, \ldots, {\bf x}_{n}^{\prime} {\bf \varphi}}\) and variance v, minimum variance unbiased estimates are given for the moments of L, where log L is normal with mean \({{\bf x}^{\prime} {\bf \varphi}}\) and variance v. These estimates converge to wrong values if the normality assumption is false. In the latter case estimates based on any M-estimate of \({{\bf \varphi}}\) are available of bias \({O\left(n^{-1}\right)}\) and \({O\left(n^{-2}\right)}\). More generally, these are given for any smooth function of \({\left({\bf \varphi}, F\right)}\), where F is the unknown distribution of the residuals. The regression functions need not be linear.     

Estimation of the regression slope by means of Gini’s cograduation index

The simple linear model \(Y_i = \alpha + \beta \, x_i + \epsilon _i\) \((i=1,2, \ldots ,N \ge 2)\) is considered, where the \(x_i\)’s are given constants and \(\epsilon _1, \epsilon _2 , \ldots , \epsilon _N\) are independent identically distributed (iid) with continuous distribution function F. An estimator \(\tilde{\beta }\) of the slope parameter is proposed, based on a stochastic process which makes use of Gini’s cograduation index. The properties of \(\tilde{\beta }\) and of the related confidence interval are studied. Some comparisons are given, in terms of asymptotic relative efficiency, with other estimators of \(\beta \) including that obtained with the method of least squares.     

Is the service sector different in size heterogeneity?

We examined whether significant differences in size heterogeneity exist between the service and the manufacturing industries by using PL exponents as the proxy for intra-industry size heterogeneity. For the purpose, we analyzed firm size distribution (FSD) and estimated the PL exponents, on the right tails of FSD, of the service and manufacturing industries in Korea for the period 2008–2012 using the Business Activity Survey dataset created by the Korean National Statistical Office As a result, we observed that the estimates of the PL exponents for the service industry are lower than those for the manufacturing industry (\(\upalpha _\mathrm{Service}<\upalpha _\mathrm{Manufacturing}\)) regardless of size variable, year, and dataset. This relationship may be related to the weaker negative relationship between the size and growth of the service industry, which made the slope of the PL distribution in the right tail of the FSD smoother. This finding implies that size heterogeneity may be more distinctive in the service industry than in the manufacturing industry. In addition, the PL exponents of sales were larger than those of assets and smaller than those of employees (\(\upalpha _\mathrm{Asset}<\upalpha _\mathrm{Sales}<\upalpha _\mathrm{Employee}\)) regardless of industry, year, and dataset. We also observed the PL exponents in the survived-firm dataset to decrease, compared to those in the all-firm dataset.     

Approximate maximum likelihood estimation for stochastic differential equations with random effects in the drift and the diffusion

Consider N independent stochastic processes \((X_i(t), t\in [0,T])\), \(i=1,\ldots , N\), defined by a stochastic differential equation with random effects where the drift term depends linearly on a random vector \(\Phi _i\) and the diffusion coefficient depends on another linear random effect \(\Psi _i\). For these effects, we consider a joint parametric distribution. We propose and study two approximate likelihoods for estimating the parameters of this joint distribution based on discrete observations of the processes on a fixed time interval. Consistent and \(\sqrt{N}\)-asymptotically Gaussian estimators are obtained when both the number of individuals and the number of observations per individual tend to infinity. The estimation methods are investigated on simulated data and show good performances.     

On consistency of the weighted least squares estimators in a semiparametric regression model

This paper is concerned with the semiparametric regression model \(y_i=x_i\beta +g(t_i)+\sigma _ie_i,~~i=1,2,\ldots ,n,\) where \(\sigma _i^2=f(u_i)\), \((x_i,t_i,u_i)\) are known fixed design points, \(\beta \) is an unknown parameter to be estimated, \(g(\cdot )\) and \(f(\cdot )\) are unknown functions, random errors \(e_i\) are widely orthant dependent random variables. The p-th (\(p>0\)) mean consistency and strong consistency for least squares estimators and weighted least squares estimators of \(\beta \) and g under some more mild conditions are investigated. A simulation study is also undertaken to assess the finite sample performance of the results that we established. The results obtained in the paper generalize and improve some corresponding ones of negatively associated random variables.     

On purely sequential estimation of an inverse Gaussian mean

The first part of this paper deals with developing a purely sequential methodology for the point estimation of the mean \(\mu \) of an inverse Gaussian distribution having an unknown scale parameter \(\lambda \). We assume a weighted squared error loss function and aim at controlling the associated risk function per unit cost by bounding it from above by a known constant \(\omega \). We also establish first-order and second-order asymptotic properties of our stopping rule. The second part of this paper deals with obtaining a purely sequential fixed accuracy confidence interval for the unknown mean \(\mu \), assuming that the scale parameter \(\lambda \) is known. First-order asymptotic efficiency and asymptotic consistency properties are also built of our proposed procedures. We then provide extensive sets of simulation studies and real data analysis using data from fatigue life analysis to show encouraging performances of our proposed stopping strategies.     

An iterative computational scheme for solving the coupled Hamilton–Jacobi–Isaacs equations in nonzero-sum differential games of affine nonlinear systems

In this paper, we present iterative or successive approximation methods for solving the coupled Hamilton–Jacobi–Isaacs equations (HJIEs) arising in nonzero-sum differential game for affine nonlinear systems. We particularly consider the ones arising in mixed \({\mathcal H}_{2}/{\mathcal H}_{\infty }\) control. However, the approach is perfectly general and can be applied to any others including those arising in the N-player case. The convergence of the method is established under fairly mild assumptions, and examples are solved to demonstrate the utility of the method. The results are also specialized to the coupled algebraic Riccati equations arising typically in mixed \({\mathcal H}_{2}/{\mathcal H}_{\infty }\) linear control. In this case, a bound within which the optimal solution lies is established. Finally, based on the iterative approach developed, a local existence result for the solution of the coupled-HJIEs is also established.     

Blocked factor aliased effect-number pattern and column rank of blocked regular designs

In factorial experiments, estimation precision of specific factor effects depends not only on design selection but also on factor assignments to columns of selected designs. Usually, different columns in a design play different roles when estimating factor effects. Zhou et al. (Can J Stat 41:540-555, 2013) introduced a factor aliased effect-number pattern (F-AENP) and proposed a column ranking scheme for all the GMC \(2^{n-m}\) designs with \(5N/16+1\le n\le N-1\), where \(N=2^{n-m}\). In this paper, we first introduce a blocked factor aliased effect-number pattern (B-F-AENP) for blocked regular designs as an extension of the F-AENP. Then, by using the B-F-AENP, we propose a column ranking scheme for all the B\(^1\)-GMC \(2^{n-m}:2^s\) designs with \(5N/16+1\le n\le N-1\), as well as an assignment strategy for important factors.     

Iterated Kalai–Smorodinsky–Nash compromise

In this paper, we introduce a new two-person bargaining solution, which we call iterated Kalai–Smorodinsky–Nash compromise (IKSNC). For its characterization, we present an axiom called \(\varGamma \)-Decomposability which is satisfied by any solution that is decomposable with respect to a given reference solution \(\varGamma \). We show that the IKSNC solution is uniquely characterized by \(\varGamma \)-Decomposability whenever \(\varGamma \) satisfies the standard axioms of Independence of Equivalent Utility Representations and Symmetry, along with three additional axioms, namely Restricted Monotonicity of Individually Best Extensions, Weak Independence of Irrelevant Alternatives, and Weak Pareto Optimality under Symmetry.     

The variance of the discrepancy distribution of rounding procedures,and sums of uniform random variables

When \(\ell \) probabilities are rounded to integer multiples of a given accuracy n, the sum of the numerators may deviate from n by a nonzero discrepancy. It is proved that, for large accuracies \(n \rightarrow \infty \), the limiting discrepancy distribution has variance \(\ell /12\). The relation to the uniform distribution over the interval \([-1/2, 1/2]\), whose variance is 1 / 12, is explored in detail.     

A new measure of association between random variables

We propose a new measure of association between two continuous random variables X and Y based on the covariance between X and the log-odds rate associated to Y. The proposed index of correlation lies in the range [\(-1\), 1]. We show that the extremes of the range, i.e., \(-1\) and 1, are attainable by the Fr\(\acute{\mathrm{e}}\)chet bivariate minimal and maximal distributions, respectively. It is also shown that if X and Y have bivariate normal distribution, the resulting measure of correlation equals the Pearson correlation coefficient \(\rho \). Some interpretations and relationships to other variability measures are presented. Among others, it is shown that for non-negative random variables the proposed association measure can be represented in terms of the mean residual and mean inactivity functions. Some illustrative examples are also provided.     

Linear estimation of location and scale parameters using partial maxima

Consider an i.i.d. sample \({X^*_{1},X^*_{2},\ldots,X^*_{n}}\) from a location-scale family, and assume that the only available observations consist of the partial maxima (or minima) sequence, \({X^*_{1:1},X^*_{2:2},\ldots,X^*_{n:n}}\), where \({X^*_{j:j}=\max\{ X^*_1, \ldots,X^*_j \}}\). This kind of truncation appears in several circumstances, including best performances in athletics events. In the case of partial maxima, the form of the BLUEs (best linear unbiased estimators) is quite similar to the form of the well-known Lloyd’s (in Biometrica 39:88–95, 1952) BLUEs, based on (the sufficient sample of) order statistics, but, in contrast to the classical case, their consistency is no longer obvious. The present paper is mainly concerned with the scale parameter, showing that the variance of the partial maxima BLUE is at most of order O(1/ log n), for a wide class of distributions.     

Proper strong-Fibonacci games

We define proper strong-Fibonacci (PSF) games as the subset of proper homogeneous weighted majority games which admit a Fibonacci representation. This is a homogeneous, type-preserving representation whose ordered sequence of type weights and winning quota is the initial string of Fibonacci numbers of the one-step delayed Fibonacci sequence. We show that for a PSF game, the Fibonacci representation coincides with the natural representation of the game. A characterization of PSF games is given in terms of their profile. This opens the way up to a straightforward formula which gives the number \(\varPsi (t)\) of such games as a function of t, number of non-dummy players’ types. It turns out that the growth rate of \(\varPsi (t)\) is exponential. The main result of our paper is that, for two consecutive t values of the same parity, the ratio \(\varPsi (t+2)/\varPsi (t)\) converges toward the golden ratio \({\varPhi }\).     

Reliability parameters estimation for parallel systems under imperfect repair

We consider in this paper a parallel system consisting of \(\eta \) identical components. Each component works independently of the others and has a Weibull distributed inter-failure time. When the system fails, we assume that the repair maintenance is imperfect according to the Arithmetic Reduction of Age models (\(ARA_{m}\)) proposed by Doyen and Gaudoin. The purpose of this paper is to generate a simulated failure data of the whole system in order to forecast the behavior of the failure process. Besides, we estimate the maintenance efficiency and the reliability parameters of an imperfect repair following \(ARA_{m}\) models using maximum likelihood estimation method. Our method is tested with several data sets available from related sources. The real data set corresponds to the time between failures of a compressor which is tested by Likelihood Ratio Test (LR). An analysis of the importance and the effect of the memory order of imperfect repair classes (\(ARA_{m}\)) will be discussed using LR test.     

Asymptotics of self-weighted M-estimators for autoregressive models

In this paper, we consider a stationary autoregressive AR(p) time series \(y_t=\phi _0+\phi _1y_{t-1}+\cdots +\phi _{p}y_{t-p}+u_t\). A self-weighted M-estimator for the AR(p) model is proposed. The asymptotic normality of this estimator is established, which includes the asymptotic properties under the innovations with finite or infinite variance. The result generalizes and improves the known one in the literature.     

Asymptotic inference for a one-dimensional simultaneous autoregressive model

A nonstationary simultaneous autoregressive model \({X^{(n)}_k=\alpha \Big(X^{(n)}_{k-1}+X^{(n)}_{k+1}\Big)+\varepsilon_k, k=1, 2, \ldots , n-1}\), is investigated, where \({X^{(n)}_0}\) and \({X^{(n)}_n}\) are given random variables. It is shown that in the unstable case α = 1/2 the least squares estimator of the autoregressive parameter converges to a functional of a standard Wiener process with a rate of convergence n ², while in the stable situation |α| < 1/2 the estimator is biased but asymptotically normal with a rate n ^1/2.     

Every member of the core is as respectful as any other

We strategically separate different core outcomes. The natural counterparts of a core allocation in a strategic environment are the α-core, the β-core and the strong equilibrium, modified by assuming that utility is transferable in a strategic context as well. Given a core allocation ω of a convex transferable utility (TU) game \(v\), we associate a strategic coalition formation game with \( \left( {v, \omega } \right) \) in which ω survives, while most other core allocations are eliminated. If the TU game is strictly convex, the core allocations respected by the TU-α-core, the TU-β-core and the TU-strong equilibrium shrink to ω only in the canonical family of coalition formation games associated with \( \left( {v, \omega } \right) \). A mechanism, which strategically separates core outcomes from noncore outcomes for each convex TU game according to the TU-strong equilibrium notion is reported.