For this example we create a small data set that has a Poisson distribution

Make/O/N=50 data1=30+poissonNoise(5)

To run bootstrap, in each iteration we draw (with replacement) 100 numbers and compute standard WaveStats on the drawn sample.

StatsResample/N=100 /ws=2/ITER=5 data1

The results for five iterations are shown in the table:

Increasing the number of samples does not change the results substantially:

StatsResample/N=10000 /ws=2/ITER=5 data1

The results for five iterations are shown in the table:

Jacknife Example

We can use the Jacknife test to analyze the standard deviation of the same data. For this application our user function can be written in the form:

Function jackFunc(inWave)
Wave inWave
WaveStats/Q inWave
return V_sdev
End

To run the Jacknife test execute:

StatsResample/jckn=jackFunc data1

The results are displayed in the table:

It is interesting to observe the differences between the standard deviation obtained by the bootstrap approach (above) to the oen obtained here.