Search

We attempt to identify the 25 most-cited statistical papers, providing some brief commentary on each paper on our list. This list consists, to a great extent, of papers that are on non-parametric methods, have applications in the life sciences, or deal with the multiple comparisons problem. We also list the most-cited papers published in 1993 or later. In contrast to the overall most-cited papers, these are predominately papers on Bayesian methods and wavelets. We briefly discuss some of the issues involved in the use of citation counts. [ABSTRACT FROM AUTHOR]

The mean residual lifetime (MRL) of a unit is its expected additional lifetime provided that it has survived until time t. The MRL estimation problem has been frequently addressed in the literature since it has wide applications in statistics, reliability and survival analysis. In this paper, we consider the problem of estimating the MRL in ranked set sampling when actual quantifications of a concomitant variable are available. To exploit the additional information of the concomitant variable, we introduce several MRL estimators based on some regression techniques. We then compare them with the standard MRL estimator in simple random sampling using Monte Carlo simulation and a real dataset from the Surveillance, Epidemiology, and End Results Program. Our results indicate the superiority of the procedures that we have developed when the quality of ranking is fairly good. [ABSTRACT FROM AUTHOR]

For checking time series models, the Ljung–Box, Li–Mak and Zhu–Wang statistics play an important role, which use the Pearson's correlation coefficient to implement (squared) residual (partial) autocorrelation tests. In this paper, we replace the Pearson's correlation coefficient with a new rank correlation coefficient and propose a new test statistic to conduct diagnostic checks for residuals in autoregressive moving average models, autoregressive conditional heteroscedasticity models and integer-valued time series models, respectively. We conduct simulations to assess the performance of the new test statistic, and compare it with existing ones, and the results show the superiority of the proposed one. We use three real examples to exhibit its usefulness. [ABSTRACT FROM AUTHOR]

In this paper, Bayesian estimates are derived for the location and scale parameters of the Laplace distribution based on complete, Type-I, and Type-II censored samples under different prior settings. Subsequently, Bayesian point and interval estimates, as well as the associated statistical inference, are discussed in detail. The developed methods are then applied to two real data sets for illustrative purposes. Moreover, a detailed Monte Carlo simulation study is carried out for evaluating the performance of the inferential methods developed here. Finally, we provide a brief discussion of the established results to demonstrate their practical utility and present some associated problems of further interest. Overall, this study fills an existing gap in the development of Bayesian inferential techniques for the parameters of the two-parameter Laplace distribution, making this research innovative and offering more investigative implications. It showcases the potential for broader methodological applications of Bayesian inference for complex real-world data sets, especially in scenarios involving different forms of censoring. This research provides a critical tool for statistical analysis in different fields such as engineering and finance, where the Laplace distribution is frequently adopted as a fundamental model. [ABSTRACT FROM AUTHOR]

The C statistic, also known as the Cash statistic, is often used in astronomy for the analysis of low-count Poisson data. The main advantage of this statistic, compared to the more commonly used χ 2 statistic, is its applicability without the need to combine data points. This feature has made the C statistic a very useful method to analyze Poisson data that have small (or even null) counts in each resolution element. One of the challenges of the C statistic is that its probability distribution, under the null hypothesis that the data follow a parent model, is not known exactly. This paper presents an effort towards improving our understanding of the C statistic by studying (a) the distribution of C statistic for a fully specified model, (b) the distribution of Cmin resulting from a maximum-likelihood fit to a simple one-parameter constant model, i.e. a model that represents the sample mean of N Poisson measurements, and (c) the distribution of the associated Δ C statistic that is used for parameter estimation. The results confirm the expectation that, in the high-count limit, both C statistic and Cmin have the same mean and variance as a χ 2 statistic with same number of degrees of freedom. It is also found that, in the low-count regime, the expectation of the C statistic and Cmin can be substantially lower than for a χ 2 distribution. The paper makes use of recent X-ray observations of the astronomical source PG 1116+215 to illustrate the application of the C statistic to Poisson data. [ABSTRACT FROM AUTHOR]

The Cash statistic, also known as the C statistic, is commonly used for the analysis of low-count Poisson data, including data with null counts for certain values of the independent variable. The use of this statistic is especially attractive for low-count data that cannot be combined, or re-binned, without loss of resolution. This paper presents a new maximum-likelihood solution for the best-fit parameters of a linear model using the Poisson-based Cash statistic. The solution presented in this paper provides a new and simple method to measure the best-fit parameters of a linear model for any Poisson-based data, including data with null counts. In particular, the method enforces the requirement that the best-fit linear model be non-negative throughout the support of the independent variable. The method is summarized in a simple algorithm to fit Poisson counting data of any size and counting rate with a linear model, by-passing entirely the use of the traditional χ 2 statistic. [ABSTRACT FROM AUTHOR]

In recent years, some interested data can be recorded only if the values fall within an interval range, and the responses are often subject to censoring. Attempting to perform effective statistical analysis with censored, especially heavy-tailed and asymmetric data, can be difficult. In this paper, we develop a novel linear regression model based on the proposed skewed generalized t distribution for censored data. The likelihood-based inference and diagnostic analysis are established using the Expectation/Conditional Maximization Either algorithm in conjunction with smoothing approximate functions. We derive relevant measures to perform global influence for this novel model and develop local influence analysis based on the conditional expectation of the complete-data log-likelihood function. Some useful perturbation schemes are discussed. We illustrate the finite sample performance and the robustness of the proposed method by simulation studies. The proposed model is compared with other procedures based on a real dataset, and a sensitivity analysis is also conducted. [ABSTRACT FROM AUTHOR]

Researchers in statistical shape analysis often analyze outlines of objects. Even though these contours are infinite-dimensional in theory, they must be discretized in practice. When discretizing, it is important to reduce the number of sampling points considerably to reduce computational costs, but to not use too few points so as to result in too much approximation error. Unfortunately, determining the minimum number of points needed to achieve sufficiently approximate the contours is computationally expensive. In this paper, we fit regression models to predict these lower bounds using characteristics of the contours that are computationally cheap as predictor variables. However, least squares regression is inadequate for this task because it treats overestimation and underestimation equally, but underestimation of lower bounds is far more serious. Instead, to fit the models, we use the LINEX loss function, which allows us to penalize underestimation at an exponential rate while penalizing overestimation only linearly. We present a novel approach to select the shape parameter of the loss function and tools for analyzing how well the model fits the data. Through validation methods, we show that the LINEX models work well for reducing the underestimation for the lower bounds. [ABSTRACT FROM AUTHOR]

Discrete data are collected in many application areas and are often characterised by highly-skewed distributions. An example of this, which is considered in this paper, is the number of visits to a specialist, often taken as a measure of demand in healthcare. A discrete Weibull regression model was recently proposed for regression problems with a discrete response and it was shown to possess desirable properties. In this paper, we propose the first Bayesian implementation of this model. We consider a general parametrization, where both parameters of the discrete Weibull distribution can be conditioned on the predictors, and show theoretically how, under a uniform non-informative prior, the posterior distribution is proper with finite moments. In addition, we consider closely the case of Laplace priors for parameter shrinkage and variable selection. Parameter estimates and their credible intervals can be readily calculated from their full posterior distribution. A simulation study and the analysis of four real datasets of medical records show promises for the wide applicability of this approach to the analysis of count data. The method is implemented in the R package BDWreg. [ABSTRACT FROM AUTHOR]

Adaptive designs for multi-armed clinical trials have become increasingly popular recently because of their potential to shorten development times and to increase patient response. However, developing response-adaptive designs that offer patient-benefit while ensuring the resulting trial provides a statistically rigorous and unbiased comparison of the different treatments included is highly challenging. In this paper, the theory of Multi-Armed Bandit Problems is used to define near optimal adaptive designs in the context of a clinical trial with a normally distributed endpoint with known variance. We report the operating characteristics (type I error, power, bias) and patient-benefit of these approaches and alternative designs using simulation studies based on an ongoing trial. These results are then compared to those recently published in the context of Bernoulli endpoints. Many limitations and advantages are similar in both cases but there are also important differences, specially with respect to type I error control. This paper proposes a simulation-based testing procedure to correct for the observed type I error inflation that bandit-based and adaptive rules can induce. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the mean change-point models based on associated sequences. Under some weak conditions, we obtain a limit distribution of CUSUM statistic which can be used to judge the mean change-mount δ n is satisfied or dissatisfied n 1 / 2 δ n = o (1). We also study the consistency of sample covariances and change-point location statistics. Based on Normality and Lognormality data, some simulations such as empirical sizes, empirical powers and convergence are presented to test our results. As an important application, we use CUSUM statistics to do the mean change-point analysis for a financial series. [ABSTRACT FROM AUTHOR]

The aim of this paper is to provide some practical aspects of point and interval estimates of the global maximum of a function using extreme value theory. Consider a real-valued function f:D→ defined on a bounded interval D such that f is either not known analytically or is known analytically but has rather a complicated analytic form. We assume that f possesses a global maximum attained, say, at u*∈D with maximal value x*=max u f(u)≐f(u*). The problem of seeking the optimum of a function which is more or less unknown to the observer has resulted in the development of a large variety of search techniques. In this paper we use the extreme-value approach as appears in Dekkers et al. [A moment estimator for the index of an extreme-value distribution, Ann. Statist. 17 (1989), pp. 1833-1855] and de Haan [Estimation of the minimum of a function using order statistics, J. Amer. Statist. Assoc. 76 (1981), pp. 467-469]. We impose some Lipschitz conditions on the functions being investigated and through repeated simulation-based samplings, we provide various practical interpretations of the parameters involved as well as point and interval estimates for x*. [ABSTRACT FROM AUTHOR]

A special source of difficulty in the statistical analysis is the possibility that some subjects may not have a complete observation of the response variable. Such incomplete observation of the response variable is called censoring. Censorship can occur for a variety of reasons, including limitations of measurement equipment, design of the experiment, and non-occurrence of the event of interest until the end of the study. In the presence of censoring, the dependence of the response variable on the explanatory variables can be explored through regression analysis. In this paper, we propose to examine the censorship problem in context of the class of asymmetric, i.e., we have proposed a linear regression model with censored responses based on skew scale mixtures of normal distributions. We develop a Monte Carlo EM (MCEM) algorithm to perform maximum likelihood inference of the parameters in the proposed linear censored regression models with skew scale mixtures of normal distributions. The MCEM algorithm has been discussed with an emphasis on the skew-normal, skew Student-t-normal, skew-slash and skew-contaminated normal distributions. To examine the performance of the proposed method, we present some simulation studies and analyze a real dataset. [ABSTRACT FROM AUTHOR]

Ebegil, Özdemir and Gökpinar [[14]] adopted the Ridge and Liu-type estimators into some shrinkage estimators using a median ranked set sample in the presence of multicollinearity. Sancar and Inan [[15]] proposed a Liu-type logistic estimator as a two-parameter estimator to overcome the multicollinearity problem in the logistic model to insist that using the two-phase approaches is not guaranteed to overcome multicollinearity. They applied those approaches to the 2017 Annual Industry and Service Statistics data and found that restricted calibration estimators were more efficient than the generalized regression estimator in estimating the variables with high variance such as turnover. For classical estimation methods, they used maximum likelihood estimators, moment's estimators, least squares estimators, and weighted least squares estimators, percentile estimators, the ordinary least squares estimators, best linear unbiased estimators, L-moments estimators, trimmed L-moments estimators, and Bain and Engelhardt [[3]] estimators. [Extracted from the article]

From the prediction viewpoint, mode regression is more attractive since it pay attention to the most probable value of response variable given regressors. On the other hand, high-dimensional data are very prevalent as the advance of the technology of collecting and storing data. Variable selection is an important strategy to deal with high-dimensional regression problem. This paper aims to propose a variable selection procedure for high-dimensional mode regression via combining nonparametric kernel estimation method with sparsity penalty tactics. We also establish the asymptotic properties under certain technical conditions. The effectiveness and flexibility of the proposed methods are further illustrated by numerical studies and the real data application. [ABSTRACT FROM AUTHOR]

Forecasting of future snow depths is useful for many applications like road safety, winter sport activities, avalanche risk assessment and hydrology. Motivated by the lack of statistical forecasts models for snow depth, in this paper we present a set of models to fill this gap. First, we present a model to do short-term forecasts when we assume that reliable weather forecasts of air temperature and precipitation are available. The covariates are included nonlinearly into the model following basic physical principles of snowfall, snow aging and melting. Due to the large set of observations with snow depth equal to zero, we use a zero-inflated gamma regression model, which is commonly used to similar applications like precipitation. We also do long-term forecasts of snow depth and much further than traditional weather forecasts for temperature and precipitation. The long-term forecasts are based on fitting models to historic time series of precipitation, temperature and snow depth. We fit the models to data from six locations in Norway with different climatic and vegetation properties. Forecasting five days into the future, the results showed that, given reliable weather forecasts of temperature and precipitation, the forecast errors in absolute value was between 3 and 7 cm for different locations in Norway. Forecasting three weeks into the future, the forecast errors were between 7 and 16 cm. [ABSTRACT FROM AUTHOR]

In this paper, we propose a two-stage functional principal component analysis method in age-period-cohort (APC) analysis. The first stage of the method considers the age-period effect with the fitted values treated as an offset; and the second stage of the method considers the residual age-cohort effect conditional on the already estimated age-period effect. An APC version of the model in functional data analysis provides an improved fit to the data, especially when the data are sparse and irregularly spaced. We demonstrate the effectiveness of the proposed method using body mass index data stratified by gender and ethnicity. [ABSTRACT FROM AUTHOR]

A key challenge in rainfall estimation is spatio-temporal variablility. Weather radars are used to estimate precipitation with high spatial and temporal resolution. Due to the inherent errors in radar estimates, spatial interpolation has been often employed to calibrate the estimates. Kriging is a simple and popular spatial interpolation method, but the method has several shortcomings. In particular, the prediction is quite unstable and often fails to be performed when sample size is small. In this paper, we proposed a flexible and efficient spatial interpolator for radar rainfall estimation, with several advantages over kriging. The method is illustrated using a real-world data set. [ABSTRACT FROM AUTHOR]

The standard regression model designed for real space is not suitable for compositional variables; it should be considered, whether the response and/or covariates are of compositional nature. There are usually three types of multiple regression model with compositional variables: Type 1 refers to the case where all the covariates are compositional data and the response is real; Type 2 is the opposite of Type 1; Type 3 relates to the model with compositional response and covariates. There have been some models for the three types. In this paper, we focus on Type 3 and propose multiple linear regression models including model in the simplex and model in isometric log-ratio (ilr) coordinates. The model in the simplex is based on matrix product, which can project a-part composition to another-part composition, and can deal with different number of parts of compositional variables. Some theorems are given to point out the relationship of parameters between the proposed models. Moreover, the inference for parameters in proposed models is also given. Real example is studied to verify the validity and usefulness of proposed models. [ABSTRACT FROM AUTHOR]

Dispersion indexes with respect to the Poisson and binomial distributions are widely used to assess the conformity of the underlying distribution from an observed sample of the count with one or the other of these theoretical distributions. Recently, the exponential variation index has been proposed as an extension to nonnegative continuous data. This paper aims to gather to study the unified definition of these indexes with respect to the relative variability of a nonnegative natural exponential family of distributions through its variance function. We establish the strong consistency of the plug-in estimators of the indexes as well as their asymptotic normalities. Since the exact distributions of the estimators are not available in closed form, we consider the test of hypothesis relying on these estimators as test statistics with their asymptotic distributions. Simulation studies globally suggest good behaviours of these tests of hypothesis procedures. Applicable examples are analysed, including the lesser-known references such as negative binomial and inverse Gaussian, and improving the very usual case of the Poisson dispersion index. Concluding remarks are made with suggestions of possible extensions. [ABSTRACT FROM AUTHOR]

Application of the exact statistical inference frequently leads to non-standard probability distributions of the considered estimators or test statistics. The exact distributions of many estimators and test statistics can be specified by their characteristic functions, as is the case for the null distribution of the Bartlett's test statistic. However, analytical inversion of the characteristic function, if possible, frequently leads to complicated expressions for computing the distribution function and the corresponding quantiles. An efficient alternative is the well-known method based on numerical inversion of the characteristic functions, which is, however, ignored in popular statistical software packages. In this paper, we present the explicit characteristic function of the corrected Bartlett's test statistic together with the computationally fast and efficient implementation of the approach based on numerical inversion of this characteristic function, suggested for evaluating the exact null distribution used for testing homogeneity of variances in several normal populations, with possibly unequal sample sizes. [ABSTRACT FROM AUTHOR]

Deciding on the best statistical method to apply when the response variable is ordinal is essential because the way the categories are ordered in the data is relevant as it could change the results of the analysis. Although the models for continuous variables have similarities to those for ordinal variables, this paper presents the advantages of the use of the ordering information on the outcomes with methods developed for modeling ordinal data such as the ordered stereotype model. The novelty of this article lies in showing the dangers of assigning equally spaced scores to ordered response categories in statistical analysis, which are illustrated with a simulation study and a case study. We propose a new way to use the score parameters, which incorporates the fitted spacing dictated by the data. Additionally, this article uses score parameter estimates in the ordered stereotype model to propose a new measure to calculate continuous medians in the raw data: the adjusted c-median. It benefits the general audience who can easily understand the median as a summary statistic. Supplementary materials for this article are available online. [ABSTRACT FROM AUTHOR]

This paper introduces a new two-parameter exponentiated discrete Lindley distribution. A wide range of its structural properties are investigated. This includes the shape of the probability mass function, hazard rate function, moments, skewness, kurtosis, stress–strength reliability, mean residual lifetime, mean past lifetime, order statistics and L-moment statistics. The hazard rate function can be increasing, decreasing, decreasing–increasing–decreasing, increasing–decreasing–increasing, unimodal, bathtub, and J-shaped depending on its parameters values. Two methods are used herein to estimate the model parameters, namely, the maximum likelihood, and the proportion. A detailed simulation study is carried out to examine the bias and mean square error of maximum likelihood and proportion estimators. The flexibility of the proposed model is explained by using four distinctive data sets. It can serve as an alternative model to other lifetime distributions in the existing statistical literature for modeling positive real data in many areas. [ABSTRACT FROM AUTHOR]

In this paper, we are employing the generalized linear model (GLM) in the form lij=χλ to decompose the symmetry model into the class of models discussed in Tomizawa (1992). In this formulation, the random component would be the observed counts fij with an underlying Poisson distribution. This approach utilizes the non-standard log-linear model and our focus in this paper therefore relates to models that are decompositions of the complete symmetry model. That is, models that are implied by the symmetry models. We develop factor and regression variables required for the implementation of these models in SAS PROC GENMOD and SPSS PROC GENLOG. We apply this methodology to analyse the three 4 × 4 contingency table, one of which is the Japanese Unaided distance vision data. Results obtained in this study are consistent with those from the numerous literature on the subject. We further extend our applications to the 6 × 6 Brazilian social mobility data. We found that both the quasi linear diagonal-parameters symmetry (QLDPS) and the quasi 2-ratios parameter symmetry (Q2RPS) models fit the Brazilian data very well. Parsimonious models being the QLDPS and the quasi-conditional symmetry (QCS) models. The SAS and SPSS programs for implementing the models discussed in this paper are presented in Appendices A, B and C. [ABSTRACT FROM AUTHOR]

This paper deals with the well-studied problem of how best to estimate the number of mutually exclusive and exhaustive classes in a population, based on a sample from it. Haas & Stokes review and provide non-parametric approaches, but there are associated difficulties especially for small sampling fractions and/or widely varying population class sizes. Sichel provided 'GIGP' methodology, for this problem and for other purposes; this paper utilizes the three-parameter GIGP distribution for this problem, and also for the estimation of the number of classes of size 1, as an alternative to the non-parametric approaches. Methodological and computational issues are considered, and examples indicate the potential for GIGP. [ABSTRACT FROM AUTHOR]

Until now, various acceptance reliability sampling plans have been developed based on different life tests of items. However, the statistical effect of the acceptance sampling tests on the reliability characteristic of the lots accepted in the test has not been appropriately addressed. In this paper, we deal with an acceptance reliability sampling plan under a ‘general framework’ and discuss the corresponding statistical effect of the acceptance sampling tests. The lifetime of the population before the acceptance test and that of population ‘conditional on the acceptance’ in the sampling test are stochastically compared. The improvement of reliability characteristics of the population conditional on the acceptance in the sampling test is precisely analyzed. [ABSTRACT FROM AUTHOR]

This paper uses data for 44 countries from 1970-1990, to investigate the relationship between economic growth and carbon dioxide emissions. Empirical results are obtained from a structural model from the empirical growth literature modified to include environmental 'bads'. Results suggest that richer countries exhibit technical progress in a way that economizes on carbon dioxide emissions but that poorer countries do not. Furthermore, there is no indication that the growth process is leading poorer countries to move towards the adoption of the same pollution-ameliorating technology as characterizes richer countries. [ABSTRACT FROM AUTHOR]

Cluster analysis is one of the most widely used method in statistical analyses, in which homogeneous subgroups are identified in a heterogeneous population. Due to the existence of the continuous and discrete mixed data in many applications, so far, some ordinary clustering methods such as, hierarchical methods, k-means and model-based methods have been extended for analysis of mixed data. However, in the available model-based clustering methods, by increasing the number of continuous variables, the number of parameters increases and identifying as well as fitting an appropriate model may be difficult. In this paper, to reduce the number of the parameters, for the model-based clustering mixed data of continuous (normal) and nominal data, a set of parsimonious models is introduced. Models in this set are extended, using the general location model approach, for modeling distribution of mixed variables and applying factor analyzer structure for covariance matrices. The ECM algorithm is used for estimating the parameters of these models. In order to show the performance of the proposed models for clustering, results from some simulation studies and analyzing two real data sets are presented. [ABSTRACT FROM AUTHOR]

This paper considers the problem of making statistical inferences about a parameter when a narrow interval centred at a given value of the parameter is considered special, which is interpreted as meaning that there is a substantial degree of prior belief that the true value of the parameter lies in this interval. A clear justification of the practical importance of this problem is provided. The main difficulty with the standard Bayesian solution to this problem is discussed and, as a result, a pseudo-Bayesian solution is put forward based on determining lower limits for the posterior probability of the parameter lying in the special interval by means of a sensitivity analysis. Since it is not assumed that prior beliefs necessarily need to be expressed in terms of prior probabilities, nor that post-data probabilities must be Bayesian posterior probabilities, hybrid methods of inference are also proposed that are based on specific ways of measuring and interpreting the classical concept of significance. The various methods that are outlined are compared and contrasted at both a foundational level, and from a practical viewpoint by applying them to real data from meta-analyses that appeared in a well-known medical article. [ABSTRACT FROM AUTHOR]

The statistical modeling of big data bases constitutes one of the most challenging issues, especially nowadays. The issue is even more critical in case of a complicated correlation structure. Variable selection plays a vital role in statistical analysis of large data bases and many methods have been proposed so far to deal with the aforementioned problem. One of such methods is the Sure Independence Screening which has been introduced to reduce dimensionality to a relatively smaller scale. This method, though simple, produces remarkable results even under both ultra high dimensionality and big scale in terms of sample size problems. In this paper we dealt with the analysis of a big real medical data set assuming a Poisson regression model. We support the analysis by conducting simulated experiments taking into consideration the correlation structure of the design matrix. [ABSTRACT FROM AUTHOR]

In many practical situations, a statistical practitioner often faces a problem of classifying an object from one of the segmented (or screened) populations where the segmentation was conducted by a set of screening variables. This paper addresses this problem, proposing and studying yet another optimal rule for classification with segmented populations. A class ofq-dimensional rectangle-screened elliptically contoured (RSEC) distributions is considered for flexibly modeling the segmented populations. Based on the properties of the RSEC distributions, a parametric procedure for the segmented classification analysis (SCA) is proposed. This includes motivation for the SCA as well as some theoretical propositions regarding its optimal rule and properties. These properties allow us to establish other important results which include an efficient estimation of the rule by the Monte Carlo expectation–conditional maximization algorithm and an optimal variable selection procedure. Two numerical examples making use of utilizing a simulation study and a real dataset application and advocating the SCA procedure are also provided. [ABSTRACT FROM PUBLISHER]

Estimation of the lifetime distribution of industrial components and systems yields very important information for manufacturers and consumers. However, obtaining reliability data is time consuming and costly. In this context, degradation tests are a useful alternative approach to lifetime and accelerated life tests in reliability studies. The approximate method is one of the most used techniques for degradation data analysis. It is very simple to understand and easy to implement numerically in any statistical software package. This paper uses time series techniques in order to propose a modified approximate method (MAM). The MAM improves the standard one in two aspects: (1) it uses previous observations in the degradation path as a Markov process for future prediction and (2) it is not necessary to specify a parametric form for the degradation path. Characteristics of interest such as mean or median time to failure and percentiles, among others, are obtained by using the modified method. A simulation study is performed in order to show the improved properties of the modified method over the standard one. Both methods are also used to estimate the failure time distribution of the fatigue-crack-growth data set. [ABSTRACT FROM PUBLISHER]

In this paper, we consider statistical diagnostic for non-parametric regression models with right-censored data based on empirical likelihood. First, the primary model is transformed to the non-parametric regression model. Then, based on empirical likelihood methodology, we define some diagnostic statistics. At last, some simulation studies show that our proposed procedure can work fairly well. [ABSTRACT FROM PUBLISHER]

The first step in statistical analysis is the parameter estimation. In multivariate analysis, one of the parameters of interest to be estimated is the mean vector. In multivariate statistical analysis, it is usually assumed that the data come from a multivariate normal distribution. In this situation, the maximum likelihood estimator (MLE), that is, the sample mean vector, is the best estimator. However, when outliers exist in the data, the use of sample mean vector will result in poor estimation. So, other estimators which are robust to the existence of outliers should be used. The most popular robust multivariate estimator for estimating the mean vector is S-estimator with desirable properties. However, computing this estimator requires the use of a robust estimate of mean vector as a starting point. Usually minimum volume ellipsoid (MVE) is used as a starting point in computing S-estimator. For high-dimensional data computing, the MVE takes too much time. In some cases, this time is so large that the existing computers cannot perform the computation. In addition to the computation time, for high-dimensional data set the MVE method is not precise. In this paper, a robust starting point for S-estimator based on robust clustering is proposed which could be used for estimating the mean vector of the high-dimensional data. The performance of the proposed estimator in the presence of outliers is studied and the results indicate that the proposed estimator performs precisely and much better than some of the existing robust estimators for high-dimensional data. [ABSTRACT FROM PUBLISHER]

Beta regression models provide an adequate approach for modeling continuous outcomes limited to the interval (0, 1). This paper deals with an extension of beta regression models that allow for explanatory variables to be measured with error. The structural approach, in which the covariates measured with error are assumed to be random variables, is employed. Three estimation methods are presented, namely maximum likelihood, maximum pseudo-likelihood and regression calibration. Monte Carlo simulations are used to evaluate the performance of the proposed estimators and the naïve estimator. Also, a residual analysis for beta regression models with measurement errors is proposed. The results are illustrated in a real data set. [ABSTRACT FROM PUBLISHER]

This paper demonstrates how to plan a contingent valuation experiment to assess the value of ecologically produced clothes. First, an appropriate statistical model (the trinomial spike model) that describes the probability that a randomly selected individual will accept any positive bid, and if so, will accept the bid A, is defined. Secondly, an optimization criterion that is a function of the variances of the parameter estimators is chosen. However, the variances of the parameter estimators in this model depend on the true parameter values. Pilot study data are therefore used to obtain estimates of the parameter values and a locally optimal design is found. Because this design is only optimal given that the estimated parameter values are correct, a design that minimizes the maximum of the criterion function over a plausable parameter region (i.e. a minimax design) is then found. [ABSTRACT FROM AUTHOR]

In this paper, we consider the notions of data depth for ordering multivariate data and propose a classification rule based on the distribution of some depth functions in . The equivalence of the proposed classification rule to optimal Bayes rule is discussed under suitable conditions. The performance of the proposed classification method is investigated in low- and high-dimensional setting using real datasets. Also, the performance of the proposed classification method is illustrated in comparison to some other depth-based classifiers using simulated data sets. [ABSTRACT FROM AUTHOR]

Since the pioneering work by Koenker and Bassett [27], quantile regression models and its applications have become increasingly popular and important for research in many areas. In this paper, a random effects ordinal quantile regression model is proposed for analysis of longitudinal data with ordinal outcome of interest. An efficient Gibbs sampling algorithm was derived for fitting the model to the data based on a location-scale mixture representation of the skewed double-exponential distribution. The proposed approach is illustrated using simulated data and a real data example. This is the first work to discuss quantile regression for analysis of longitudinal data with ordinal outcome. [ABSTRACT FROM AUTHOR]

CD4 and viral load play important roles in HIV/AIDS studies, and the study of their relationship has received much attention with well-known results. However, AIDS datasets are often highly complex in the sense that they typically contain outliers, measurement errors, and missing data. These data complications can greatly affect statistical analysis results, but much of the literature fail to address these issues in data analysis. In this paper, we re-visit the important relationship between CD4 and viral load and propose methods which simultaneously address outliers, measurement errors, and missing data. We find that the strength of the relationship may be severely mis-estimated if measurement errors and outliers are ignored. The proposed methods are general and can be used in other settings, where jointly modelling several different types of longitudinal data is required in the presence of data complications. [ABSTRACT FROM PUBLISHER]

This paper treats the problem of estimating the Mahalanobis distance for the purpose of detecting outliers in high-dimensional data. Three ridge-type estimators are proposed and risk functions for deciding an appropriate value of the ridge coefficient are developed. It is argued that one of the ridge estimator has particularly tractable properties, which is demonstrated through outlier analysis of real and simulated data. [ABSTRACT FROM AUTHOR]

This paper describes a statistical method for estimating data envelopment analysis (DEA) score confidence intervals for individual organizations or other entities. This method applies statistical panel data analysis, which provides proven and powerful methodologies for diagnostic testing and for estimation of confidence intervals. DEA scores are tested for violations of the standard statistical assumptions including contemporaneous correlation, serial correlation, heteroskedasticity and the absence of a normal distribution. Generalized least squares statistical models are used to adjust for violations that are present and to estimate valid confidence intervals within which the true efficiency of each individual decision-making unit occurs. This method is illustrated with two sets of panel data, one from large US urban transit systems and the other from a group of US hospital pharmacies. [ABSTRACT FROM AUTHOR]

Parameter estimation is the first step in constructing control charts. One of these parameters is the process mean. The classical estimators of the process mean are sensitive to the presence of outlying data and subgroups which contaminate the whole data. In existing robust estimators for the process mean, the effects of the presence of the individual outliers are being considered, while, in this paper, a robust estimator is being proposed to reduce the effect of outlying subgroups as well as the individual outliers within a subgroup. The proposed estimator was compared with some classical and robust estimators of the process mean. Although, its relative efficiency is fourth among the estimators tested, its robustness and efficiency are large when the outlying subgroups are present. Evaluation of the results indicated that the proposed estimator is less sensitive to the presence of outliers and the process mean performs well when there are no individual outliers or outlying subgroups. [ABSTRACT FROM AUTHOR]

Measurement error (ME) is the difference between the true unknown value of a variable and the data assigned to that variable during the measuring process. The multiple correlation coefficient quantifies the strength of the relationship between the dependent and independent variable(s) in regression modeling. In this paper, we show that ME in the dependent variable results in a negative bias in the multiple correlation coefficient, making the relationship appear weaker than it should. The adjusted R2 provides regression modelers an unbiased estimate of the multiple correlation coefficient. However, due to the ME induced bias in the multiple correlation coefficient, the otherwise unbiased adjusted R2 under-estimates the variance explained by a regression model. This paper proposes two statistics for estimating the multiple correlation coefficient, both of which take into account the ME in the dependent variable. The first statistic uses all unbiased estimators, but may produce values outside the [0,1] interval. The second statistic requires modeling a single data set, created by including descriptive variables on the subjects used in a gage study. Based on sums of squares, the statistic has the properties of an R2: it measures the proportion of variance explained; has values restricted to the [0,1] interval; and the endpoints indicate no variance explained and all variance explained respectively. We demonstrate the methodology using data from a study of cervical spine range of motion in children. [ABSTRACT FROM AUTHOR]

The first two stages in modelling times series are hypothesis testing and estimation. For long memory time series, the second stage was studied in the paper published in [M. Boutahar et al., Estimation methods of the long memory parameter: monte Carlo analysis and application, J. Appl. Statist. 34(3), pp. 261-301.] in which we have presented some estimation methods of the long memory parameter. The present paper is intended for the first stage, and hence completes the former, by exploring some tests for detecting long memory in time series. We consider two kinds of tests: the non-parametric class and the semi-parametric one. We precise the limiting distribution of the non-parametric tests under the null of short memory and we show that they are consistent against the alternative of long memory. We perform also some Monte Carlo simulations to analyse the size distortion and the power of all proposed tests. We conclude that for large sample size, the two classes are equivalent but for small sample size the non-parametric class is better than the semi-parametric one. [ABSTRACT FROM AUTHOR]

The generalized standard two-sided power (GTSP) distribution was mentioned only in passing by Kotz and van Dorp Beyond Beta, Other Continuous Families of Distributions with Bounded Support and Applications, World Scientific Press, Singapore, 2004. In this paper, we shall further investigate this three-parameter distribution by presenting some novel properties and use its more general form to contrast the chronology of developments of various authors on the two-parameter TSP distribution since its initial introduction. GTSP distributions allow for J-shaped forms of its pdf, whereas TSP distributions are limited to U-shaped and unimodal forms. Hence, GTSP distributions possess the same three distributional shapes as the classical beta distributions. A novel method and algorithm for the indirect elicitation of the two-power parameters of the GTSP distribution is developed. We present a Project Evaluation Review Technique example that utilizes this algorithm and demonstrates the benefit of separate powers for the two branches of activity GTSP distributions for project completion time uncertainty estimation. [ABSTRACT FROM AUTHOR]

The area of marked-point processes is well developed but simulation is still a challenging problem when mark correlations are to be included. In this paper we propose the use of simulated annealing to incorporate the spatial mark correlation into the simulations of correlated marked-point processes. Such a simulation has wide applications in areas such as inference and goodness-of-fit investigations of proposed models. The technique is applied to a forest dataset for which the results are extremely encouraging. [ABSTRACT FROM AUTHOR]

In this paper we use a penalized likelihood approach to image warping in the context of discrimination and averaging. The choice of average image is formulated statistically by minimizing a penalized likelihood, where the likelihood measures the similarity between images after warping and the penalty is a measure of distortion of a warping. The notions of measures of similarity are given in terms of normalized image information. The measures of distortion are landmark based. Thus we use a combination of landmark and normalized image information. The average defined in the paper is also extended by allowing random perturbation of the landmarks. This strategy improves averages for discrimination purposes. We give here real applications from medical and biological areas. [ABSTRACT FROM AUTHOR]

Consider the regression model . In a variety of situations, an estimate of VAR ( Y | X )?=?? ( X ) is needed. The paper compares the small-sample accuracy of five estimators of ? ( X ). The results suggest that the optimal estimator is a somewhat complex function of the underlying distributions. In terms of mean squared error, one of the estimators, which is based in part on a non-robust version of Cleveland's smoother, performed about as well as a bagged version of the so-called running interval smoother, but the running interval smoother was found to be preferable in terms of bias. A modification of Cleveland's smoother, stemming from Ruppert et al. (1997), achieves its intended goal of reducing bias when the error term is homoscedastic, but under heteroscedasticity, bias can be high, and in terms of mean squared error it does not compete well with the kernel method considered in the paper. When ? has a heavy-tailed distribution, a robust version of Cleveland's smoother performed particularly well except in some situations where X has a heavy-tailed distribution as well. A negative feature of using Cleveland's robust smoother is relatively high bias, and when there is heteroscedasticity and X has a heavy-tailed distribution, a kernel-type method and the running interval smoother give superior results in terms of both mean squared error and bias. [ABSTRACT FROM AUTHOR]

The propensity for military service (PMS) of young Americans is an important issue for our Armed Forces. Since the 1990s, the PMS of young Americans has steadily declined. Overtime, a declining PMS may cause military mission degradation, lowering of military recruitment standards, base closures, and reinstatement of the unpopular military draft system. This paper investigates the moderator effect of prior military service on the Generalized Exchange–PMS relationship. Generalized exchange is when indirect benefits such as preserving freedom and the American way of life accrue to the larger society because of an individual's military service. This paper uses a structural equation modelling approach to analyse the moderating effect of prior military exposure on prospective recruits regarding their PMS. Findings indicate that the group of prospective recruits with prior military exposure had higher levels of PMS than the group without such exposure, that is, the young people with prior military exposure are more likely to enlist in the military than the young Americans with no prior military exposure. [ABSTRACT FROM AUTHOR]

Winters are a difficult period for the National Health Service (NHS) in the United Kingdom (UK), due to the combination of cold weather and the increased likelihood of respiratory infections, especially influenza. In this article we present a proper statistical time series approach for modelling and analysing weekly hospital admissions in the West Midlands in the UK during the period week 15/1990 to week 14/1999. We consider three variables, namely, hospital admissions, general practitioner consultants, and minimum temperature. The autocorrelations of each series are shown to decay hyperbolically. The correlations of hospital admission and the lag of other series also decay hyperbolically but with different speed and directions. One of the main objectives of this paper is to show that each of the three series can be represented by a Fractional Differenced Autoregressive integrated moving average model, (FDA). Further, the hospital admission winter and summer residuals shows significant interdependency, which may be interpreted as hidden periodicities within the last 10-years time interval. The short-range (8 weeks) forecasting of hospital admission of the FDA model and a fourth-order AutoRegressive AR(4) model are quite similar. However, our results reveal that the long-range forecasting of FDA is more realistic. This implies that, using the FDA approach, the respective authority can plan for winter pressure properly. [ABSTRACT FROM AUTHOR]