Breeding males, but not females, of Fukomys mole-rats use stronger bites to defend reproductive monopoly. Fukomys

This dataset supports the analyses and conclusions presented in the above-titled manuscript. It includes in vivo bite force measurements, morphological data, and associated sample metadata from four cooperatively breeding species of African mole-rats (genus Fukomys). These data provide the empirical foundation for testing hypotheses about sex- and status-specific differences in bite force and its morphological correlates in relation to intraspecific competition and reproductive roles.The supplementary materials include two structured datasets and one annotated R script. These resources allow full replication of the analyses and figures presented in the manuscript and ensure compliance with the Nature Portfolio’s data transparency standards.1. Dataset: Dataset_biteforce.txtDescription:This tab-delimited text file contains the primary dataset used to analyze in vivo bite force (BF) in 404 adult individuals of four species of African mole-rats from the genus Fukomys: F. anselli (treated taxonomically as F. micklemi), F. darlingi, F. mechowii, and F. damarensis. Each row corresponds to one individual.Key variables include:Column Name DescriptionIDchip Unique chip or specimen IDanimal_id Animal identifier codespecies Nominal species (e.g., F. anselli, F. darlingi)species1 Taxonomic grouping used in analysis (e.g., F. micklemi includes F. anselli)status Breeding status: B (breeder) or N (non-breeder)sex Sex of the individual: M or Fstatussex Combined breeding status and sex (e.g., BM, BF, NM, NF)weight Body mass in gramsHead_Width Head width (mm), measured between zygomatic archesIncisor_width Upper incisor width (mm), measured at widest pointbfmaxn Maximum in vivo bite force (N)The file also includes binary classification columns used in modeling:FNonbreeders, FBreeders, MNonbreeders, MBreeders, FMBreeders: dummy variables for model specification.2. Dataset: Dataset_headincisorwidth.txtDescription:This tab-delimited file is an extended version of the above dataset. It includes the same variables as Dataset_biteforce_review.txt, but is specifically structured and filtered for detailed analysis of head width and incisor width variation across breeding status and sex.It contains slightly more rows due to inclusion of additional individuals who may have been excluded from bite force analyses (e.g., due to refusal to bite) but were still measured morphologically.Use case:Used to analyze structural traits—head width and upper incisor width—and their scaling with body mass across sex and breeding status categories.3. Script: BiteForce_rscrip.txtDescription:This R script contains the complete analysis workflow used in the manuscript. It is written in R and uses the brms package for Bayesian multilevel modeling, along with ggplot2 and related visualization packages.The script includes:Data loading and preprocessing (e.g., log-transformation, standardization)Calculation of species-specific summaries for bite force, head width, and incisor widthVisualization of trait–body mass relationships across species and groupsBayesian model fitting for:Bite force (BF)Head width (HW)Incisor width (IW)Model comparisons (WAIC and LOO-CV)Posterior predictive checksBack-transformation of standardized estimates to original measurement scalesVisualization of conditional effects and pairwise contrastsSpecies-specific models (e.g., F. micklemi subset analysis)Dependencies:brms, ggplot2, marginaleffects, emmeans, tidybayes, dplyr, patchwork, bayesplot, ggdist, loo, ggpubr, ggeffectsReproducibility Note:The script is organized for full reproducibility. All figures and results reported in the manuscript (Figures 1–4, Tables 1–3, S1–S7) can be reproduced from this script using the provided datasets. The modeling follows a clearly documented structure using well-defined priors and diagnostics.