Reliability characteristics of Farlie–Gumbel–Morgenstern family of bivariate distributions

Abstract

In this paper, we study the Farlie–Gumbel–Morgenstern family of bivariate distributions from a reliability point of view. The properties of this family of distributions and the association between the two variables are investigated by studying the local dependence function and the association measure defined by Clayton (1978 Clayton, D.G. (1978). A model for association in bivariate life tables and its applications in epidemiological studies of familial tendency in chronic disease incidence. Biometrika 65:141–151.[Crossref], [Web of Science ®] , [Google Scholar]). We also study the effect of the association parameter on the hazard components, the failure rate of the series system, and the regression mean residual life of a parallel system. Stochastic comparisons with respect to the association parameter are also studied. Some examples are provided to illustrate the results.     

Moments of order statistics of Topp–Leone distribution

We derive explicit algebraic expressions for both of the single and product moments of order statistics from Topp–Leone distribution. We also give an identity about single moments of order statistics. These expressions will be useful for computational purposes.     

Measures of biomarker dependence using a copula-based multivariate epsilon–skew–normal family of distributions

In this note we develop a new multivariate copula model based on epsilon–skew–normal marginal densities for the purpose of examining biomarker dependency structures. We illustrate the flexibility and utility of this model via a variety of graphical tools and a data analysis example pertaining to salivary biomarker. The multivariate normal model is a sub-model of the multivariate epsilon–skew–normal distribution.     

Characteristics of some classes of space–time covariance functions

Although a wide list of classes of space–time covariance functions is now available, selecting an appropriate class of models for a variable under study is still difficult and it represents a priority problem with respect to the choice of a particular model of a specified class. Then, knowing the characteristics of various classes of covariances, and their auxiliary functions, and matching those with the characteristics of the empirical space–time covariance surface might be helpful in the selection of a suitable class. In this paper some characteristics, such as behavior at the origin, asymptotic behavior, nonseparability and anisotropy aspects, are studied for some well known classes of covariance models of stationary space–time random fields. Moreover, some important issues related to modeling choices are described and a case study is presented.     

A modification of the Whittaker–Henderson method of graduation

This paper presents a modified Whittaker–Henderson (WH) Method of Graduation. After giving a closed-form solution, we show that it is of practical use because it provides not only a smoothed series identical to that of the WH graduation, but also an extrapolation beyond the sample limit of current data. In addition, we introduce two other penalized least squares problems and show that they provide the same results as those of the modified WH graduation.     

Estimation of the Parameters of the Birnbaum–Saunders Distribution

Some alternative estimators to the maximum likelihood estimators of the two parameters of the Birnbaum–Saunders distribution are proposed. Most have high efficiencies as measured by root mean square error and are robust to departure from the model as well as to outliers. In addition, the proposed estimators are easy to compute. Both complete and right-censored data are discussed. Simulation studies are provided to compare the performance of the estimators.     

Additional Tables for Steel–Dwass–Critchlow–Fligner Distribution-Free Multiple Comparisons of Three Treatments

ABSTRACT

Additional critical points are presented for the Steel–Dwass–Critchlow–Fligner distribution-free multiple comparison procedure for comparing all pairs of three population medians in the one-way layout. A computational technique developed by van de Wiel is used to find critical points yielding an experimentwise error rate of approximately 0.01, 0.05, and 0.10 for a total sample size of at most 30, with individual sample sizes from 4 to 10 and a maximum sample size of at least 8, and for equal sample sizes from 8 to 14. Additional discussion is given regarding step-down testing methods and the dangers of using the Steel–Dwass–Critchlow–Fligner procedure with unequal sample sizes if two of the sample sizes are very small.     

On the effect of item nonresponse on the estimation of a two–panel–waves wage equation

Summary: In this paper, we present results of the estimation of a two–panel–waves wage equation based on completely observed units and on a multiply imputed data set. In addition to the survey information, reliable income data is available from the register. These external data are used to assess the reliability of wage regressions that suffer from item nonresponse. The findings reveal marked differences between the complete case analyses and both versions of multiple imputation analyses. We argue that the results based on the multiply imputed data sets are more reliable than those based on the complete case analysis.* We would like to thank Statistics Finland for providing the data. We are also very grateful to Susanna Sandström and Marjo Pyy–Martikainen for their helpful advice using the Finnish data. Helpful comments from Joachim Winter and participants of the Workshop on Item Nonresponse and Data Quality in Large Social Surveys, Basel, October, 2003, on an earlier version of the paper are greatfully acknowledged. Further, we would like to thank three anonymous referees and the editor for helpful comments and suggestions.     

A Monte Carlo comparison of Jarque–Bera type tests and Henze–Zirkler test of multivariate normality

In the paper, tests for multivariate normality (MVN) of Jarque-Bera type, based on skewness and kurtosis, have been considered. Tests proposed by Mardia and Srivastava, and the combined tests based on skewness and kurtosis defined by Jarque and Bera have been taken into account. In the Monte Carlo simulations, for each combination of p = 2, 3, 4, 5 number of traits and n = 10(5)50(10)100 sample sizes 10,000 runs have been done to calculate empirical Type I errors of tests under consideration, and empirical power against different alternative distributions. Simulation results have been compared to the Henze–Zirkler’s test. It should be stressed that no test yet proposed is uniformly better than all the others in every combination of conditions examined.     

A Bayesian analysis of the Conway–Maxwell–Poisson cure rate model

The purpose of this paper is to develop a Bayesian analysis for the right-censored survival data when immune or cured individuals may be present in the population from which the data is taken. In our approach the number of competing causes of the event of interest follows the Conway–Maxwell–Poisson distribution which generalizes the Poisson distribution. Markov chain Monte Carlo (MCMC) methods are used to develop a Bayesian procedure for the proposed model. Also, some discussions on the model selection and an illustration with a real data set are considered.     

Rate of convergence for asymptotic variance of the Horvitz–Thompson estimator

Drawing distinct units without replacement and with unequal probabilities from a population is a problem often considered in the literature (e.g. Hanif and Brewer, 1980, Int. Statist. Rev. 48, 317–355). In such a case, the sample mean is a biased estimator of the population mean. For this reason, we use the unbiased Horvitz–Thompson estimator (1951). In this work, we focus our interest on the variance of this estimator. The variance is cumbersome to compute because it requires the calculation of a large number of second-order inclusion probabilities. It would be helpful to use an approximation that does not need heavy calculations. The Hájek (1964) variance approximation provides this advantage as it is free of second-order inclusion probabilities. Hájek (1964) proved that this approximation is valid under restrictive conditions that are usually not fulfilled in practice. In this paper, we give more general conditions and we show that this approximation remains acceptable for most practical problems.     

A class of saturated row–column designs

A new class of row–column designs is proposed. These designs are saturated in terms of eliminating two-way heterogeneity with an additive model. The (m,s)-criterion is used to select optimal designs. It turns out that all (m,s)-optimal designs are binary. Square (m,s)-optimal designs are constructed and they are treatment-connected. Thus, all treatment contrasts are estimable regardless of the row and column effects.     

Estimation of reliability of multicomponent stress–strength for a Kumaraswamy distribution

This article deals with the Bayesian and non Bayesian estimation of multicomponent stress–strength reliability by assuming the Kumaraswamy distribution. Both stress and strength are assumed to have a Kumaraswamy distribution with common and known shape parameter. The reliability of such a system is obtained by the methods of maximum likelihood and Bayesian approach and the results are compared using Markov Chain Monte Carlo (MCMC) technique for both small and large samples. Finally, two data sets are analyzed for illustrative purposes.     

Adaptive optimal scaling of Metropolis–Hastings algorithms using the Robbins–Monro process

ABSTRACT

We present an adaptive method for the automatic scaling of random-walk Metropolis–Hastings algorithms, which quickly and robustly identifies the scaling factor that yields a specified overall sampler acceptance probability. Our method relies on the use of the Robbins–Monro search process, whose performance is determined by an unknown steplength constant. Based on theoretical considerations we give a simple estimator of this constant for Gaussian proposal distributions. The effectiveness of our method is demonstrated with both simulated and real data examples.     

Characterization of th estandard normal probability density within the class of self–adjoint (self–reciprocal) densities

We consider symmetric densities p statisfying the well–known integral equation[ULM0001] The solution proves to be standard normal provided that[ILM0001] for some integer     

Default Bayesian analysis of the Behrens–Fisher problem

In the Bayesian approach, the Behrens–Fisher problem has been posed as one of estimation for the difference of two means. No Bayesian solution to the Behrens–Fisher testing problem has yet been given due, perhaps, to the fact that the conventional priors used are improper. While default Bayesian analysis can be carried out for estimation purposes, it poses difficulties for testing problems. This paper generates sensible intrinsic and fractional prior distributions for the Behrens–Fisher testing problem from the improper priors commonly used for estimation. It allows us to compute the Bayes factor to compare the null and the alternative hypotheses. This default procedure of model selection is compared with a frequentist test and the Bayesian information criterion. We find discrepancy in the sense that frequentist and Bayesian information criterion reject the null hypothesis for data, that the Bayes factor for intrinsic or fractional priors do not.