Blog Archives

Bayesian analysis of sensory profiling data, part 2

February 9, 2014
By
Bayesian analysis of sensory profiling data, part 2

Last week I made the core of a Bayesian model for sensory profiling data. This week the extras need to be added. That is, there are a bunch of extra interactions and the error is dependent on panelists and descriptors.Note that where last week I pointe...

Read more »

Bayesian analysis of sensory profiling data

February 2, 2014
By
Bayesian analysis of sensory profiling data

I looked at Bayesian analysis of sensory profiling data in May and June 2012. I do remember not being totally happy with the result and computations taking a bit more time than I wanted. But now it is 2014, I can use STAN and I have been thinking about...

Read more »

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 »