Diagnostics for generalized Poisson regression models with errors in variables

In this paper, we develop diagnostic methods for generalized Poisson regression (GPR) models with errors in variables based on the corrected likelihood. The one-step approximations of the estimates in the case-deletion model are given and case-deletion and local influence measures are presented. Meanwhile, based on a corrected score function, the testing statistics for the significance of dispersion parameters in GPR models with measurement errors are investigated. Finally, illustration of our methodology is given through numerical examples.     

Computations of Bayesian t-values for multiple comparisons

Results from a simulation study of the power of eight statistics for testing that a sample is form a uniform distribution on the unit interval are reported. Power is given for each statistic against four classes if alternatives. The statistics studied include the discrete Pearson chi-square with ten and twenty cells, X² ₁₀ and X² ₂₀; Kolmogorov-smirov, D; Cramer-Von Mises, W²; Watson, U²; Anderson-Darling, A; Greenwood. G;and a new statistic called O A modified form of each of these statistic is also studied by first transforming the sample using a transformation given by Durbin. On the basis of the results observed in this study, the Watson U² statistic is recommended as a general test for uniformity.     

Small-sample comparison of the exact and asymptotic upper tail probabilities of chi-squared goodness-of-fit statistics: the binomila and the mixture binomial

The exact and asymptotic upper tail probabilities (α = .10, .05, .01, .001) of the three chi-squared goodness-of-fit statistics Pearson's X ², likelihood ratioG ², and powerdivergence statisticD ²(λ), with λ= 2/3 are compared by complete enumeration for the binomial and the mixture binomial. For the two-component mixture binomial, three cases have been distinguished. 1. Both success probabilities and the mixing weights are unknwon. 2. One of the two success probabilities is known. And 3., the mixing weights are known. The binomial was investigated for the number of cellsk, being between 3 and 6 with sample sizes between 5 and 100, for k = 7 with sample sizes between 5 and 45, and for k = 10 with sample sizes ranging from 5 to 20. For the mixture binomial, solely k = 5 cells were considered with sample sizes from 5 to 100 and k = 8 cells with sample sizes between 4 and 20. Rating the relative accuracy of the chi-squared approximation in terms of ±10% and ±20% intervals around α led to the following conclusions for the binomial: 1. Using G2 is not recommendable. 2. At the significance levels α=.10 and α=.05X ² should be preferred over D ²; D ² is the best choice at α = .01. 3. Cochran's (1954; Biometrics, 10, 417-451) rule for the minimum expectation when using X ² seems to generalize to the binomial for G ² and D ² ; as a compromise, it gives a rather strong lower limit for the expected cell frequencies in some circumstances, but a rather liberal in others. To draw similar conclusions concerning the mixture binomial was not possible, because in that case, the accuracy of the chi-squared approximation is not only a function of the chosen test statistic and of the significance level, but also heavily depends on the numerical value of theinvolved unknown parameters and on the hypothesis to be tested. Thereto, the present study may give rise only to warnings against the application of mixture models to small samples.     

A joint model of binary and longitudinal data with non-ignorable missingness,with application to marital stress and late-life major depression in women

Understanding how long-term marital stress affects major depressive disorder (MDD) in older women has clinical implications for the treatment of women at risk. In this paper, we consider the problem of predicting MDD in older women (mean age 60) from a marital stress scale administered four times during the preceding 20-year period, with a greater dropout by women experiencing marital stress or MDD. To analyze these data, we propose a Bayesian joint model consisting of: (1) a linear mixed effects model for the longitudinal measurements, (2) a generalized linear model for the binary primary endpoint, and (3) a shared parameter model for the missing data mechanism. Our analysis indicates that MDD in older women is significantly associated with higher levels of prior marital stress and increasing marital stress over time, although there is a generally decreasing trend in marital stress. This is the first study to propose a joint model for incompletely observed longitudinal measurements, a binary primary endpoint, and non-ignorable missing data; a comparison shows that the joint model yields better predictive accuracy than a two-stage model. These findings suggest that women who experience marital stress in mid-life need treatment to help prevent late-life MDD, which has serious consequences for older persons.     

Diagnostic tools for approximate Bayesian computation using the coverage property

Approximate Bayesian computation (ABC) is an approach to sampling from an approximate posterior distribution in the presence of a computationally intractable likelihood function. A common implementation is based on simulating model, parameter and dataset triples from the prior, and then accepting as samples from the approximate posterior, those model and parameter pairs for which the corresponding dataset, or a summary of that dataset, is ‘close’ to the observed data. Closeness is typically determined though a distance measure and a kernel scale parameter. Appropriate choice of that parameter is important in producing a good quality approximation. This paper proposes diagnostic tools for the choice of the kernel scale parameter based on assessing the coverage property, which asserts that credible intervals have the correct coverage levels in appropriately designed simulation settings. We provide theoretical results on coverage for both model and parameter inference, and adapt these into diagnostics for the ABC context. We re‐analyse a study on human demographic history to determine whether the adopted posterior approximation was appropriate. Code implementing the proposed methodology is freely available in the R package abctools .     

Higher Moments and Prediction‐Based Estimation for the COGARCH(1,1) Model

COGARCH models are continuous time versions of the well‐known GARCH models of financial returns. The first aim of this paper is to show how the method of prediction‐based estimating functions can be applied to draw statistical inference from observations of a COGARCH(1,1) model if the higher‐order structure of the process is clarified. A second aim of the paper is to provide recursive expressions for the joint moments of any fixed order of the process. Asymptotic results are given, and a simulation study shows that the method of prediction‐based estimating function outperforms the other available estimation methods.     

A note on inference for P(X < Y) for right truncated exponentially distributed data

In this paper, a likelihood based analysis is developed and applied to obtain confidence intervals and p values for the stress-strength reliability R  =  P(X  <  Y) with right truncated exponentially distributed data. The proposed method is based on theory given in Fraser et al. (Biometrika 86:249–264, 1999) which involves implicit but appropriate conditioning and marginalization. Monte Carlo simulations are used to illustrate the accuracy of the proposed method.     

A family of shrinkage estimators for Weibull shape parameter in censored sampling

In this paper a new class of shrinkage estimators has been introduced for the shape parameter in an independently identically distributed two-parameterWeibull model under censored sampling. The main idea is to incorporate the prior guessed value by correcting the standard estimator, which is essentially an unbiased estimator, with optimally weighted ratios of the guessed value and the standard estimator, instead of considering a convex combination of the standard estimator and the difference of the guessed value and the standard estimator. The resulting estimator dominates the standard estimator in a surprisingly large neighborhood of the guessed value. The suggested estimator has also been compared with the minimum mean squared error estimator and a class of estimators suggested by Singh and Shukla in IAPQR Trans 25(2), 107–118, 2000. It is found that the suggested class of estimators has lesser bias as well as lesser mean squared error than its competitors subject to certain conditions.      

Designing completely randomized experiments to detect ordered treatment effects: an empirical study

In designing experiments the researcher frequently must decide as to how to allocate fixed resources among k factor levels (Cox (1958)). This study investigates the effects on the power of a test caused by changes in: the sample size (n); the number of factor levels (k); the allocation of fixed total observations (N) among k and n: the shift parameter (ø); the type of parent population sampled; and, the type of ordered location alternative involved. Using Monte Carlo methods the powers of eight test procedures specifically devised to detect ordered treatment effects under completely randomized designs were evaluated along with those of the more general one-way F test. The results are of interest to researchers in all fields of application.     

A new variable selection method for uniform designs

As an important class of space-filling designs, uniform designs (UDs) choose a set of points over a certain domain such that these points are uniformly scattered, under a specific discrepancy measure. They have been applied successfully in many industrial and scientific experiments since they appeared in 1980. A noteworthy and practical advantage is their ability to investigate a large number of high-level factors simultaneously with a fairly economical set of experimental runs. As a result, UDs can be properly used as experimental plans that are intended to derive the significant factors from a list of many potential ones. To this end, a new screening procedure is introduced via penalized least squares. A simulation study is conducted to support the proposed method, which reveals that it can be considered quite promising and expedient, as judged in terms of Type I and Type II error rates.