Blog Archives

Convering plots to data

January 18, 2014
By
Convering plots to data

It is a problem which occurs ever so often in applied work, you have a plot, but you want the data. There are at least two programs which can help you there; PlotDigitizer and Engauge Digitizer. I got both on my openSuse machine. Both are available for...

Read more »

Converting a JAGS model to STAN

January 11, 2014
By

For my first experience with STAN I wanted to convert my last JAGS program into STAN. This was a bit more difficult than I expected. The JAGS program was Fe concentration in rainwater including values below detection level.DataData has been explained b...

Read more »

Le Monde puzzle 847 in Julia

January 4, 2014
By

This week I wanted to play around with Julia and exporting the results. I found http://xianblog.wordpress.com/2013/12/29/le-monde-puzzle-847/ to be just the right size to play around with. CodeA function to check if a triplet has the desired property In : function lemonde847(xx) a=,2] b=,8] c=,4] sum(kron(c,kron(a,b)))endOut:lemonde847 (generic function with 1 method)Just a check -...

Read more »

Sinterklaas and Santa Claus gave presents

December 27, 2013
By
Sinterklaas and Santa Claus gave presents

In December my nice little netbook acquired a shiny new openSUSE 13.1, including all the goodies I might want to use in 2014; R, Stan, Julia and Jags. So here is what you might get when you set up a computer for statistics today.JagsI have been using J...

Read more »

Merry Christmas

December 20, 2013
By
Merry Christmas

Merry Christmas for all readers. Please enjoy my tree. codepart <- list(x0=0,y0=0,x1=0,y1=1,    branch1=NULL,branch2=NULL,extend=NULL,    lwd=1,depth=0,col='springgreen')par(mfrow=c(1,1),mar=c(5, 4, 4, 2) + 0.1)segplot <-...

Read more »

plotting y and log(y) in one figure

December 15, 2013
By
plotting y and log(y) in one figure

Sometimes I have the desire to plot both on the linear and on the log scale. To save space just two figures is not my solution. I want to reuse the x-axis, legend, title. This post examines possibilities to do so with standard plot tools, lattice and g...

Read more »

Fukushima region radiation decrease

December 8, 2013
By
Fukushima region radiation decrease

Fukushima is still a topic which gets headlines, and somewhere in a comment was a link to actual and historical radiation data: Fukushima prefecture radioactivity measurement map.  It takes a bit of time to load but then you have a map of the Fukushima region with current radiation levels at measurement stations. Somewhere on that website there is...

Read more »

JAGS model Fe concentration in rainwater including values below detection level

December 1, 2013
By
JAGS model Fe concentration in rainwater including values below detection level

In my previous post I ignored the fact that some data are below the detection level. I would not know how to handle those in a mixed model from lme4 or nlme. However, JAGS can handle these values. Next to that I kept the usual independent variables, su...

Read more »

Dutch Rainwater Composition 1992-2011

November 24, 2013
By
Dutch Rainwater Composition 1992-2011

Last week I examined rainwater composition 1992 to 2005. There is additional data, but National Institute for Public Health has changed equipment in 2005. This week I will add those data.DataData is in a number of spreadsheets. The scrip...

Read more »

Dutch Rainwater Composition 1992-2005.

November 17, 2013
By
Dutch Rainwater Composition 1992-2005.

After reading Blog About Stats' Open Data Index Blog Post I decided to browse a bit in the Open Data Index. Choosing Netherlands and following Emission of Pollutants I ended on a page from National Institute for Public Health. The page&n...

Read more »