Bayesian inference on extreme value distribution using upper record values

In this paper we address estimation and prediction problems for extreme value distributions under the assumption that the only available data are the record values. We provide some properties and pivotal quantities, and derive unbiased estimators for the location and rate parameters based on these properties and pivotal quantities. In addition, we discuss mean-squared errors of the proposed estimators and exact confidence intervals for the rate parameter. In Bayesian inference, we develop objective Bayesian analysis by deriving non informative priors such as the Jeffrey, reference, and probability matching priors for the location and rate parameters. We examine the validity of the proposed methods through Monte Carlo simulations for various record values of size and present a real data set for illustration purposes.     

Bayesian analysis for the two-parameter Pareto distribution based on record values and times

Doostparast and Balakrishnan (Pareto record-based analysis, Statistics, under review) recently developed optimal confidence intervals as well as uniformly most powerful tests for one- and two-sided hypotheses concerning shape and scale parameters, for the two-parameter Pareto distribution based on record data. In this paper, on the basis of record values and inter-record times from the two-parameter Pareto distribution, maximum-likelihood and Bayes estimators as well as credible regions are developed for the two parameters of the Pareto distribution. For illustrative purposes, a data set on annual wages of a sample of production-line workers in a large industrial firm is analysed using the proposed procedures.     

Bayesian inference and prediction of the Rayleigh distribution based on ordered ranked set sampling

Based on ordered ranked set sample, Bayesian estimation of the model parameter as well as prediction of the unobserved data from Rayleigh distribution are studied. The Bayes estimates of the parameter involved are obtained using both squared error and asymmetric loss functions. The Bayesian prediction approach is considered for predicting the unobserved lifetimes based on a two-sample prediction problem. A real life dataset and simulation study are used to illustrate our procedures.     

Statistical inference about the shape parameter of the Weibull distribution by upper record values

Summary In this paper, we provide some pivotal quantities to test and establish confidence interval of the shape parameter on the basis of the firstn observed upper record values. Finally, we give some examples and the Monte Carlo simulation to assess the behaviors (including higher power and more shorter length of confidence interval) of these pivotal quantities for testing null hypotheses and establishing confidence interval concerning the shape parameter under the given significance level and the given confidence coefficient, respectively.     

The E-Bayesian and hierarchical Bayesian estimations for the proportional reversed hazard rate model based on record values

In this paper, E-Bayesian and hierarchical Bayesian estimations of the shape parameter, when the underlying distribution belongs to the proportional reversed hazard rate model, are considered. Maximum likelihood, Bayesian and E-Bayesian estimates of the unknown parameter and reliability function are obtained based on record values. The Bayesian estimates are derived based on squared error and linear–exponential loss functions. It is pointed out that some previously obtained order relations of E-Bayesian estimates are inadequate and these results are improved. The relationship between E-Bayesian and hierarchical Bayesian estimations is obtained under the same loss functions. The comparison of the derived estimates is carried out by using Monte Carlo simulations. A real data set is analysed for an illustration of the findings.     

Inference for the generalized Rayleigh distribution based on progressively censored data

In this paper, and based on a progressive type-II censored sample from the generalized Rayleigh (GR) distribution, we consider the problem of estimating the model parameters and predicting the unobserved removed data. Maximum likelihood and Bayesian approaches are used to estimate the scale and shape parameters. The Gibbs and Metropolis samplers are used to predict the life lengths of the removed units in multiple stages of the progressively censored sample. Artificial and real data analyses have been performed for illustrative purposes.     

Bayesian Prediction for Progressively Type-II Censored Data from the Rayleigh Model

ABSTRACT

The Rayleigh distribution is proposed to be the underlying model from which observables are to be predicted by using Bayesian approach. Progressively Type-II censored data from the Rayleigh distribution is considered and the two-sample prediction technique is used. Numerical computations and a simulation are given to illustrate the performance of the procedures.     

Bayesian estimation and prediction for some life distributions based on record values

Some statistical data are most easily accessed in terms of record values. Examples include meteorology, hydrology and athletic events. Also, there are a number of industrial situations where experimental outcomes are a sequence of record-breaking observations. In this paper, Bayesian estimation for the two parameters of some life distributions, including Exponential, Weibull, Pareto and Burr type XII, are obtained based on upper record values. Prediction, either point or interval, for future upper record values is also presented from a Bayesian view point. Some of the non-Bayesian results can be achieved as limiting cases from our results. Numerical computations are given to illustrate the results.     

Shortest tolerance intervals controlling both tails of the exponential distribution based on record values

This paper deals with computing shortest width tolerance intervals controlling both tails of the exponential distribution on the basis of record values. Equal-tailed and shortest tolerance factors are derived. The expected widths of these tolerance intervals are evaluated via a Monte Carlo simulation study. Finally, two illustrative examples are also included.     

Lower bounds on the expectations of upper record values

For infinite sequences of independent random variables with identical continuous distributions, we establish optimal lower bounds on the deviations of the expectations of record values from population means in units generated by the central absolute moments of various orders. The bounds are non-negative for the classic record values, and non-positive for the other kth records with k?2. We also provide analogous bounds for the record increments.     

Bayesian estimation and prediction based on lognormal record values

In this paper we consider the problems of estimation and prediction when observed data from a lognormal distribution are based on lower record values and lower record values with inter-record times. We compute maximum likelihood estimates and asymptotic confidence intervals for model parameters. We also obtain Bayes estimates and the highest posterior density (HPD) intervals using noninformative and informative priors under square error and LINEX loss functions. Furthermore, for the problem of Bayesian prediction under one-sample and two-sample framework, we obtain predictive estimates and the associated predictive equal-tail and HPD intervals. Finally for illustration purpose a real data set is analyzed and simulation study is conducted to compare the methods of estimation and prediction.     

Bayesian process monitoring schemes for the two-parameter exponential distribution

In this paper a Bayesian procedure is applied to obtain control limits for the location and scale parameters, as well as for a one-sided upper tolerance limit in the case of the two-parameter exponential distribution. An advantage of the upper tolerance limit is that it monitors the location and scale parameter at the same time. By using Jeffreys’ non-informative prior, the predictive distributions of future maximum likelihood estimators of the location and scale parameters are derived analytically. The predictive distributions are used to determine the distribution of the “run-length” and expected “run-length”. A dataset given in Krishnamoorthy and Mathew (2009 Krishnamoorthy, K., and T. Mathew. 2009. Statistical Tolerance Regions: Theory, Applications and Computation. Wiley Series in Probability and Statistics.[Crossref] , [Google Scholar]) are used for illustrative purposes. The data are the mileages for some military personnel carriers that failed in service. The paper illustrates the flexibility and unique features of the Bayesian simulation method.     

Estimation with the generalized exponential distribution based on record values and inter-record times

The maximum likelihood and Bayesian approaches have been considered for the two-parameter generalized exponential distribution based on record values with the number of trials following the record values (inter-record times). The maximum likelihood estimates are obtained under the inverse sampling and the random sampling schemes. It is shown that the maximum likelihood estimator of the shape parameter converges in mean square to the true value when the scale parameter is known. The Bayes estimates of the parameters have been developed by using Lindley's approximation and the Markov Chain Monte Carlo methods due to the lack of explicit forms under the squared error and the linear-exponential loss functions. The confidence intervals for the parameters are constructed based on asymptotic and Bayesian methods. The Bayes and the maximum likelihood estimators are compared in terms of the estimated risk by the Monte Carlo simulations. The comparison of the estimators based on the record values and the record values with their corresponding inter-record times are performed by using Monte Carlo simulations.     

Bayesian shrinkage estimation based on Rayleigh type-II censored data

In this paper, we construct a Bayes shrinkage estimator for the Rayleigh scale parameter based on censored data under the squared log error loss function. Risk-unbiased estimator is derived and its risk is computed. A Bayes shrinkage estimator is obtained when a prior point guess value is available for the scale parameter. Risk-bias of the Bayes shrinkage estimator is considered. A comparison between the proposed Bayes shrinkage estimator and the risk-unbiased estimator is provided using calculation of the relative efficiency. A numerical example is presented for illustrative and comparative purposes.     

Bayesian prediction of order statistics with fixed and random sample sizes based on k-record values from Pareto distribution

In this paper, the two-parameter Pareto distribution is considered and the problem of prediction of order statistics from a future sample and that of its geometric mean are discussed. The Bayesian approach is applied to construct predictors based on observed k-record values for the cases when the future sample size is fixed and when it is random. Several Bayesian prediction intervals are derived. Finally, the results of a simulation study and a numerical example are presented for illustrating all the inferential procedures developed here.     

Comparisons of several Pareto distributions based on record values

In order to avoid wrong conclusions in any further analysis, it is of importance to conduct a formal comparison for characteristic quantities of the distributions. These characteristic quantities we are familiar with include mean, quantity and reliability function, and so on. In this paper, we consider two tests aiming at the comparisons for function of parameters in Pareto distribution based on record values. They are generalized p-value-based test and parametric bootstrap-based test, respectively. The resulting procedures are easy to compute and are applicable to small samples. A simulation study is conducted to investigate and compare the performance of the proposed tests. A phenomenon we note is that generalized p-value-based test almost uniformly outperforms the parametric bootstrap-based test.     

An objective Bayesian analysis of the two-parameter half-logistic distribution based on progressively type-II censored samples

This paper develops an objective Bayesian analysis method for estimating unknown parameters of the half-logistic distribution when a sample is available from the progressively Type-II censoring scheme. Noninformative priors such as Jeffreys and reference priors are derived. In addition, derived priors are checked to determine whether they satisfy probability-matching criteria. The Metropolis–Hasting algorithm is applied to generate Markov chain Monte Carlo samples from these posterior density functions because marginal posterior density functions of each parameter cannot be expressed in an explicit form. Monte Carlo simulations are conducted to investigate frequentist properties of estimated models under noninformative priors. For illustration purposes, a real data set is presented, and the quality of models under noninformative priors is evaluated through posterior predictive checking.     

Estimation and prediction of the Kumaraswamy distribution based on record values and inter-record times

ABSTRACT

The maximum likelihood and Bayesian approaches for estimating the parameters and the prediction of future record values for the Kumaraswamy distribution has been considered when the lower record values along with the number of observations following the record values (inter-record-times) have been observed. The Bayes estimates are obtained based on a joint bivariate prior for the shape parameters. In this case, Bayes estimates of the parameters have been developed by using Lindley's approximation and the Markov Chain Monte Carlo (MCMC) method due to the lack of explicit forms under the squared error and the linear-exponential loss functions. The MCMC method has been also used to construct the highest posterior density credible intervals. The Bayes and the maximum likelihood estimates are compared by using the estimated risk through Monte Carlo simulations. We further consider the non-Bayesian and Bayesian prediction for future lower record values arising from the Kumaraswamy distribution based on record values with their corresponding inter-record times and only record values. The comparison of the derived predictors are carried out by using Monte Carlo simulations. Real data are analysed for an illustration of the findings.