In this readme is explained how to run the experiments needed for thesignal-to-noise ratio plot in:H. Valpola, T. 謘tman and J. Karhunen"Nonlinear independent factor analysis by hierarchical models"ICA2003You need the Bayes Blocks library fromhttp://www.cis.hut.fi/projects/bayes/and have the environment variable PYTHONPATH to point to theinstallation for running these experiments.To get a list of possible parameters to use with runhnfa type%  python runhnfa.py --helpTo run one experiment (each experiment takes several hours to run)% mkdir experiments% cd experiments% python ../runhnfa.py -p nl1 -n8 -i10000 -s0 -l22All the experiment used for the signal-to-noise ratio plot are listedin file 'all_runs'.To calculate signal-to-noise ratio to all results in experiments directory:% cd experiments% python ../calculate_corrs.pyThis will save the results to cost_snr.asc. First column is the valueof the cost function and the second is the signal-to-noise ratio. Thefollowing matlab commands makes the signal-to-noise ratio plot.load cost_snr.ascfonts=16;markers=16;plot(cost_snr(:,1)/1000,cost_snr(:,2),'.','MarkerSize',markers);xlabel('Cost / 1000','FontSize',fonts);ylabel('SNR / dB','FontSize',fonts);set(gca,'FontSize',fonts);