Conditional logistic regression (CLR) is widely used to analyze habitat selection and movement of animals when resource availability changes over space and time. Observations used for these analyses are typically autocorrelated, which biases model-based variance estimation of CLR parameters. This bias can be corrected using generalized estimating equations (GEE), an approach that requires partitioning the data into independent clusters. Here we establish the link between clustering rules in GEE and their effectiveness to remove statistical biases in variance estimation of CLR parameters.The current lack of guidelines is such that broad variation in clustering rules can be found among studies (e.g., 14–450 clusters) with unknown consequences on the robustness of statistical inference. We simulated datasets reflecting conditions typical of field studies. Longitudinal data were generated based on several parameters of habitat selection with varying strength of autocorrelation and some individuals having more observations than others. We then evaluated how changing the number of clusters impacted the effectiveness of variance estimators. Simulations revealed that 30 clusters were sufficient to get unbiased and relatively precise estimates of variance of parameter estimates. The use of destructive sampling to increase the number of independent clusters was successful at removing statistical bias, but only when observations were temporally autocorrelated and the strength of inter-individual heterogeneity was weak. GEE also provided robust estimates of variance for different magnitudes of unbalanced datasets. Our simulations demonstrate that GEE should be estimated by assigning each individual to a cluster when at least 30 animals are followed, or by using destructive sampling for studies with fewer individuals having intermediate level of behavioural plasticity in selection and temporally autocorrelated observations. The simulations provide valuable information to build reliable habitat selection and movement models that allow for robustness of statistical inference without removing excessive amounts of ecological information.

Primary production can determine the outcome of management actions on ecosystem properties, thereby defining sustainable management. Yet human agencies commonly overlook spatio-temporal variations in productivity by recommending fixed resource extraction thresholds. We studied the influence of forest productivity on habitat disturbance levels that boreal caribou – a threatened, late-seral ungulate under top-down control – should be able to withstand. Based on 10 years of boreal caribou monitoring, we found that adult survival and recruitment to populations decreased with landscape disturbance, but increased with forest productivity. This benefit of productivity reflected the net outcome of an increase in resources for apparent competitors and predators of caribou, and a more rapid return to the safety of mature conifer forests. We estimated 3-fold differences in forest harvesting levels that caribou populations could withstand due to variations in forest productivity. The adjustment of ecosystem provisioning services to local forest productivity should provide strong conservation and socio-economic advantages.

No abstract