I’ve been performing some simulation studies comparing a Bayesian to a more traditional frequentist estimation approach in a particular problem. To do this I’ve been using the excellent JAGS package, calling it from R. One of the issues I’ve faced is the question of how long to run the MCMC sampler in the Bayesian approach. Use too few iterations, and the chains will not have converged to their stationary distribution, such that the samples will not be from the posterior distribution of the model parameters. In regular data analysis situations, one can make use of the extensive diagnostic toolkit which has been developed over the years. The most popular of these is I believe to examine trace plots from multiple chains, started with dispersed initial values, and also Gelman and Rubin’s Rhat measure.
Multiarm trials – should we allow for multiplicity?
Last week I listened to a great presentation about new trial designs by Mahesh Parmar, director of the Medical Research Council Clinical Trial Unit in London. Among the topics he touched on were multi-arm trials (and extensions), as an attractive alternative to the classic two arm trial. There seem to be a number of advantages to such a trial design, in which in the simplest case, the trial randomizes patients to either control or one of a number of experimental treatments.
R squared/correlation depends on variance of predictor
I’ve written about R squared a few times before. In a discussion I was involved with today the question was raised as to how/whether the R squared in a linear regression model with a single continuous predictor depends on the variance of the predictor variable. The answer to the question is of course yes.