To adjust or not to adjust for baseline when analyzing repeated binary responses? The case of complete data when treatment comparison at study end is of interest

The benefits of adjusting for baseline covariates are not as straightforward with repeated binary responses as with continuous response variables. Therefore, in this study, we compared different methods for analyzing repeated binary data through simulations when the outcome at the study endpoint is of interest. Methods compared included chi‐square, Fisher's exact test, covariate adjusted/unadjusted logistic regression (Adj.logit/Unadj.logit), covariate adjusted/unadjusted generalized estimating equations (Adj.GEE/Unadj.GEE), covariate adjusted/unadjusted generalized linear mixed model (Adj.GLMM/Unadj.GLMM). All these methods preserved the type I error close to the nominal level. Covariate adjusted methods improved power compared with the unadjusted methods because of the increased treatment effect estimates, especially when the correlation between the baseline and outcome was strong, even though there was an apparent increase in standard errors. Results of the Chi‐squared test were identical to those for the unadjusted logistic regression. Fisher's exact test was the most conservative test regarding the type I error rate and also with the lowest power. Without missing data, there was no gain in using a repeated measures approach over a simple logistic regression at the final time point. Analysis of results from five phase III diabetes trials of the same compound was consistent with the simulation findings. Therefore, covariate adjusted analysis is recommended for repeated binary data when the study endpoint is of interest. Copyright © 2015 John Wiley & Sons, Ltd.     

Bayesian comparative study on binary time series

In this paper, we consider the Bayesian analysis of binary time series with different priors, namely normal, Students' t, and Jeffreys prior, and compare the results with the frequentist methods through some simulation experiments and one real data on daily rainfall in inches at Mount Washington, NH. Among Bayesian methods, our results show that the Jeffreys prior perform better in most of the situations for both the simulation and the rainfall data. Furthermore, among weakly informative priors considered, Student's t prior with 7 degrees of freedom fits the data most adequately.     

Computing the conditional stationary distribution in Markov chains of level-dependent M/G/1-type

This paper considers the computation of the conditional stationary distribution in Markov chains of level-dependent M/G/1-type, given that the level is not greater than a predefined threshold. This problem has been studied recently and a computational algorithm is proposed under the assumption that matrices representing downward jumps are nonsingular. We first show that this assumption can be eliminated in a general setting of Markov chains of level-dependent G/G/1-type. Next we develop a computational algorithm for the conditional stationary distribution in Markov chains of level-dependent M/G/1-type, by modifying the above-mentioned algorithm slightly. In principle, our algorithm is applicable to any Markov chain of level-dependent M/G/1-type, if the Markov chain is irreducible and positive-recurrent. Furthermore, as an input to the algorithm, we can set an error bound for the computed conditional distribution, which is a notable feature of our algorithm. Some numerical examples are also provided.     

Sample size determination for testing equality in frequency data under an incomplete block crossover design

When there are more than two treatments under comparison, we may consider the use of the incomplete block crossover design (IBCD) to save the number of patients needed for a parallel groups design and reduce the duration of a crossover trial. We develop an asymptotic procedure for simultaneously testing equality of two treatments versus a control treatment (or placebo) in frequency data under the IBCD with two periods. We derive a sample size calculation procedure for the desired power of detecting the given treatment effects at a nominal-level and suggest a simple ad hoc adjustment procedure to improve the accuracy of the sample size determination when the resulting minimum required number of patients is not large. We employ Monte Carlo simulation to evaluate the finite-sample performance of the proposed test, the accuracy of the sample size calculation procedure, and that with the simple ad hoc adjustment suggested here. We use the data taken as a part of a crossover trial comparing the number of exacerbations between using salbutamol or salmeterol and a placebo in asthma patients to illustrate the sample size calculation procedure.     

Estimation of coefficient of dispersion using auxiliary information

In this note, we consider the problem of estimation of coefficient of dispersion of a study variable by making use of known value of the coefficient of dispersion of an auxiliary variable. We propose ratio and regression-type estimators. We derive expressions of bias and variance of the proposed estimators to the first order of approximation. The relative efficiencies of the ratio and regression-type estimators with respect to the naïve estimator are investigated through a simulation study.     

Bayesian inference of three-parameter bathtub-shaped lifetime distribution

We investigate a Bayesian inference in the three-parameter bathtub-shaped lifetime distribution which is obtained by adding a power parameter to the two-parameter bathtub-shaped lifetime distribution suggested by Chen (2000). The Bayes estimators under the balanced squared error loss function are derived for three parameters. Then, we have used Lindley's and Tierney–Kadane approximations (see Lindley 1980; Tierney and Kadane 1986) for computing these Bayes estimators. In particular, we propose the explicit form of Lindley's approximation for the model with three parameters. We also give applications with a simulated data set and two real data sets to show the use of discussed computing methods. Finally, concluding remarks are mentioned.     

Nonparametric kernel estimation of the impact of tax policy on the demand for private health insurance in Australia

This paper is motivated by our attempt to answer a policy question: how is private health insurance take‐up in Australia affected by the income threshold at which the Medicare Levy Surcharge (MLS) kicks in? We propose a new difference deconvolution kernel estimator for the location and size of regression discontinuities. We also propose a bootstrapping procedure for estimating the confidence interval for the estimated discontinuity. Performance of the estimator is evaluated by Monte Carlo simulations before it is applied to estimating the effect of the income threshold of MLS on the take‐up of private health insurance in Australia, using contaminated data.     

Optimal Bandwidth Selection in Heteroskedasticity–Autocorrelation Robust Testing

This paper considers studentized tests in time series regressions with nonparametrically autocorrelated errors. The studentization is based on robust standard errors with truncation lag M=bT for some constant b∈(0, 1] and sample size T. It is shown that the nonstandard fixed‐b limit distributions of such nonparametrically studentized tests provide more accurate approximations to the finite sample distributions than the standard small‐b limit distribution. We further show that, for typical economic time series, the optimal bandwidth that minimizes a weighted average of type I and type II errors is larger by an order of magnitude than the bandwidth that minimizes the asymptotic mean squared error of the corresponding long‐run variance estimator. A plug‐in procedure for implementing this optimal bandwidth is suggested and simulations (not reported here) confirm that the new plug‐in procedure works well in finite samples.     

An Efficient Class of Estimators for the Population Mean Using Auxiliary Information in Stratified Random Sampling

This article addresses the problem of estimating the population mean in stratified random sampling using the information of an auxiliary variable. A class of estimators for population mean is defined with its properties under large sample approximation. In particular, various classes of estimators are identified as particular member of the suggested class. It has been shown that the proposed class of estimators is better than usual unbiased estimator, usual combined ratio estimator, usual product estimator, usual regression estimator and Koyuncu and Kadilar (2009 Koyuncu, N., Kadilar, C. (2009). Ratio and product estimators in stratified random sampling. J. Statist. Plan. Infere. 139:2552–2558.[Crossref], [Web of Science ®] , [Google Scholar]) class of estimators. The results have been illustrated through an empirical study.