On the construction of hypercubic design

Hypercubic design was introduced by Shah (1958) and Kusumoto(1965) in which the t=v^m treatments are represented by n-plets (x₁, x₂, ..., x_m),where 1≤ x₁, x₂, ..., x_m ≤v, and two treatments are said to be i-th associates if they differ in exactly i components. This paper deals with the construction of some hypercubic designs and gives their application to confounding plans for symmetrical factorial experiments. The proposed methods prove to be quite flexible in terms of choice of possible block sizes and are easy to use.     

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.     

Limiting distributions of Kolmogorov-Lévy-type statistics under the alternative

Let X_i, 1 ≤ i ≤ n, be independent identically distributed random variables with a common distribution function F, and let G be a smooth distribution function. We derive the limit distribution of □{ρ_α(F_n, G) - _α(F, G)}, where F_n is the empirical distribution function based on X₁,…,X_n and _α is a Kolmogorov-Lévy-type metric between distribution functions. For α ≤ 0 and two distribution functions F and G the metric p_α is given by p_α(F, G) = inf {? ≤ 0: G(x - α?) - ? F(x) ≤ G(x + α?) + ? for all x ?}.     

Nonparametric Estimation of a Regression Function: Limiting Distribution2

Consider the regression model Y_i= g(x_i) + e_i, i = 1,…, n, where g is an unknown function defined on [0, 1], 0 = x₀ < x₁ < … < x_n≤ 1 are chosen so that max_1≤i≤n(x_i-x_{i- 1}) = 0(n^-1), and where {e_i} are i.i.d. with Ee₁= 0 and Var e₁ - s?². In a previous paper, Cheng & Lin (1979) study three estimators of g, namely, g_1n of Cheng & Lin (1979), g_2n of Clark (1977), and g_3n of Priestley & Chao (1972). Consistency results are established and rates of strong uniform convergence are obtained. In the current investigation the limiting distribution of &_in, i = 1, 2, 3, and that of the isotonic estimator g**_n are considered.     

Bayesian one-sample prediction of future observations under Pareto distribution

Suppose that the length of time in years for which a business operates until failure has a Pareto distribution. Let t ₁?t ₂?t _r denote the survival lifetimes of the first r of a random sample of n businesses. Bayesian predictions are to be made on the ordered failure times of the remaining (n???r) businesses, using the conditional probability function. Numerical examples are given to illustrate our results.     

Lois limites pour les statistiques d'ordre dans le cas non identiquement distribue

Let X₁, X₂,… be a sequence of independent random variables with distribution functions F₁, where 1 ≤ i ≤ n, and for each n ≥ 1 let X_1,n ≤… ≤ X_n,n denote the order statistics of the first n random variables. Under suitable hypotheses about the F₁, we characterize the limit distribution functions H(x) for which P(X_k,n ? a_nx + b_n) → H(x), where a_n > 0 and b_n are real constants. We consider the cases where κ = κ(n) satisfies √n {κ(n)/n — λ} → 0 and √n {κ(n)/n — λ} → ∞ separately.     

ASYMPTOTIC AND SMALL SAMPLE STATISTICAL PROPERTIES OF RANDOM FRAILTY VARIANCE ESTIMATES FOR SHARED GAMMA FRAILTY MODELS

This paper concerns maximum likelihood estimation for the semiparametric shared gamma frailty model; that is the Cox proportional hazards model with the hazard function multiplied by a gamma random variable with mean 1 and variance θ. A hybrid ML-EM algorithm is applied to 26 400 simulated samples of 400 to 8000 observations with Weibull hazards. The hybrid algorithm is much faster than the standard EM algorithm, faster than standard direct maximum likelihood (ML, Newton Raphson) for large samples, and gives almost identical results to the penalised likelihood method in S-PLUS 2000. When the true value θ₀ of θ is zero, the estimates of θ are asymptotically distributed as a 50–50 mixture between a point mass at zero and a normal random variable on the positive axis. When θ₀ > 0, the asymptotic distribution is normal. However, for small samples, simulations suggest that the estimates of θ are approximately distributed as an x ? (100 ? x)% mixture, 0 ≤ x ≤ 50, between a point mass at zero and a normal random variable on the positive axis even for θ₀ > 0. In light of this, p-values and confidence intervals need to be adjusted accordingly. We indicate an approximate method for carrying out the adjustment.     

A Note on the Calculation of Pr{X 1 < X 2 < ··· < Xk }

Suppose that X_i are independent random variables, and that X_i has cdf F_i (x), 1 ≤ i ≤ k. Many statistical problems involve the probability Pr{X ₁ < X ₂ < ··· < X_k }. In this note a numerical method is proposed for computing this probability.     

Simultaneous confidence bands for expectile functions

Expectile regression, as a general M smoother, is used to capture the tail behaviour of a distribution. Let (X ₁,Y ₁),…,(X _n ,Y _n ) be i.i.d. rvs. Denote by v(x) the unknown τ-expectile regression curve of Y conditional on X, and by v _n (x) its kernel smoothing estimator. In this paper, we prove the strong uniform consistency rate of v _n (x) under general conditions. Moreover, using strong approximations of the empirical process and extreme value theory, we consider the asymptotic maximal deviation sup_0≤x≤1|v _n (x)?v(x)|. According to the asymptotic theory, we construct simultaneous confidence bands around the estimated expectile function. Furthermore, we apply this confidence band to temperature analysis. Taking Berlin and Taipei as an example, we investigate the temperature risk drivers to these two cities.     

FIRST-EXIT TIMES FOR POISSON SHOT NOISE

For a shot-noise process X(t) with Poisson arrival times and exponentially diminishing shocks of i.i.d. sizes, we consider the first time T _b at which a given level b > 0 is exceeded. An integral equation for the joint density of T _b and X(T _b ) is derived and, for the case of exponential jumps, solved explicitly in terms of Laplace transforms (LTs). In the general case we determine the ordinary LT of the function ? P(T _b > t) in terms of certain LTs derived from the distribution function H(x; t) = P(X(t) ≤ x), considered as a function of both variables x and t. Moreover, for G(t, u) = P(T _b > t, X(t) < u), that is the joint distribution function of sup_{0 ≤ s ≤ t} X(s) and X(t), an integro-differential equation is presented, whose unique solution is G(t, u).     

A class of tests for testing ‘new is better than used’

A class of tests is proposed for testing H₀ F?(x) = e^?λx, λ > 0, x≥0 vs. H₁ F?(x + y) ≤ F?(x)F?(y), x, y≥0, with strict inequality for some x, y ≥ 0 (F = new is better than used). Efficiency comparisons of some tests within the class are made and a new test is proposed on the basis of these comparisons. Consistency and the asymptotic normality of the class of tests is proved under fairly broad conditions on the underlying entities.     

A Class of Multivariate Distributions

In this paper we introduce a class of multivariate distributions, known as the generalized Liouville distribution and defined by the functional form (0 ≤x_i < ∞, αi > 0, β_i > 0, q_i > 0). It is shown that such distributions can be used to derive both the Dirichlet distribution and the beta distribution of the second kind.     

Partial and inverse autocorrelations in portmanteau-type tests for time series

We present a decomposition of the correlation coefficient between x_t and x_t?k into three terms that include the partial and inverse autocorrelations. The first term accounts for the portion of the autocorrelation that is explained by the inner variables {x_t?1 , x_t?2 , …, x _{t? k+1}}, the second one measures the portion explained by the outer variables {x _t+1, x _t+2, } ∪ {x _t?k?1, x _t?k?2,…} and the third term measures the correlation between x _t and x_t?k given all other variables. These terms, squared and summed, can form the basis of three portmanteau-type tests that are able to detect both deviation from white noise and lack of fit of an entertained model. Quantiles of their asymptotic sample distributions are complicated to derive at an adequate level of accuracy, so they are approximated using the Monte Carlo method. A simulation experiment is carried out to investigate significance levels and power of each test, and compare them to the portmanteau test.     

Admissibility intervals for linear correlation coefficients

The correlation coefficient is widely used to quantify the degree of association between two quantitative variables. By resorting to the geometric representation of the linear correlation coefficient, it is possible to calculate the upper and lower bounds of the correlation coefficient between two variables x ₁,x ₂ when the correlation coefficients with a third variable x ₃ are available. Implications in observational studies, where x ₃ could be a proxy of a target variable x ₂, whose direct measurement is too expensive or impractical, are discussed.     

Characterizations based on conditional expectations

For given real functionsg andh, first we give necessary and sufficient conditions such that there exists a random variableX satisfying thatE(g(X)|X≥y)=h(y)r _x (y),∀y ∈ C _x , whereC _x andT _X are the support and the failure rate function ofX, respectively. These extend the results of Ruiz and Navarro (1994) and Ghitany et al. (1995). Next we investigate necessary and sufficient conditions such thath(y)=E(g(X)|X≥y), for a given functionh. Support for this research was provided in part by the National Science Council of the Republic of China, Grant No. NSC 86-2115-M-110-014 and NSC 88-2118-M-110-001     

Local asymptotic normality of intermediate and central order statistics

Consider n independent random variables Z_i,…, Z_n on R with common distribution function F, whose upper tail belongs to a parametric family F(t) = F_θ(t),t ≥ x₀, where θ ∈ ? ? R ^d. A necessary and sufficient condition for the family F_θ, θ ∈ ?, is established such that the k-th largest order statistic Z_n?k+1:n alone constitutes the central sequence yielding local asymptotic normality ( LAN ) of the loglikelihood ratio of the vector (Z_n?i+1:n)¹ _i=kof the k largest order statistics. This is achieved for k = k(n)→_n→∞∞ with k/n→_n→∞ 0.

In the case of vectors of central order statistics ( Z_r:n, Z_r+1:n,…, Z_s:n ), with r/n and s/n both converging to q ∈ ( 0,1 ), it turns out that under fairly general conditions any order statistic Z_m:n with r ≤ m ≤s builds the central sequence in a pertaining LAN expansion.These results lead to asymptotically optimal tests and estimators of the underlying parameter, which depend on single order statistics only     

The one-way analysis of variance with ordered alternatives: A modification of bartholomew's Ē2 Test

Tests of homogeneity of normal means with the alternative restricted by an ordering on the means are considered. The simply ordered case, μ₁ ≤ μ₂ ≤ ··· ≤ μ_k, and the simple tree ordering, μ₁ ≤ μ_j, for; j= 2, 3,…, k, are emphasized. A modification of the likelihood-ratio test is proposed which is asymptotically equivalent to it but is more robust to violations of the hypothesized orderings. The new test has power at the points satisfying the hypothesized ordering which is similar to that of the likelihood-ratio test provided the degrees of freedom are not too small. The modified test is shown to be unbiased and consistent.     

Estimation of Ordered Location Parameters: The Exponential Distribution

The problem of component wise estimation of ordered location parameters θ ₁ (θ ₁ ≤ θ ₂) of two independent exponential distributions is investigated. The scale parameters are assumed to be unequal but known. Independent random samples of unequal sample sizes are drawn from two populations and the estimators admissible among the mixed estimators of θ ₁ and θ ₂ are obtained. It is shown that the minimum risk estimators (MREs) of θ ₁ and θ ₂ without assuming θ ₁ ≤ θ ₂ are inadmissible when one does assume that θ ₁ ≤ θ ₂. The efficiencies of mixed estimators relative to MREs (without assuming θ ₁ ≤ θ ₂) are tabulated for equal sample sizes and equal scale parameters.