# Statistical signals of copying are robust to time- and space-averaging

Mason Youngblood, Helena Miton, Olivier Morin

► To cite this version:

Mason Youngblood, Helena Miton, Olivier Morin. Statistical signals of copying are robust to time- and space-averaging. Evolutionary Human Sciences, 2023, 5, pp.e10. 10.1017/ehs.2023.5. hal-04307756

HAL Id: hal-04307756

https://hal.science/hal-04307756v1

Submitted on 26 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

Statistical signals of copying are robust to time- and space-averaging

Mason Youngblood¹* ID, Helena Miton¹,² and Olivier Morin¹,³

¹Minds and Traditions Group, Max Planck Institute for Geoanthropology, Jena, Germany, ²Santa Fe Institute, Santa Fe, New Mexico, USA and ³Institut Jean Nicod, ENS, EHESS, PSL University, CNRS, Paris, France

*Corresponding author. E-mail: masonyoungblood@gmail.com

(Received 29 August 2022; revised 23 March 2023; accepted 4 April 2023)

**Abstract**

Cattle brands (ownership marks left on animals) are subject to forces influencing other graphic codes: the copying of constituent parts, pressure for distinctiveness and pressure for complexity. The historical record of cattle brands in some US states is complete owing to legal registration, providing a unique opportunity to assess how sampling processes leading to time- and space-averaging influence our ability to make inferences from limited datasets in fields like archaeology. In this preregistered study, we used a dataset of ~81,000 Kansas cattle brands (1990–2016) to explore two aspects: (1) the relative influence of copying, pressure for distinctiveness and pressure for complexity on the creation and diffusion of brand components; and (2) the effects of time- and space-averaging on statistical signals. By conducting generative inference with an agent-based model, we found that the patterns in our data are consistent with copying and pressure for intermediate complexity. In addition, by comparing mixed and structured datasets, we found that these statistical signals of copying are robust to, and possibly boosted by, time- and space-averaging.

**Keywords:** copying; complexity; time-averaging; space-averaging; generative inference

**Social media summary:** Cattle brands provide clues about how time- and space-averaging might affect statistical signals in archaeological data.

**1. Introduction**

Cattle branding, or the use of hot irons to leave marks on animals, is a common practice that dates back to at least ~1500 BCE (Insoll et al., 2015; van der Moezel, 2016; Wolfenstine, 1970). Cattle brands are emblems – simple graphic codes that do not encode language (Morin et al., 2020). In North America, cattle brands were introduced by the Spanish in the sixteenth century (Brand, 1961) and are applied to one part of the animals' body to denote ownership and deter theft. Today, cattle branding in most parts of the USA is regulated by state organisations that register new brands and ensure that they are unique. Cattle brands in the USA are usually composed of two to four components – constituent parts that include letters, numbers and symbols (Figure 1). Letters often correspond to initials or ranch names, but sometimes the reverse is true and ranches (and even towns) are named after cattle brands (Lombard & Du Plessis, 2016, 2019). The symbols in brands are more arbitrary and include everything from punctuation marks to state outlines and animals.

In any symbolic system used for communication, including simple emblems like cattle brands, there needs to be a way to distinguish symbols from one another. If more distinctive symbols are less likely

© The Author(s), 2023. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Examples of two cattle brands registered in the US state of Kansas in 1884 (https://www.kansasmemory.org/item/309140/). Both brands are composed of two components: on the left a stylised ‘N’ and ‘P’, and on the right a rotated ‘(’ and ‘S’).</figcaption></figure>

to be confused with other symbols, then pressure for distinctiveness should cause symbolic systems to
become more effective over time. In the evolution of writing and speech, for example, pressure for
distinctiveness is thought to have been fundamental in maximising the signal space for better expres-
sivity (de Boer, 2005, 2015; Flemming, 2002; Liljencrants & Lindblom, 1972; Smith & Kirby, 2008).
Distinctiveness is probably an important consideration when ranchers create new cattle brands,
since differentiating animals from one another is their primary function (Gracy II, 2016).
Distinctiveness is interesting to us chiefly because it cannot be maximised without countering two
important mechanisms that also shape emblems: copying and pressure for reduced complexity.

When graphic elements are copied the resulting emblems are less distinctive than they would have
been if they were independently created. Successful graphic forms, when they spread to new users, tend
to become less informative as a result. The designs that ancient Greek city-states minted on their coins
were often adopted by entirely different cities. When this happened, the amount of information that
coins carried about the cities that minted them declined: the same emblem became a poorer indicator
of its city of origin (Pavlek et al., 2019). Likewise, the copying of heraldic coats of arms in medieval
Europe decreased the distinctiveness of copied arms (Morin & Miton, 2018).

Pressure for reduced complexity, which occurs because simpler symbols are easier to perceive, learn and produce (Miton & Morin, 2021; Smith, 2020), can also reduce the distinctiveness of emblems. One possible definition of complexity is the number of ways that a shape could vary from its actual form (van der Helm, 2014). Intuitively, there are relatively few ways to transform the letter ‘O’ into another letter without adding to it, whereas a character like ☰ has many more degrees of freedom. Complex shapes are more likely to be distinctive because they can vary along a greater number of dimensions (Miton & Morin, 2021). This notion predicts (among other things) that written letters should be more complex in larger scripts (e.g. Devanagari vs. Greek), a phenomenon that has been shown using a variety of metrics for visual complexity (Chang et al., 2018; Miton & Morin, 2021). Graphic codes generally evolve to have reduced complexity over time (Kelly et al., 2021; Tamariz & Kirby, 2015), although this tendency is not universal (Miton & Morin, 2019; Tran et al., 2021).

1.1. First aim: distinctiveness, copying and complexity

Variation in cattle brand designs is also likely to be shaped by pressure for distinctiveness, copying and complexity, but the relative importance of these forces is unknown. We make a distinction between distinctiveness (outcome at the level of brands) and *pressure for distinctiveness* (process at the level of components). In our model, pressure for distinctiveness occurs when the rarest components have the highest probability of being included in brand designs. This component-level process is one of several factors (including complexity, as described above) that affect the resulting distinctiveness of brands.

State regulations forbid the exact copying of another rancher's brand: each ranch must register a
brand that is unique in one way or another – meaning that its component symbols, the way they
are rotated or arranged, their location on the animal, or all of the above, are unique.<fnref n="1" /> However, they do not forbid the copying of components or even sets of components, as long as the result is not an identical replication of another brand. Ranchers may copy components out of convenience from published brand books, use common components because of standardised design or reading conventions (e.g. from branding manuals), or modify existing brands with new components (e.g. as observed among pastoralist groups and cattle thieves) (Arnold, 1965; Landais, 2001; Thomas, 1967; Wolfenstine, 1970). Among East African pastoralists, the copying of brands within male lineages is thought to explain the high level of sharing between groups, and possibly even the occurrence of contemporary brands in rock art in the area (Lynch & Robbins, 1977). Other aspects of cattle branding, such as the branding rituals themselves, are also known to be copied (Bird Rondeau, 2013).

As for complexity, operationalised here as the number of components in a brand, the peculiar technology of cattle branding imposes constraints of its own. In the context of cattle branding, designs that are too complex may be more difficult to physically blacksmith, more difficult to read from a distance, and more painful for the animal (Newman, 2007; Stamp, 2013). On the other hand, pressure for reduced complexity in cattle brands is probably limited by the fact that simpler brands are more easily modified by cattle thieves (Adams, 1906; Gracy II, 2016). Cattle theft remains a significant problem (Keen, 2013), and while other technologies such as radio-frequency identification and ear-tagging are less painful for the animals, they are also less effective in deterring theft (American Association of Bovine Practitioners, 2020).

The first aim of this study was to infer the relative influence of pressure for distinctiveness, copying and complexity in the cultural evolution of cattle brands in the US state of Kansas between 1990 and 2016. To do this, we conducted generative inference by fitting an agent-based model (ABM) of cattle brand invention, in which simulated ranchers create new brands under varying levels of copying or pressure for distinctiveness, and pressure for complexity, to the spatiotemporal dynamics of brands and their constituent components. We used a simplified version of the model that we originally pre-registered and made no specific hypotheses. The original preregistered analysis and hypotheses can be found in the first version of the preprint,<fnref n="2" /> and a justification for departure from preregistration is included in the SI.

## 1.2. Second aim: time- and space-averaging

The second aim of this study was to use the high temporal and spatial resolution of our data to assess the influence of time- and space-averaging on the predictability of data. Cultural transmission processes are most easily differentiated from one another by their spatiotemporal dynamics (Barrett, 2019; Wilder & Kandler, 2015). Time- and space-averaging, or the mixing of temporal or spatial structure in data (Perreault, 2019), disrupts such dynamics and makes it much more difficult to accurately infer processes from patterns in data (Garvey, 2018). Averaging, also referred to as mixing, is usually an organic result of how artefacts are deposited in the archaeological record, but it is also sometimes imposed by researchers for analytical convenience (Lyman, 2003; Perreault, 2018). In some cases, time averaging may dilute evidence for other processes and increase the fit of data to neutral models (Brantingham, 2003), as empirically observed in fossil assemblages (Tomašových & Kidwell, 2010). However, simulation-based evidence suggests that time-averaging can sometimes make it more difficult to correctly detect random copying depending on the methods used (Madsen, 2012; Porčić, 2015; Premo, 2014). The effects of space-averaging are less well understood, but there is some empirical evidence that time- and space-averaging affects statistical signals in similar ways (Lyman, 2003).

We investigated the effect of time- and space-averaging on signals of copying at the level of brand components (i.e. subunits). In cultural systems where variants are composed of subunits, copying should leave statistical signals analogous to linkage disequilibrium (i.e. genetic alleles being non-randomly associated; Lewontin, 1995) or to pointwise mutual information (i.e. one item carrying

<footnote n="1">¹https://agriculture.ks.gov/divisions-programs/division-of-animal-health/brands-program</footnote>

<footnote n="2">²https://psyarxiv.com/h5t7b/</footnote>
information about the presence of another; Shannon, 1951). For example, if components are some-
times copied wholesale, or together in groups, then the co-occurrence of components in brands
will be higher than expected by chance (Morin & Miton, 2018). The shuffling model is a modified
derivation of pointwise mutual information that predicts the prevalence of variants assuming that
their subunits contain no information about one another (Morin & Miton, 2018). In other words,
the shuffling model predicts what the prevalence of brands would be if all components in the popu-
lation were randomly copied proportional to their current frequencies, with no wholesale copying. In
their study of heraldic coats of arms, Morin and Miton found that their shuffling model managed to
predict the popularity of coats of arms to a substantial degree. They also found that the presence of
heraldic coats of arms is better predicted by a random copying model when it is applied to space-
averaged data from all of Europe rather than from individual provinces (Morin & Miton, 2018). In
this case, space-mixing might disrupt local clusters of design elements that result from wholesale copy-
ing of coats of arms within families, thus improving the predictions of a model that assumes random
copying of individual design elements.

Do time- and space-mixing boost signals of random copying? To answer this question, we com-
pared the performance of an adapted version of Morin and Miton's shuffling model (Morin &
Miton, 2018) on structured, time-mixed and space-mixed datasets. For each of these datasets, we
first used the shuffling model to calculate the predicted prevalence of every possible combination of
components. Then, we used linear mixed models to compare the accuracy of these predictions after
accounting for factors like complexity. We hypothesised that the shuffling model would perform better
on the time- and space-mixed datasets compared with structured datasets. The structured datasets
built for this goal were also used to test whether brands that are closer to each other in time tend
to be more similar.

## 2. Methods

### 2.1 Data

The data for this study come from the Brands Program of the Kansas Department of Agriculture, which has periodically published books containing all registered cattle brands in the state since 1941.<fnref n="3" /> In the early 1990s Kansas' Brands Program began to use a 13-digit coding system to catalogue all of their cattle brands, which has been used in the seven most recent brand books: 1990, 1998, 2003, 2008, 2014, 2015 and 2016. Kansas' coding system includes information about brands' components, angles of rotation and location on the body. Of the brand books that include this coding system, only the latter four are available in original PDF formats that can be accurately transcribed with optical character recognition (OCR). The three earlier books are only available as scans of the printed books, and thus have noise and artefacts that interfere with OCR. We used a custom neural network in Tesseract OCR to automatically extract the cattle brand codes and zip codes from 2008, 2014, 2015 and 2016, and manually transcribed the brand book from 1990 to maximise the temporal depth of our data. 1998 and 2003 were excluded owing to time and labour constraints. The estimated error rate in transcription, based on a manual check of 1500 brands, was 0% for brand codes and 0.4% for locations. Details about the automated and manual transcription can be seen in the data GitHub.<fnref n="4" />

Each brand code has 13 digits, composed of four three-digit component codes and one single-digit
location code. The three-digit component codes contain the abbreviation for the component and a
number denoting its angle of rotation. Kansas' coding system has some redundancies, such askre and /'kre having separate abbreviations despite being rotated versions of the same component. We collapsed and converted redundant abbreviations into single categories with corrected angles of rotation.
Kansas' Brands Program only allows duplicated cattle brands if they are located on different parts of
the animal. Duplicated brand codes (as opposed to the cattle brands themselves) can also occur when

<footnote n="3">³https://cdm16884/contentdm.oclc.org/digital/collection/p16884coll5/id/317/rec/2</footnote>

<footnote n="4">⁴https://github.com/masonyoungblood/cattle_brand_data</footnote>
the same components are combined in a different arrangement. We chose to remove duplicate brand codes from within the same zip code, which usually occur when a family or ranch has registered a single brand multiple times for different locations, but we chose to keep duplicate brand codes from different zip codes, which usually correspond to different arrangements of the same components. In total, we have 81,063 brand codes from 1990 to 2016, comprised 103 unique components (details in the code GitHub<fnref n="5" />).

## 2.2. Generative inference (first aim)

The ABM simulates ranchers creating new brands every year based on the components present in the brands used by ranchers in the entire state, as well as the dual constraint of simplicity/complexity. Prior simulations showed that including angles of rotation in the model yielded results incompatible with the real data, so we chose to exclude them from our analysis (Figure S4).

New brands are created by sampling components from existing brands within the entire state. All of the components are compiled into a frequency table that is used for weighted random sampling. The probability $P(x)$ that a rancher uses a particular component $x$ is based on the frequency of $x$ raised to an exponent $C$, normalised against the probability of using other components in the population:

$$P(x) = \frac{F_x^C}{\sum_{y=1}^n F_y^C}$$

Negative values of $C$ cause the rarest components to have the highest probability of adoption, thus leading to more distinctive brands. We operationalise this situation ($C < 1$), where the probability of adoption is negatively related to frequency, as pressure for distinctiveness. $C = 0$ produces completely random brands independent of component frequencies. Positive values of $C$ cause the common components to have the highest probability of adoption and introduces redundancy in brand designs. We operationalise this situation ($C > 1$), where the probability of adoption is positively correlated with frequency, as copying. A value of $C = 1$ leads to random (i.e. unbiased) copying, where components are used proportional to their current frequencies. The parameterisation of $C$ in this model is directly based on standard models of anticonformity and conformity in cultural transmission (Crema et al., 2014, 2016; Lachlan et al., 2018; Youngblood et al., 2021; Youngblood & Lahti, 2022), and a full comparison can be seen in Figure S3. To simulate pressure for simplicity/complexity, the number of components in each new brand is drawn from a Poisson distribution where $\lambda$ is the parameter of interest (henceforth complexity; see Figure S1 for examples). At lower values of $\lambda$ (e.g. 0.5) one-components brands are much more common, and at higher values of $\lambda$ (e.g. 15) four-components brands are much more common. An intermediate value of $\lambda$ (e.g. 3) leads to an equal probability of two- and three-component brands and a low probability of one- and four-component brands. No new components were added to Kansas' standardised set during our study period, so we do not simulate the innovation of new component types.

At the beginning of each year in the ABM a set of $N_{new}$ brands is created, where $N_{new}$ is the average number of new brands that appear each year in the observed data. Each new brand is assigned a zip code, where the probability of each zip code is proportional to its frequency in all years of the observed data. After new brands are created a set of $N_{old}$ random brands is removed, where $N_{old}$ is the average number of brands that disappear each year in the observed data.

The ABM is initialised with the data from 1990 and runs through every year until 2016. In order to fit the parameters of the ABM to the observed data we calculated the following set of summary statistics from 2008, 2014, 2015 and 2016:

1. proportion of components that are the most common type;

2. proportion of components that are the most rare type;

<footnote n="5">⁵https://github.com/masonyoungblood/CattleBrandABM</footnote>
3. Shannon's diversity index of components;

4. Simpson's diversity index of components;

5. Jaccard index of components (zip codes);

6. Morisita-Horn index of components (zip codes);

7. Jaccard index of components (counties);

8. Morisita-Horn index of components (counties);

9. mean Levenshtein distance between brands (from random 10%).

Note that these summary statistics are calculated using all of the brands that are present in a given year (e.g. 2014 includes all brands created before and during 2014 that have not yet been removed from the population). This is process is conducted separately for the real and simulated data so that the two can be compared. These summary statistics were chosen because they capture the overall diversity of components with an emphasis on both common (1 and 4) and rare (2 and 3) variants, the spatial diversity of components at the level of both zip codes (5 and 6) and counties (7 and 8), and a proxy for the average perceptual distance between brands (9). For all diversity metrics we calculated their Hill number counterparts, because they are measured on the same scale (Roswell et al., 2021) and better account for relative abundance (Chao et al., 2014). Shannon's diversity index emphasises rarer types whereas Simpson's diversity index emphasises more common types. The Jaccard and Morisita-Horn indices were similarly chosen for their complementarity. The Morisita-Horn index is a commonly used abundance-based beta diversity index (Chao et al., 2006), whereas the Jaccard index is the most robust of the incidence-based beta diversity indices to sampling error (Schroeder & Jenkins, 2018). We calculated beta diversity at both the zip code and county level to assess spatial diversity at two different resolutions. The mean Levenshtein distance (a measure of edit distance), or the minimum number of insertions, deletions and substitutions required to convert one sequence to another, was calculated from a random 10% of brands.

The same summary statistics were calculated from the real data from 2008, 2014, 2015 and 2016. Then, the random forest version of approximate Bayesian computation (ABC) was conducted using the *abcrf* package in R (Raynal et al., 2019), with the following steps:

1. 500,000 iterations of the ABM were run to generate simulated summary statistics for different values of $\lambda$ and C.

2. The output of these simulations were combined into a reference table with the simulated summary statistics as predictor variables, and the parameter values as outcome variables. $\lambda$ was log-transformed because it is non-negative, as recommended for non-linear regression-based ABC (Blum & François, 2010; Sisson et al., 2018).

3. A random forest of 1000 regression trees was constructed for each parameter using bootstrap samples of 80% of the data from the reference table. Random forest parameters were tuned using a random 10% of the data with the *tuneRanger* package in R ($\lambda$, split features = 6 and minimum node size = 3; C, split features = 17 and minimum node size = 3) (Probst et al., 2018).

4. Each trained random forest was provided with the observed summary statistics, and each regression tree was used to generate posterior distributions for both parameters.

Calculating all nine summary statistics for each of the four years means that we conducted ABC with 36 summary statistics in total. Random forest ABC is robust to the number of summary statistics (Pudlo et al., 2016) so we did not have to reduce the dimensionality prior to inference (Blum et al., 2013). As an additional robustness check, we fit the agent-based model to simulated data from two combinations of known parameter values: (1) $\lambda = 3$ and $C = -1$ (intermediate complexity and pressure for distinctiveness); and (2) $\lambda = 3$ and $C = 1$ (intermediate complexity and random copying). In both cases, the random forest ABC recovered the correct input parameters (Figures S6 and S7). Visualisations of model output for the SI were conducted with principal component analysis and uniform manifold approximation and projection (McInnes et al., 2018).
C was assigned a normal prior distribution with a mean of 0 and standard deviation of 2. λ was assigned a gamma prior distribution with a shape parameter of 0.9 and a rate parameter of 0.2. This distribution was chosen because, based on 200,000 samples, it includes a relatively even spread across the component probabilities shown in Figure S1 while allowing for some extreme values (minimum of ~0, maximum of 75).

## 2.3. Shuffling model (second aim)

The shuffling model takes a set of brands and calculates a predicted prevalence for every possible combination of the components within those brands. When it performs well, the shuffling model assigns a higher predicted prevalence to brands that actually exist in the data. We predict that this will be more true for time- and space-mixed datasets than for the original, structured datasets. According to the shuffling model, the predicted prevalence of a brand is simply the product of the frequency of its components in a specific subset of data. For example, the predicted prevalence (S) of a brand with three components A, B and C would be:

$$S_{ABC} = F_A \times F_B \times F_C \times N_{other}$$

where $F_A$ is the frequency of A among brands that do not include B or C, $F_B$ is the frequency of B among brands that do not include A or C, $F_C$ is the frequency of C among brands that do not include A or B and $N_{other}$ is the total number of brands excluding the focal brand.

To apply this model, we first needed to generate structured and mixed datasets. The structured datasets were split into groups according to time, space and complexity. For time, brands were grouped into two age groups: the ‘old’ group (O) included all brands that appear in the 1990 brand book, and the ‘young’ group (Y) included all brands that appear in later brand books. For space, brands were grouped into four rectangular areas of roughly equal size following county borders: northwest (NW), southwest (SW), northeast (NE) and southeast (SE) (see Figure S2). For complexity, brands were grouped into two categories: two-component and three-component brands. One-component brands were excluded because they lack the combinatorics required for the shuffling model, and four-component brands were excluded because they are relatively rare in the dataset (~2%). This structure was used to split the data into 16 subsets, each of which represents a unique combination of time, space and complexity (e.g. O-NE-2 is old northeastern brands with two components).

Once we had structured datasets we generated time- and space-mixed datasets to test our hypotheses about the effects of time- and space- averaging. Each mixed dataset has as many brands as the structured dataset that it imitates. For time-mixed datasets, brands were randomly drawn from the same area independently of time. The time-mixed O-NE-2 dataset was constructed by randomly sampling n two-component brands from northeast Kansas across all time periods, where n is the number of brands in the structured subset. For space-mixed datasets, brands were randomly drawn from the same time dataset independently of space. The space-mixed O-NE-2 dataset was constructed by randomly sampling old two-component brands from all of Kansas.

For each structured and mixed dataset, we inventoried all of the components present in the brands, listed all of their possible combinations, and used the shuffling model to assign a predicted prevalence to each combination. For each possible combination we noted whether it was actually present in the dataset or not.

To compare the accuracy of the shuffling model when applied to the structured and mixed datasets, we built two linear mixed effects models with predicted prevalence as the outcome variable: one with the structured and time-mixed datasets, and the other with the structured and space-mixed datasets. Whether the dataset was mixed (MIXED: 0 = N, Y = 1), whether the combination of components was actually present in the data (ACTUAL: 0 = N, Y = 1) and complexity (COMPLEXITY: 2 = 0, 3 = 1) were potential fixed effects. Space, time and the specific combination of components were potential varying intercepts. The shuffling model assigns a predicted prevalence of 0 when a component occurs
multiple times in a target brand or only occurs in a single brand in the dataset, so all combinations
with a predicted prevalence of 0 were dropped from the modelling. The remaining predicted prevalence values fit a lognormal distribution so we log-transformed them prior to modelling. Model choice was conducted using Akaike's information criterion (AIC) calculated from frequentist models fit with the *lme4* package in R (Bates et al., 2015). The best fitting versions of each model were run as Bayesian models in Stan using the *brms* package in R (Bürkner, 2017), using mean field approximation with 5000 samples from 50,000 iterations.

Our first prediction was that ACTUAL = 1 would have a positive effect on the predicted prevalence of a combination of components, which would indicate that the shuffling model works better than chance. Our second prediction was that there would be an informative and significantly positive interaction term ACTUAL*MIXED, or in other words that the predicted prevalence differential between ACTUAL = 1 and ACTUAL = 0 is more important when MIXED = 1. This would indicate that the performance of the shuffling model is higher in mixed datasets. Finally, complexity was included as a control variable because three-component brands have an additional multiplied proportion that probably reduces their prevalence values.

## 2.4. Temporal distance (supplemental analysis)

The structured datasets used for the shuffling model were also used to test whether brands that are closer to each other in time tend to be more similar. First, we pooled together both the old and young brands from each spatial subset (e.g. O-NE-2 + Y-NE-2 = NE-2) and calculated the Levenshtein distance (LD) between each pair of brands. As a reminder, old brands are those that appear in the 1990 brand book and young brands are those that appear in later brand books but not in the 1990 issue. Then, we ran a linear mixed effects model with LD as the outcome variable, whether the brands are both old (OO), both young (YY) or one of each (OY) as the predictor variable, and the identities of the compared brands as random effects (i.e. one intercept for the first brand and one for the second). We predicted that pairs of brands that are OY would have a higher LD relative to OO and YY, or, in other words, that brands from different time periods would be more distinct from one another.

# 3. Results

## 3.1. Generative inference (first aim)

The posterior distributions and estimates of the two parameters can be seen in Figure 2 and Table 1, respectively. The posterior for complexity ($\lambda$) has a median of 1.55 and a relatively tight 95% confidence interval (CI), suggesting that the data are consistent with pressure for intermediate complexity. As a reminder, $\lambda$ controls the shape of the Poisson distribution from which the number of components in new brands is drawn, so a value of 1.55 gives two components a slightly higher probability than three components (Figure S1). The posterior for copying strength (C) has a median of 0.92 with a relatively tight 95% CI, suggesting that the data is consistent with copying rather than pressure for distinctiveness, at a level slightly below random copying (C = 1). Posterior simulations demonstrate that the fitted parameter values produce output close to the observed data (Figure S5).

## 3.2. Shuffling model (second aim)

Brand appears to be the only grouping variable that explains a high level of variance in the data (time-mixed, ICC<sub>brand</sub> = 0.926, ICC<sub>space</sub> = 0.021, ICC<sub>time</sub> = 0.009; space-mixed, ICC<sub>brand</sub> = 0.923, ICC<sub>space</sub> = 0.003, ICC<sub>time</sub> = 0.002), so we included it as a varying intercept in both models. Adding complexity as a control variable improved model fit (time-mixed, ΔAIC = 5120; space-mixed, ΔAIC = 4852). Adding ACTUAL as a fixed effect improved model fit (time-mixed, ΔAIC = 1377; space-mixed: ΔAIC = 2388) and adding ACTUAL*MIXED as an interaction further improved model fit (time-mixed, ΔAIC = 23;
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The posterior distributions for the two parameters of the ABM computed with random forest ABC, plotted against the priors (dotted).</figcaption></figure>

<table><caption>**Table 1.** The median estimates, 95% quantile-based credible intervals, and root mean square log/logit-transformed errors for each parameter. $\lambda$ controls the shape of the Poisson distribution from which the number of components in new brands is drawn, and $C$ is the exponent to which the components frequencies are raised.</caption><thead><tr><th></th><th>M</th><th>95% CI</th><th>RMSLE</th></tr></thead><tbody><tr><td>λ: complexity</td><td>1.55</td><td>[0.56, 1.82]</td><td>0.13</td></tr><tr><td>C: strength</td><td>0.92</td><td>[0.44, 0.96]</td><td>0.25</td></tr></tbody></table>

space-time, $\Delta$AIC = 403). The best fitting models for both the time-mixed and space-mixed datasets had the following specification:

$$ \log(S) \approx \text{Normal}(\mu_i, \sigma) $$

$$ \mu_i = \alpha_{\text{BRAND}[i]} + (\beta_C \times \text{COMPLEXITY}) + (\beta_A \times \text{ACTUAL}) + (\beta_M \times \text{MIXED}) + (\beta_{AM} \times \text{ACTUAL} \times \text{MIXED}) $$

$$ \alpha_j \approx \text{Normal}(0, 1) $$

$$ \beta_{C, A, M, AM} \approx \text{Normal}(0, 0.5) $$

$$ \sigma \approx \text{Exponential}(3) $$

The results of the best fitting models can be seen in Table 2.

As predicted, we found that both the effect of ACTUAL and the interaction between ACTUAL and MIXED are positive, indicating that the shuffling model works better than chance and performs better in both time- and space-mixed datasets. Figure 3 demonstrates the accuracy of the shuffling model's predictions when applied to a structured, time-mixed and space-mixed dataset (i.e. Y-SE-2, which has the highest sample size). The effect is subtle, but when applied to mixed datasets (right panel) the shuffling model has more true and fewer false positives, apparently because more of the possible combinations of common components are represented.

We found that complexity has the expected negative effect on predicted prevalence. The R-hat values are all equal to 1, and the effective sample sizes are all greater than 1000, showing that both models have converged. The results of a frequentist version of the best fitting model can be seen in Table S2.
<table><caption>**Table 2.** The mean estimates and 95% credible intervals for each parameter in the best fitting models for the time-mixed (left) and space-mixed (right) subsets</caption><thead><tr><th rowspan="2"></th><th colspan="2">Time-mixed</th><th colspan="2">Space-mixed</th></tr><tr><th>Estimate</th><th>95% CI</th><th>Estimate</th><th>95% CI</th></tr></thead><tbody><tr><td>α: intercept</td><td>0.75</td><td>[0.75, 0.76]</td><td>0.74</td><td>[0.74, 0.74]</td></tr><tr><td>β<sub>C</sub>: complexity</td><td>-1.28</td><td>[-1.28, -1.28]</td><td>-1.21</td><td>[-1.21, -1.21]</td></tr><tr><td>β<sub>A</sub>: actual</td><td>0.13</td><td>[0.12, 0.13]</td><td>0.16</td><td>[0.15, 0.16]</td></tr><tr><td>β<sub>M</sub>: mixed</td><td>-0.007</td><td>[-0.009, -0.005]</td><td>-0.02</td><td>[-0.02, -0.01]</td></tr><tr><td>β<sub>AM</sub>: actual*mixed</td><td>0.02</td><td>[0.02, 0.03]</td><td>0.05</td><td>[0.04, 0.06]</td></tr><tr><td>σ: standard deviation</td><td>0.29</td><td>[0.29, 0.29]</td><td>0.29</td><td>[0.29, 0.29]</td></tr></tbody></table>

As a follow-up, we repeated our generative inference analysis after randomly mixing the years and locations of the real brands (not preregistered) and found that the detected signals of copying were robust to time- and space-mixing (Figure S8).

### 3.3. Temporal distance (supplementary aim)

The temporal distance model had the following specification:

$$LD \approx \text{Binomial} (3, p_i)$$

$$\text{logit}(p_i) = \alpha_{\text{BRAND}_A[i]} + \delta_{\text{BRAND}_B[i]} + (\beta_{\text{OO}} \times \text{OO}) + (\beta_{\text{YY}} \times \text{YY})$$

$$\alpha_j \approx \text{Normal} (0, 1)$$

$$\delta_k \approx \text{Normal} (0, 1)$$

$$\beta_{\text{OO,YY}} \approx \text{Normal} (0, 0.5)$$

The results of the temporal distance model can be seen in Table 3. As predicted, we found that the effect of compared brands being from the same time period (e.g. OO or YY) on their Levenshtein distance was negative. In other words, we found that brands from different time periods are more distinct from one another. The R-hat values are all equal to 1, and the effective sample sizes are all greater than 1000, showing that the model has converged.

## 4. Discussion

Based on the results of generative inference, the spatiotemporal dynamics of Kansas cattle brands between 1990 and 2016 are consistent with copying of components across the entire state. Ranchers do not maximise the distinctiveness of their brands beyond the basic constraint of uniqueness, but they do appear to create brands of intermediate complexity. In support of our hypothesis, we found that time- and space-averaging improves the ability of the shuffling model, which assumes random copying at the level of components, to predict the presence and absence of brands in our data. This finding is supported by the fact that our generative inference results remain consistent with copying after the years and locations of brands have been shuffled. Finally, brands from the same time period are more similar to one another than brands from different time periods, as would be expected when copying is present.
<figure><img src="image_5.png" /><figcaption>**Figure 3.** The accuracy of the shuffling model's predictions when applied to a structured (left panel), time-mixed and space-mixed (right panel) version of Y-SE-2, the dataset with the highest sample size. Each point above the diagonal is a unique combination of components, and components on both axes are ranked by their commonness in the entire dataset. Green is a true positive (exists and $S \ge 1$), blue is a true negative (does not exist and $S < 1$), orange is a false positive (does not exist and $S \ge 1$), and yellow is a false negative (exists and $S < 1$).</figcaption></figure>

<table><caption>**Table 3.** The mean estimates and 95% credible intervals for each parameter in the best fitting model.</caption><thead><tr><td></td><td>Estimate</td><td>95% CI</td></tr></thead><tbody><tr><td>α<sub>j</sub>: intercept a</td><td>0.60</td><td>[0.59, 0.60]</td></tr><tr><td>δ<sub>k</sub>: intercept b</td><td>1.05</td><td>[1.05, 1.05]</td></tr><tr><td>β<sub>OO</sub>: both old</td><td>-0.20</td><td>[-0.20, -0.19]</td></tr><tr><td>β<sub>YY</sub>: both young</td><td>-0.30</td><td>[-0.30, -0.30]</td></tr></tbody></table>

Our finding that signals of copying are robust to, and possibly boosted by, time- and space-averaging is significant in light of the previous literature on this topic. Several simulation studies have found that time-averaging makes it more difficult to correctly detect random copying using classic neutrality tests like Slatkin's exact test (Madsen, 2012; Porčić, 2015; Premo, 2014). Such tests are very sensitive to evenness, and thus the long tail of distributions (Hill, 1997) and departures from neutrality induced by time-averaging appear to be due to the overrepresentation of rare and fleeting variants (Premo, 2014). Previous simulation studies have assumed that an infinite number of new variants
<figure><img src="image_6.png" /><figcaption>**Figure 4.** (a) A basic simulation of cultural transmission over three timesteps, with a fully-connected population of 30 agents transmitting two variants (A and B) with a conformity bias ($C=2$). A and B are equally represented in the first timestep. The bar plot in the top right shows the frequencies of A and B in $t_3$, and the bar plot in the top left shows the frequencies of A and B in a time-mixed subsample of the same size from $t_1$ to $t_3$. The dashed lines show the expected frequencies under random copying. (b) The simulated component frequency distributions from the cattle brand ABM under several conditions (five iterations each). Orange and blue are the frequency distributions from 2016 when $C=0$ (frequency information does not matter) and $C=10$ (extreme conformity), respectively. Yellow and green are the frequency distributions from time-mixed subsamples of the same size collected across all years when $C=0$ and $C=10$, respectively. The black dashed line is the frequency distribution expected when $C=1$ (random copying).</figcaption></figure>

are possible, though, and time-averaging may not have this effect in systems like cattle brands where variation is bounded because new variants are produced from a finite set of constituent parts. Additionally, Premo (2014) found that if the entire frequency distribution of variants is considered, then signals of random copying appear to be robust to time-averaging (or even slightly boosted by it, even at low sample sizes). This is consistent with our findings, as well as empirical evidence that time-averaging boosts signals of neutral evolution in fossil assemblages (Tomašových & Kidwell, 2010).

Figure 4 provides an example of the logic of averaging boosting signals of random copying, with a basic simulation of cultural transmission over three timesteps. In $t_1$ the two cultural variants A and B are equally represented in the population. In $t_2$ and $t_3$ a conformity bias ($C=2$) causes B to increase in frequency, as shown in the plot of the full population in $t_3$ in the top right. If we only have access to a subsample from $t_1$ to $t_3$, though, the frequency of the rare variant (A) is higher and the frequency of the common variant (B) is lower (see bottom right), bringing them closer to the expectation under random copying (dashed line). If we simulate this with our cattle brand ABM, using the same starting point and temporal dynamics as our main analysis, we see a similar pattern (Figure 4). The orange and blue lines show the simulated frequency distributions of components in 2016 when $C=0$ (frequency information does not matter) and when $C=10$ (extreme conformity). If we take random subsamples of the same size from all years, shown in yellow and green, it brings the frequency distributions closer to the expectation under random copying, shown in black. When $C=0$ the time-mixing boosts the frequencies of common components at the expense of rare ones, and when $C=10$ it boosts the frequencies of rare components at the expense of common ones.

There are several important differences between cattle brands and typical archaeological datasets that should be highlighted. First, the relative time-scale of averaging in our study is much smaller than that typically seen in archaeological assemblages. The effects of time-averaging generally increase with the aggregation window (Perreault, 2018; Wilder & Kandler, 2015), with the onset of effects occurring when the window is about the same as the average lifetime of variants (Madsen, 2012). Cattle brands are often inherited within families and can have relatively long use-lives, and we suspect that the ~25-year window of our study is close to the average lifetime of a cattle brand. Our study shows that time-averaging can have a significant impact on the detectability of copying with a much narrower time window than previously assumed. Second, time-averaging has bigger effects when innovation rates are high (Madsen, 2012), and no new component types were innovated during
the course of our study period. Combined with the relatively small averaging window, we think that this explains the subtlety of the observed effects of time- and space-averaging (see Table 2 and Figure 3). Finally, the time- and space-averaging in our study occurred at the level of brands, while copying and pressure for distinctiveness were modelled at the level of components. In previous studies, models have assumed that cultural variants are copied discretely as types at the same level of analysis at which time- and space-averaging occurs (Madsen, 2012; Porčić, 2015; Premo, 2014). If we use pottery traditions as an example, a ‘type’ is actually a complex combination of different techniques or design elements. Copying probably occurs at the level of these techniques and elements (or combinations of them), whereas time- and space-averaging occur at the level of the artefacts. We think that continuous cultural variation may contain additional information that is affected differently by averaging or missing data (Premo, 2021).

**Supplementary material.** To view supplementary material for this article, please visit https://doi.org/10.1017/ehs.2023.5

**Acknowledgments.** We would like to thank Lilia Kurnosova, Elia Weber and Agnes Dietrich for their assistance with manual coding of cattle brands, and members of our lab, the Minds and Traditions Research Group at the Max Planck Institute for Geoanthropology, for feedback on earlier versions of this manuscript.

**Author contributions.** M.Y. transcribed the data with optical character recognition and manual coding with research assistants. M.Y. designed and conducted the agent-based modelling with feedback from H.M. and O.M. O.M. and M.Y. designed the time- and space-averaging analysis and M.Y. conducted it. M.Y., H.M., and O.M. wrote the preregistration and manuscript.

**Financial support.** This work has received funding from the Max Planck Society and from the ‘Frontiers in Cognition’ EUR grant, ANR-17-EURE-0017 EUR.

**Conflicts of interest.** We declare no conflicts of interest.

**Data and code availability.** The preregistration document is available on OSF (https://osf.io/d4qv3), the data and transcription methodology are available on our data GitHub (https://github.com/masonyoungblood/cattle_brand_data), and the agent-based model is available on our code GitHub (https://github.com/masonyoungblood/CattleBrandABM).

## References

Adams, A. (1906). *Cattle brands: A collection of western camp-fire stories*. https://www.gutenberg.org/ebooks/12281

American Association of Bovine Practitioners (2020). *Summary: AABP Animal Welfare Committee Branding Working Group activities*.

Arnold, O. (1965). *Irons in the fire: Cattle brand lore*. Abelard-Schuman.

Barrett, B. J. (2019). Equifinality in empirical studies of cultural transmission. *Behavioural Processes*, *161*, 129–138. https://doi.org/10.1016/j.beproc.2018.01.011

Bates, D., Mächler, M., Bolker, B., & Walker, S. (2015). Fitting linear mixed-effects models using lme4. *Journal of Statistical Software*, *67*(1), 1–48. https://doi.org/10.18637/jss.v067.i01

Bird Rondeau, J. J. (2013). *Communicating social identity: The sensory performance of cattle branding*. University of Calgary.

Blum, M. G. B., & François, O. (2010). Non-linear regression models for approximate Bayesian computation. *Statistical Computing*, *20*, 63–73. https://doi.org/10.1007/s11222-009-9116-0

Blum, M. G. B., Nunes, M. A., Prangle, D., & Sisson, S. A. (2013). A comparative review of dimension reduction methods in approximate Bayesian computation. *Statistical Science*, *28*(2), 189–208.

Brand, D. D. (1961). The early history of the range cattle industry in northern Mexico. *Agricultural History*, *35*(3), 132–139.

Brantingham, P. J. (2003). A neutral model of stone raw material procurement. *American Antiquity*, *68*(3), 487–509. https://doi.org/10.2307/3557105

Bürkner, P. C. (2017). brms: An R package for Bayesian multilevel models using Stan. *Journal of Statistical Software*, *80*(1). https://doi.org/10.18637/jss.v080.i01

Chang, L. Y., Chen, Y. C., & Perfetti, C. A. (2018). GraphCom: A multidimensional measure of graphic complexity applied to 131 written languages. *Behavior Research Methods*, *50*(1), 427–449. https://doi.org/10.3758/s13428-017-0881-y

Chao, A., Chazdon, R. L., Colwell, R. K., & Shen, T. J. (2006). Abundance-based similarity indices and their estimation when there are unseen species in samples. *Biometrics*, *62*(2), 361–371. https://doi.org/10.1111/j.1541-0420.2005.00489.x

Chao, A., Gotelli, N. J., Hsieh, T. C., Sander, E. L., Ma, K. H., Colwell, R. K., & Ellison, A. M. (2014). Rarefaction and extrapolation with Hill numbers: A framework for sampling and estimation in species diversity studies. *Ecological Monographs*, *84*(1), 45–67. https://doi.org/10.1890/13-0133.1
Crema, E. R., Edinborough, K., Kerig, T., & Shennan, S. J. (2014). An approximate Bayesian computation approach for infer-
ring patterns of cultural evolutionary change. *Journal of Archaeological Science, 50*, 160–170. https://doi.org/10.1016/j.jas.
2014.07.014

Crema, E. R., Kandler, A., & Shennan, S. (2016). Revealing patterns of cultural transmission from frequency data: Equilibrium and non-equilibrium assumptions. *Scientific Reports*, 6(December), 1–10. https://doi.org/10.1038/srep39122

de Boer, B. (2005). Evolution of speech and its acquisition. *Adaptive Behavior, 13*(4), 281–292. https://doi.org/10.1177/
105971230501300405

de Boer, B. (2015). Biology, culture, evolution and the cognitive nature of sound systems. *Journal of Phonetics*, 53, 79–87.
https://doi.org/10.1016/j.wocn.2015.07.001

Flemming, E. S. (2002). *Auditory representations in phonology*. Routledge.

Garvey, R. (2018). Current and potential roles of archaeology in the development of cultural evolutionary theory.
*Philosophical Transactions of the Royal Society B: Biological Sciences*, 373(1743). https://doi.org/10.1098/rstb.2017.0057

Gracy II, D. B. (2016). A cowman’s-eye view of the information ecology of the Texas cattle industry from the Civil War to World War I. *Information & Culture, 51*(2), 164–191. https://doi.org/10.7560/IC51202

Hill, M. O. (1997). An evenness statistic based on the abundance-weighted variance of species proportions. *Oikos*, 79(2), 413–416.
https://doi.org/10.2307/3546027

Insoll, T., Clack, T., & Rege, O. (2015). Mursi ox modification in the Lower Omo Valley and the interpretation of cattle rock art in Ethiopia. *Antiquity*, 89(343), 91–105. https://doi.org/10.15184/aqy.2014.31

Keen, J. (2013, April). Cattle thefts spur new interest in livestock branding. USA Today. https://eu.usatoday.com/story/news/
nation/2013/04/13/cattle-theft-branding/2079323/

Kelly, P., Winters, J., Miton, H., & Morin, O. (2021). The predictable evolution of letter shapes: An emergent script of West Africa recapitulates historical change in writing systems. *Current Anthropology*, 62(6). https://doi.org/10.1086/717779

Lachlan, R. F., Ratmann, O., & Nowicki, S. (2018). Cultural conformity generates extremely stable traditions in bird song.
*Nature Communications*, 9. https://doi.org/10.1038/s41467-018-04728-1

Landais, E. (2001). The marking of livestock in traditional pastoral societies. *Revue Scientifique et Technique*, 20(2), 463–479. http://dx.doi.org/10.20506/rst.20.2.1286

Lewontin, R. C. (1995). The detection of linkage disequilibrium in molecular sequence data. *Genetics*, *140*(1), 377–388.
https://doi.org/10.1093/genetics/140.1.377

Liljencrants, J., & Lindblom, B. (1972). Numerical simulation of vowel quality systems: The role of perceptual contrast.
*Language*, 48(4), 839–862. https://doi.org/10.2307/411991

Lombard, C. G., & Du Plessis, T. (2016). Beyond the branding iron: Cattle brands as heritage place names in the state of Montana. *Names*, 64(4), 224–233. https://doi.org/10.1080/00277738.2016.1223119

Lombard, C. G., & Du Plessis, T. (2019). The socio-onomastic features of American cattle brands. *Names*. https://doi.org/10.1080/00277738.2018.1490517

Lyman, L. R. (2003). The influence of time averaging and space averaging on the application of foraging theory in zooarchaeology. *Journal of Archaeological Science*, 30(5), 595–610. https://doi.org/10.1016/S0305-4403(02)00236-4

Lynch, M., & Robbins, L. H. (1977). Animal brands and the interpretation of rock art in East Africa. *Current Anthropology*, 18
(3), 538–539.

Madsen, M. E. (2012). Unbiased cultural transmission in time-averaged archaeological assemblages. SSRN Electronic Journal.
https://doi.org/10.2139/ssrn.2037622

McInnes, L., Healy, J., & Melville, J. (2018). UMAP: Uniform manifold approximation and projection for dimension reduction. *ArXiv*. https://arxiv.org/abs/1802.03426

Miton, H., & Morin, O. (2019). When iconicity stands in the way of abbreviation: No Zipfian effect for figurative signals. *PLoS ONE*, 14(8), 1–19. https://doi.org/10.1371/journal.pone.0220793

Miton, H., & Morin, O. (2021). Graphic complexity in writing systems. *Cognition*, 214, 104771. https://doi.org/10.1016/j.cognition.2021.104771

Morin, O., Kelly, P., & Winters, J. (2020). Writing, graphic codes, and asynchronous communication. *Topics in Cognitive Science*, 12(2), 727–743. https://doi.org/10.1111/tops.12386

Morin, O., & Miton, H. (2018). Detecting wholesale copying in cultural evolution. *Evolution and Human Behavior*, 39(4),
392–401. https://doi.org/10.1016/j.evolhumbehav.2018.03.004

Newman, R. (2007). A guide to best practice husbandry in beef cattle: Branding, castrating and dehorning. Meat & Livestock Australia.

Pavlek, B., Winters, J., & Morin, O. (2019). Ancient coin designs encoded increasing amounts of economic information over centuries. *Journal of Anthropological Archaeology*, 56, 101103. https://doi.org/10.1016/j.jaa.2019.101103

Perreault, C. (2018). Time-averaging slows down rates of change in the archaeological record. *Journal of Archaeological Method and Theory*, 25(3), 953–964. https://doi.org/10.1007/s10816-018-9364-4

Perreault, C. (2019). *The Quality of the Archaeological Record*. The University of Chicago Press.

Porčić, M. (2015). Exploring the effects of assemblage accumulation on diversity and innovation rate estimates in neutral, conformist, and anti-conformist models of cultural transmission. *Journal of Archaeological Method and Theory*, 22(4), 1071–1092.
Premo, L. S. (2014). Cultural transmission and diversity in time-averaged assemblages. *Current Anthropology*, *55*(1), 105–114.
https://doi.org/10.1086/674873

Premo, L. S. (2021). Population size limits the coefficient of variation in continuous traits affected by proportional copying error (and why this matters for studying cultural transmission). *Journal of Archaeological Method and Theory*, *28*(2), 512–534.
https://doi.org/10.1007/s10816-020-09464-9

Probst, P., Wright, M., & Boulesteix, A.-L. (2018). Hyperparameters and tuning strategies for random forest. *Wiley Interdisciplinary Reviews: Data Mining and Knowledge Discovery*. https://doi.org/10.1002/widm.1301

Pudlo, P., Marin, J. M., Estoup, A., Cornuet, J. M., Gautier, M., & Robert, C. P. (2016). Reliable ABC model choice via random forests. *Bioinformatics*, *32*(6), 859–866.
https://doi.org/10.1093/bioinformatics/btv684

Raynal, L., Marin, J.-M., Pudlo, P., Ribatet, M., Robert, C. P., & Estoup, A. (2019). ABC random forests for Bayesian parameter inference. *Bioinformatics*, *35*(10), 1720–1728.
https://doi.org/10.1093/bioinformatics/bty867

Roswell, M., Dushoff, J., & Winfree, R. (2021). A conceptual guide to measuring species diversity. *Oikos*, *130*(3), 321–338.
https://doi.org/10.1111/oik.07202

Schroeder, P. J., & Jenkins, D. G. (2018). How robust are popular beta diversity indices to sampling error. *Ecosphere*, *9*(2).
https://doi.org/10.1002/ecs2.2100

Shannon, C. E. (1951). Prediction and entropy of printed English. *Bell System Technical Journal*, *30*(1), 50–64.
https://doi.org/10.1002/j.1538-7305.1951.tb01366.x

Sisson, S. A., Fan, Y., & Beaumont, M. A. (2018). *Handbook of approximate Bayesian computation*. CRC Press.
https://doi.org/10.1201/9781315117195

Smith, K. (2020). How culture and biology interact to shape language and the language faculty. *Topics in Cognitive Science*, *12*
(2), 690–712.
https://doi.org/10.1111/tops.12377

Smith, K., & Kirby, S. (2008). Cultural evolution: Implications for understanding the human language faculty and its evolution. *Philosophical Transactions of the Royal Society B*, *363*(1509), 3591–3603.
https://doi.org/10.1098/rstb.2008.0145

Stamp, J. (2013, April). Decoding the range: The secret language of cattle branding. *Smithsonian Magazine*. https://www.smithsonianmag.com/arts-culture/decoding-the-range-the-secret-language-of-cattle-branding-45246620/

Tamariz, M., & Kirby, S. (2015). Culture: Copying, compression, and conventionality. *Cognitive Science*, *39*(1), 171–183.
https://doi.org/10.1111/cogs.12144

Thomas, F. J. (1967). *Mission cattle brands*. Ten Finger Press.

Tomašových, A., & Kidwell, S. M. (2010). The effects of temporal resolution on species turnover and on testing metacommunity models. *American Naturalist*, *175*(5), 587–606.
https://doi.org/10.1086/651661

Tran, N.-H., Waring, T., Atmaca, S., & Beheim, B. A. (2021). Entropy trade-offs in artistic design: A case study of Tamil kolam. *Evolutionary Human Sciences*, *3*(e23), 1–14.
https://doi.org/10.1017/ehs.2021.14

van der Helm, P. A. (2014). *Simplicity in vision: A multidisciplinary account of perceptual organization*. Cambridge University Press.

van der Moezel, K. V. J. (2016). *Of marks and meaning: A palaeographic, semiotic–cognitive, and comparative analysis of the identity marks from Deir El-Medina*. Universiteit Leiden.

Wilder, B., & Kandler, A. (2015). Inference of cultural transmission modes based on incomplete information. *Human Biology*, *87*(3), 193–204.

Wolfenstein, M. R. (1970). *Manual of Brands and marks*. University of Oklahoma Press.

Youngblood, M., & Lahti, D. (2022). Content bias in the cultural evolution of house finch song. *Animal Behaviour*, *185*, 37–48.
https://doi.org/10.1016/j.anbehav.2021.12.012

Youngblood, M., Stubbersfield, J. M., Morin, O., Glassman, R., & Acerbi, A. (2021). Negativity bias in the spread of voter fraud conspiracy theory tweets during the 2020 US election. *PsyArXiv*. https://doi.org/10.31234/osf.io/2jksg

<footnote>Cite this article: Youngblood M, Miton H, Morin O (2023). Statistical signals of copying are robust to time- and space-averaging. *Evolutionary Human Sciences* 5, e10, 1–15. https://doi.org/10.1017/ehs.2023.5</footnote>