occuComm function

Fit a community occupancy model with species-level random effects