Blog Archives

Bridge hand distribution: simulation vs exact calculation

December 8, 2012
By
Bridge hand distribution: simulation vs exact calculation

Recently I played bridge with my friends. Being frustrated with several consecutive poor hand distributions we asked ourselves a question what is the probability of having a hand good enough for a small slam. A well known rule of thumb is that you need...

Read more »

ESSA2013 Conference

November 24, 2012
By
ESSA2013 Conference

It has been just announced that during ESSA2013 conference I am planning to organize a special track on "Statistical analysis of simulation models". I hope to get some presentations using GNU R to promote it in social simulation community.It is obvious that GNU R excels in analysis of simulation data. However, very often it can be neatly...

Read more »

Possible error in Bayesian bootstrap

November 8, 2012
By
Possible error in Bayesian bootstrap

After my last post on Bayesian bootstrap I got a question why the sample from Dirichlet distribution is taken as weights for calculating mean in the procedure and not as weights used for sampling from the original data set. Actually this mistake i...

Read more »

Simple Bayesian bootstrap

November 2, 2012
By

Bootstrapping is a very popular statistical technique. However, its Bayesian analogue proposed by Rubin (1981) is not very common. I was looking for an example of its implementation in GNU R and could not find one so I decided to write a snippet presen...

Read more »

Plotting Watts-Strogatz model

September 21, 2012
By
Plotting Watts-Strogatz model

Recently I wanted to reproduce Figure 2 from Watts and Strogatz (1998). The task using igraph is simple but an interesting task was annotation of the resulting plot.Watts-Strogatz model generates graphs that have so called small-world ne...

Read more »

Simulation metamodeling with constraints

September 14, 2012
By
Simulation metamodeling with constraints

Last week I have posted about using simulation metamodeling to verify results of analytical solution of the model. After posting it I realized that the solution presented there can be improved by using knowledge of simulation model struc...

Read more »

Simulation metamodeling with GNU R

September 7, 2012
By
Simulation metamodeling with GNU R

I am one of the organizers of ESSA2013 conference that will take place in September 2013 in Warsaw, Poland. The conference scope is social simulation and in particular methods of statistical analysis of simulation output (metamodeling). As we have just issued Call for Papers for the conference so I decided to post a simple example of a metamodel.Recently I had...

Read more »

Exporting ctree object to Asymptote

August 25, 2012
By
Exporting ctree object to Asymptote

When producing regression or classification trees (standard rpart or ctree from party package) in GNU R I am often unsatisfied with the default plots they produce. One of many possible solutions is to export a tree plot to Asymptote.The ...

Read more »

An example of OOP in GNU R using S4 Classes

August 18, 2012
By
An example of OOP in GNU R using S4 Classes

Recently I have discussed with my friend from WLOG Solutions an implementation of banking cash management engine in GNU R. The code made a nice use of S4 classes so I thought it would be worth showing as an example.The problemEvery commercial bank need...

Read more »

Animation basics for a vacation

August 5, 2012
By
Animation basics for a vacation

Since I have a vacation this time I decided to implement some entertaining graphics. I have chosen to animate a Cassini oval.The task is can be accomplished using polar equation:The implementation of the animation is given by the following code:library...

Read more »