These data include all the code needed to generate the figures for the accompanying paper.  . The code needed to generate the figures is provided nearby,  but simulation output needs to be read in from the auxiliary files (listed below) to reproduce the exact output.  Alternatively, new simulations can be generated using the notebook.

"Pedigrees N=30 neutral II"
"Initial population N=30 II"
"p, zA, zD, N=30 neutral III"
"p, zA, zD, N=30 neutral IV"
"p, zA, zD, N=30 selected k=5";
"Neutral identities 16th Jan"
"selected reps N=30 selected k=1,2 300 reps III"

The first  section of the main notebook summarises how to calculate identity coefficients from a pedigree; see the ‘notebook Algorithm …. “ for a fuller explanation.  

Some code (immediately above figures) is hidden: Use "Cell Properties>Open" to reveal.

In the Mathematica notebook, ancillary runs are hidden (double-click brackets on the right to see them).

A free Mathematica reader is available from wolfram.com