occuCOP function

Fit the occupancy model using count data