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.     

Dispersive ordering of fail-safe systems with heterogeneous exponential components

Let X ₁, . . . , X _n be independent exponential random variables with respective hazard rates λ₁, . . . , λ _n , and Y ₁, . . . , Y _n be independent and identically distributed random variables from an exponential distribution with hazard rate λ. Then, we prove that X _2:n , the second order statistic from X ₁, . . . , X _n , is larger than Y _2:n , the second order statistic from Y ₁, . . . , Y _n , in terms of the dispersive order if and only if

$\lambda\geq \sqrt{\frac{1}{{n\choose 2}}\sum_{1\leq i < j\leq n}\lambda_i\lambda_j}.$

We also show that X _2:n is smaller than Y _2:n in terms of the dispersive order if and only if

$ \lambda\le\frac{\sum^{n}_{i=1} \lambda_i-{\rm max}_{1\leq i\leq n} \lambda_i}{n-1}. $

Moreover, we extend the above two results to the proportional hazard rates model. These two results established here form nice extensions of the corresponding results on hazard rate, likelihood ratio, and MRL orderings established recently by Pǎltǎnea (J Stat Plan Inference 138:1993–1997, 2008), Zhao et al. (J Multivar Anal 100:952–962, 2009), and Zhao and Balakrishnan (J Stat Plan Inference 139:3027–3037, 2009), respectively.

     

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.     

A Cramér-type large deviation theorem for sums of functions of higher order non-overlapping spacings

Let U ₁, U ₂, . . . , U _n–1 be an ordered sample from a Uniform [0,1] distribution. The non-overlapping uniform spacings of order s are defined as \({G_{i}^{(s)} =U_{is} -U_{(i-1)s}, i=1,2,\ldots,N^\prime, G_{N^\prime+1}^{(s)} =1-U_{N^\prime s}}\) with notation U ₀ = 0, U _n = 1, where \({N^\prime=\left\lfloor n/s\right\rfloor}\) is the integer part of n/s. Let \({ N=\left\lceil n/s\right\rceil}\) be the smallest integer greater than or equal to n/s, f _m (u), m = 1, 2, . . . , N, be a sequence of real-valued Borel-measurable functions. In this article a Cramér type large deviation theorem for the statistic \({f_{1,n} (nG_{1}^{(s)})+\cdots+f_{N,n} (nG_{N}^{(s)} )}\) is proved.     

Almost sure central limit theorem for the products of <Emphasis Type="Italic">U</Emphasis>-statistics

Let (X _n ) be a sequence of i.i.d random variables and U _n a U-statistic corresponding to a symmetric kernel function h, where h ₁(x ₁) = Eh(x ₁, X ₂, X ₃, . . . , X _m ), μ = E(h(X ₁, X ₂, . . . , X _m )) and ? ₁ = Var(h ₁(X ₁)). Denote \({\gamma=\sqrt{\varsigma_{1}}/\mu}\), the coefficient of variation. Assume that P(h(X ₁, X ₂, . . . , X _m ) > 0) = 1, ? ₁ > 0 and E|h(X ₁, X ₂, . . . , X _m )|³ < ∞. We give herein the conditions under which

$\lim_{N\rightarrow\infty}\frac{1}{\log N}\sum_{n=1}^{N}\frac{1}{n}g\left(\left(\prod_{k=m}^{n}\frac{U_{k}}{\mu}\right)^{\frac{1}{m\gamma\sqrt{n}}}\right) =\int\limits_{-\infty}^{\infty}g(x)dF(x)\quad {\rm a.s.}$

for a certain family of unbounded measurable functions g, where F(·) is the distribution function of the random variable \({\exp(\sqrt{2} \xi)}\) and ξ is a standard normal random variable.

     

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.     

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.     

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.     

Asymptotic non-deficiency of the Bayes sequential estimation in a family of transformed Chi-square distributions

Bayes sequential estimation in a family of transformed Chi-square distributions using a linex loss function and a cost c > 0 for each observation is considered in this paper. It is shown that an asymptotic pointwise optimal rule (A.P.O.) is asymptotically non-deficient, i.e., the difference between the Bayes risk of the A.P.O. rule and the Bayes risk of the optimal procedure is of smaller order of magnitude than c, the cost of single observation, as c → 0.     

Consistency of a nonparametric estimation of a density functional

LetX be a random variable with distribution functionF and density functionf. Let ? and ψ be known measurable functions defined on the real lineR and the closed interval [0, 1], respectively. This paper proposes a smooth nonparametric estimate of the density functional \(\theta = \int\limits_R \phi (x) \psi \left[ {F (x)} \right]f^2 (x) dx\) based on a random sampleX ₁, ...,X _n fromF using a kernel functionk. The proposed estimate is given by \(\hat \theta = (n^2 a_n )^{ - 1} \mathop \sum \limits_{i = 1}^n \mathop \sum \limits_{j = 1}^n \phi (X_i ) \psi \left[ {\hat F (X_i )} \right]k\left[ {(X_i - X_j )/a_n } \right]\) , where \(\hat F(x) = n^{ - 1} \mathop \sum \limits_{i = 1}^n K\left[ {(x - X_i )/a_n } \right]\) with \(K (w) = \int\limits_{ - \infty }^w {k (u) } du\) . The estimate \(\hat \theta \) is shown to be consistent both in the weak and strong sense and is used to estimate the asymptotic relative efficiency of various nonparametric tests, with particular reference to those using the Chernoff-Savage statistic.     

Hermite ranks and U-statistics

We focus on the asymptotic behavior of $U$ -statistics of the type $$\begin{aligned} \sum _{1\le i\ne j\le n} h(X_i,X_j)\\ \end{aligned}$$ in the long-range dependence setting, where $(X_i)_{i\ge 1}$ is a stationary mean-zero Gaussian process. Since $(X_i)_{i\ge 1}$ is Gaussian, $h$ can be decomposed in Hermite polynomials. The goal of this paper is to compare the different notions of Hermite rank and to provide conditions for the remainder term in the decomposition to be asymptotically negligeable.     

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.     

The concept of weak exchangeability and its applications

A finite sequence of binary random variables is called a weak exchangeable sequence of order m if the sequence consists of m random vectors such that the elements within each random vector are exchangeable in the usual sense and the different random vectors are dependent. The exact and asymptotic joint distributions of the m-dimensional random vector whose elements include the number of successes in each exchangeable sequence are derived. Potential applications of the concept of weak exchangeability are discussed with illustrative examples.     

On generalized transforms of distribution functions

In this note we discuss the following problem. LetX andY to be two real valued independent r.v.'s with d.f.'sF and ?. Consider the d.f.F*? of the r.v.X oY, being o a binary operation among real numbers. We deal with the following equation: $$\mathcal{G}^1 (F * \phi ,s) = \mathcal{G}^2 (F,s)\square \mathcal{G}^3 (\phi ,s)\forall s \in S$$ where \(\mathcal{G}^1 ,\mathcal{G}^2 ,\mathcal{G}^3 \) are real or complex functionals, т another binary operation ands a parameter. We give a solution, that under stronger assumptions (Aczél 1966), is the only one, of the problem. Such a solution is obtained in two steps. First of all we give a solution in the very special case in whichX andY are degenerate r.v.'s. Secondly we extend the result to the general case under the following additional assumption: $$\begin{gathered} \mathcal{G}^1 (\alpha F + (1 - \alpha )\phi ,s) = H[\mathcal{G}^i (F,s),\mathcal{G}^i (\phi ,s);\alpha ] \hfill \\ \forall \alpha \in [0,1]i = 1,2,3 \hfill \\ \end{gathered} $$ .     

R-optimal designs for multi-factor models with heteroscedastic errors

In this paper, we consider the R-optimal design problem for multi-factor regression models with heteroscedastic errors. It is shown that a R-optimal design for the heteroscedastic Kronecker product model is given by the product of the R-optimal designs for the marginal one-factor models. However, R-optimal designs for the additive models can be constructed from R-optimal designs for the one-factor models only if sufficient conditions are satisfied. Several examples are presented to illustrate and check optimal designs based on R-optimality criterion.     

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.     

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.     

{\varvec{P}}-value model selection criteria for exponential families of increasing dimension

Let $\mathcal{M }_{\underline{i}}$ be an exponential family of densities on $[0,1]$ pertaining to a vector of orthonormal functions $b_{\underline{i}}=(b_{i_1}(x),\ldots ,b_{i_p}(x))^\mathbf{T}$ and consider a problem of estimating a density $f$ belonging to such family for unknown set ${\underline{i}}\subset \{1,2,\ldots ,m\}$ , based on a random sample $X_1,\ldots ,X_n$ . Pokarowski and Mielniczuk (2011) introduced model selection criteria in a general setting based on p-values of likelihood ratio statistic for $H_0: f\in \mathcal{M }_0$ versus $H_1: f\in \mathcal{M }_{\underline{i}}\setminus \mathcal{M }_0$ , where $\mathcal{M }_0$ is the minimal model. In the paper we study consistency of these model selection criteria when the number of the models is allowed to increase with a sample size and $f$ ultimately belongs to one of them. The results are then generalized to the case when the logarithm of $f$ has infinite expansion with respect to $(b_i(\cdot ))_1^\infty $ . Moreover, it is shown how the results can be applied to study convergence rates of ensuing post-model-selection estimators of the density with respect to Kullback–Leibler distance. We also present results of simulation study comparing small sample performance of the discussed selection criteria and the post-model-selection estimators with analogous entities based on Schwarz’s rule as well as their greedy counterparts.