Blog Archives

Converting plots to data II

January 26, 2014
By
Converting plots to data II

There were quite some reactions on lasts week's post converting plots to data. Two additional programs were mentioned; WebPlotDigitzier and DataThief. Note that google or wikipedia has a bunch of alternatives. Datathief is shareware, a java program, so...

Read more »

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 »