Information about Cultural Propagation

Published on May 31, 2008

Author: epokh

Source: slideshare.net

M. N. KUPERMAN PHYSICAL REVIEW E 73, 046139 2006 foods, or musical style, readings, sports, etc. In turn, a fur- 1 ther subclassiﬁcation of each feature into categories will = i 4 2 i serve to denote the different preferences or traits. The sim- plest consideration to achieve such subdivision is to consider cannot be claimed to have a monotonic behavior when the that any cultural feature may take on any of q different val- system, as deﬁned, evolves in time. Suppose that as a result ues, the same for any k. In principle the values are only of the interaction between i and j at time t there is a change labels, so despite the fact that we use numbers for the clas- in the value of ik. We will call m the neighborhood of i i siﬁcation, any set of symbols would work as well. Thus, an m hm such that for any of the individuals hm i , k t individual i is culturally characterized by a cultural vector of = k t ; and i the set of neighbors hn such that hn t i n k F components ik, each one adopting values ranging between = ik t + 1 . The rest of the neighbors will be included in the 1 and q. The way to culturally compare two individuals i and set li. If at each time step only one change is allowed, we j is to measure their cultural overlap, i,j as follows: can calculate = t + 1 − t = i t + 1 − i t by con- sidering that F = i j . 1 i t+1 = i,j t + i,j t −1 i,j k k l m k=1 j i j i The individuals are situated on the vertices of a graph. + i,j t +1 Two of them are neighbors when linked by an edge. In the j n i original model, the underlying network was a bidimensional regular lattice 12 . The topology of the network was later = i,j t − 1+ 1. 5 generalized in 14 , considering amongst others, small world j m n i j i j i SW networks 2 . Here, we also consider SW networks. As the individuals interact with the set of their neighbors, Thus it is useful to deﬁne at this point the local overlap of a given t+1 = t + − . 6 i i individual i as The change in is thus = − , which is not necessarily equal or greater than zero. By introducing a modiﬁcation in i = i,j , 2 the dynamics of the original model we can assure that this j i condition will be fulﬁlled and thus can think of a Lyapunov functional for the system. where i is the set of neighbors of i. We will also deﬁne the We will consider different types of dynamics, each one quadratic local overlap as associated with a corresponding Lyapunov function but at the same time with a clear social interpretation of the behav- ior of the individuals. The underlying network will be built 2 i = i,j . 3 up following the procedure described in Ref. 22 . In the j i original model of SW networks, a single parameter p, run- Starting from an initial distribution of cultural vectors, the ning from 0 to 1, characterizes the degree of disorder of the individuals evolve by analyzing and interacting with their network, respectively, ranging from a regular lattice to a environment, adapting their cultural preferences according to completely random graph. The construction of these net- the tendencies of the neighborhood. The original numerical works starts from a regular, one-dimensional, periodic lattice simulations proceed as follows in Refs. 12,14,15 . At time of N elements linked to 2K nearest neighbors. Then each of step t, a randomly chosen individual i and one of its neigh- the sites is visited, rewiring K of its links with probability p. bors j are evaluated. Their cultural overlap is calculated to Values of p within the interval 0 , 1 produce a continuous decide whether they will interact or not. The interaction takes spectrum of small world networks. place with a probability i,j t / F, in which case one of the CULTURAL EXCHANGE DYNAMICS features ik such that ik j j k is set equal to k. Though it is evident that i,j t Case 1: Restricted cultural afﬁnity i,j t − 1 , the interaction may affect as well the overlaps between i and the rest of the neighbor- The original model proposed by Axelrod considered a hood and thus the change on i cannot be anticipated. very special case of biased dynamics for the interactions of Interesting results were obtained in Ref. 15 by consid- the individuals. Despite the fact that individuals are im- ering the whole process of cultural dissemination as an opti- mersed in their neighborhood, this was ignored by requiring mization problem. In that work, by analyzing a one dimen- that the individual interact with only one of its neighbors. sional system with interaction amongst the ﬁrst neighbors, Taking this fact into account, a ﬁrst adaptation of the original the authors found a Lyapúnov potential that allowed them to model consists in deciding whether to change or not the analyze the stability of the states at which the system re- value of the chosen feature by weighting the decision with a mained frozen after some evolution time. We are interested further evaluation of the inﬂuence of the neighborhood. in extending those results to more general situations, namely, Given that the individual i interacts with j, the possibility of SW and other complex networks. The global overlap adopting k for ik will depend on the result of an evaluation j 046139-2

CULTURAL PROPAGATION ON SOCIAL NETWORKS PHYSICAL REVIEW E 73, 046139 2006 following a sort of majority rule. If by accepting the change n m n m i t+1 = i t + i t+1 − i t+1 + i t − i t , of the value of ik i will share the value taken by ik with a bigger group than if by rejecting the change, then i accepts = t +2 n t+1 − m t + − . i i i the change. With a probability 1/2 the change is accepted in case of equality. This is translated into the following situa- On the other hand, we have Eq. 2 . Finally tion. The change is accepted whenever , and with prob- n m ability 1/2 when = . Under this condition, t+1 − t £2 t + 1 − £2 t = − 2 i t+1 − i t . 0. So, we will take £1 t = − t as the Lyapunov function The condition to be fulﬁlled is of this dynamics, and the system will evolve to reach a local n m or absolute minimum. i t+1 − i t 0 Case 2: Complete cultural afﬁnity that corresponds to the imposed constraint. It is important to note that in all the cases, the monocul- The former rule assures that the global cultural back- tural state corresponds to minimum value of the Lyapunov ground grows or at least is maintained constant, while the function £iM . We can use this value for normalization, such individuals knows that accepting the change will warrant be- that Li = £i / £iM . ing in a bigger group regarding the changing feature. But basing the decision of the individual on the comparison of only one feature seems quite myopic. We can propose an- NUMERICAL RESULTS other condition, making the individual base the decision on a further evaluation of the local partial overlaps m In what follows we will include results corresponding to i n = j m i,j and i = j n i,j . Despite that one feature the cases 1 and 2 as well as those corresponding to the Ax- i i elrod’ s original model, for which we have not deﬁned a was chosen to be changed, the individual decides whether to Lyapunov function. We have performed extensive numerical adopt the new value or not by weighting the whole cultural simulations of the described model, considering different dy- overlap with its neighborhood and not by analyzing what namics. The networks have N = 104 vertices and connectivity happens with the speciﬁc feature to be changed. This is K = 2. A typical realization starts with the generation of the equivalent to saying that the individual will favor a majority random network and the initialization of the state of the el- weighted by deeper cultural afﬁnity. Now, we can no longer ements. After a transient period, the duration of which de- say that t + 1 − t 0. We must look for another quan- pends on the parameters of the particular simulation, a mac- tity £2 t , such that £2 t + 1 − £2 t 0. In what follows we roscopic stationary state is achieved. The computations are 1 show that £ t 2 = − t + t , with t = 2 i i t , satis- then repeated for several thousand time steps to perform sta- ﬁes the required condition. tistical averages. We consider that the system has achieved a Let us consider that in the proposed interaction between i macroscopic stationary state when the corresponding and j, ik will be change by k. There are three classes of j Lyapunov function of the systems reaches a stationary value. individuals among the neighbors of i, those belonging to m, i We will see that this does not imply that the system is steady those belonging to n and the rest, that will be grouped in li. i in a particular microscopical state. Indeed, the conﬁguration The local partial overlap li = j l i,j will no be affected, of the system ﬂuctuates among states associated to equal i regardless of whether or not the interaction takes place. On values of the Lyapunov function. There are several aspects the contrary, if the interaction occurs at time t, m t + 1 i characterizing the asymptotic evolution of the system to a = m t − , and n t + 1 = n t + . i i i stationary value of the Lyapunov function. We also analyze If no interaction is allowed, £2 t + 1 = £2 t . If on the con- the behavior of the system when governed by the original trary, the change is accepted, we have dynamics, in which case the steady state is not characterized by a Lyapunov function and achieves a microscopically fro- £2 t + 1 − £2 t = i t+1 − i t + i t+1 − i t . zen state. We can expand the right-hand side of the former equation by In all the calculations, we took F = 10 and several values considering sums over m, n, and li. On one side we have of q, ranging from 2 to 80. At each time step only one change i i was allowed, the system was updated asynchronically. We i t+1 = i,j t 2 + i,j t −1 2 considered that one unity of time corresponded to N time j l j m steps. i i For each of the dynamics described above, we have ana- 2 + i,j t +1 lyzed several aspects of the evolution of the system. First we j n i have calculated the proportion of overlaps between indi- viduals corresponding to three cases, a 0 when i,j = 0, 2 = i,j t + 1−2 i,j t b F when i,j = F, and c a when 0 F. Cases m i,j j i j i a and b correspond to the situation when no change in the + 1+2 t . 7 system is possible because the interaction of two individuals: i,j n in case a because no interaction will occur when the cul- j i tural overlap is equal to zero, in the case b because indi- Expanding and regrouping terms we get viduals are already culturally identical. The only active links 046139-3

M. N. KUPERMAN PHYSICAL REVIEW E 73, 046139 2006 FIG. 1. Proportion of active links a and of complete overlap FIG. 2. Proportion of changes vs time, for different values of q links F vs time, with q = 2 full , q = 5 dashed , q = 10 dotted , q and p as in Fig. 1. Axelrod’s case. = 15 dotted-dashed . Each plot correspond to a different value of p: A p = 0, B p = 0.01, C p = 0.5, D p = 0.9. Axelrod’s case. als. Each change corresponds to a component of any cultural vector that changed its value. We show the amount of are those corresponding to case c . Then we have calculated changes in a unit of time normalized to the maximum value the corresponding Lyapunov function when deﬁned to allowed N, the number of proposed changes. show how its value evolves monotonically to a steady one. Figures 1 a and 2 a correspond to an ordered underlying Though this does not provide any information about the in- network. The system goes to a state where only non active ner structure of the system, or about the existence of clusters, links survive, that is, a → 0. At the same time, while the it helps us to have an idea of the amount of cultural differ- system reaches a steady state, associated with the number of entiation that is present. To show that though the Lyapunov changes approaching 0, the system achieves a monocultural function reaches a steady value, but the system is not in a state when q F but goes to a multicultural state for higher steady state, we calculated the amount of changes that occur values of q. in each time step. This was also useful to show that in the When some disorder is introduced into the network, the Axelrod case, the system attained a frozen state, with no behavior of the system is more complex. By looking at Fig. changes. 2 we can see that there are two different behaviors for or- dered and very disordered networks while the intermediate Axelrod’s case case, p = 0.01 shows a mixture of both. The system, in or- This case corresponds to the original model 12,14,15 dered networks evolves rather fast to a state of low multicul- where individuals interact with only one of their neighbors at turality or monoculturality. When the disorder is increased, each time step. The interaction is mediated by the cultural the initial disorder survives for longer times. At the end, the afﬁnity, deﬁned through the cultural overlap i,j . The stron- system ends in a monocultural state except when p = 0 and ger the afﬁnity is, the greater the possibility of interaction q F. Figure 1, showing the number of changes in time con- between two subjects. In the present work, the individuals ﬁrms what was mentioned before. We have not observed are located on networks with different degrees of disorder. sharp transitions while varying the p value, indeed, the be- The ordered case p = 0, corresponds to a one-dimensional havior of the system undergoes a smooth change as the spa- lattice with interactions between the ﬁrst and second neigh- tial disorder is increased. bors. As stated before, we did not ﬁnd a Lyapunov function for Case 1: Restricted cultural afﬁnity this case, and we restrict the displayed results to the time In the following cases the calculation of the Lyapunov dependence of the proportion of overlaps F and a, and of function will provide us additional information about the the proportion of changes in the individuals’ cultural proﬁles. system behavior. As in the previous case, Fig. 3 displays the We recall that 0 = 1 − F + a . Figure 1 displays the time time evolution, averaged over 1000 realizations, of the val- evolution, averaged over 1000 realizations, of the values F ues F and a, evaluating these quantities on networks of and a, corresponding to the amount of overlaps i,j = F and varying disorder. Figure 5 shows the evolution in time of the 0 i,j F, normalized to the total number of links KN. We proportion of changes. The evolution of the normalized evaluate these quantities on networks with different degree Lyapunov function corresponding to this case, L1 is plotted of disorder, namely, p = 0 , 0.01, 0.5, 0.9. in Fig. 4. Starting from the ordered case, Fig. 3 a , we ob- A deeper insight into what is happening is obtained by serve that results do not differ so much from what we have analyzing the data contained in Fig. 2. There we show the previously seen. Again, the system goes to a state where only proportion of changes in the cultural vectors of the individu- non active links survive, reaching a steady state, and achiev- 046139-4

CULTURAL PROPAGATION ON SOCIAL NETWORKS PHYSICAL REVIEW E 73, 046139 2006 FIG. 5. Proportion of changes vs time, for different values of q FIG. 3. Proportion of active links a and of complete overlap and p as in Fig. 1. Case 1. links F vs time, for different values of q and p as in Fig. 1. Case 1. ing monoculturality when q F and a certain degree of mul- number of changes remains above zero in all cases. Again, ticulturality for higher values of q. the mean value of changes behaves in a nontrivial way when This time we can recur to the Lyapunov function to see q or p change. that the absolute minimum is reached when q F, but the Perhaps the most interesting feature is the interplay be- system remains in a frozen state of multiculturality when q tween the effect of the spatial disorder and the values of q. F. It is interesting to observe that the number of changes, This can be better observed by analyzing the behavior of the Fig. 5 a , goes to zero. Lyapunov function. In some cases the disorder introduced by When disorder is included on the network, the behavior of the network prevents the system from achieving the previ- the system displays nontrivial effects as can be can observed ously reached monocultural state, but on the other hand, the in Figs. 3 b –3 d . The number of active links is different ﬁnal degree of multiculturality depends in a very interesting from zero, even when a steady value for the Lyapunov func- way from both parameters. An interesting nonmonotonic be- tion is reached. Though the monocultural state is the absolute havior of L1 in terms of p can be observed depicted in Fig. 6. minimum, it is note attained by the system, who ﬁnishes trapped in local minimum. In Figs. 4 b –4 d we observe that Case 2: Complete cultural afﬁnity the Lyapunov function decreases monotonically to attain a The ﬁrst aspect that we can observe for this case is that steady state but not to the absolute minimum. On the other independently of the degree of disorder of the network, the hand, the steady values depend non monotonically on the state of monoculturality is never achieved, as shown in Fig. disorder of the network. Despite the fact that £1 remains 7. We can again verify the interplay between the parameters steady, the state of the system is not frozen. This afﬁrmation q and p and their effect on the behavior of the system. An- comes from the observation of Fig. 5, where we ﬁnd that the FIG. 4. Normalized Lyapunov function L1 vs time, for different FIG. 6. Steady value of −L1 vs p, with q = 2 full , q = 5 values of q and p as in Fig. 1. dashed , q = 10 dotted , q = 15 dotted-dashed . 046139-5

M. N. KUPERMAN PHYSICAL REVIEW E 73, 046139 2006 FIG. 7. Proportion of active links a and of complete overlap FIG. 9. Proportion of changes vs time, for different values of q links F vs time, for different values of q and p as in Fig. 1. Case 3. and p as in Fig. 1. Case 2. other issue to be observed is the time scale. The evolution system. Further analysis 13 of the relative size of the larg- towards a steady value of the Lyapunov function is much est cultural domain revealed an order disorder transition with faster than before, as observed in Fig. 8. At the same time we q, the number of different traits, playing the role of the con- observe that by increasing the disorder the amount of active trol parameter. Under a threshold value qc F , the system links grows. converges to a monocultural uniform state. Above qc F the This is associated to the fact observed in Fig. 9, where we system freezes in a multicultural state, that can be associated can see how the amount of changes in the ﬁnal state also to polarization. The stability of the multicultural states was increases with p. The behavior of the Lyapunov function analyzed in Ref. 14 by perturbing the system when frozen simply veriﬁes that the system reaches a steady value and in a multicultural state and showing the further convergence that this value is far from being the absolute minimum. In all to the monocultural state. Perturbations were associated to the cases, the steady value decreases with q. cultural drift. In this work we proposed a different sort of generalization of Axelrod’ s model. We modiﬁed the model to include in- CONCLUSIONS teractions among several individuals within a neighborhood Axelrod’s model shows how a microscopical local pro- cess of interaction, leading to convergence provokes the emergence of global polarization. In previous works, the model was used to analyze the effect of the number of cul- tural aspects and traits on the steady conﬁguration of the FIG. 10. Asymptotic proportion of inactive links 0 solid and FIG. 8. Normalized Lyapunov function L2 vs time, for different F dashed for different values of q. With bold line: p = 0.9 and thin values of q and p as in Fig. 1. line p = 0, solid line. a Case 1, b case 2. 046139-6

CULTURAL PROPAGATION ON SOCIAL NETWORKS PHYSICAL REVIEW E 73, 046139 2006 or to let each individual evaluate the changes in its cultural turality to multiculturality at different values of q. This can preferences by analyzing those of its neighbors. The cultural be explain by recalling that in a disordered network the clus- inﬂuence of the environment was already studied in Ref. terization of the system is lower and thus, the existence of 16 . clusters of culture reﬂected in a polarized situation is no Different ways of considering this extended interaction longer achieved. When the slightest disorder is added to the were shown. For each, an associated Lyapunov function was network, the number of links with overlap equal to zero de- found, letting us analyze the convergence of the system to- cays. In case 2, the transition to the multicultural state occurs wards an absolute or local minima. The disorder of the sys- tem was not reduced to that introduced by the initial condi- at lower values of q when compared with previous results. It tion by increasing the value of q, but also included in the is important to recall that multiculturality presents here a spatial distribution of the agents. For this purpose we ana- different character. The change of the rules of interaction lyzed the effect of the disorder of the underlying network introduces a new interesting behavior. Not only does the considering small world networks of varying disorder. The amount of active links not go to zero, with the exception results linked to this aspects can be compared with previous when the underlying lattice is ordered and q = 2, but also the results and thus unveil the effect of the newly deﬁned inter- system reaches a situation when the Lyapunov function action of each individual with the whole neighborhood. As adopts a steady value but the system is not frozen. The con- already known, increasing the value of q leads the system to ﬁguration of the system changes in time, as can be observed undergo a transition from monoculturality to multiculturality. from the ﬁgures displaying the number of changes in time. However, when the dynamics of the system corresponds to The results presented here complement what was already the case 1, the effect of spatial disorder attempts against this found in the analysis of the model ﬁrst proposed by Axelrod. effect. Figure 10 shows the asymptotic values of 0, a, and The interesting feature is that the system, despite reaching a F for different values of q and dynamics. Figure 10 a cor- steady situation, does not remains static. Some aspects still responding to the case 1 and Fig. 10 b to the case 2. In Fig. deserve further analysis. Among them we will consider in a 10 a it is possible to observe the transition from monocul- future work the inclusion of noise. 1 W. Weidlich, Phys. Rep. 204, 1 1991 . 13 C. Castellano, M. Marsili, and A. Vespignani, Phys. Rev. Lett. 2 D. J. Watts, Small Worlds Princeton University Press, Prince- 85, 3536 2000 . ton, 1999 . 14 K. Klemm, V. M. Eguíluz, R. Toral, and M. San Miguel, Phys. 3 R. Axelrod, The Complexity of Cooperation Princeton Univer- Rev. E 67, 026120 2003 . sity Press, Princeton, 1997 . 15 K. Klemm, V. M. Eguíluz, R. Toral, and M. San Miguel, Phys. 4 K. Sznajd-Weron and J. Sznajd, Int. J. Mod. Phys. C 11, 1157 Rev. E 67, 045101 R 2003 . 2000 . 16 Y. Shibanai, S. Yasuno, and I. Ishiguro, J. Conﬂict Resolut. 45, 5 D. H. Zanette, Phys. Rev. E 65, 041908 2002 . 6 M. Kuperman and G. Abramson, Phys. Rev. Lett. 86, 2909 80 2001 . 2001 . 17 J. C. Gonzalez-Avella, M. G. Cosenza, and K. Tucci, Phys. 7 R. Pastor-Satorras and A. Vespignani, Phys. Rev. Lett. 86, Rev. E 72, 065102 R 2005 ; nlin.AO/0511013. 3200 2001 . 18 C. Klemm, V. M. Eguluz, R. Toral, and M. San Miguel, J. 8 P. Holme and A. Grönlund, J. Artif. Soc. Soc. Simul. 8 2005 . Econ. Dyn. Control 29, 321 2005 . 9 M. Kuperman and D. Zanette, Eur. Phys. J. B 26, 387 2002 . 19 G. Weisbuch, Eur. Phys. J. B 38, 339 2004 . 10 H. Yizhaq, B. Portnov, and E. Meron, Envir. Plan. A 36, 149 20 G. Weisbuch, G. Deffuant, and F. Amblard, Physica A 353, 2004 . 555 2005 . 11 A. T. Bernardes, D. Stauffer, and J. Kertesz, Eur. Phys. J. B 21 G. Fath and M. Sarvary, Physica A 348, 611 2005 . 25, 123 2002 . 22 D. J. Watts and S. H. Strogatz, Nature London 393, 440 12 R. Axelrod, J. Conﬂict Resolut. 41, 203 1997 . 1998 . 046139-7

## Add a comment