The influence of social structure on the propagation of social ...

14 downloads 5286 Views 719KB Size Report
Feb 9, 2008 - The most famous examples for such traditions in animals .... we mapped the social networks as graphs where the vertices represent individuals ...
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution and sharing with colleagues. Other uses, including reproduction and distribution, or selling or licensing copies, or posting to personal, institutional or third party websites are prohibited. In most cases authors are permitted to post their version of the article (e.g. in Word or Tex form) to their personal website or institutional repository. Authors requiring further information regarding Elsevier’s archiving and manuscript policies are encouraged to visit: http://www.elsevier.com/copyright

ARTICLE IN PRESS

Journal of Theoretical Biology 252 (2008) 77–86 www.elsevier.com/locate/yjtbi

The influence of social structure on the propagation of social information in artificial primate groups: A graph-based simulation approach Bernhard Voelkl, Ronald Noe¨ De´partement Ecologie, Physiologie & Ethologie, IPHC, Centre National pour la Recherche Scientifique, Universite´ Louis Pasteur, 23 rue Becquerel, Strasbourg 67087, France Received 16 August 2007; received in revised form 7 December 2007; accepted 2 February 2008 Available online 9 February 2008

Abstract Observations of primate groups have shown that social learning can lead to the development of temporal stable traditions or even proto-culture. The social structure of primate groups is highly diverse and it has been proposed that differences in the group structure shall influence the patterns of social information transmission. While empirical studies have mainly focused on the psychological mechanisms of social learning in individuals, the phenomenon of information propagation within the group has received relatively little attention. This might be due to the fact that formal theories that allow actual testing have not been formulated, or were kept too simple, ignoring the social dynamics of multi-agent societies. We want to propose a network approach to social information transmission that (1) preserves the complexity of the social structure of primate groups and (2) allows direct application to empirical data. Results from simulation experiments with artificial group structures confirm that association patterns of group-members influence the expected speed of information transmission during the propagation process. Introducing a forgetting rate shows that under certain conditions the proportion of informed individuals will reach a stable rate in some systems while it will drop to zero in others. This suggests that the likelihood to observe temporal stable traditions shall differ between social systems with different structure. r 2008 Elsevier Ltd. All rights reserved. Keywords: Information transmission; Social system; Weighted graph; Social network

1. Introduction In the early 1950s, the Japanese primatologist Kinji Imanishi used an imaginary dialogue between an evolutionist, a layman, a monkey and a wasp to demonstrate that culture—defined as the non-genetic transmission of habits—was entirely possible, and even likely, for animals other than humans (de Waal, 1999). This proposition was of course vigorously opposed by those who wanted the term culture to be reserved for humans alone. Whether or not some animal species have a ‘‘full blown’’ culture was, and still is, a matter of fierce debate. However, plenty of evidence has been collected that animals are capable of Corresponding author. Tel.: +33 3 8810 7450; fax: +33 3 8810 6906.

E-mail addresses: [email protected] (B. Voelkl), [email protected] (R. Noe¨). 0022-5193/$ - see front matter r 2008 Elsevier Ltd. All rights reserved. doi:10.1016/j.jtbi.2008.02.002

transmitting information socially and that such information transmission can lead to temporally stable traditions. The most famous examples for such traditions in animals include birdsong dialects (Kroodsma and Baylis, 1982), milk bottle opening in British tits (Hinde and Fisher, 1951), potato washing in Japanese macaques (Kawai, 1965), termite fishing (Goodall, 1968; Lohnsdorf, 2006), nut cracking (Biro et al., 2003; Sugiyama and Koman, 1979) and the use of medical plants (Huffmam and Seifu, 1989; Huffman and Hirata, 2004) in chimpanzees. Experimental studies with guppies demonstrated that fish are able to learn new routes and foraging tasks socially both in the laboratory and in the wild (Laland and vanBergen, 2003; Reader et al., 2003). So far the majority of experimental studies on social information transmission have focused on the psychological mechanisms responsible for social learning (see Heyes

ARTICLE IN PRESS 78

B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

and Galef, 1996; Zentall and Galef, 1988; and chapters therein) and neglected social processes as potential influences on both the likelihood of transmission and the type of information that can be acquired socially (Coussi-Korbel and Fragaszy, 1995). The traditional experimental approach to social learning where individuals are tested in isolation after they have observed a physically separated demonstrator was repeatedly criticised for two reasons. First, it was argued that the understanding of social learning from ecological and comparative perspectives would be aided by experiments on social learning in groups rather than in individuals (Lepoivre and Pallaud, 1985). Second, it was criticised that this type of experiments says little about the diffusion processes of socially learned behaviour through populations (Laland et al., 1996) and about the influence of the social dynamics on the transmission process (Coussi-Korbel and Fragaszy, 1995). Studies of diffusion processes in pigeons have demonstrated that information transmission can be inhibited by the individuals’ propensity to adopt certain roles as either producers or scroungers (Giraldeau and Lefebvre, 1986, 1987; Laland and Plotkin, 1990). In the laboratory transmission, chain studies (Curio et al., 1978; Horner et al., 2006; Laland and Plotkin, 1990) have tried to simulate passage of information transmission along a line of several individuals, by iterating the classical social learning paradigm and using the observer of one test as a demonstrator for the next test. This paradigm was taken a step further by Galef and Allen (1995) who adopted a method by Jacobs and Campbell (1961) for the study of the stability of arbitrary traditions in humans. In their experiment founder colonies of rats were taught an arbitrary food preference. Individual members of the founder colony were then slowly replaced with naı¨ ve subjects and even several generations after the last founder had been removed, the arbitrary food preference was still evident. The same effect could be demonstrated in guppies where preferences for routes to a food source were transmitted over several generations (Laland and Williams, 1997). While these studies have convincingly demonstrated that information can actually spread in a population by means of social learning, they have still not addressed the question how information spreads within natural populations. The lack of studies in this direction might be primarily due to missing predictions how information shall spread in such populations, or that existing theories are kept too general to allow actual testing (Laland and vanBergen, 2003; Reader, 2004). Field observations—mainly of primate groups—revealed a long list of seemingly socially transmitted habits. However, due to the observational character of these studies it is usually not possible to exclude asocial learning processes as candidate explanations for these behaviours. To support the claim that these behaviours were in fact socially transmitted, it was suggested that the spreading pattern of the behaviour itself can reveal something about

the learning process involved (Galef, 1991). Lefebvre (1995) re-analysed the data of 21 cases of presumably culturally transmitted feeding innovations in primate groups. Social learning processes are argued to result in accelerating diffusion curves, such as the logistic, the exponential, or the hyperbolic sine (Boyd and Richerson, 1985; Cavalli-Sforza and Feldman, 1981; Laland et al., 1996) as the number of demonstrators of the behaviour increases over time. In contrast, asocial learning should result in a linear diffusion curve. According to Lefebvre, the primate data supported the assumption that the innovations were transmitted socially. However, recently (Laland and vanBergen, 2003) demonstrated that even asocial learning processes can result in accelerating diffusion curves if one incorporates variation in the individuals’ learning rates—a truly justified assumption— into the model. Reader (2004) summarises over 20 datasets where diffusion curves have been used to distinguish social from asocial learning. He reported that the shapes of the curves were quite variable and often did not fit the predictions sufficiently. As a likely reason for this discrepancy he suggested that the models were so far too simple, neglecting the effects of the population’s social structure. Recently, Kendal et al. (2007) have investigated social learning processes in small groups of callitrichid monkeys. They simulated the effects of three learning processes—stimulus enhancement, observational learning and asocial learning—and selected the modelparameters that provided the best fit to the data. Logarithmic and inverse functions provided the best fit for most diffusion curves. They found no evidence for observational learning; however that might be due to the nature of the foraging task. Socio-spatial parameters might play only a marginal role in such egalitarian species studied under spatially confined conditions and were not included in the model. To investigate effects of the social structure on the propagation dynamics, we developed a stochastic model that allows us to make predictions about the expected patterns of information transmission in closed social units. This model shall include varying transmission rates for different individuals according to the social structure of the system and it should be formulated in a general way that does not restrict its applicability to only a small class of social systems. 2. Method 2.1. Graph-representation of social systems Social systems can be described in many different ways. For the study of social dynamics it is, however, crucial that the complexity of the system is not lessened by an overly reductionist descriptor. A method to avoid this is to represent the group structure by a socio-matrix (Moreno, 1946). For the study of information transmission, the critical variable is the likelihood with which

ARTICLE IN PRESS B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

information is passed from one specific individual to any other. In some cases—depending on the type of information—this likelihood might be basically determined by the proximity of the individuals (Coussi-Korbel and Fragaszy, 1995). If this is the case, some measure of proximity such as the association index (Clutton-Brock et al., 1982), might be the best descriptor of the social structure. In other cases, where more detailed information has to be obtained or where information can easily be concealed from unwanted onlookers, other measures like nearest-neighbour or grooming data might be more appropriate. To investigate information flow in small-scale societies we mapped the social networks as graphs where the vertices represent individuals and edges between the vertices represent the ‘‘connections’’ between the individuals (see Scott, 1991; Wasserman and Faust, 1994 for extensive treatment of this method). Such an approach has been used to study network dynamics in human (e.g. Moreno, 1946; Nyblom et al., 2003; Watts, 1999; Watts and Strogatz, 1998) and animal (e.g. Croft et al., 2004; Fewell, 2003; Flack et al., 2006; Lusseau and Newman, 2004) social systems. For a recent review on the network approach in the behavioural sciences see (Krause et al., 2007). In our case, the edges represent spatiotemporal associations of a kind that enable the transmission of a specific information unit. The weight of an edge represents the likelihood that such an association is formed between the two individuals of a specific dyad. Henceforth we will refer to the formation of such an association as an interaction. In our examples, data were either derived from proximity matrices where entries in the matrix represent the numbers of instances where two individuals were observed close to each other or from grooming matrices where entries represent the numbers of grooming interactions between two individuals. As we define a group as a set of individuals that interact with each other on a regular basis, the group is represented as a complete graph with ðnk Þ edges. Social groups can differ in the total amount of interactions during a certain time period and the distribution of these interactions between the individuals. As it is trivial that information will be transmitted faster when individuals interact more frequently (all other parameters being constant), we investigated only the question how the group structure, i.e. the distribution of interaction likelihoods might affect information transmission. For this reason interaction likelihoods for dyads were always standardised by dividing through the sum of all matrix elements. For groups of equal size, a higher sum of all matrix elements would mean that individuals were interacting at a higher rate or spent more time in close proximity. However, this conclusion is only valid if data were collected with the same sampling protocol. As this is often not the case when comparing data from different field studies, standardisation will not only eliminate effects due to group size but also due to different sampling procedures.

79

2.2. Simulation To simulate information flow in these networks, we set up the following rules: (1) We assume that all individuals of the group are naı¨ ve with respect to a specific trait (or some sort of information). We represent the knowledge status of the individuals by assigning weights to the vertices. At the onset of the simulation, all individuals are naı¨ ve: oðvÞ ¼ 0

8v 2 V ðG; wÞ.

(1.1)

(2) Thereafter, we randomly choose one individual as inventor by assigning one bit of information to that individual, oðvrÞ ¼ 1. This individual is henceforth called an ‘‘informed individual’’. We further assume that this invention is made only once, as this is suggested to be the typical case for inventions in animal groups (Laland and Reader, 1999; Laland and vanBergen, 2003).1 (3) In each round, two individuals are randomly selected to interact with each other, with the likelihood that this specific dyad is selected being determined by the weight of the edge connecting the corresponding vertices of the graph. If both individuals are naı¨ ve, they will stay naı¨ ve. If one individual is naı¨ ve and the other individual is informed, then the naı¨ ve individual will learn the trait and thereafter both individuals are informed. If both individual are informed, nothing changes: ( oðvÞ ¼

0; 1;

if oðuÞ ¼ oðvÞ ¼ 0; if oðuÞ ¼ 1 _ oðvÞ ¼ 1

8fu; vg 2 V ðG; wÞ.

(1.2) (4) The simulation is terminated as soon as all individuals have acquired the information. This scenario is basically equivalent to the SI-model in epidemic modelling, where there are only two classes of individuals, susceptibles and infected, and once an individual got infected it will not change its status thereafter. The simulation was implemented as a routine for the Mathematica 6.0 software package from Wolfram Research Inc. 3. Experiment 1 The first simulation experiment was carried out in order to determine whether the group structure can influence transmission speed and spreading pattern of socially 1 Questions how multiple inventions would influence propagation dynamics are—though interesting—beyond the scope of this study.

ARTICLE IN PRESS B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

80

propagated information. For this purpose we constructed five artificial group structures (Fig. 1), termed equal, chain, star, linearA and linearB. By using artificial structures we can apply the same construction algorithm to any required group size. Furthermore, this approach allows us to investigate a broader range of networks than when we were just relying on existing data of real social groups. However, to demonstrate how real world data from observational studies fit into this modelling approach we included a socio-matrix of a group of 16 rhesus monkeys (Macaca mulatta) studied by (Sade, 1972), termed rhesus16. The simulation was run 1000 times for each of the five systems for the group sizes of 6–60 in increments of 6 and 10,000 times for the rhesus16 system.

Equal

Chain

Star

3.1. Results: transmission times

LinearA

LinearB

Rhesus16

To compare average transmission times for the five systems, we plotted the mean time (measured in rounds) for the first to the nth transmission in each group of group size n (Table 1). We found clear differences between the investigated systems in both the mean time until the individuals acquired the information and the shape of the transmission curve (Fig. 2). Out of the five investigated systems, the system equal seems to be the most efficient system (i.e. fastest in terms of transmission speed). The system star shows lower transmission times in the intermediate phase of the transmission process, but only for group sizes larger than 12 and it has still higher transmission times in the early and late phases. Not surprisingly (except, perhaps, for some military officers), the system chain where information can only be passed down along one line performed worst. The systems linearA and linearB were designed to be more similar to natural societies where some individuals are more in the centre in the sense that they enjoy more interactions with others or

Table 1 Evaluated parameters Average transmission times

Fig. 1. Representation of the five artificial systems used in this study. In the graph representations vertices are individuals and edges are interactions between individuals. The thickness of the edges is proportional to the interaction rates between the individuals. Vertices are arranged in a circle and the distance between vertices has no specific meaning. Examples are given for group size n ¼ 6. The systems are (a) ‘‘equal’’ with equal interaction rates between all individuals; (b) ‘‘chain’’ which can be thought of as a row of individuals were each one interacts just with its two neighbours, except for the outmost individuals which have only one neighbour; (c) ‘‘star’’ where all individuals interact only with one individual in the centre of the group; (d) ‘‘linearA’’ and (e) ‘‘linearB’’ where interaction rates vary systematically with some individuals having overall high interaction rates and others having overall low interaction rates. The equations describe the construction algorithms for group size n.

Characteristic path length

Average path lengtha

No. of informed individuals a

Mean number of rounds in which the ith (i ¼ 1; :::; n  1) transmission took place, averaged over all simulations for each system and group size Mean of the means of the shortest path lengths connecting each vertex v 2 V ðGÞ of the graph G derived from the socio-matrix to all other vertices of G Mean of the means of the shortest path lengths connecting each vertex v 2 V ðTÞ of the transmission tree T derived from the simulation to all other vertices of T At a given round of a simulation the number of vertices with a vertex weight of 1, jfvjov ¼ 1gj

While some authors use the terms characteristic path length and average path length interchangeably, we use average path length for the transmission tree that resulted from a simulation and characteristic path length for the graph representing the group structure.

ARTICLE IN PRESS B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

400

Round

Gr ou

3000 1500

200 6

6 18

ps iz

0

18 30

ion iss

e

81

30

45 42

m ans

54

Tr

45

42

30

30

15

15

54

100

90 1200

0 6

-100

18 30

45 42

30 54

15

40

600 6

6

0

18 30 42

-10

18 30

45

42

30 54

45 30

15

15

54

Fig. 2. (a) Transmission time for group sizes of 6–60 and the first to the nth transmission for system equal. For the systems (b) chain, (c) star, (d) linearA, and (e) linearB the relative transmission time (i.e. number of rounds needed more or less than for equal) is given. Values are means over 1000 simulations for each group size and each system. Notice unequal scales for the number of rounds.

have higher association indices, while other individuals are more in the periphery in the sense that they have fewer interaction with others or lower association indices. These systems, too, transmit information slower than the equal system, even though linearB is in its performance still quite close to system equal.

L

8

We calculated the average path length for the transmission tree of each simulation as apl ¼ 1=ðn  1ÞSdðviÞ, where n ¼ jV ðTÞj and dðviÞ is the distance of vertex vi from the root r of tree T (Fig. 3). The average path length is a measure for the closeness of the vertices of the graph. Median values for the average path length over all simulations for each structure are plotted against group size in (Fig. 4). The two systems chain and star are exceptional with respect that they are not complete graphs. While these systems do not conform with our original definitions of a group, we still included them in our study because they represent two extreme cases. The system star is the system with the shortest characteristic path length which is given by ð2n  2Þ=n, and the system chain is the system with the longest characteristic path length given by ðn þ 1Þ=3. These systems deliver the lower and upper benchmark (Lovejoy and Loch, 2003) for the expected average path length. All other systems are expected to have average path lengths between these two extremes. The systems linearA and linearB did not differ from system equal in their average path length and even the average path length of the rhesus16 system with several zero weights did not deviate substantially from these systems. As the path length might have considerable consequences for the

R 6

11

5

17

10 K

O

13 Q

4 G

M

E

3

7

16 J

I

B

A 9

F

D

12

H

14

2

C

3.2. Results: path length

1

N

15 P

Fig. 3. Example for a transmission tree for group size 18. Boxes represent individuals ‘‘A’’ to ‘‘R’’, arrows show the direction of the transmission with numbers indicating the temporal order. In the first round information was assigned to a randomly chosen individual (in this example ‘‘L’’) and spread thereafter according to the interaction probabilities derived from the proximity matrix.

expected variation of a trait in a population (see Section 6), we set up a second simulation experiment to investigate this phenomenon more closely. 4. Experiment 2 While in experiment 1 the average path length of the constructed systems equal, linearA and linearB was nearly identical and also the average path length of the system rhesus16 did not deviate much from the former, it is

ARTICLE IN PRESS equal chain star linearA linearB rhesus16

5 4 3

MES

Average Pathlength

6

PDF

B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

82

1 MS

*

2 *

6

2

24

42

60

Groupsize Fig. 4. Means for the average path length for the five artificial systems for group sizes 6–60 and the system rhesus16. For the system star the average path length is increasing approximately linear (y  0:35 þ 0:3x), thus in order to keep the figure in a reasonable proportion, only the first three values are plotted for this system.

possible to construct systems that will yield much higher (or lower) average path lengths. For example, one can take the chain system and fill up the missing connections by connections with weights very close to zero. This system would then fulfil all requirements of a regular group but would still yield an average path length close to the one of the chain system. The same procedure can, of course, be used for the star system. By gradually increasing the weights of these groups, one can generate systems for any expected average path length between the lower and upper benchmark. However, systems with extreme average path lengths might be both highly arbitrary and highly unlikely. While the question of arbitrariness awaits empirical surveys, we tried to answer the question of likelihood by sampling systems with randomly allocated interaction rates. To estimate the expected variance in the average path length we constructed 1000 systems by randomly allocating edge weights (representing interaction likelihoods) between 0 and 1 to each edge of a complete graph for group size 16. Additionally, we constructed 200 systems where random interaction rates s were assigned to the individual agents and edge weights were given by the product of the interaction rates of the two agents oij ¼ si sj . For each of these 1200 systems, we ran the simulation as described above with 1000 repetitions.

3 ME Path length

4

5 MC

Fig. 5. Probability density function (PDF) and mean estimates (MES) for the average path length for the system equal (black line, ME), 1000 systems with randomly assigned interaction rates for all dyads (light grey), and 200 systems with randomly assigned interaction rates for the individuals (dark grey). The expected means for the systems star (MS) and chain (MC) are indicated by solid black lines.

length was slightly higher with 2.4670.51, but still very close to the mean for the system equal. In total 95% of the estimates for the mean average path length of the random systems were between 2.38 and 2.44, thus covering only 1.5% of the bandwidth of possible mean values given by the expected means for the system star (1.87) and chain (5.67). 5. Experiment 3 So far we have assumed that once an individual acquired the information it will keep it until the end of the simulation. Now we want to ask, how a propensity to forget the information again would influence the performance of the systems. To do this we ran another set of simulations with an additional rule: After each round a randomly chosen individual will forget the information with a certain probability (forgetting rate Pðoðui Þ ! 0Þ) and become a naı¨ ve individual again. This scenario is equivalent to the SIS-model in epidemic modelling. Furthermore, we had to introduce a new stopping rule: instead of terminating the simulation as soon as all individuals are informed (a criterion that might not be reached when the forgetting rate is reasonably high), the simulation was terminated after 2000 rounds. This simulation series was carried out with all five systems for a group size of 16 and the rhesus16 system. We simulated transmission with forgetting rates of 0.1%, 0.5%, 1%, 2%, 3% and 5%.

4.1. Results 5.1. Results For a group size of n ¼ 16 the mean average path length for the random systems was 2.4170.46 (mean7S.D.; based on 1000 random systems with 1000 simulations, each). These values were very close to the estimates for the mean average path length of system equal, 2.3970.46 (Fig. 5). For the 200 systems with randomly allocated interaction rates for the individuals, the mean average path

Transmission curves for single simulations of this experiment differed from the previous ones in that they are not necessarily monotonically increasing. If in one round no naı¨ ve individual adopts the information, but one individual forgets it, then the transmission curve is actually decreasing (Fig. 6a). If the number of informed individuals

ARTICLE IN PRESS B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

16 Informed individuals

16 Informed individuals

83

12 8 4

12 8 4 0

0 100

200 300 Rounds

400

100

500

200 300 Rounds

400

500

Proportion of simulations with extinct information

Fig. 6. Examples for transmission curves of single simulations for system equal and group size n ¼ 16 after introduction of a forgetting rate. With a low forgetting rate of (a) 0.01 the number of informed individuals will reach n and will fluctuate only slightly thereafter. With a high forgetting rate of (b) 0.05 the number of informed individuals will not reach n, and will fluctuate more. If the number of informed individuals reaches at one point zero, it cannot be transmitted thereafter and the information is said to have died out.

1

1

1

0.8

0.8

0.8

0.6

0.6

0.4

C

C

E,S,A,B,R

E,S,A,B,R 1000

1000

2000

1

C

0.8

R

0.6

0.4 0.2 1000

2000

1000 C,R

2000

1 C,E,S,A,B,R

A,B 0.8 E S 0.6

0.8

A,B 0.4 E S 0.2

A,B,R E,S

0.2

2000

1

0.6

0.6 0.4

0.4 0.2

0.2

C

0.4 0.2 1000 Rounds

2000

1000

2000

Fig. 7. Proportion of 100 simulations where the information was extinct in the respective round for forgetting rates of (a) 0.001, (b) 0.005, (c) 0.01, (d) 0.02, (e) 0.03, and (f) 0.05 for the systems equal (E), chain (C), star (S), linearA (A), linearB (B), and rhesus16 (R), respectively.

becomes at one point of the simulation zero, the information cannot be transmitted any further and we say that the information ‘‘died out’’ (Fig. 6b). For low forgetting rates, the information died out in a small proportion of the simulations within the first hundred rounds but this proportion reached a stable phase thereafter for all six systems. Though, the time until this equilibrium state was reached differed between the systems (Fig. 7a and b). Almost trivial, for high forgetting rates the information died out in all systems (Fig. 7f). However, for intermediate forgetting rates we found that the information died out only in some systems, while it remained stable over the 2000 rounds period in others (Fig. 7c–e).

6. Discussion As information propagation in social groups cannot be understood without considering the underlying social network structure (Croft et al., 2005), any formal or stochastic model that is put up to produce testable hypotheses has to incorporate the specifics of the social system. Here we presented a network approach to study the propagation of information in small-scale social groups. Simulating information-flow in artificial as well as natural group structures we found that the speed of information transmission increased as a function of the proportion of informed individuals in a group. Such a ‘‘bush-fire effect’’

ARTICLE IN PRESS 84

B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

is a general characteristic of social trait propagation and well documented in the literature (Cavalli-Sforza and Feldman, 1981; Giraldeau and Caraco, 2000). However, we could also show that the strength of this effect varied markedly between the investigated systems. This is in accordance with the suggestions made by Coussi-Korbel and Fragaszy (1995), that the structure of an animal group shall influence the efficiency of social learning. Empirical studies on primates (Bonnie and de Waal, 2006; Nahallage and Huffman, 2007) and fish (Reader and Laland, 2000) delivered evidence that information spreads faster within classes of closely related animals or individuals of the same sex. As members of such classes often spend more time with each other than with members of other classes, these findings are in accordance with our model. By introducing a forgetting rate (Section 5) we could show, that such differences in propagation speed can have considerable impact on the longevity of cultural traditions. Interestingly, the systems’ sensitivity to forgetting is not equivalent with their efficiency to transmit information: the most efficient system in terms of propagation speed (system equal) was not the most robust against forgetting. This forgetting rate can either be understood as the forgetting of the information by single individuals or that informed individuals die or leave the group and are replaced by naı¨ ve individuals. While the former case will be of interest primarily when modelling horizontal propagation within individuals of one generation, the latter interpretation enables us to model vertical transmission over many generations. If we assume that the transmission process from one individual to the other is a source for ‘‘transcription errors’’, then the length of a transmission chain shall influence the extent to which the behaviour is modified. Therefore, the average path length of the information chain can provide us with an estimate which variability of a trait we might expect in a population with a given structure. Except for the highly arbitrary systems chain and star, we could not find clear differences in the expected path lengths of the investigated systems. This was also the case for 1200 systems with randomly allocated interaction rates. Thus, it seems that contrary to our expectations for a broad range of systems group structure does not influence the expected average path length, and hence the expected variability of a trait. In this respect, the systems resemble ‘‘small worlds’’ (Watts and Strogatz, 1998) where the mean path length is proportional to the logarithm of the group size. The dynamics of information transmission within populations equal those of disease transmission which have been studied extensively in the past (see Daley and Gani, 1999; Hethcote, 2000). If populations are fully mixed, that is, any individual has the same fixed probability per unit of time of contracting the disease from any infected individual in the population, the number of infected individuals can be modelled by a system of ordinary differential equations. While the basic SI model allows only for two classes of individuals, susceptibles and infected, the SIR model assumes that infected become

removed probabilistically or after a certain time period and the SIS model considers infected that become susceptible again (Hethcote, 2000). However, modelling the number of infected and susceptible individuals as continuous variables gives only poor approximations in small populations. In such cases, the system can be modelled as a Markov process with discrete population and time, where the probabilities of transmission from ðs; iÞ to ðs  1; i þ 1Þ and ðs; i  1Þ are given by bsi DT and gi DT, respectively (Witten and Poulter, 2007). Recently, network approaches to disease (Eubank et al., 2004; Liljeros et al., 2001; May and Lloyd, 2001; Newman, 2003; Nowak and Bangham, 1996; Verdasca et al., 2004; Witten and Poulter, 2007) and information (Kincaid, 2000; Kossinets and Watts, 2006; Valente, 2005; Watts and Strogatz, 1998) transmission have proven to allow more detailed predictions as they release the assumption of fully-mixed populations. Primate groups are small—they usually consist of less than hundred animals—and despite the fact that every animal interacts with any other they are not fully-mixed as interaction rates differ significantly between dyads. To account for these characteristics we represented the groups as complete weighted graphs. As several frequently used measures like characteristic path length, connectivity or betweenness are only of limited or no use to analyse such graphs (Barrat et al., 2004; Lovejoy and Loch, 2003), we opted for a simulation approach, where information flow was modelled as a stochastic process. The outcome of each simulation was a propagation tree that allowed us to measure transmission speed and average path length. Because we hypothesised that relative interaction frequencies are basically determined by spatial proximity we used a measure of proximity as predictor for information transfer. As proximity matrices are by their very nature symmetric (whenever A is close to B, B is also close to A), the systems could be represented as simple graphs. However, if there is reason to believe, that the likelihood of information-flow between two individuals is not symmetric (e.g. subordinates are more likely to learn from dominants), then our approach can be easily extended using directed multigraphs where each pair of vertices is connected by two edges, one in each direction. We did not attempt to develop a framework that requires to classify social systems or to rank them along a single dimension. While a simple classification scheme would facilitate generalisations about system behaviour, we believe that the representation of social systems by sociomatrices has several advantages. (1) It is less reductionist than any classification system, as it preserves the complexity of the system. (2) It is more flexible and allows the investigation of any possible system. There is no uncertainty about class-membership and there are no exceptional structures that do not fit into the framework. (3) Is directly applicable to empirical data, as socio-matrices are well-established tools in the social and behavioural sciences and as these socio-matrices can be mapped directly onto the adjacency matrix of the graph.

ARTICLE IN PRESS B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

With this study we could demonstrate the potential influence of group structure on the propagation of information. The next steps will be to extend the bandwidth by investigating properties of a broad range of existing social systems, and to test predictions stemming from these models experimentally in captive and wild primate groups. Such tests will consist of four components: establishing the social structure of the study group, making predictions for the most likely spreading patterns based on the social data, introducing a novel behaviour and recording its spreading pattern within the group. The suggestion that mechanisms of cultural learning in non-human primates might be homologous to those in humans, stimulated research on cultural traditions especially in primates. However, this modelling approach shall not be restricted to primates, but it can be applied to any species living in social groups.

Acknowledgements We want to thank C. Kasper, M. Steurer, L. Huber for fruitful discussions and C. Sueur and one anonymous referee for valuable comments on the manuscript. This research was funded by the 6th Framework Program of the European Union as part of the NEST Project GEBACO— The Genetic Basis of Cooperation, Con. 28696.

References Barrat, A., Barthelemy, M., Pastor-Satorras, R., Vespignani, A., 2004. The architecture of complex weighted networks. Proc. Natl. Acad. Sci. USA 101, 3747–3752. Biro, D., Inoue-Nakamura, N., Tonooga, R., Yamakoshi, G., Sousa, C., Matsuzawa, T., 2003. Cultural innovation and transmission of tool use in wild chimpanzees: evidence from field experiments. Anim. Cognit. 6, 213–223. Bonnie, K.E., de Waal, F.B.M., 2006. Affiliation promotes the transmission of a social custom: handclasp grooming among captive chimpanzees. Primates 47, 27–34. Boyd, R., Richerson, P.J., 1985. Culture and the Evolutionary Process. University of Chicago Press, Chicago. Cavalli-Sforza, L., Feldman, M.W., 1981. Cultural Transmission and Evolution: A Quantitative Approach. Princeton University Press, Princeton. Clutton-Brock, T.H., Guinness, F.E., Albon, S.D., 1982. Red Deer: Behavior and Ecology of Two Sexes. University of Chicago Press, Chicago. Coussi-Korbel, S., Fragaszy, D.M., 1995. On the relation between social dynamics and social learning. Anim. Behav. 50, 1441–1453. Croft, D.P., Krause, J., James, R., 2004. Social networks in the guppy (Poecilia reticulata). Proc. R. Soc. Lond. B (Suppl.) 271, 516–519. Croft, D.P., James, R., Ward, A.J.W., Botham, M.S., Mawdsley, D., Krause, J., 2005. Assortative interactions and social networks in fish. Oecologia 143, 211–219. Curio, E., Ernst, U., Vieht, W., 1978. Cultural transmission of enemy recognition: one function of mobbing. Science 202, 899–901. Daley, D.J., Gani, J., 1999. Epidemic Modelling. Cambridge University Press, Cambridge. de Waal, F.B.M., 1999. Cultural primatology comes of age. Nature 399, 635 doi:10.1038/21310.

85

Eubank, S., Guclu, H., Kumar, V.S., Marathe, M.V., Srinivasan, A., Toroczkai, Z., Wang, N., 2004. Modelling disease outbreaks in realistic urban social networks. Nature 429, 180–184. Fewell, J.H., 2003. Social insect networks. Science 301, 1867–1870. Flack, J.C., Girvan, M., de Waal, F.B.M., Krakauer, D.C., 2006. Policing stabilizes construction of social niches in primates. Nature 439, 426–429. Galef, B.G., 1991. Traditions in animals: field observations and laboratory analyses. In: Bekoff, M., Jamieson, D. (Eds.), Interpretation and Explanation in the Study of Behavior, Vol 1: Interpretation, Intentionality and Communication. Westview Press, Boulder. Galef, B.G., Allen, C., 1995. A new model system for studying behavioural traditions in animals. Anim. Behav. 50, 705–717. Giraldeau, L.A., Lefebvre, L., 1986. Exchangeable producer and scrounger roles in a captive flock of feral pigeons: a case for the skill pool effect. Anim. Behav. 34, 797–803. Giraldeau, L.A., Lefebvre, L., 1987. Scrounging prevents cultural transmission of food-finding behaviour in pigeons. Anim. Behav. 35, 387–394. Giraldeau, L.A., Caraco, T., 2000. Social Foraging Theory. Princeton University Press, Princeton, NJ. Goodall, J., 1968. The behaviour of free-ranging chimpanzees in the Gombe Stream Reserve. Anim. Behav. Monogr. 1, 161–311. Hethcote, H.W., 2000. The mathematics of infectious diseases. SIAM Rev. 42, 599–653. Heyes, C.M., Galef, B.G. (Eds.), 1996. Social Learning in Animals: The Roots of Culture. Academic Press, San Diego. Hinde, R.A., Fisher, J., 1951. Further observations on the opening of milk bottles by British birds. British Birds 44, 393–396. Horner, V., Whiten, A., Flynn, E., de Waal, F.B.M., 2006. Faithful copying of foraging techniques along cultural transmission chains by chimpanzees and children. Proc. Natl. Acad. Sci. 103, 13878–13883. Huffman, M.A., Hirata, S., 2004. An experimental study of leaf swallowing in captive chimpanzees: insights into the origin of a selfmedicative behavior and the role of social learning. Primates 45, 113–118. Huffmam, M.A., Seifu, M., 1989. Observations on the illness and consumption of a possibly medical plant, Veronica amygdalina (Del.), by a wild chimpanzee in the Mahale Mountains National Park, Tanzania. Primates 30, 51–63. Jacobs, R.C., Campbell, D.T., 1961. The perpetuation of an arbitrary tradition through several generations of a laboratory microculture. J. Abnormal Social Psychol. 62, 649–658. Kawai, M., 1965. Newly-acquired pre-cultural behavior of the natural troop of Japanese monkeys on Koshima islet. Primates 6 (1), 1–30. Kendal, J.R., Kendal, R.L., Laland, K.N., 2007. Quantifying and modelling social learning processes in monkey populations. Int. J. Psychol. Psychol. Therapy 7, 123–138. Kincaid, D.L., 2000. Social networks, ideation, and contraceptive behavior in Bangladesh: a longitudinal analysis. Soc. Sci. Med. 50, 215–231. Kossinets, G., Watts, D.J., 2006. Empirical analysis of an evolving social network. Science 311, 88–90. Krause, J., Croft, D.P., James, R., 2007. Social network theory in the behavioural sciences: potential applications. Behav. Ecol. Sociobiol. 62, 15–27. Kroodsma, D.E., Baylis, J.R., 1982. A world survey for evidence of vocal learning in birds. In: Kroodsma, D.E., Miller, E.H. (Eds.), Acoustic Communication in Birds, vol. 2. Academis Press, New York, pp. 311–337. Laland, K.N., Plotkin, H.C., 1990. Social learning and social transmission of foraging information in Norway rats (Rattus norvegicus). Anim. Learn. Behav. 18, 246–251. Laland, K.M., Reader, S.M., 1999. Foraging innovation in the guppy. Anim. Behav. 57, 331–340. Laland, K.N., vanBergen, Y., 2003. Experimental studies of innovation in the guppy. In: Reader, S.M., Laland, K.N. (Eds.), Animal Innovation. Oxford University Press, Oxford, pp. 155–173.

ARTICLE IN PRESS 86

B. Voelkl, R. Noe¨ / Journal of Theoretical Biology 252 (2008) 77–86

Laland, K.N., Williams, K., 1997. Shoaling generates social learning of foraging information in guppies. Anim. Behav. 53, 1161–1169. Laland, K.N., Richerson, P.J., Boyd, R., 1996. Developing a theory of animal social learning. In: Heyes, C.M., Galef, B.G. (Eds.), Social Learning in Animals: The Roots of Culture. Academic Press, San Diego, pp. 129–154. Lefebvre, L., 1995. Culturally-transmitted feeding behaviour in primates: evidence for accelerating learning rates. Primates 36, 227–239. Lepoivre, H., Pallaud, B., 1985. Social facilitation in a troop of Guinea baboons (Papio papio) living in an enclosure. Behav. Process. 11, 405–418. Liljeros, F., Edling, C.R., Amaral, L.A.N., Stanley, H.E., Aberg, Y., 2001. The web of human sexual contacts. Nature 411, 907–908 doi:10.1038/ 35082140. Lohnsdorf, E.V., 2006. What is the role of mothers in the acquisition of termite-fishing behaviors in wild chimpanzees (Pan troglodytes schweinfurthii)? Anim. Cognit. 9, 36–46. Lovejoy, W.S., Loch, C.H., 2003. Minimal and maximal characteristic path length in connected sociomatrices. Social Netw. 25, 333–347. Lusseau, D., Newman, M.E.J., 2004. Identifying the role that individual animals play in their social network. Proc. R. Soc. Lond. B 271, 477–481. May, R.M., Lloyd, A.L., 2001. Infection dynamics on scale-free networks. Phys. Rev. E 64, 066112. Moreno, J.L., 1946. Sociogram and sociomatrix: a note to the paper by Forsyth and Katz. Sociometry 9, 348–349. Nahallage, C.A.D., Huffman, M.A., 2007. Age-specific functions of stone handling, a solitary-object play behavior, in Japanese macaques (Macaca fuscata). Am. J. Primatol. 69, 1–15. Newman, M.E.J., 2003. Spread of epidemic disease on networks. Phys. Rev. E 66, 016128. Nowak, M.A., Bangham, C.R.M., 1996. Population dynamics of immune responses to persistent viruses. Science 272, 74–79.

Nyblom, J., Borgatti, S., Roslakka, J., Salo, M.A., 2003. Statistical analysis of network data—an application to diffusion of innovation. Social Netw. 25, 175–195. Reader, S.M., 2004. Distinguishing social and asocial learning using diffusion dynamics. Learn. Behav. 32, 90–104. Reader, S.M., Laland, K.M., 2000. Diffusion of foraging innovations in the guppy. Anim. Behav. 60, 175–180. Reader, S.M., Kendal, J.R., Laland, K.M., 2003. Social learning of foraging sites and escape routes in wild Trinidadian guppies. Anim. Behav. 66, 729–739. Sade, S.D., 1972. Sociometrics of Macaca mulatta I. Linkages and cliques in grooming matrices. Folia Primatol. 18, 196–223. Scott, J., 1991. Social Network Analysis: A Handbook. Sage Publications, Newbury Park. Sugiyama, Y., Koman, J., 1979. Tool-using and making behavior in wild chimpanzees at Bossou, Guinea. Primates 20, 513–524. Valente, T.W., 2005. Network models and methods for studying the diffusion of innovations. In: Carrington, P.J., et al. (Eds.), Models and Methods in Social Network Analysis. Cambridge University Press, Cambridge, pp. 98–116. Verdasca, J.A., Telo da Gama, M.M., Nunes, A., Bernardino, N.R., Pacheco, J.M., Gomes, M.C., 2004. Recurrent epidemics in small world networks. arXiv:q-bio.PE, 0408002 v1. Wasserman, S., Faust, K., 1994. Social Network Analysis: Methods and Applications. Cambridge University Press, Cambridge. Watts, D.J., 1999. Small Worlds: The Dynamics of Networks between Order and Randomness. Princeton University Press, Princeton, NJ. Watts, D.J., Strogatz, S.H., 1998. Collective dynamics of ‘small-world’ networks. Nature 393, 440–442 doi:10.1038/30918. Witten, G., Poulter, G., 2007. Simulations of infectious diseases on networks. Comput. Biol. Med. 37, 195–205. Zentall, T.R., Galef, B.G. (Eds.), 1988. Social Learning. Psychological and Biological Perspectives. Lawrence Erlbaum Associates Publishers, Hillsdale, NJ.