A Bayesian method of sample size determination with practical applications

Summary.  The problem motivating the paper is the determination of sample size in clinical trials under normal likelihoods and at the substantive testing stage of a financial audit where normality is not an appropriate assumption. A combination of analytical and simulation-based techniques within the Bayesian framework is proposed. The framework accommodates two different prior distributions: one is the general purpose fitting prior distribution that is used in Bayesian analysis and the other is the expert subjective prior distribution, the sampling prior which is believed to generate the parameter values which in turn generate the data. We obtain many theoretical results and one key result is that typical non-informative prior distributions lead to very small sample sizes. In contrast, a very informative prior distribution may either lead to a very small or a very large sample size depending on the location of the centre of the prior distribution and the hypothesized value of the parameter. The methods that are developed are quite general and can be applied to other sample size determination problems. Some numerical illustrations which bring out many other aspects of the optimum sample size are given.     

Some large trial properties of minimum likelihood allocation

The method of minimum likelihood allocation (MLA) for allocating subjects to treatments in a clinical trial amounts to checking at each stage which allocation would lead an outside observer to find the least evidence of a relationship between treatment and factors of prognostic significance, assuming that the observer would use a linear exponential model. One advantage of MLA is that results from game theory and likelihood theory can be used to prove it has desirable long run properties. Two of these demonstrated here are (1) ‘consistency’, in the sense that the average likelihood ratio which measures design imbalance tends to zero, and (2) ‘efficiency’ in the sense that the variance estimates of treatment effects will tend to be minimized in the long run.     

Some computational problems related to multinomial trials

Efficient numerical algorithms are developed to evaluate several probabilities related to multinomial trials.In the first part of the paper, the probability distribution of the number of trials until the alternatives j, j = 1,… m, have occurred at least i_j times is computed. The multinomial trials involve the m alternatives l,…, m, with positive probabilities P_l-P_m of occurrence. In the second part, several aspects of a multinomial subset selection problem, discussed by S. S. Gupta and K. Nagel, are investigated.     

Probability computations for an extension of the geometric distribution

An extension of the stochastic process associated with the geometric distribution is presented. Combinatorial arguments are used to derive probabilities for various events of interest. Probabilities are approximated by evaluating truncated series. Bounds on the errors of approximation are developed. An example is presented and some additional applications are noted.     

Robust inference from multiple test statistics via permutations: a better alternative to the single test statistic approach for randomized trials

Formal inference in randomized clinical trials is based on controlling the type I error rate associated with a single pre‐specified statistic. The deficiency of using just one method of analysis is that it depends on assumptions that may not be met. For robust inference, we propose pre‐specifying multiple test statistics and relying on the minimum p‐value for testing the null hypothesis of no treatment effect. The null hypothesis associated with the various test statistics is that the treatment groups are indistinguishable. The critical value for hypothesis testing comes from permutation distributions. Rejection of the null hypothesis when the smallest p‐value is less than the critical value controls the type I error rate at its designated value. Even if one of the candidate test statistics has low power, the adverse effect on the power of the minimum p‐value statistic is not much. Its use is illustrated with examples. We conclude that it is better to rely on the minimum p‐value rather than a single statistic particularly when that single statistic is the logrank test, because of the cost and complexity of many survival trials. Copyright © 2013 John Wiley & Sons, Ltd.     

More Powerful Tests from Confidence Interval p Values

In this article the problem of comparing two independent binomial populations is considered. It is shown that the test based on the confidence interval p value of Berger and Boos often is uniformly more powerful than the standard unconditional test. This test also requires less computational time.     

论庭审“合作”——基于语用学视角

庭审是一种高度程序化、机构化的活动,庭审的各个参与主体都有明确的角色定位与权力认知,其间又存在复杂的目的关系。尽管庭审中充斥着各种各样的不遵守合作原则的情况,表现出明显的不合作现象,但是,掌控话语权的法官、公诉人和律师等会利用自己的权力,使用各种策略和手段来迫使当事人合作,使庭审重新回到正常轨道,从而在宏观层面实现庭审的合作。     

Evaluation of single and mixed verbal operant arrangements for teaching mands and tacts

The purpose of this series of experiments was to evaluate the effects of mixed mand-tact arrangements on the acquisition of mands and tacts in preschool-aged children. In Experiment 1, the effects of three training arrangements (mand-only training, tact-only training, and mand-tact training) were investigated with 3 typically developing children. Rates of acquisition in single (mand-only and tact-only) versus mixed (mand-tact) presentation were comparable, in contrast to earlier investigations. Experiment 2 attempted to clarify the equivocal findings of Experiment 1 by directly replicating the Carroll and Hesse (1987) investigation with 2 typically developing children. Results again demonstrated no clear benefit of mixed verbal operant training on tact acquisition. In Experiment 3, these same arrangements were evaluated with a boy with autism and included assessments to determine that a relevant establishing operation was in effect prior to each mand training session. Experiment 3 again failed to demonstrate the facilitative effects of mand-tact training on the acquisition of mands or tacts. Taken together, the data from these three experiments fail to support the improved efficiency of mand-tact training suggested by prior studies. Findings are discussed in the context of future research investigating mixed verbal operant arrangements.     

The potential for bias in reporting of industry-sponsored clinical trials

Concerns about potentially misleading reporting of pharmaceutical industry research have surfaced many times. The potential for duality (and thereby conflict) of interest is only too clear when you consider the sums of money required for the discovery, development and commercialization of new medicines. As the ability of major, mid-size and small pharmaceutical companies to innovate has waned, as evidenced by the seemingly relentless decline in the numbers of new medicines approved by Food and Drug Administration and European Medicines Agency year-on-year, not only has the cost per new approved medicine risen: so too has the public and media concern about the extent to which the pharmaceutical industry is open and honest about the efficacy, safety and quality of the drugs we manufacture and sell. In 2005 an Editorial in Journal of the American Medical Association made clear that, so great was their concern about misleading reporting of industry-sponsored studies, henceforth no article would be published that was not also guaranteed by independent statistical analysis. We examine the precursors to this Editorial, as well as its immediate and lasting effects for statisticians, for the manner in which statistical analysis is carried out, and for the industry more generally.     

Adaptive Percolation Using Subjective Likelihoods

A phenomenon that I call “adaptive percolation” commonly arises in biology, business, economics, defense, finance, manufacturing, and the social sciences. Here one wishes to select a handful of entities from a large pool of entities via a process of screening through a hierarchy of sieves. The process is not unlike the percolation of a liquid through a porous medium. The probability model developed here is based on a nested and adaptive Bayesian approach that results in the product of beta-binomial distributions with common parameters. The common parameters happen to be the observed data. I call this the percolated beta-binomial distribution . The model turns out to be a slight generalization of the probabilistic model used in percolation theory. The generalization is a consequence of using a subjectively specified likelihood function to construct a probability model. The notion of using likelihoods for constructing probability models is not a part of the conventional toolkit of applied probabilists. To the best of my knowledge, a use of the product of beta-binomial distributions as a probability model for Bernoulli trials appears to be new. The development of the material of this article is illustrated via data from the 2009 astronaut selection program, which motivated this work.