RF code from test_stats_sp.py into simple_hrf_dataset so we could create "realistic" ones for testing

This is probably not all what should/could be done, but could be a beginning point, so I thought to get it merged.

All of that madness was manually cherry-picked from the content of "not-to-be-merged" https://github.com/PyMVPA/PyMVPA/pull/191