Estimable functions in the nonlinear model

When the method of least squares is used to estimate the parameters in a general model and the generated system of normal equations is linearly dependent, the estimate of the vector of parameters which satisfies the criterion is not unique. However, there exist certain functions of the estimated vector of parameters which are invariant to the least squares solution obtained from the normal equations. We define those invariant functions to be estimable, and present a technique to determine the functions of the parameters which are estimable for the general model. The method results in solving either a linear first order partial differential equation or a system of linear first order partial differential equations corresponding, respectively, to a single or multiple dependency between columns of the Jacobian matrix of the mean of the model. The usual results concerning estimability for linear models are a special case of the general results developed.     

Simulation of the p-generalized Gaussian distribution

In this paper, we introduce the p-generalized polar methods for the simulation of the p-generalized Gaussian distribution. On the basis of geometric measure representations, the well-known Box–Muller method and the Marsaglia–Bray rejecting polar method for the simulation of the Gaussian distribution are generalized to simulate the p-generalized Gaussian distribution, which fits much more flexibly to data than the Gaussian distribution and has already been applied in various fields of modern sciences. To prove the correctness of the p-generalized polar methods, we give stochastic representations, and to demonstrate their adequacy, we perform a comparison of six simulation techniques w.r.t. the goodness of fit and the complexity. The competing methods include adapted general methods and another special method. Furthermore, we prove stochastic representations for all the adapted methods.     

Income vulnerability of rural households in Bangladesh: a comparison between Bayesian and classical methods

The geographical location and the monsoon climate render Bangladesh highly vulnerable to natural hazards, deteriorating the country's socio-economic stability. This study is based on 500 randomly chosen rural households from the Household Income and Expenditure Survey [Bangladesh Bureau of Statistics, Planning Division, Ministry of Planning, Government of the People's Republic of Bangladesh, Dhaka, 2006]. The objectives are to estimate the income vulnerability of rural households and to check whether the Bayesian approaches (natural conjugate prior and non-informative prior estimates) have any superiority over the classical (feasible generalized least square (FGLS)) method. The poverty level, measured from the data, is 24%; whereas the vulnerability estimates, using FGLS, natural conjugate prior and non-informative prior are 31%, 69% and 82%, respectively. Vulnerability estimates by the Bayesian natural conjugate prior approach is found to have greater efficiency compared with FGLS and non-informative prior approaches.     

An EM algorithm for the estimation of parameters of bivariate generalized exponential distribution under random left censoring

In this paper, we consider the four-parameter bivariate generalized exponential distribution proposed by Kundu and Gupta [Bivariate generalized exponential distribution, J. Multivariate Anal. 100 (2009), pp. 581–593] and propose an expectation–maximization algorithm to find the maximum-likelihood estimators of the four parameters under random left censoring. A numerical experiment is carried out to discuss the properties of the estimators obtained iteratively.     

Approximate Bayes estimation of the parameters and reliability function of a mixture of two inverse Weibull distributions under Type-2 censoring

The main goal of this paper is to develop the approximate Bayes estimation of the five-dimensional vector of the parameters and reliability function of a mixture of two inverse Weibull distributions (MTIWD) under Type-2 censoring. Usually, the posterior distribution is complicated under the scheme of Type-2 censoring and the integrals that are involved cannot be obtained in a simple explicit form. In this study, we use Lindley's [Approximate Bayesian method, Trabajos Estadist. 31 (1980), pp. 223–237] approximate form of Bayes estimation in the case of an MTIWD under Type-2 censoring. Later, we calculate the estimated risks (ERs) of the Bayes estimates and compare them with the corresponding ERs of the maximum-likelihood estimates through Monte Carlo simulation. Finally, we analyse a real data set using the findings.     

Comparison of ANOVA-F and ANOM tests with regard to type I error rate and test power

A Monte Carlo simulation was conducted to compare the type I error rate and test power of the analysis of means (ANOM) test to the one-way analysis of variance F-test (ANOVA-F). Simulation results showed that as long as the homogeneity of the variance assumption was satisfied, regardless of the shape of the distribution, number of group and the combination of observations, both ANOVA-F and ANOM test have displayed similar type I error rates. However, both tests have been negatively affected from the heterogeneity of the variances. This case became more obvious when the variance ratios increased. The test power values of both tests changed with respect to the effect size (Δ), variance ratio and sample size combinations. As long as the variances are homogeneous, ANOVA-F and ANOM test have similar powers except unbalanced cases. Under unbalanced conditions, the ANOVA-F was observed to be powerful than the ANOM-test. On the other hand, an increase in total number of observations caused the power values of ANOVA-F and ANOM test approach to each other. The relations between effect size (Δ) and the variance ratios affected the test power, especially when the sample sizes are not equal. As ANOVA-F has become to be superior in some of the experimental conditions being considered, ANOM is superior in the others. However, generally, when the populations with large mean have larger variances as well, ANOM test has been seen to be superior. On the other hand, when the populations with large mean have small variances, generally, ANOVA-F has observed to be superior. The situation became clearer when the number of the groups is 4 or 5.     

Asymptotic properties of a rank estimate in linear regression with symmetric non-identically distributed errors

In this article, a simple linear regression model with independent and symmetric but non-identically distributed errors is considered. Asymptotic properties of the rank regression estimate defined in Jaeckel [Estimating regression coefficients by minimizing the dispersion of the residuals, Ann. Math. Statist. 43 (1972), pp. 1449–1458] are studied. We show that the studied estimator is consistent and asymptotically normally distributed. The cases of bounded and unbounded score functions are examined separately. The regularity conditions of the article are exemplified for finite mixture distributions.     

The umvue and bayes estimate of reliability of mixed failure time distribution

We study the reliability estimates of the non-standard mixture of degenerate (degenerated at zero) and exponential distributions. The Uniformly Minimum Variance Unbiased Estimator (UMVUE) and Bayes estimator of the reliability for some selective prior when the mixing proportion is known and unknown are derived. The Bayes risk is computed for each Bayes estimator of the reliability. A simulated study is carried out to assess the performance of the estimators alongwith the true and Maximum Likelihood Estimate (MLE) of the reliability. An example from Vannman (1991) is also discussed at the end of the paper.     

Comparison of Confidence Intervals on Between Group Variance in Unbalanced Heteroscedastic One-Way Random Models

A computer algorithm for computing the alternative distributions of the Wilcoxon signed rank statistic under shift alternatives is discussed. An explicit error bound is derived for the numeric integration approximation to these distributions.

A nonparametric process control procedure in which the standard CUSUM procedure is applied to the Wilcoxon signed rank statistic is discussed. In order to implement this procedure, the distribution of the Wilcoxon statistic under shift of the underlying distribution from its point of symmetry needs to be computed. The average run length of the nonparametric and parametric CUSUM are compared.