The Lamé class of Lorenz curves

In this paper, the class of Lamé Lorenz curves is studied. This family has the advantage of modeling inequality with a single parameter. The family has a double motivation: it can be obtained from an economic model and from simple transformations of classical Lorenz curves. The underlying cumulative distribution functions have a simple closed form, and correspond to the Singh–Maddala and Dagum distributions, which are well known in the economic literature. The Lorenz order is studied and several inequality and polarization measures are obtained, including Gini, Donaldson–Weymark–Kakwani, Pietra, and Wolfson indices. Some extensions of the Lamé family are obtained. Fitting and estimation methods under two different data configurations are proposed. Empirical applications with real data are given. Finally, some relationships with other curves are included.     

Bayesian Analysis of a Generalized Poisson Distribution

A generalized form of the Poisson Distribution with two parameters will be estimated by the Bayesian technique. When one of the parameters is known, several important parametric functions will be estimated and a numerical comparison with estimates obtained by the methods of maximum likelihood and unbiased minimum variance will be drawn. The simplicity of the posterior distribution of the unknown parameter enables us to construct exact probability intervals, and to devise a statistic to test the homogeneity of several populations. When the two parameters are unknown, dependent priors are being considered. Although the posterior distributions are sensitive to the choice of the prior, the posterior estimates are very stable and we use the Pearson system of curves to construct approximate posterior confidence limits for the parameters.     

Generating Ordered Families of Lorenz Curves by Strongly Unimodal Distributions

From any strongly unimodal density on the real line, it is possible to generate a one-parameter family of Lorenz curves. The resulting families of Lorenz curves are Lorenz ordered with respect to the indexing parameter. Symmetry of the unimodal density results in the generation of symmetric Lorenz curves. A related characterization of the normal distribution is presented.     

Influence functions of two families of robust estimators under proportional scatter matrices

In this paper, under a proportional model, two families of robust estimates for the proportionality constants, the common principal axes and their size are discussed. The first approach is obtained by plugging robust scatter matrices on the maximum likelihood equations for normal data. A projection- pursuit and a modified projection-pursuit approach, adapted to the proportional setting, are also considered. For all families of estimates, partial influence functions are obtained and asymptotic variances are derived from them. The performance of the estimates is compared through a Monte Carlo study.     

Statistical Comparison of the Goodness of Fit Delivered by Five Families of Distributions Used in Distribution Fitting

A basic assumption in distribution fitting is that a single family of distributions may deliver useful representation to the universe of available distributions. To date, little study has been conducted to compare the relative effectiveness of these families. In this article, five families are compared by fitting them to a sample of 20 distributions, using 2 fitting objectives: minimization of the L ₂ norm and four-moment matching. Values of L ₂ norm associated with the fitted families are used as input data to test for significant differences. The Pearson family and the RMM (Response Modeling Methodology) family significantly outperforms all other families.     

Ratio and product estimators in stratified random sampling

Khoshnevisan et al. [2007. A general family of estimators for estimating population mean using known value of some population parameter(s). Far East Journal of Theoretical Statistics 22, 181–191] have introduced a family of estimators using auxiliary information in simple random sampling. They have showed that these estimators are more efficient than the classical ratio estimator and that the minimum value of the mean square error (MSE) of this family is equal to the value of MSE of regression estimator. In this article, we adapt the estimators in this family to the stratified random sampling and motivated by the estimator in Searls [1964. Utilization of known coefficient of kurtosis in the estimation procedure of variance. Journal of the American Statistical Association 59, 1225–1226], we also propose a new family of estimators for the stratified random sampling. The expressions of bias and MSE of the adapted and proposed families are derived in a general form. Besides, considering the minimum cases of these MSE equations, the efficient conditions between the adapted and proposed families are obtained. Moreover, these theoretical findings are supported by a numerical example with original data.     

On Minimum Variance Unbiased Estimators

The problem of finding minimum variance unbiased estimators of various parameters for parametric distributions is an important one in statistics. This article gives analytical formulas for the minimum variance unbiased estimators of parametric functions, which are usually used in a classroom, for two types of densities. The first type is the one-parameter regular exponential family, and the second is a two-parameter family of a continuous random variable whose range depends on the unknown parameters.     

Functional linear discriminant analysis for irregularly sampled curves

We introduce a technique for extending the classical method of linear discriminant analysis (LDA) to data sets where the predictor variables are curves or functions. This procedure, which we call functional linear discriminant analysis ( FLDA ), is particularly useful when only fragments of the curves are observed. All the techniques associated with LDA can be extended for use with FLDA. In particular FLDA can be used to produce classifications on new (test) curves, give an estimate of the discriminant function between classes and provide a one- or two-dimensional pictorial representation of a set of curves. We also extend this procedure to provide generalizations of quadratic and regularized discriminant analysis.     

Local principal curves

Principal components are a well established tool in dimension reduction. The extension to principal curves allows for general smooth curves which pass through the middle of a multidimensional data cloud. In this paper local principal curves are introduced, which are based on the localization of principal component analysis. The proposed algorithm is able to identify closed curves as well as multiple curves which may or may not be connected. For the evaluation of the performance of principal curves as tool for data reduction a measure of coverage is suggested. By use of simulated and real data sets the approach is compared to various alternative concepts of principal curves.     

基于B-样条基底展开的曲线聚类方法

随着大数据时代的来临,近年来函数型数据分析方法成为研究的热点问题,针对曲线的聚类分析方法引起了学界的关注.给出一种曲线聚类的方法:以L2距离作为亲疏程度的度量,在B样条基底函数展开表述下,将曲线本身信息、曲线变化信息引入聚类算法构建,并实现了曲线聚类与传统多元统计聚类方法的对接.作为应用,以城乡收入函数聚类实例验证了该曲线聚类方法,结果表明,在引入曲线变化信息的情况下,比仅考虑曲线本身信息能够取得更好的聚类效果.     

Comparison of Separable Components in Different Samples

Abstract.  Imagine we have two different samples and are interested in doing semi- or non-parametric regression analysis in each of them, possibly on the same model. In this paper, we consider the problem of testing whether a specific covariate has different impacts on the regression curve in these two samples. We compare the regression curves of different samples but are interested in specific differences instead of testing for equality of the whole regression function. Our procedure does allow for random designs, different sample sizes, different variance functions, different sets of regressors with different impact functions, etc. As we use the marginal integration approach, this method can be applied to any strong, weak or latent separable model as well as to additive interaction models to compare the lower dimensional separable components between the different samples. Thus, in the case of having separable models, our procedure includes the possibility of comparing the whole regression curves, thereby avoiding the curse of dimensionality. It is shown that bootstrap fails in theory and practice. Therefore, we propose a subsampling procedure with automatic choice of subsample size. We present a complete asymptotic theory and an extensive simulation study.     

An alternative method for global and partial comparison of two diagnostic systems based on ROC curves

In this paper, an alternative method for the comparison of two diagnostic systems based on receiver operating characteristic (ROC) curves is presented. ROC curve analysis is often used as a statistical tool for the evaluation of diagnostic systems. However, in general, the comparison of ROC curves is not straightforward, in particular, when they cross each other. A similar difficulty is also observed in the multi-objective optimization field where sets of solutions defining fronts must be compared with a multi-dimensional space. Thus, the proposed methodology is based on a procedure used to compare the performance of distinct multi-objective optimization algorithms. In general, methods based on the area under the ROC curves are not sensitive to the existence of crossing points between the curves. The new approach can deal with this situation and also allows the comparison of partial portions of ROC curves according to particular values of sensitivity and specificity of practical interest. Simulations results are presented. For illustration purposes, considering real data from newborns with very low birthweight, the new method was applied in order to discriminate the better index for evaluating the risk of death.     

The Exponentiated G Poisson Model

We study a new family of distributions defined by the minimum of the Poisson random number of independent identically distributed random variables having a general exponentiated G distribution. Some mathematical properties of the new family including ordinary and incomplete moments, quantile and generating functions, mean deviations, order statistics and their moments, reliability, and Shannon entropy are derived. Maximum likelihood estimation of the model parameters is investigated. Two special models of the new family are discussed. We perform an application to a real data set to show the potentiality of the proposed family.     

Instability of Regression Coefficients Illustrated

In the two-factor factorial fixed effects ANOVA, no interaction between factors may be demonstrated graphically by connecting population cell means to form a set of parallel piecewise linear curves. The three-factor case is considerably more complicated, because both three-factor and two-factor interactions may be present. However, if no three-factor interaction is assumed, or, in addition, if one or more of the two-factor interactions is also assumed to be zero, then the resulting constrained model may be characterized in a manner similar to the two-factor case. Details are presented for three such constrained models.     

Local Power Analyses of Goodness-of-fit Tests for Copulas

Abstract.  The asymptotic behaviour of several goodness-of-fit statistics for copula families is obtained under contiguous alternatives. Many comparisons between a Cramér–von Mises functional of the empirical copula process and new moment-based goodness-of-fit statistics are made by considering their associated asymptotic local power curves. It is shown that the choice of the estimator for the unknown parameter can have a significant influence on the power of the Cramér–von Mises test and that some of the moment-based statistics can provide simple and efficient goodness-of-fit methods.     

Bayesian spatio‐temporal models based on discrete convolutions

Abstract: The authors consider a class of models for spatio‐temporal processes based on convolving independent processes with a discrete kernel that is represented by a lower triangular matrix. They study two families of models. In the first one, spatial Gaussian processes with isotropic correlations are convoluted with a kernel that provides temporal dependencies. In the second family, AR(p) processes are convoluted with a kernel providing spatial interactions. The covariance structures associated with these two families are quite rich. Their covariance functions that are stationary and separable in space and time as well as time dependent nonseparable and nonisotropic ones.     

A Small Sample Test for the Equality of Intraclass Correlation Coefficients Under Unequal Family Sizes for Several Populations

The set of all distinct blocks of a BIBD(v,b,r,k,λ) is referred to as the support of the design. In this paper, the family of BIB designs with v=9 and k=3 is studied from the view of possible support sizes, b*'s. A table is constructed of designs with support sizes belonging to {12,18,20,21,…,84}, for minimum possible b in each case and for any larger admissible b. In constructing this table the methods of trade-off and composition of designs are utilized     

An identity for the joint distribution of order statistics and its applications

It is shown that the joint distribution function of several order statistics can be expressed in terms of the distribution functions of the maximum (or the minimum) order statistics of suitable subsamples. The result is used to derive explicit expressions for the expectations of functions of order statistics from certain families of distributions which include the exponential distribution and the power function distribution. The results generalize earlier work by the authors for a single order statistic.     

Minimum variance unbiased estimation in natural exponential families generated by stable distributions

The class of nature exponential families generated by stable distributions has been introduced in different contexts by several authors. Tweedie (1984) and Jorgensen (1987) studied this class in the context of generalized liner models and exponential dispersion models. Bar-Lev and Enis (1986) introduced this class in the context of the property of reproducibility in natural exponential families and Hougaard (1986) found the distributions in this class to be natural candidates for applications as survival distributions in life tables for heterogeneous populations. In this note, we consider such a class in the context of minimum variance unbiased estimation. For each family in this class, we obtain an explicit expression for the uniformly minimum variance unbiased estimator for the r-th cumlant, the density function, and the reliability function.     

MINIMUM QUADRATIC DISTANCE ESTIMATION FOR A PARAMETRIC FAMILY OF DISCRETE DISTRIBUTIONS DEFINED RECURSIVELY

The distribution function of a random sum can easily be computed iteratively when the distribution of the number of independent identically distributed elements in the sum is itself defined recursively. Classical estimation procedures for such recursive parametric families often require specific distributional assumptions (e.g. Poisson, Negative Binomial). The minimum distance estimator proposed here is an estimator within a larger parametric family. The estimator is consistent, efficient when the parametric family is truncated, and can be made either robust or asymptotically efficient when the parametric family has infinite range. Its asymptotic distribution is derived. A brief illustration with Automobile Insurance data is included.