sanssouci-org / sanssouci.python

Post hoc inference via multiple testing
GNU General Public License v3.0
6 stars 3 forks source link

Full vectorisation of one- and two-sample tests #34

Closed pneuvial closed 2 years ago

pneuvial commented 2 years ago

Currently the functions get_permuted_p_values and get_permuted_p_values_one_sample perform a loop on the permutation/sign flips.

This operation can be vectorized for classical tests, as done in R: https://github.com/pneuvial/sanssouci/blob/develop/R/rowZTests.R

This should first be done for one-sample tests.

pneuvial commented 2 years ago

Decision: we won't implement this because this leads to storing matrices of size m x B, which is a bad idea if m = 400,000 and B = 10,000