How can I get p-values for “I” or “D” similarity statistic using nicheOverlap function of dismo R package?

Question

I am working with ensemble rasters of different species of butterflies in R and I want to compare the pattern of maps. Basically, I want to verify how similar the distributions