Publications by Wingfeet
When should I change to snow tires in Netherlands
The Royal Netherlands Meteorological Institute has weather information by day for a number of Dutch stations. In this post I want to use those data for a practical problem: when should I switch to winter tires? (or is that snow tires? In any case nails or metal studs are forbidden in Netherlands). Netherlands does not have laws prescr...
3085 sym 2 img
Change in temperature in Netherlands over the last century
I read a post ‘race for the warmest year’ at sargasso.nl. They used a plot, originating from Ed Hawkins to see how 2014 progressed to be warmest year. Obviously I wanted to make the same plot using R. In addition, I wondered which parts of the year had most increased in temperature.DataSimilar to last week, data are acquired from ...
5148 sym 6 img
SAS PROC MCMC in R: Nonlinear Poisson Regression Models
In exercise 61.1 the problem is that the model has bad mixing. In the SAS manual the mixing is demonstrated after which a modified distribution is used to fix the model.In this post the same problem is tackled in R; MCMCpack, RJags, RStan and LaplaceDemon. MCMCpack has quite some mixing problems, RStan seems to do best.DataTo quote th...
5125 sym 16 img
Monthly Weather in Netherlands
When I downloaded the KNMI meteorological data, the intention was to do something which takes more than just the computers memory. While it is clearly not big data, at the very least 100 years of daily data is not small either. So I took along a load of extra variables to see what trouble I would run into. I did not run out of memory,...
2877 sym 10 img
Merry Christmas
Based on The DO loop, since I wanted a fractal Christmas tree and there is no point in inventing what has been made already. Besides, this is not the first time this year that I used R to do what has been done in SAS.Code# http://blogs.sas.com/content/iml/2012/12/14/a-fractal-christmas-tree/# Each row is a 2×2 linear transformation ...
1968 sym 2 img
A time series contest attempt
I saw the post a time series contest on Rob J Hyndman’s blog. Since I am still wanting to play around with some bigger data sets, so I went to the source website https://drive.google.com/folderview?id=0BxmzB6Xm7Ga1MGxsdlMxbGllZnM&usp=sharing and got myself the data. One warning, if you are reading this to know how to get close to ...
6908 sym 6 img
Capital in the Netherlands, 2006-2013
I guess most people have heard of Piketty and his book capital in the twenty-First century. I don’t have that book, but the media attention has made me wonder if I could see any trends in Dutch public data. As I progressed with this post, I have concluded that these data could not tell me much about longer term trends, however, one ...
15023 sym 12 img
Multivariate analysis of death rate on the map of Europe
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 translated in RGB colors and placed in the map of Europe.DataDeath rates are from Euro...
3693 sym 4 img
SAS PROC MCMC example in R: Logistic Regression Random-Effects Model
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 factorial layout, with two types of seeds, O. aegyptiaca 75 and O. aegyptiaca 73, and two ro...
9671 sym
Odds, odds ratio, language and intuition
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 background. My mother tongue is Dutch, which does not have a word for odds. Hence ...
1129 sym