Last week I created a JAGS model combining data from two paper helicopter datasets. This week, I will use the model to find the longest flying one.PredictingThe JAGS/RJAGS system has no predict() function that I know of. What I therefore did is adapt t... [Read more...]

The paper helicopter is one of the devices to explain about design of experiments. The aim is to create the longest flying paper helicopter by means of experimental design.Paper helicopters are a nice example, because they are cheap to make, easy to te... [Read more...]

This is the time when we celebrate the end of the second world war in The Netherlands, so I thought to do somwthig with data from that era. One of the things I enjoyed were books on the sea warfare, such as 'The Cruel Sea'by Nicholas Monsarrat. In that... [Read more...]

A while ago I had a post: 'Should I use premium Diesel? Setup. Since that time the data has been acquired. This post describes the results.DataData is registered by me in 2014 and 2015. 2014 has standard Diesel, while 2015 has premium. Both are fr... [Read more...]

During the Christmas break I worked on some code to show unemployment by NUTS 2 region. At that point no 2014 data was available. When I noticed the 214 was available I dug up the code and plotted again.Data and CodeAs written, the code was made beginn... [Read more...]

This week I wanted to repeat the Theoph PK analysis of two weeks ago in Stan. It also suddenly dawned on me. For a non-linear model in classical statistics we think it normal to provide good initializations. In contrast, for those MCMC samples I rememb... [Read more...]

After last week's post on Theoph I noticed the MEMSS library contained more PK data sets. The Cefamandole data is an IV data set with six subjects. It is somewhat challenging, since there seem to be several elimination rates.DataData are displayed... [Read more...]

In this post I am running the Theoph dataset from MEMSS (Data sets and sample analyses from Pinheiro and Bates, "Mixed-effects Models in S and S-PLUS" (Springer, 2000)) in a JAGS model. To quote the MEMSS manual:'Boeckmann, Sheiner and Beal (1994... [Read more...]

Continuing from last weeks post, I am now looking at space launch sites.DataData are from the main table. In addition, this sites table was manually browsed for interpretation of abbreviations.List of most important sitesJust by running counts the most... [Read more...]

I saw a link to a list with all rocket launches into space the other day. This post contains some plots concerning day of launch made from that.DataData is a fixed format file with eleven columns. Reading fixed format is not very difficult, however, it... [Read more...]

While browsing through the internet I ran into the medicines under evaluation page of the European Medicines Agency. There was a .pdf for each month and I thought it would be interesting to look if I could analyze those data. Each .pdf contains various... [Read more...]

I am slowly working my way through the PROC MCMCexamples. Regarding these data, the SAS manual says: 'This example uses the pump failure data of Gaver and O’Muircheartaigh (1987) to illustrate how to fit a multilevel random-effects model with PROC MCMC. The number of failures and the time of operation ... [Read more...]

Since I drive quite a lot, I have some interest in getting the most km out every Euro spent on fuel. One thing to change is the fuel. The oil companies have a premium fuel, which is supposed to be better for both engine and fuel consumption. On the oth... [Read more...]

In Bayesian Approaches to Clinical Trials and Health-Care Evaluation (David J. Spiegelhalter, Keith R. Abrams, Jonathan P. Myles) they mention that an non informative prior should be uniform over the range of interest. They combine this with the d... [Read more...]

Studying on in Bayesian Approaches to Clinical Trials and Health-Care Evaluation (David J. Spiegelhalter, Keith R. Abrams, Jonathan P. Myles) they state that if you have one success in a trial, that a design is needed to know what that means. For examp... [Read more...]

I am working through Bayesian Approaches to Clinical Trials and Health-Care Evaluation (David J. Spiegelhalter, Keith R. Abrams, Jonathan P. Myles) (referred to as BACTHCE from here on). In chapter three I saw an example (3.13) where I wanted to d... [Read more...]

Beginning 2013 I made a post plotting unemployment in Europe. Last year I did the same. Now that the unemployment numbers of December are on Eurostat, I am making them again. The plots shown are unemployment and its first derivative, both smoothed.You ... [Read more...]

I was reading a statistics book the other day. I am at the beginning. In this section I read '(we) report results as odds ratios, which is more intuitive'. I must have read sentences stating this any number of times. But I don't agree.It may be my back... [Read more...]

In this post I will run SAS example Logistic Regression Random-Effects Model in four R based solutions; Jags, STAN, MCMCpack and LaplacesDemon. To quote the SAS manual: 'The data are taken from Crowder (1978). The Seeds data set is a 2 x 2 fa... [Read more...]

Eurostat has information on death rates by cause and NUTS 2 region. I am trying to get this visually displayed on the map. To get there I map all causes to three dimensions via a principal components analysis. These three dimensions are subsequently tr... [Read more...]