I am using the indicspecies R package to run an indicator species analysis on a large dataset to see what species are indicators for 3 different regions over the years 2001 to 2021. My matrix is ~550 columns (one for each species in the dataset) and has 21 rows for each year. I am currently running the analysis at nperm=999 and it has been running for over 24 hours.