A comparison of the Hosmer–Lemeshow,Pigeon–Heyse,and Tsiatis goodness-of-fit tests for binary logistic regression under two grouping methods

Algebraic relationships between Hosmer–Lemeshow (HL), Pigeon–Heyse (J²), and Tsiatis (T) goodness-of-fit statistics for binary logistic regression models with continuous covariates were investigated, and their distributional properties and performances studied using simulations. Groups were formed under deciles-of-risk (DOR) and partition-covariate-space (PCS) methods. Under DOR, HL and T followed reported null distributions, while J² did not. Under PCS, only T followed its reported null distribution, with HL and J² dependent on model covariate number and partitioning. Generally, all had similar power. Of the three, T performed best, maintaining Type-I error rates and having a distribution invariant to covariate characteristics, number, and partitioning.     

Bayesian inference for joint location and scale nonlinear models with skew-normal errors

A regression model with skew-normal errors provides a useful extension for ordinary normal regression models when the dataset under consideration involves asymmetric outcomes. In this article, we explore the use of Markov Chain Monte Carlo (MCMC) methods to develop a Bayesian analysis for joint location and scale nonlinear models with skew-normal errors, which relax the normality assumption and include the normal one as a special case. The main advantage of these class of distributions is that they have a nice hierarchical representation that allows the implementation of MCMC methods to simulate samples from the joint posterior distribution. Finally, simulation studies and a real example are used to illustrate the proposed methodology.     

Tempered stable Ornstein– Uhlenbeck processes: A practical view

We study the one-dimensional Ornstein–Uhlenbeck (OU) processes with marginal law given by tempered stable and tempered infinitely divisible distributions. We investigate the transition law between consecutive observations of these processes and evaluate the characteristic function of integrated tempered OU processes with a view toward practical applications. We then analyze how to draw a random sample from this class of processes by considering both the classical inverse transform algorithm and an acceptance–rejection method based on simulating a stable random sample. Using a maximum likelihood estimation method based on the fast Fourier transform, we empirically assess the simulation algorithm performance.     

Gibbs sampling method for the Bayesian adaptive elastic net

This article considers the adaptive elastic net estimator for regularized mean regression from a Bayesian perspective. Representing the Laplace distribution as a mixture of Bartlett–Fejer kernels with a Gamma mixing density, a Gibbs sampling algorithm for the adaptive elastic net is developed. By introducing slice variables, it is shown that the mixture representation provides a Gibbs sampler that can be accomplished by sampling from either truncated normal or truncated Gamma distribution. The proposed method is illustrated using several simulation studies and analyzing a real dataset. Both simulation studies and real data analysis indicate that the proposed approach performs well.     

Approximation of the Fisher price index by using the Lloyd–Moulton index: Simulation study

The Lloyd–Moulton price index does not make use of current-period expenditure data and, as it is commonly known, it allows us to approximate superlative indices, in particular the Fisher price index. This is a very important property for the inflation measurement and the Consumer Price Index bias calculations. In this article, we verify the utility of the Lloyd–Moulton price index in the Fisher price index approximation. We propose a simple modification of that index which reduces the variation of the estimator of an unknown parameter in this index formula. We also examine the influence of the price volatility on the quality of estimation of the parameter from the Lloyd–Moulton formula.     

On the goodness-of-fit test based on the ratio of cumulative hazard functions with the Type II censored data

For comparing two cumulative hazard functions, we consider an extension of the Kullback–Leibler information to the cumulative hazard function, which is concerning the ratio of cumulative hazard functions. Then we consider its estimate as a goodness-of-fit test with the Type II censored data. For an exponential null distribution, the proposed test statistic is shown to outperform other test statistics based on the empirical distribution function in the heavy censoring case against the increasing hazard alternatives.     

Comparison of confidence intervals for the Behrens-Fisher problem

The Behrens–Fisher problem concerns the inferences for the difference between means of two independent normal populations without the assumption of equality of variances. In this article, we compare three approximate confidence intervals and a generalized confidence interval for the Behrens–Fisher problem. We also show how to obtain simultaneous confidence intervals for the three population case (analysis of variance, ANOVA) by the Bonferroni correction factor. We conduct an extensive simulation study to evaluate these methods in respect to their type I error rate, power, expected confidence interval width, and coverage probability. Finally, the considered methods are applied to two real dataset.     

Estimation of the exponential Pareto II distribution parameters

In this article, we estimate the parameters of exponential Pareto II distribution by two new methods. The first one is based on the principle of maximum entropy (POME) and the second is by Kullback–Leibler divergence of survival function (KLS). Monte Carlo simulated data are used to evaluate these methods and compare them with the maximum likelihood method. Finally, we fit this distribution to a set of real data by estimation procedures.     

Estimating and contextualizing the attenuation of odds ratios due to non collapsibility

The odds ratio is a measure commonly used for expressing the association between an exposure and a binary outcome. A feature of the odds ratio is that its value depends on the choice of the distribution over which the probabilities in the odds ratio are evaluated. In particular, this means that an odds ratio conditional on a covariate may have a different value from an odds ratio marginal on the covariate, even if the covariate is not associated with the exposure (not a confounder). We define the individual odds ratio (IORs) and population odds ratios (PORs) as the ratio of the odds of the outcome for a unit increase in the exposure, respectively, for an individual in the population and for the whole population, in which case the odds are averaged across the population. The attenuation of conditional odds ratio, marginal odds ratio, and PORs from the IOR is demonstrated in a realistic simulation exercise. The degree of attenuation differs in the whole population and in a case–control sample, and the property of invariance to outcome-dependent sampling is only true for the IOR. The relevance of the non collapsibility of odds ratios in a range of methodological areas is discussed.