Semi-empirical likelihood ratio confidence intervals for the difference of two sample means

We all know that we can use the likelihood ratio statistic to test hypotheses and construct confidence intervals in full parametric models. Recently, Owen (1988,Biometrika,75, 237–249; 1990,Ann. Statist.,18, 90–120) has introduced the empirical likelihood method in nonparametric models. In this paper, we combine these two likelihoods together and use the likelihood ratio to construct confidence intervals in a semiparametric problem, in which one model is parametric, and the other is nonparametric. A version of Wilks's theorem is developed.     

A proof of independent Bartlett correctability of nested likelihood ratio tests

It is well known that likelihood ratio statistic is Bartlett correctable. We consider decomposition of a likelihood ratio statistic into 1 degree of freedom components based on sequence of nested hypotheses. We give a proof of the fact that the component likelihood ratio statistics are distributed mutually independently up to the order O(1/n) and each component is independently Bartlett correctable. This was implicit in Lawley (1956, Biometrika, 43, 295–303) and proved in Bickel and Ghosh (1990, Ann. Statist., 18, 1070–1090) using a Bayes method. We present a more direct frequentist proof.     

On the joint distribution of studentized order statistics

In this paper, the joint distribution of some special linear combinations of the (internally) studentized order statistics are derived for both normal and exponential populations; the exact relationship between their pdf's is also obtained. The exact sampling distributions of studentized extreme deviation statistic, which has been proposed by Pearson and Chandra Sekar (1936,Biometrika,28, 308–320), are derived for these two populations. An application to the most powerful location and scale invariant test is discussed briefly.     

Multinomial logistic regression algorithm

The lower bound principle (introduced in Böhning and Lindsay 1988, Ann. Inst. Statist. Math., 40, 641–663), Böhning (1989, Biometrika, 76, 375–383) consists of replacing the second derivative matrix by a global lower bound in the Loewner ordering. This bound is used in the Newton-Raphson iteration instead of the Hessian matrix leading to a monotonically converging sequence of iterates. Here, we apply this principle to the multinomial logistic regression model, where it becomes specifically attractive.Supplement to Monotonicity of quadratic-approximation algorithms by Böhning and Lindsay (1988). Ann. Inst. Statist. Math., 40, 641–663.This research was supported by the German Research Foundation.     

A Statistic for Testing the Null Hypothesis of Elliptical Symmetry

We present and study a procedure for testing the null hypothesis of multivariate elliptical symmetry. The procedure is based on the averages of some spherical harmonics over the projections of the scaled residual (1978, N. J. H. Small, Biometrika65, 657–658) of the d-dimensional data on the unit sphere of ^d. We find, under mild hypothesis, the limiting null distribution of the statistic presented, showing that, for an appropriate choice of the spherical harmonics included in the statistic, this distribution does not depend on the parameters that characterize the underlying elliptically symmetric law. We describe a bivariate simulation study that shows that the finite sample quantiles of our statistic converge fairly rapidly, with sample size, to the theoretical limiting quantiles and that our procedure enjoys good power against several alternatives.     

Some autoregressive moving average processes with generalized Poisson marginal distributions

Some simple models are introduced which may be used for modelling or generating sequences of dependent discrete random variables with generalized Poisson marginal distribution. Our approach for building these models is similar to that of the Poisson ARMA processes considered by Al-Osh and Alzaid (1987,J. Time Ser. Anal.,8, 261–275; 1988,Statist. Hefte,29, 281–300) and McKenzie (1988,Adv. in Appl. Probab.,20, 822–835). The models have the same autocorrelation structure as their counterparts of standard ARMA models. Various properties, such as joint distribution, time reversibility and regression behavior, for each model are investigated.     

Estimation in Linear Models with Random Effects and Errors-in-Variables

The independent variables of linear mixed models are subject to measurement errors in practice. In this paper, we present a unified method for the estimation in linear mixed models with errors-in-variables, based upon the corrected score function of Nakamura (1990, Biometrika, 77, 127–137). Asymptotic normality properties of the estimators are obtained. The estimators are shown to be consistent and convergent at the order of n ^–1/2. The performance of the proposed method is studied via simulation and the analysis of a data set on hedonic housing prices.     

Recurrence relations for single and product moments of progressive Type-II right censored order statistics from exponential and truncated exponential distributions

In this paper, we establish several recurrence relations satisfied by the single and product moments of progressive Type-II right censored order statistics from an exponential distribution. These relations may then be used, for example, to compute all the means, variances and covariances of exponential progressive Type-II right censored order statistics for all sample sizes n and all censoring schemes (R ₁, R ₂, ..., R _m ), mn. The results presented in the paper generalize the results given by Joshi (1978, Sankhy Ser. B, 39, 362–371; 1982, J. Statist. Plann. Inference, 6, 13–16) for the single moments and product moments of order statistics from the exponential distribution.To further generalize these results, we consider also the right truncated exponential distribution. Recurrence relations for the single and product moments are established for progressive Type-II right censored order statistics from the right truncated exponential distribution.     

Relationships for moments of order statistics from the right-truncated generalized half logistic distribution

In this paper, we establish several recurrence relations satisfied by the single and the product moments for order statistics from the right-truncated generalized half logistic distribution. These relationships may be used in a simple recursive manner in order to compute the single and the product moments of all order statistics for all sample sizes and for any choice of the truncation parameter P. These generalize the corresponding results for the generalized half logistic distribution derived recently by Balakrishnan and Sandhu (1995, J. Statist. Comput. Simulation, 52, 385–398).Earlier went by the name R. A. Sandhu.     

Rank estimates in a class of semiparametric two-sample models

We consider a two-sample semiparametric model involving a real parameter and a nuisance parameter F which is a distribution function. This model includes the proportional hazard, proportional odds, linear transformation and Harrington-Fleming models (1982, Biometrika, 69, 533–546). We propose two types of estimates based on ranks. The first is a rank approximation to Huber's M-estimates (1981, Robust Statistics, Wiley) and the second is a Hodges-Lehmann type rank inversion estimate (1963, Ann. Math. Statist., 34, 598–611). We obtain asymptotic normality and efficiency results. The estimates are consistent and asymptotically normal generally but fully efficient only for special cases.Research partially supported by National Science Foundation Grant DMS-86-02083 and National Institute of General Medical Sciences Grant SSS-Y1RO1GM35416-01     

Recurrence relations among moments of order statistics from two related sets of independent and non-identically distributed random variables

Some recurrence relations among moments of order statistics from two related sets of variables are quite well-known in the i.i.d. case and are due to Govindarajulu (1963a, Technometrics, 5, 514–518 and 1966, J. Amer. Statist. Assoc., 61, 248–258). In this paper, we generalize these results to the case when the order statistics arise from two related sets of independent and non-identically distributed random variables. These relations can be employed to simplify the evaluation of the moments of order statistics in an outlier model for symmetrically distributed random variables.     

A Notion of α-Monotonicity with Generalized Multiplications

The multiplications of van Harn et al. (1982, Z. Wahrsch. Verw. Gebiete, 61, 97–118) are used to generalize the definition of -monotonicity of Olshen and Savage (1970, J. Appl. Probab., 7, 21–34) and Steutel (1988, Statist. Neerlandica, 42, 137–140) for distributions with support in Z ₊ and R ₊. Several characterizations are offered and a convolution property is established. Some relevant stability equations are solved and a relationship with the important concept of self-decomposability is noted. Poisson mixtures are used to deduce results for the R ₊-case from those for the Z ₊-case.     

Variable location and scale kernel density estimation

Variable (bandwidth) kernel density estimation (Abramson (1982,Ann. Statist.,10, 1217–1223)) and a kernel estimator with varying locations (Samiuddin and El-Sayyad (1990,Biometrika,77, 865–874)) are complementary ideas which essentially both afford bias of orderh ⁴ as the overall smoothing parameterh 0, sufficient differentiability of the density permitting. These ideas are put in a more general framework in this paper. This enables us to describe a variety of ways in which scale and location variation may be extended and/or combined to good theoretical effect. This particularly includes extending the basic ideas to provide new kernel estimators with bias of orderh ⁶. Technical difficulties associated with potentially overly large variations are fully accounted for in our theory.     

Laurent–Padé Approximants to Four Kinds of Chebyshev Polynomial Expansions. Part I. Maehly Type Approximants

Laurent Padé–Chebyshev rational approximants, A _m (z,z ^–1)/B _n (z,z ^–1), whose Laurent series expansions match that of a given function f(z,z ^–1) up to as high a degree in z,z ^–1 as possible, were introduced for first kind Chebyshev polynomials by Clenshaw and Lord [2] and, using Laurent series, by Gragg and Johnson [4]. Further real and complex extensions, based mainly on trigonometric expansions, were discussed by Chisholm and Common [1]. All of these methods require knowledge of Chebyshev coefficients of f up to degree m+n. Earlier, Maehly [5] introduced Padé approximants of the same form, which matched expansions between f(z,z ^–1)B _n (z,z ^–1) and A _m (z,z ^–1). The derivation was relatively simple but required knowledge of Chebyshev coefficients of f up to degree m+2n. In the present paper, Padé–Chebyshev approximants are developed not only to first, but also to second, third and fourth kind Chebyshev polynomial series, based throughout on Laurent series representations of the Maehly type. The procedures for developing the Padé–Chebyshev coefficients are similar to that for a traditional Padé approximant based on power series [8] but with essential modifications. By equating series coefficients and combining equations appropriately, a linear system of equations is successfully developed into two sub-systems, one for determining the denominator coefficients only and one for explicitly defining the numerator coefficients in terms of the denominator coefficients. In all cases, a type (m,n) Padé–Chebyshev approximant, of degree m in the numerator and n in the denominator, is matched to the Chebyshev series up to terms of degree m+n, based on knowledge of the Chebyshev coefficients up to degree m+2n. Numerical tests are carried out on all four Padé–Chebyshev approximants, and results are outstanding, with some formidable improvements being achieved over partial sums of Laurent–Chebyshev series on a variety of functions. In part II of this paper [7] Padé–Chebyshev approximants of Clenshaw–Lord type will be developed for the four kinds of Chebyshev series and compared with those of the Maehly type.     

Score and Wald tests for the multivariate growth curve model with missing data

We present the score and Wald test analogues to Srivastava's (1985, Comm. Statist. A—Theory Methods, 14, 775–792) likelihood ratio tests for the multivariate growth curve model with missing data, and illustrate their use with data from an immunotherapy experiment (Fukushima et al. (1982, Int. J. Cancer, 29, 107–112, 113–117)).     

Influence Analysis for Linear Measurement Error Models

In this paper, we present a unified diagnostic method for linear measurement error models based upon the corrected likelihood of Nakamura (1990, Biometrika, 77, 127–137). Both global influence and local influence are discussed. The case-deletion model and mean-shift outlier model are considered, and they are shown to be approximately equivalent. Several diagnostic measures are derived and discussed. It is found that they can be written in terms of the residual and leverage measure. Some existing results are improved. Numerical example illustrates that our method is useful for diagnosing influential observations.     

Truncated selection procedures for the most probable event and the least probable event

This paper proposes two sequential procedures for selecting respectively the multinomial cell with the largest cell probability and the multinomial cell with the smallest cell probability. The stopping rule for both procedures uses truncation of the procedure studied by Ramey and Alam (1979, Biometrika, 66, 171–173). A property of the least favorable configuration of the proposed procedures is proved, which partially solves a conjecture given in Ramey and Alam (1979). The proposed procedures are compared with other procedures which have been considered in the literature and are found to be better in certain respects.     

Nonparametric test of restricted interchangeability

Exact and large sample distributions of the rank order test under the null hypothesis of restricted interchangeability are obtained. Under given regularity conditions and under Pitman's shift in location alternative, the asymptotic relative efficiency of this nonparametric test in comparison with Votaw's (1948, Ann. Math. Statist., 19, 447–473) likelihood ratio test is given.     

Empirical likelihood for partial linear models

In this paper the empirical likelihood method due to Owen (1988,Biometrika,75, 237–249) is applied to partial linear random models. A nonparametric version of Wilks' theorem is derived. The theorem is then used to construct confidence regions of the parameter vector in the partial linear models, which has correct asymptotic coverage. A simulation study is conducted to compare the empirical likelihood and normal approximation based method. Research supported by NNSF of China and a grant to the first author for his excellent Ph.D. dissertation work in China. Research supported by Hong Kong RGC CERG No. HKUST6162/97P.     

Bayesian Prediction Analysis for Growth Curve Model Using Noninformative Priors

We apply a Bayesian approach to the problem of prediction in an unbalanced growth curve model using noninformative priors. Due to the complexity of the model, no analytic forms of the predictive densities are available. We propose both approximations and a prediction-oriented Metropolis-Hastings sampling algorithm for two types of prediction, namely the prediction of future observations for a new subject and the prediction of future values for a partially observed subject. They are illustrated and compared through real data and simulation studies. Two of the approximations compare favorably with the approximation in Fearn (1975, Biometrika, 62, 89–100) and are very comparable to the more accurate Rao-Blackwellization from Metropolis-Hastings sampling algorithm.