occuRN function

Fit the occupancy model of Royle and Nichols (2003)