Let's use a state of the art machine learning library to build an old-fashioned Gabor wavelet pyramid model

Simulating how voxel selection in fMRI caused regression to the mean in a recently-retracted paper