Data from: Field evaluation of Nzi and modified Nzi traps for sampling of Stomoxys calcitrans at dairy farms in Uruguay
Data files
Mar 26, 2026 version files 27.50 KB
-
README.md
5.14 KB
-
statistic_code.txt
5.17 KB
-
trap_analisis.xlsx
17.19 KB
Abstract
Stomoxys calcitrans (Linnaeus, 1758), the stable fly, is a hematophagous pest that affects animal livestock worldwide by lowering productivity and mechanically transmitting pathogens. This study evaluated the effectiveness of the Nzi trap and a modified version (with and without bait) for capturing stable flies on two dairy farms in southwestern Uruguay, under conditions of natural infestation by these flies. The study used a 4×4 Latin square design with five replicates conducted independently at each dairy farm. Traps at each dairy were emptied daily to identify and count stable flies before being rotated and reset at the next planned location. Among all configurations, standard Nzi traps captured the highest number of S. calcitrans compared to the other trap evaluated, regardless of whether an attractant was used. This study represents the first evaluation in South America of Nzi traps for S. calcitrans used in combination with a commercial bait. Findings indicate that Nzi traps are a potentially effective tool for capturing S. calcitrans at dairy farms.
Dataset DOI: 10.5061/dryad.xsj3tx9vs
Description of the data and file structure
This dataset was collected to evaluate the capture efficiency of two types of traps—Nzi and a locally developed "modified Nzi" trap—for sampling the stable fly (Stomoxys calcitrans) in dairy farms in Uruguay. The study was conducted at two different farms (Farm 1 and Farm 2), where traps were rotated daily among four specific positions at varying distances from the dairy facilities. The experimental efforts included testing the effect of attractants (with/without), trap design, and the influence of spatial distance on fly abundance. Data collection was carried out daily over several weeks in early 2023, recording counts of S. calcitrans, horse flies (Tabanidae), and other dipterans.
Files and variables
File: statistic_code.txt
Description: R script used for data processing and statistical analysis. It includes the code for model selection (Poisson vs. Negative Binomial), Generalized Linear Mixed Models (GLMM), pairwise comparisons (post-hoc tests), and the analysis of the effect of distance on fly abundance.
Note: This file contains all the packages used and the statistical processes executed. The script follows the logic of the Navure software interface for execution in R.
File: trap_analisis.xlsx
Description: Main dataset containing the daily capture records of Stomoxys calcitrans and other dipterans. The data includes experimental variables such as trap design, use of attractants, farm location, and spatial positioning.
Variables:
- date: Date of sampling (YYYY-MM-DD).
- farm: Identification of the dairy farm (1 = Farm 1; 2 = Farm 2).
- week: Consecutive week of the study (1 to 5).
- day: Consecutive day within the sampling period.
- position: Specific location ID (1, 2, 3, or 4) where the trap was placed.
- trap: Categorical variable combining trap type and attractant (e.g., NZImod_atray).
- Trap_1: Trap design type ("Nzi" for standard; "Nzimod" for modified design).
- Attract: Chemical attractant status ("with" or "No").
- Stomoxys_calcitrans: Number of stable flies captured (Count/Unit).
- horse_fly: Number of horse flies (Tabanidae) captured (Count/Unit).
- other_dipteous: Number of other non-target dipterans captured (Count/Unit).
File: Supplementary_Parodi_cv.docx
Description: Supplementary technical documentation. It provides:
- Technical Diagrams: Measurements and components of the modified Nzi trap.
- Fabric Analysis: Reflectance spectra (FORS) of the blue vinyl used compared to standard phthalogen blue cotton.
- Site Maps: Spatial distribution of traps in Farm 1 and Farm 2 relative to the dairy facilities
Code/software
The data analysis for this study was performed using R Statistical Software (version 4.0) through the NAVURE interface (2023). The analysis follows a workflow designed to compare trap efficiency and spatial effects on fly abundance using Generalized Linear Mixed Models (GLMM).
Required R Packages: To execute the provided script (statistic_code.txtThe following libraries must be installed and loaded:
- Model Fitting:
lme4(for GLMM),MASS(for Negative Binomial models). - Post-hoc Analysis:
emmeans,multcomp,multcompView, andcar(for Type III ANOVA tables). - Data Manipulation & Visualization:
tidyverse(includingdplyr,ggplot2),readxl, andstats.
Workflow and File Relationship:
- Data Loading: The script imports the raw data from
trap_analisis.csv(originally an .xlsx file). - Preprocessing: Variables such as Trap_1, Attract, farm, week, day, and position are converted into factors for statistical modeling.
- Model Selection: The script compares a Poisson model against a Negative Binomial (NB) model using AIC/BIC criteria to account for overdispersion in fly counts.
- Statistical Inference:
- The final NB model evaluates the effects of trap design and attractants.
- Pairwise comparisons and significance letters (Compact Letter Display) are generated for interpretation.
- A secondary analysis incorporates spatial distances (in meters) to evaluate the interaction between distance and sampling date.
- Model Stability: A simplified model is compared against the full interaction model to ensure result stability as per peer-review recommendations.
Access information
Other publicly accessible locations of the data:
- None. This is the primary and original repository for these data. Any subsequent publications resulting from this dataset will refer to this Dryad DOI.
Data was derived from the following sources:
- None. The data are original and were collected firsthand by the authors through field experiments conducted at dairy farms in Uruguay. No third-party datasets or external sources were used to derive this information.
This study was conducted during the southern hemisphere summer (January and February 2023) at two dairy farms 3.1 km apart (Farm 1: -34.369338º S -57.879860º W, Farm 2: -34.345090º S -57.895231º W, Colonia, Uruguay, Supplementary Figure S1). The sampling period was based on previous work by Pérez (2021), which showed a seasonal peak between December and February. These dairies were small-scale farms, with 60 (Farm 1) and 80 (Farm 2) lactating Holstein cows. In previous years, farm owners reported severe stable fly infestations, which led to decreased weight gain, lower milk yields, and skin lesions in their cattle.
Traps, bait, and experimental design
Two traps were evaluated: a typical Nzi trap (Mihok, 2002) and a locally assembled modified version (dataset D1 and Figure S2 in the supplementary material). Both were constructed with a local blue vinyl fabric; the modified trap had the same true surface area as the Nzi but incorporated three side entrances. Each trap was tested with or without a new commercial bait for S. calcitrans (Adios Moscas- Stomoxys calcitransTM, Products P651-Lure, Batch: 175.731, Chemtica International S.A, Costa Rica). It contains two ingredients in a plastic bag: 1) 1-octen-3-ol (formulated as a solid tablet-like matrix) and 2) a proprietary mixture of other chemicals in powder form (not yet disclosed). The bait was reconstituted in water at room temperature (300-350 mL/per container). The attractant package was affixed to the upper section of each trap.
The experimental design followed a 4×4 Latin square arrangement with five replicates conducted independently at each farm. Four different trap types were evaluated: the standard Nzi trap, a modified Nzi trap, and both versions equipped with bait. At each farm, four fixed sampling points were established. Sampling points were located near the dairy, with sites separated by more than 50 meters. Traps were rotated daily among the four points in such a way that each trap occupied each position once during a 4-day cycle. This rotation scheme was repeated five times (i.e., five full Latin square replicates), ensuring that each trap was tested at each location a total of five times. The same protocol was implemented identically at both farms, allowing for replication and comparison under different environmental conditions. At the end of the experiment, a total of 40 collections were obtained for each trap type (20 at each farm), resulting in 160 trap collections overall (4 trap types × 40 collections per trap).
Species identification
Identification of stable flies was carried out using a stereomicroscope (Zeiss, model Stemi 508, Jena, Germany) according to Zumpt (1973). For each capture, we classified and quantified the specimens.
Characterization of the blue fabric
All traps were made from an opaque, shiny vinyl as a local alternative to phthalogen blue IF3GM cotton. A FORS (Fiber Optic Reflectance Spectroscopy) spectrum for this blue vinyl indicated that the main colorant was similar but not identical to the spectrum of copper phthalocyanine (CuPc) or Pigment Blue 15 (Poldi and Caglio, 2013), the exception being a lower peak in the blue at 450 nm. Spectra were taken with an Ocean Optics USB-4000 spectrophotometer using a deuterium-halogen light source (Florida, USA) and a Spectralon WS-1-SL white standard (Labsphere, USA). Full details are provided in supplementary material D2 and Figure S3.
Statistical analysis
All analyses were performed using Generalized Linear Mixed Models with a negative binomial distribution using Navure software (NAVURE, 2023). To select the best model, data normality and AIC/BIC fits were compared. Trap performance was evaluated based on the number of flies captured, with trap type (Nzi and modNzi), bait (with or without), and their interaction included as fixed factors, while farm, trap position, day, and week were included as random effects. Differences were considered statistically significant at p < 0.05. To address the distance hypothesis, the exact distance (m) from the dairy was included as a fixed quantitative factor in this model. To assess the stability of capture rates over time, an interaction term between distance and date was tested using a likelihood ratio (LR) test.
