Estimation of parameters in logistic and log-logistic distribution with grouped data

In many situations, instead of a complete sample, data are available only in grouped form. For example, grouped failure time data occur in studies in which subjects are monitored periodically to determine whether failure has occurred in the predetermined intervals. Here the model under consideration is the log-logistic distribution. This paper demonstrates the existence and uniqueness of the MLEs of the parameters of the logistic distribution under mild conditions with grouped data. The times with the maximum failure rate and the mode of the p.d.f. of the log-logistic distribution are also estimated based on the MLEs. The methodology is further studied with simulations and exemplified with a data set with artificially introduced grouping from a locomotive life test study.     

Orthogonal projections and the geometry of estimating functions

In this paper, a notion of generalized inner product spaces is introduced to study optimal estimating functions. The basic technique involves an idea of orthogonal projection first introduced by Small and McLeish (1988, 1989, 1991, 1992, 1994). A characterization of orthogonal projections in generalized inner product spaces is given. It is shown that the orthogonal projection of the score function into a linear subspace of estimating functions is optimal in that subspace, and a characterization of optimal estimating functions is given. As special cases of the main results of this paper, we derive the results of Godambe (1985) on the foundation of estimation in stochastic processes, the result of Godambe and Thompson (1989) on the extension of quasi-likelihood, and the generalized estimating equations for multivariate data due to Liang and Zeger (1986). Also we have derived optimal estimating functions in the Bayesian framework.     

Estimating moving average parameters in the presence of measurement error

Suppose that a moving average time series X_t is not observed, but instead Y_t = X_t + ?_t is observed, where ?_t, is measurement error. Estimation of the parameters of X_t has previously been considered under the assumption that X_t and ?_t are uncorrelated. The case where X_t and ?_t have known cross covariances is considered here, and a method is described for estimating the parameters of X_t. A simulation compares four estimators for a MA(1) series parameter in the presence of measurement error.     

Test of equality of normal means in the absence of independent estimator of variance

Let X_l,…,X_n be normally and independently distributed with means θ_l,…,θ_nand a cornmorl variance. Thus there are n observations and n+i unknwon parameters. A test of the null hypothesis that, the θ_i's are all zero and the alternative that the vector (θ_l,…,θ_n) lies in a convex cone with its vertex a.t the origin is connsidered in this paper. It is shown that under a mild condition the likelihood ratio test is possible. The ordinary one sided t - test belongs to the class of tests considered in this paper. The hypothesis of equality of means against the simple order alternative can be tested in certain cases .     

Parametric and nonparametric confidence intervals for estimating the difference of means of two skewed populations

In this paper, we have reviewed and proposed several interval estimators for estimating the difference of means of two skewed populations. Estimators include the ordinary-t, two versions proposed by Welch [17] and Satterthwaite [15], three versions proposed by Zhou and Dinh [18], Johnson [9], Hall [8], empirical likelihood (EL), bootstrap version of EL, median t proposed by Baklizi and Kibria [2] and bootstrap version of median t. A Monte Carlo simulation study has been conducted to compare the performance of the proposed interval estimators. Some real life health related data have been considered to illustrate the application of the paper. Based on our findings, some possible good interval estimators for estimating the mean difference of two populations have been recommended for the researchers.     

Self-updating clustering algorithm for estimating the parameters in mixtures of von Mises distributions

The EM algorithm is the standard method for estimating the parameters in finite mixture models. Yang and Pan [25] proposed a generalized classification maximum likelihood procedure, called the fuzzy c-directions (FCD) clustering algorithm, for estimating the parameters in mixtures of von Mises distributions. Two main drawbacks of the EM algorithm are its slow convergence and the dependence of the solution on the initial value used. The choice of initial values is of great importance in the algorithm-based literature as it can heavily influence the speed of convergence of the algorithm and its ability to locate the global maximum. On the other hand, the algorithmic frameworks of EM and FCD are closely related. Therefore, the drawbacks of FCD are the same as those of the EM algorithm. To resolve these problems, this paper proposes another clustering algorithm, which can self-organize local optimal cluster numbers without using cluster validity functions. These numerical results clearly indicate that the proposed algorithm is superior in performance of EM and FCD algorithms. Finally, we apply the proposed algorithm to two real data sets.     

An efficient and fast algorithm for estimating the parameters of two-dimensional sinusoidal signals

In this paper we propose a computationally efficient algorithm to estimate the parameters of a 2-D sinusoidal model in the presence of stationary noise. The estimators obtained by the proposed algorithm are consistent and asymptotically equivalent to the least squares estimators. Monte Carlo simulations are performed for different sample sizes and it is observed that the performances of the proposed method are quite satisfactory and they are equivalent to the least squares estimators. The main advantage of the proposed method is that the estimators can be obtained using only finite number of iterations. In fact it is shown that starting from the average of periodogram estimators, the proposed algorithm converges in three steps only. One synthesized texture data and one original texture data have been analyzed using the proposed algorithm for illustrative purpose.     

An emplrical bayes approach in estimating the parameters of the or more geometric distributions

The estimation of the parameters of two or more geometric distribuionsis considered by usinq an empirical Bayesian approach. Robbins (1983) gave empirical Bayes estimates if the number of distributions N is large, buthere we consider the cascwhore N is small. The parameters of the prior distribution areest imated by looking at maximum like lihood and momentest imation methods.     

Decomposing the predictive performance of the moving average trading rule of technical analysis: the contribution of linear and non-linear dependencies in stock returns

The main purpose of this work is to decompose the predictive performance of the moving average (MA) trading rule and find out the portion that could be attributed to the possible exploitation of linear and non-linear dependencies in stock returns. Data from the General Index of the Athens Stock Exchange, from the Standard and Poor-500 Index of the New York Stock Exchange and from the Austrian Traded Index of the Vienna Stock Exchange are filtered by linear filters so as the resulting simulated ‘returns’ exhibit no serial correlation. Applying MA trading rules to both the original and the simulated indices and using a new statistical testing procedure that takes into account the sensitivity of the performance of the trading rule as a function of the length of the MA it is found that the predictive performance of the trading rule is clearly weakened when applied to the simulated indices indicating that a substantial part of the rule's predictive performance is due to the exploitation of linear dependencies in stock returns. This weakening is uneven; in general the shorter the MA length the more pronounced the attenuation.     

Invariant tests based on M-estimators,estimating functions,and the generalized method of moments

We study the invariance properties of various test criteria which have been proposed for hypothesis testing in the context of incompletely specified models, such as models which are formulated in terms of estimating functions (Godambe, 1960) or moment conditions and are estimated by generalized method of moments (GMM) procedures (Hansen, 1982), and models estimated by pseudo-likelihood (Gouriéroux, Monfort, and Trognon, 1984b,c) and M-estimation methods. The invariance properties considered include invariance to (possibly nonlinear) hypothesis reformulations and reparameterizations. The test statistics examined include Wald-type, LR-type, LM-type, score-type, and C(α)?type criteria. Extending the approach used in Dagenais and Dufour (1991), we show first that all these test statistics except the Wald-type ones are invariant to equivalent hypothesis reformulations (under usual regularity conditions), but all five of them are not generally invariant to model reparameterizations, including measurement unit changes in nonlinear models. In other words, testing two equivalent hypotheses in the context of equivalent models may lead to completely different inferences. For example, this may occur after an apparently innocuous rescaling of some model variables. Then, in view of avoiding such undesirable properties, we study restrictions that can be imposed on the objective functions used for pseudo-likelihood (or M-estimation) as well as the structure of the test criteria used with estimating functions and generalized method of moments (GMM) procedures to obtain invariant tests. In particular, we show that using linear exponential pseudo-likelihood functions allows one to obtain invariant score-type and C(α)?type test criteria, while in the context of estimating function (or GMM) procedures it is possible to modify a LR-type statistic proposed by Newey and West (1987) to obtain a test statistic that is invariant to general reparameterizations. The invariance associated with linear exponential pseudo-likelihood functions is interpreted as a strong argument for using such pseudo-likelihood functions in empirical work.     

Comparison of adjusted and unadjusted marginal probabilities for estimating binomial parameters in observational studies

In this study, we demonstrate how generalized propensity score estimators (Imbens’ weighted estimator, the propensity score weighted estimator and the generalized doubly robust estimator) can be used to calculate the adjusted marginal probabilities for estimating the three common binomial parameters: the risk difference (RD), the relative risk (RR), and the odds ratio (OR). We further conduct a simulation study to compare the estimated RD, RR, and OR using the adjusted and the unadjusted marginal probabilities in terms of the bias and mean-squared error (MSE). Although there is no clear winner in terms of the MSE for estimating RD, RR, and OR, simulation results surprisingly show thatthe unadjusted marginal probabilities produce the smallest bias compared with the adjusted marginal probabilities in most of the estimates. Hence, in conclusion, we recommend using the unadjusted marginal probabilities to estimate RD, RR, and OR, in practice.     

A comparison of various methods for estimating the parameters in mixtures of von mises distributions

In this paper we compare five methods for estimating the unknown parameters in a mixture of two von Mises distributions. We propose a new method based on the characteristic function and compare it with the classical methods based on maximum likelihood and moments. Thus far these methods have been successfully applied only to linear data. Here we show that the application to circular data is reasonably straightforward and that convergence to the final estimates is fairly rapid. For various simulated known mixtures the results obtained are satisfactory. Finally, we introduce a modification of the method of moments which is considerably faster in CPU time than any of the other methods used and gives good results.     

贫困测量中恩格尔系数的失效及分析

 恩格尔系数在我国的适用性问题是学术界重要议题之一。本文利用对哈尔滨市和沈阳市贫困人口的调查数据,分析了恩格尔系数在我国城市贫困测量中失效的问题及原因。研究者尝试对恩格尔系数的分子、分母进行了重新界定和适当的修正,研究检验表明,从宏观测量和微观测量两个方面,修正恩格尔系数都能更准确的反映我国城市贫困家庭的真实生活水平。     

Properties of the elasticity of a continuous random variable. A special look at its behavior and speed of change

Belzunce et al. (1995 Belzunce, F., Candel, J., Ruiz, J.M. (1995). Ordering of truncated distributions through concentration curves. Sankhya 57:375–383. [Google Scholar]) define the elasticity for non negative random variables as the reversed proportional failure rate (RPFR). Veres-Ferrer and Pavía (2012 Veres-Ferrer, E.J., Pavía, J.M. (2012). La elasticidad: una nueva herramienta para caracterizar distribuciones de probabilidad. Rect@ 13:145–158. [Google Scholar], 2014b Veres-Ferrer, E.J., Pavía, J.M. (2014b). On the relationship between the reversed hazard rate and elasticity. Stat. Pap. 55:275–284.[Crossref], [Web of Science ®] , [Google Scholar]) interpret it in economic terms, extending its definition to variables that can also take negative values, and briefly present the role of elasticity in characterizing probability distributions. This paper highlights a set of properties demonstrated by elasticity, which shows many similar properties to the reverse hazard function. This paper pays particular attention to studying the increase/decrease and the speed of change of the elasticity function. These are important properties because of the characterizing role of elasticity, which makes it possible to introduce our hypotheses and knowledge about the random process in a more meaningful and intuitive way. As a general rule, it is observed the need for distinguishing between positive and negative areas of the support.     

The power functions of the likelihood ratio tests for a simply ordered trend in normal means

Likelihood ratio tests for the homogeneity of k normal means with the alternative restricted by an increasing trend are considered as well as the likelihood ratio tests of the null hypothesis that the means satisfy the trend. While the work is primarily a survey of results concerning the power functions of these tests, the extensions of some results to the case of not necessarily equal sample sizes are presented. For the case of known or unknown population variances, exact expressions are given for the power functions for k=3,4, and approximations are discussed for larger k. The topics of consistency, bias and monotonicity of the power functions are included. Also, Bartholomew's conjectures concerning minimal and maximal powers are investigated, with results of a new numerical study given.     

Subroutines for computing minimum variance unbiased estimators of functions of the parameters in the normal and gamma distributions

In this paper subroutines are given which calculate the uniformly minimum variance unbiased estimators (UMVUE’ s) of a broad class of functions of the parameters of the normal and gamma distributions. These subroutines employ the new expressions for the UMVUE’ s given recently by Gray, Watkins, and Schucany (1973), Woodward and Gray (1975), and Gray, Schucany, and Woodward (1976). In order to employ the subroutines here the user need only be able to provide a FORTRAN function subprogram to calculate derivatives of the function, either analytically or numerically.     

A method for estimating parameters and quantiles of the three-parameter inverse Gaussian distribution based on statistics invariant to unknown location

The inverse Gaussian (IG) distribution, also known as the Wald distribution, is a long-tailed positively skewed distribution and a well-known lifetime distribution. In this paper, we propose an efficient method of estimation for the parameters and quantiles of the three-parameter IG distribution, which is based on statistics invariant to unknown location. Through a Monte Carlo simulation study, we then show that the proposed method performs well compared with other prominent methods in terms of bias and variance. Finally, we present two illustrative examples.