Environmental and Farm Management Factors Associated with Tuberculosis on Cattle Farms in Northeastern Michigan

Objective—To identify major environmental and farm management factors associated with the occurrence of tuberculosis (TB) on cattle farms in northeastern Michigan. Design—Case-control study. Sample Population—17 cattle farms with infected cattle and 51 control farms. Procedure—Each case farm (laboratory confirmed diagnosis of Mycobacterium bovis infection) was matched with 2 to 4 control farms (negative whole-herd test results within previous 12 months) on the basis of type of farm (dairy or beef) and location. Cattle farm data were collected from in-person interviews and mailed questionnaires. Wildlife TB data were gathered through state wildlife surveillance. Environmental data were gathered from a satellite image-based geographic information system. Multivariable conditional logistic regression for matched analysis was performed. Results—Major factors associated with increased farm risk of TB were higher TB prevalence among wild deer and cattle farms in the area, herd size, and ponds or creeks in cattle housing areas. Factors associated with reduced farm risk of TB were greater amounts of natural open lands in the surrounding area and reducing deer access to cattle housing areas by housing cattle in barns, barnyards, or feedlots and use of electrified wire or barbed wire for livestock fencing. Conclusions and Clinical Relevance—Results suggest that certain environmental and management factors may be associated with risk of TB on cattle farms.

Space–time zero-inflated count models of Harbor seals

Environmental data are spatial, temporal, and often come with many zeros. In this paper, we included space–time random effects in zero-inflated Poisson (ZIP) and ‘hurdle’ models to investigate haulout patterns of harbor seals on glacial ice. The data consisted of counts, for 18 dates on a lattice grid of samples, of harbor seals hauled out on glacial ice in Disenchantment Bay, near Yakutat, Alaska. A hurdle model is similar to a ZIP model except it does not mix zeros from the binary and count processes. Both models can be used for zero-inflated data, and we compared space–time ZIP and hurdle models in a Bayesian hierarchical model. Space–time ZIP and hurdle models were constructed by using spatial conditional autoregressive (CAR) models and temporal first-order autoregressive (AR(1)) models as random effects in ZIP and hurdle regression models. We created maps of smoothed predictions for harbor seal counts based on ice density, other covariates, and spatio-temporal random effects. For both models predictions around the edges appeared to be positively biased. The linex loss function is an asymmetric loss function that penalizes overprediction more than underprediction, and we used it to correct for prediction bias to get the best map for space–time ZIP and hurdle models.