provide better options for multiple test correction

Bonferroni is overly conservative because voxels are not independent

We should look into FDR-based / random field methods (?)

The issue names no files, tests, or entry points. Start by locating the existing Bonferroni correction and the code that processes voxel-level results, then research suitable FDR-based and random-field methods. Done means the project offers a clearly defined, tested alternative to Bonferroni with documented behavior.

Written by the indexing model from the issue text.