On the invertibility conditions for moving average processes

Conditions for the general Moving Average process, of order q, to be invertible or borderline non-invertible are deduced. These are termed the acceptability conditions. It turns out that they depend on the magnitude of the final moving average parameter, θ _q . If ‖θ _q ‖ >1, the process is not acceptable. Should ‖θ _q ‖ = 1, the conditions, for any particular q, follow simply - if use is made of the remainder theorem. When ‖θ_q‖< 1, an appeal is made to ROUCH* E'S theorem, to establish the conditions. Analogous stationarity results immediately follow for autoregressive processes.     

On the error term in bahadur's representation of an order statistic

Bahadur (1966) presented a representation of an order statistic, giving its asymptotic distribution and the rate of convergence, under weak assumptions on the density function of the parent distribution. In this paper we consider the mean(squared) deviation of the error term in Bahadur’s approximation of the q th sample quantile (q_n ). We derive a uniform bound on the mean (squared) deviation of q_n , not depending on the value of q. An application of the given result provides the corresponding result for a kernel type estimator of the q th quantile.     

Euler(p,q) Processes and Their Application to Non Stationary Time Series with Time Varying Frequencies

We introduce Euler(p, q) processes as an extension of the Euler(p) processes for purposes of obtaining more parsimonious models for non stationary processes whose periodic behavior changes approximately linearly in time. The discrete Euler(p, q) models are a class of multiplicative stationary (M-stationary) processes and basic properties are derived. The relationship between continuous and discrete mixed Euler processes is shown. Fundamental to the theory and application of Euler(p, q) processes is a dual relationship between discrete Euler(p, q) processes and ARMA processes, which is established. The usefulness of Euler(p, q) processes is examined by comparing spectral estimation with that obtained by existing methods using both simulated and real data.     

Heine process as a q-analog of the Poisson process—waiting and interarrival times

In this study, we introduce the Heine process, {X_q(t), t > 0}, 0 < q < 1, where the random variable X_q(t), for every t > 0, represents the number of events (occurrences or arrivals) during a time interval (0, t]. The Heine process is introduced as a q-analog of the basic Poisson process. Also, in this study, we prove that the distribution of the waiting time W_{ν, q}, ν ? 1, up to the νth arrival, is a q-Erlang distribution and the interarrival times T_{k, q} = W_{k, q} ? W_{k ? 1, q},?k = 1, 2, …, ν with W_{0, q} = 0 are independent and equidistributed with a q-Exponential distribution.     

An exact decomposition for a class of exponential dispersion models: Application to sequential testing

We define the Wishart distribution on the cone of positive definite matrices and an exponential distribution on the Lorentz cone as exponential dispersion models. We show that these two distributions possess a property of exact decomposition, and we use this property to solve the following problem: given q samples (y_il,… y_iNj), i = l,…,q, from a N(μ_i,Σ_i,) distribution, test H:Σ₁ = Σ₂ = … = σ_q. Using the exact decomposition property, the classical test statistic for H, involving q parameters p_i = (N_i, - l)/2, i = 1,…,q, is replaced by a sequence of q - l test statistics for the sequence of tests H_i,:σ₁ =σ₂ = … =σ_i given that H_i-1 is true, i = 2,…,q. Each one of these test statistics involves two parameters only, p._i-1 = p₁ + … + p_i-1 and p_i. We also use the exact decomposition property to test equality of the “direction parameters” for q sample points from the exponential distribution on the Lorentz cone. We give a table of critical values for the distribution on the three-dimensional Lorentz cone. Tables of critical values in higher dimensions can easily be computed following the same method as in dimension three.     

On the asymptotics of randomness statistics

Kiefer (1959) studied the asymptotics of q-sample Cramér-von Mises nonparametric statistics when q is fixed and the sample sizes tend to infinity. Here we prove the asymptotic normality of such statistics when the sample sizes stay fixed or small while the number of samples, q, becomes large.     

Orthogonal polynomials generated by random vectors

Every random q-vector with finite moments generates a set of orthonormal polynomials. These are generated from the basis functions xn = xn11…xnqq using Gram–Schmidt orthogonalization. One can cycle through these basis functions using any number of ways. Here, we give results using minimum cycling. The polynomials look simpler when centered about the mean of X, and still simpler form when X is symmetric about zero. This leads to an extension of the multivariate Hermite polynomial for a general random vector symmetric about zero. As an example, the results are applied to the multivariate normal distribution.     

Maximum likelihood estimation for arma models in the presence of ARMA errors

An ARMA(p, q) process observed with an ARMA(c, d) error has an ARMA (p + c, k) representation with k = max(c + q, p + d) whose parameters satisfy some nonlinear constraints. Identification of the model is discussed. We develop Newton-Raphson estimators for the ARMA(p + c, k) process which take advantage of the information contained in the nonlinear restrictions. Explicit expressions for the derivatives of the restrictions are derived.     

Comparison of two life distributions on the basis of their percentile residual life functions

Let F and G be life distributions with respective failure rate functions r_F and r_G and respective 100α-percentile (0 < α < 1) residual life functions q_{α, F}, and q_{α, G}. Distribution-free two-sample tests are proposed for testing H₀: F = G against H_1,α,: q_{α, F} ≥ q_{α, G} and H_{2 α}: q_{β, F} ≥ q_β,G for all β ≥ α. This class of tests includes as a special case the test of Kochar (1981) for the alternative H*₂: r_F ≤ r_G. A theorem of Govindarajulu (1976) is extended in order to obtain asymptotic normality of the test statistics. The condition q_{α, F} ≥ q_{α, G} is implied by r_F ≤ r_G and is unrelated to the stochastic ordering F≤ G; if F and G are Weibull distributions with respective shape parameters c₁ and c₂ such that 1 ≤ C₁ < C₂, then q_α,F ≥ q_{α, G} for all α larger than a function of the parameters.     

The existence of almost difference families

Almost difference families (ADFs) were introduced by Ding and Yin as a useful generalization of almost difference sets (ADSs), and a number of infinite classes of almost difference families had been constructed. Suppose q is a prime power. To construct combinatorial designs in GF(q), one often needs to find an element xGF(q){0}, such that some polynomials in GF(q)[x] of degree one or two satisfying certain conditions. Weil's theorem on character sum estimates is very useful to do this. In this paper, a general bound for finding such x is given. By using this bound and computer searching, some known results on almost difference families by Ding and Yin are improved.     

A generalized step-up-down multiple test procedure

A generalization of step-up and step-down multiple test procedures is proposed. This step-up-down procedure is useful when the objective is to reject a specified minimum number, q, out of a family of k hypotheses. If this basic objective is met at the first step, then it proceeds in a step-down manner to see if more than q hypotheses can be rejected. Otherwise it proceeds in a step-up manner to see if some number less than q hypotheses can be rejected. The usual step-down procedure is the special case where q = 1, and the usual step-up procedure is the special case where q = k. Analytical and numerical comparisons between the powers of the step-up-down procedures with different choices of q are made to see how these powers depend on the actual number of false hypotheses. Examples of application include comparing the efficacy of a treatment to a control for multiple endpoints and testing the sensitivity of a clinical trial for comparing the efficacy of a new treatment with a set of standard treatments.     

Bayesian and Robust Bayesian analysis under a general class of balanced loss functions

For estimating an unknown parameter θ, we introduce and motivate the use of balanced loss functions of the form L_{r, w, d0}(q, d)=wr(d₀, d)+ (1-w) r(q, d){L_{\rho, \omega, \delta_0}(\theta, \delta)=\omega \rho(\delta_0, \delta)+ (1-\omega) \rho(\theta, \delta)}, as well as the weighted version q(q) L_{r, w, d0}(q, d){q(\theta) L_{\rho, \omega, \delta_0}(\theta, \delta)}, where ρ(θ, δ) is an arbitrary loss function, δ ₀ is a chosen a priori “target” estimator of q, w ? [0,1){\theta, \omega \in[0,1)}, and q(·) is a positive weight function. we develop Bayesian estimators under L_{r, w, d0}{L_{\rho, \omega, \delta_0}} with ω > 0 by relating such estimators to Bayesian solutions under L_{r, w, d0}{L_{\rho, \omega, \delta_0}} with ω = 0. Illustrations are given for various choices of ρ, such as absolute value, entropy, linex, and squared error type losses. Finally, under various robust Bayesian analysis criteria including posterior regret gamma-minimaxity, conditional gamma-minimaxity, and most stable, we establish explicit connections between optimal actions derived under balanced and unbalanced losses.     

D-Optimal Axial Designs for Quadratic and Cubic Additive Mixture Models

The paper discusses D-optimal axial designs for the additive quadratic and cubic mixture models σ_1≤i≤q(β_ix_i + β_iix²_i) and σ_1≤i≤q(β_ix_i + β_iix²_i + β_iiix³_i), where x_i≥ 0, x₁ + . . . + x_q = 1. For the quadratic model, a saturated symmetric axial design is used, in which support points are of the form (x₁, . . . , x_q) = [1 ? (q?1)δ_i, δ_i, . . . , δ_i], where i = 1, 2 and 0 ≤δ₂ <_{δ1 ≤} 1/(q ?1). It is proved that when 3 ≤q≤ 6, the above design is D-optimal if δ₂ = 0 and δ₁ = 1/(q?1), and when q≥ 7 it is D-optimal if δ₂ = 0 and δ₁ = [5q?1 ? (9q²?10q + 1)^1/2]/(4q²). Similar results exist for the cubic model, with support points of the form (x₁, . . . , x_q) = [1 ? (q?1)δ_i, δ_i, . . . , δ_i], where i = 1, 2, 3 and 0 = δ₃ <_δ2 < δ_{1 ≤}1/(q?1). The saturated D-optimal axial design and D-optimal design for the quadratic model are compared in terms of their efficiency and uniformity.     

Lp (1 < p ⩽ 2) solutions of one-dimensional BSDEs whose generator is weakly monotonic in y and non-Lipschitz in z

ABSTRACT

In this paper, we start with establishing the existence of a minimal (maximal) L^p (1 < p ? 2) solution to a one-dimensional backward stochastic differential equation (BSDE), where the generator g satisfies a p-order weak monotonicity condition together with a general growth condition in y and a linear growth condition in z. Then, we propose and prove a comparison theorem of L^p (1 < p ? 2) solutions to one-dimensional BSDEs with q-order (1 ? q < p) weak monotonicity and uniform continuity generators. As a consequence, an existence and uniqueness result of L^p (1 < p ? 2) solutions is also given for BSDEs whose generator g is q-order (1 ? q < p) weakly monotonic with a general growth in y and uniformly continuous in z.     

On the q-Weibull Distribution and Its Applications

The q-Weibull distribution is a stretched model for Weibull distribution, obtained by introducing a new pathway parameter q, which facilitates a slow transition to the Weibull as q → 1. In this article, we make a detailed study of the properties of the q-Weibull distribution and we apply it to a data on cancer remission times for which this distribution is a better fit than Weibull. Results relating to reliability properties, estimation of parameters, and applications in stress-strength analysis are also obtained.     

A Statistical Method for the Determination of the Appropriate Order in a General Class of Time Series Models

Abstract

We propose a method to determine the order q of a model in a general class of time series models. For the subset of linear moving average models (MA(q)), our method is compared with that of the sample autocorrelations. Since the sample autocorrelation is meant to detect a linear structure of dependence between random variables, it turns out to be more suitable for the linear case. However, our method presents a competitive option in that case, and for nonlinear models (NLMA(q)) it is shown to work better. The main advantages of our approach are that it does not make assumptions on the existence of moments and on the distribution of the noise involved in the moving average models. We also include an example with real data corresponding to the daily returns of the exchange rate process of mexican pesos and american dollars.     

Cramér‐von Mises regression

Consider a linear regression model with unknown regression parameters β₀ and independent errors of unknown distribution. Block the observations into q groups whose independent variables have a common value and measure the homogeneity of the blocks of residuals by a Cramér‐von Mises q‐sample statistic T_q(β). This statistic is designed so that its expected value as a function of the chosen regression parameter β has a minimum value of zero precisely at the true value β₀. The minimizer β of T_q(β) over all β is shown to be a consistent estimate of β₀. It is also shown that the bootstrap distribution of T_q(β₀) can be used to do a lack of fit test of the regression model and to construct a confidence region for β₀     

Estimating the number of clusters

Hartigan (1975) defines the number q of clusters in a d ‐variate statistical population as the number of connected components of the set {f > c}, where f denotes the underlying density function on R^d and c is a given constant. Some usual cluster algorithms treat q as an input which must be given in advance. The authors propose a method for estimating this parameter which is based on the computation of the number of connected components of an estimate of {f > c}. This set estimator is constructed as a union of balls with centres at an appropriate subsample which is selected via a nonparametric density estimator of f. The asymptotic behaviour of the proposed method is analyzed. A simulation study and an example with real data are also included.     

A Note on the Canonical Correlations From Contingency Tables

Explicit formulae are obtained for the asymptotic variances and covariances of canonical correlations which correspond to non-zero theoretical correlations in (p+ 1) x (q+1) contingency tables, with p≤q. The moments of the roots of a central Wishart matrix distributed as Wp(q; I ) are also given in general, with means, variances and covariances tabulated for p= 2, 3, 4: these may apply to canonical correlations corresponding to zeros.     

Model identification using 27 runs three level orthogonal arrays

In this paper we examine all the combinatorial non-isomorphic OA(27, q, 3, t), with 3≤q≤13 three level quantitative factors, with respect to model identification, estimation capacity and efficiency. We use the popular D-efficiency criterion to evaluate the ability of each design considered in estimating the parameters of a second-order model with adequate efficiency. The prior selection of the ‘middle’ level of factors plays an important role in the results.