
Cute Gibbs sampling for rounded observations
I was attending a course of Bayesian Statistics where this problem showed up: There is a number of individuals, say 12, who take a pass/fail test 15 times. For each individual we have recorded the number of passes, which can go from 0 to 15. Because of confidentiality issues, we are presented with roundedtotheclosestmultipleof3 data…

Analyzing a simple experiment with heterogeneous variances using asreml, MCMCglmm and SAS
I was working with a small experiment which includes families from two Eucalyptus species and thought it would be nice to code a first analysis using alternative approaches. The experiment is a randomized complete block design, with species as fixed effect and family and block as a random effects, while the response variable is growth…

INLA: Bayes goes to Norway
INLA is not the Norwegian answer to ABBA; that would probably be aha. INLA is the answer to ‘Why do I have enough time to cook a threecourse meal while running MCMC analyses?”. Integrated Nested Laplace Approximations (INLA) is based on direct numerical integration (rather than simulation as in MCMC) which, according to people ‘in…

R, Julia and genome wide selection
— “You are a pussy” emailed my friend. — “Sensu cat?” I replied. — “No. Sensu chicken” blurbed my now exfriend. What was this about? He read my post on R, Julia and the shiny new thing, which prompted him to assume that I was the proverbial old dog unwilling (or was it unable?) to…

MidJanuary flotsam: teaching edition
I was thinking about new material that I will use for teaching this coming semester (starting the third week of February) and suddenly compiled the following list of links: