Bootstrapping Reclassified Layers of Two Species in QGIS for Suitable Area Overlap

I have a very surface-level knowledge and familiarity with QGIS.

I currently have 30 reclassified QGIS layers of suitable area based off of 30 MaxEnt individual runs for 2 different species. All layers are reclassified to binary, with a 1 being suitable and a 0 being unsuitable area. The binary values are based on a logistic threshold value, and every layer has its own threshold value.

I am trying to find a way to use bootstrapping to get a distribution of suitable area overlap between the two species, with these reclassified layers as input. For example, having a random selection of any of the 30 reclassified layers of species 1 being paired with any of the 30 reclassified layers of species 2, then (because of the bootstrap) this would be repeated in order to generate multiple estimates of area of overlap (an area generated for each combination of the bootstrap). I’m looking to randomly draw one MaxEnt run for one species (that has now been reclassified in QGIS) and find the area of overlap with a MaxEnt run from the other species (also now reclassified), and then do that a bunch of times: the bootstrap. I would then use all of those results from the bootstrap to generate a distribution curve.

I have not seen any way to do this in QGIS (or bootstrap at all).

Do you have any ideas or help for me?

I am using QGIS 3.16.1-Hannover on MacOS.