Bayesian inference for circular distributions with unknown normalising constants

Very often, the likelihoods for circular data sets are of quite complicated forms, and the functional forms of the normalising constants, which depend upon the unknown parameters, are unknown. This latter problem generally precludes rigorous, exact inference (both classical and Bayesian) for circular data.Noting the paucity of literature on Bayesian circular data analysis, and also because realistic data analysis is naturally permitted by the Bayesian paradigm, we address the above problem taking a Bayesian perspective. In particular, we propose a methodology that combines importance sampling and Markov chain Monte Carlo (MCMC) in a very effective manner to sample from the posterior distribution of the parameters, given the circular data. With simulation study and real data analysis, we demonstrate the considerable reliability and flexibility of our proposed methodology in analysing circular data.     

Bayesian inference for a flexible class of bivariate beta distributions

Several bivariate beta distributions have been proposed in the literature. In particular, Olkin and Liu [A bivariate beta distribution. Statist Probab Lett. 2003;62(4):407–412] proposed a 3 parameter bivariate beta model which Arnold and Ng [Flexible bivariate beta distributions. J Multivariate Anal. 2011;102(8):1194–1202] extend to 5 and 8 parameter models. The 3 parameter model allows for only positive correlation, while the latter models can accommodate both positive and negative correlation. However, these come at the expense of a density that is mathematically intractable. The focus of this research is on Bayesian estimation for the 5 and 8 parameter models. Since the likelihood does not exist in closed form, we apply approximate Bayesian computation, a likelihood free approach. Simulation studies have been carried out for the 5 and 8 parameter cases under various priors and tolerance levels. We apply the 5 parameter model to a real data set by allowing the model to serve as a prior to correlated proportions of a bivariate beta binomial model. Results and comparisons are then discussed.     

A simple approach to the estimation of Tukey's gh distribution

ABSTRACT

The Tukey's gh distribution is widely used in situations where skewness and elongation are important features of the data. As the distribution is defined through a quantile transformation of the normal, the likelihood function cannot be written in closed form and exact maximum likelihood estimation is unfeasible. In this paper we exploit a novel approach based on a frequentist reinterpretation of Approximate Bayesian Computation for approximating the maximum likelihood estimates of the gh distribution. This method is appealing because it only requires the ability to sample the distribution. We discuss the choice of the input parameters by means of simulation experiments and provide evidence of superior performance in terms of Root-Mean-Square-Error with respect to the standard quantile estimator. Finally, we give an application to operational risk measurement.     

On the test of liu and chow's procedure for assessing equivalence in variability of bioavailability

Similar to Schuirmann's two one-sided tests procedure for assessment of bioequivalence in average bioavailability (Schuirmann,), Liu and Chow proposed a two one-sided tests procedure for assessment of equivalence of variability of bioavailability. Their procedure is derived based on the correlation between crossover differences and subject totals. In this paper, we examined the performance of their test procedure in terms of its test size and power for various situations where the intersubject variability and the intrasubject variability of the test drug product are relatively larger, similar, and smaller than that of the intrasubject variability of the reference drug product.     

A new policy for designing acceptance sampling plan based on Bayesian inference in the presence of inspection errors

The purpose of this article is to present a new policy for designing an acceptance sampling plan based on the minimum proportion of the lot that should be inspected in the presence of inspection errors. It is assumed that inspection is not perfect and every defective item cannot be detected with complete certainty. The Bayesian method is used for obtaining the probability distribution function of the number of defective items in the lot. To design this model, two constraints of producer risk and consumer risk are considered during the inspection process by using two specified points on operating characteristic curve. In order to illustrate the application of the proposed model, an example is presented. In addition, a sensitivity analysis is performed to analyze the model performance under different scenarios of process parameters and the results are elaborated. Finally, the efficiency of the proposed model is compared with the sampling method of Spencer and Kevan de Lopez (2017) at the same conditions.     

On the robustness of empirical likelihood ratio confidence intervals for location

The authors examine the robustness of empirical likelihood ratio (ELR) confidence intervals for the mean and M‐estimate of location. They show that the ELR interval for the mean has an asymptotic breakdown point of zero. They also give a formula for computing the breakdown point of the ELR interval for M‐estimate. Through a numerical study, they further examine the relative advantages of the ELR interval to the commonly used confidence intervals based on the asymptotic distribution of the M‐estimate.     

On estimating parameters of a progressively censored lognormal distribution

We consider the problem of making statistical inference on unknown parameters of a lognormal distribution under the assumption that samples are progressively censored. The maximum likelihood estimates (MLEs) are obtained by using the expectation-maximization algorithm. The observed and expected Fisher information matrices are provided as well. Approximate MLEs of unknown parameters are also obtained. Bayes and generalized estimates are derived under squared error loss function. We compute these estimates using Lindley's method as well as importance sampling method. Highest posterior density interval and asymptotic interval estimates are constructed for unknown parameters. A simulation study is conducted to compare proposed estimates. Further, a data set is analysed for illustrative purposes. Finally, optimal progressive censoring plans are discussed under different optimality criteria and results are presented.     

On the exact distribution of the likelihood ratio criterion for testing

The exact null distribution of the likelihood ratio criter- 2 ion for testing the hypothesis H: y = y~; z = a I, a unknown and UQ a given known vector against the alternative A =f H in a p-vari- ate normal population N (y,z) has been derived in the form of Meijer's G-function using mellin integral transform and also in a chisquare series form. Asymptotic behavior of the distribution of -2 log L has also been discussed. Percentage points for p=2(l)10for various level of significance and various degrees of freedom have been computed, but only selected tables have been presented in this paper.     

Parameter estimation and prediction of order statistics for the Burr Type XII distribution with Type II censoring

This article deals with the statistical inference and prediction on Burr Type XII parameters based on Type II censored sample. It is observed that the maximum likelihood estimators (MLEs) cannot be obtained in closed form. We use the expectation-maximization algorithm to compute the MLEs. We also obtain the Bayes estimators under symmetric and asymmetric loss functions such as squared error and Linex By applying Lindley's approximation and Markov chain Monte Carlo (MCMC) technique. Further, MCMC samples are used to calculate the highest posterior density credible intervals. Monte Carlo simulation study and two real-life data-sets are presented to illustrate all of the methods developed here. Furthermore, we obtain a prediction of future order statistics based on the observed ordered because of its important application in different fields such as medical and engineering sciences. A numerical example carried out to illustrate the procedures obtained for prediction of future order statistics.