Blog Archives

More ordinal data display

March 30, 2013
By
More ordinal data display

The past two weeks I made a post regarding analyzing ordinal data with R and JAGS. The calculations in the second part made me realize I could actually get top two box intervals out of R. This demonstrated here. For that I needed the inv...

Read more »

Ordinal data with JAGS

March 25, 2013
By

Last week is had a look at the standard R routines for estimating models for ordinal data. This week, I want to have a look at JAGS for examining the same data. To be honest, most of it is taking an example (inhaler) and removing code. To my surpr...

Read more »

Ordinal Data

March 17, 2013
By
Ordinal Data

I expect to be getting some ordinal data, from 5 or 9 point rating scales, pretty soon, so I am having a look ahead how to treat those. Often ANOVA is used, even though it is well known not to be ideal fro a statistical point of view, so that is the st...

Read more »

More sequential testing for triangle tests

March 10, 2013
By
More sequential testing for triangle tests

I looked before at triangle tests and at sequential testing in triangle tests (blog entry). In the latter post it was demonstrated that a sequential test is possible, without costs in desired error of the first kind. The latter because t...

Read more »

Occurrence of Alzheimer deaths in Netherlands

March 3, 2013
By
Occurrence of Alzheimer deaths in Netherlands

I wanted to know a bit more about the number of people suffering from Alzheimer. It has been getting more attention and it runs in my family. From CBS (Statistics Netherlands) I pulled the data regarding the number of people dying from it. I am af...

Read more »

Earthquakes in Netherlands

February 24, 2013
By
Earthquakes in Netherlands

In the Netherlands we have Natural Gas. Unfortunately winning this gas seems to cause some quakes. As quakes go, they are not strong. However, our buildings are not made to resist quakes, before 1986 they were unheard of, so there is some damage. It is now predicted they could get stronger and more frequent. This caused a bit of a...

Read more »

A look at strucchange and segmented

February 17, 2013
By
A look at strucchange and segmented

After last week's post it was commented that strucchange and segmented would be more suitable for my purpose. I had a look at both. Strucchange can find a jump in a time series, which was what I was looking for. In contrast segmented is more suitable f...

Read more »

Finding a jump in data

February 9, 2013
By
Finding a jump in data

The other day I saw a plot with some data with a clear jump in it. So I wondered if it were possible to determine the position of that jump. Data with a jump is easy enough made:library(ggplot2)n <- 100mydata <- data.frame(x=runif(n))mydata$...

Read more »

For descriptive statistics, values below LLOQ set to …

February 3, 2013
By
For descriptive statistics, values below LLOQ set to …

That is what I read the other day. For calculation of descriptive statistics, values below the LLOQ (lower limit of quantification)  were set to.... Then I wondered, wasn't there a trick in JAGS to incorporate the presence of missing data while es...

Read more »

European Fishing

January 27, 2013
By
European Fishing

I am playing around with Eurostat data and ggplot2 a bit more. As I progress it seems the plotting gets more easy, the data pre-processing a bit more simple and the surprises on the data stay.Eurostat dataThe data used are fish_fleet (number of ships) and fish_pr (production=catch+aquaculture). After a bit of year selection, 1992 and later, I decided to...

Read more »