Blog Archives

Unemployment of Europe in 2014 by NUTS 2 region

April 25, 2015
By
Unemployment of Europe in 2014 by NUTS 2 region

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 »

PK analysis: Theoph again

April 19, 2015
By
PK analysis: Theoph again

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 »

Cefamandole: PK after IV with six subjects in JAGS

April 12, 2015
By
Cefamandole: PK after IV with six subjects in JAGS

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 »

Hierarchical Two Compartimental PK Model

April 5, 2015
By

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 »

Space Launch Sites over Time

March 29, 2015
By
Space Launch Sites over Time

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 »

Launch into space

March 22, 2015
By
Launch into space

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 »

Medicines under evaluation

March 15, 2015
By

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 »

SAS PROC MCMC example in R: Nonlinear Poisson Regression Multilevel Random-Effects Model

March 8, 2015
By

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 »

Should I use premium Diesel? Setup

March 1, 2015
By
Should I use premium Diesel? Setup

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 »

Priors on odds and probability of success

February 22, 2015
By
Priors on odds and probability of success

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 »