This is a general question. I have a global dataset of ~60,000 protected areas (PA) distributed globally. I am modeling change in species richness change_alpha for each protected area between two time periods.