Next Article in Journal
Continuous Variables Graph States Shaped as Complex Networks: Optimization and Manipulation
Next Article in Special Issue
Co-Evolution of Complex Network Public Goods Game under the Edges Rules
Previous Article in Journal
Low Cycle Fatigue Life Prediction Using Unified Mechanics Theory in Ti-6Al-4V Alloys
Previous Article in Special Issue
Modeling the Disorder of Closed System by Multi-Agent Based Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymmetric Contrarians in Opinion Dynamics

1
CEVIPOF—Centre for Political Research, Sciences Po and CNRS, 98 rue de l’Université, 75007 Paris, France
2
Laboratory of Physics, Kochi University of Technology, Tosa Yamada, Kochi 782-8502, Japan
*
Author to whom correspondence should be addressed.
Entropy 2020, 22(1), 25; https://doi.org/10.3390/e22010025
Submission received: 28 November 2019 / Revised: 18 December 2019 / Accepted: 20 December 2019 / Published: 24 December 2019

Abstract

:
Asymmetry in contrarian behavior is investigated within the Galam model of opinion dynamics using update groups of size 3 with two competing opinions A and B. Denoting x and y the respective proportions of A and B contrarians, four schemes of implementations are studied. The first scheme activates contrarians after each series of updates with probabilities x and y for agents holding respectively opinion A and B. Second scheme activates contrarians within the update groups only against global majority with probability x when A is the majority and y when B is the majority. The third scheme considers in-group contrarians acting prior to the local majority update against both local majority and minority opinions. The last scheme activates in-group contrarians prior to the local majority update but only against the local majority. The main result is the loss of the fifty–fifty attractor produced by symmetric contrarians. Producing a bit less contrarians on its own side than the other side becomes the key to win a public debate, which in turn can guarantee an election victory. The associated phase diagram of opinion dynamics is found to exhibit a rich variety of counterintuitive results.

1. Introduction

The modeling of opinion dynamics is one of the major topics of sociophysics studies [1,2,3,4,5,6,7,8,9,10,11,12]. Most discrete models consider the dynamics of two competing choices among identical agents, who all follow the same update rule while involved in a discussion with others. For any given initial distribution of choices among the agents, applying repeatedly the update leads to a homogenization of opinions with one of the two competing opinions getting shared by more and more agents. After some updates, all agents eventually adopt the same opinion.
With respect to entropy, the initial random distribution of opinions yields a non zero value of entropy with some associated disorder. However, repeated updates decrease the disorder gradually to end up with a totally ordered state characterized by a zero entropy.
In physical systems, above situation corresponds to a kind of spin system with an initial distribution of up and down spins connected by interactions at zero temperature. Minimizing the energy leads then to a fully ordered ground state with all spins either up or down characterized by a zero entropy. A non-zero temperature increases the disorder with a non zero entropy making the free-energy to be minimized instead of the internal energy.
While the temperature is an essential physical quantity, up to now, it has not been possible to identify a social equivalent. In parallel, the concept of contrarian agents, introduced in the modeling of opinion dynamics within the Galam model [13,14,15,16,17], does create a disorder in opinion dynamics, which share features with temperature effects [18,19]. In particular, contrarian agents produce a non zero entropy making attractors to stabilize mixtures of opinions with fixed proportions of them [20,21]. The concept of contrarian behavior has generated a great deal of works among physicists [22,23,24,25,26,27,28,29] within the active field of scoiophysics [30,31,32,33,34].
The contrarian feature is defined as an additional feature to “normal” behavior of agents following local majority rules for updating their respective opinions once they are discussing in small groups. Contrary to the normal agent, the contrarian behaves in two successive steps. First, it acts as a normal agent by expressing its choice. Then, once a majority has crystallized within the discussing group, while all normal agents adopt it, the contrarian chooses to adopt the opposite one. It may be preserving its initial opinions.
The contrarian shift is independent of either choice being the majority. It is also worth to notice that activation of the contrarian behavior is probabilistic. Given an agent in an updating group, it behaves normally with a probability ( 1 c ) and turns contrarian with probability c. Accordingly, the same agent may behave normally during an update and be a contrarian during another update. However, at every sequence of simultaneous updates, a proportion c of all agents are contrarians while the complementary proportion ( 1 c ) behaves normally.
In this paper, we extend the concept of contrarian by investigating the effects of making contrarian behavior asymmetric. Instead of one proportion c of contrarians, we consider two different proportions x and y to denote the respective proportions of A and B contrarians [35].
This hypothesis implies that the contrarian feature is now related to the choice itself rather than to the agent independently of its current holding choice. It is worth emphasizing that asymmetric contrarian behavior departs from physics where asymmetric temperature does not exist. At least, it has not been defined yet.
Our main result is that any non zero value of x y destroys the hung outcome with the attractor being shifted apart from the fifty–fifty simultaneously restoring a winner over a loser. The corresponding wining strategy is now to have a bit less contrarians on its side than on the other side.
In addition, contrarian asymmetry splits the contrarian activation between two different schemes. Implementing the contrarian shift within the update groups or after the updates before reshuffling yields different outcomes, contrary to the symmetric case where both outcomes are identical.
The rest of the paper is organized as follows. Section 2 reviews the effect of symmetric contrarians using the Galam model of opinion dynamics. Section 3 introduces asymmetric contrarians with four different schemes to activate the contrarian behavior. These schemes are respectively (i) post-update local asymmetric contrarians; (ii) post-update global asymmetric contrarians; (iii) in-group asymmetric contrarians with majority/minority contrarians; (iv) in-group asymmetric contrarians with in-group majority contrarians. The four different updates equations are discussed in Section 4 while Section 5 concludes the paper.

2. Symmetric Contrarians

2.1. The Core Model

The contrarian effect has been introduced using the core Galam model of opinion dynamics. Two opinions A and B are competing within a population of agents denoted floaters. A floater does hold an opinion it has chosen. When a public debate is launched, a floater tries to convince other agents to adopt its own choice. However, a floater agent does listen to the arguments supporting the other choice and may be convinced to shift its initial choice. The complicated and complex process of individual making opinion is mimicked by using a local majority rule with every discussing agent having one vote in favor of either one of the two choices.
The underlying dynamic of opinion deploys in three steps. First, agents are randomly distributed in local groups of size r. Then, a local majority rule is applied locally within each group to update individual opinions. It makes each group ending up unanimous in supporting either A or B. In the third step, agents are reshuffled before going back to step 1. Each series of the three steps produces a variation of support for both opinions. Starting from some initial proportions, p 0 and ( 1 p 0 ) for respectively A and B, the dynamics yields with p 0 p 1 p 2 .
The updating is iterated till an attractor of the dynamics is reached with p 0 p 1 p n p n + 1 p n . Defining the local majority function P ( p ) allows identifying the associated dynamics by solving the fixed points Equation,
P ( p ) = p ,
which yields the attractors p A , p B and tipping point p t driving the dynamics. Above p t ( p 0 > p t ) the support for A increases with p 0 < p 1 < p 2 < p n p A , while below p t ( p 0 < p t ) the support decreases as p 0 > p 1 > p 2 > p n p B .
However, in real situations, the vote takes place at some update l before reaching the attractor. The dynamics may also stop by itself with people shifting to discuss another issue after some number m of updates.
In case of groups of size r = 3 , 8 local configurations are obtained as exhibited in Table 1. The respective probability of occurrence, the corresponding local majority updates and the resulting weight contributing to opinion A are also given. The resulting proportion p t + 1 at t + 1 is,
p t + 1 P ( p t ) = p t 3 + 3 p t 2 ( 1 p t ) = 2 p t 3 + 3 p t 2 ,
where p t is the A proportion at t. The associated proportion for B is ( 1 p t + 1 ) with 0 p t 1 .
From Equation (1) two attractors are p A = 1 and p B = 0 separated by a tipping point p t = 1 2 . Therefore the dynamics makes the initial majority to win the debate.

2.2. Adding Contrarian Behavior

Including the psychological feature of contrarian behavior among floaters can have substantial effects on the dynamics [18]. A contrarian departs from a floater by shifting its opinion to oppose the local majority once it has revealed itself. The contrarian shift is independent of the majority opinion is either A or B. A contrarian behavior is activated randomly with some probability c, which accounts for the proportion of contrarian behavior within a social group. Contrarian behavior may be constant for some agents and temporary for others. However, at any given time, a proportion c of floaters are contrarian.
The value c is a given fixed external parameter, which is independent of the dynamics and satisfies 0 c 1 . Most real cases have c < 1 2 with c > 1 2 producing a systematically collective flips.
For size r = 3 , incorporating contrarian shifts once the local majority has been implemented, results in applying contrarian shifts on unanimous configurations A A A and B B B. Since contrarian shift is probabilistic, 0–3 contrarians may be present in configurations A A A and B B B. Table 2 shows the various cases with the new configurations, the weight contributions to opinion A and the respective probabilities of occurrence, which are ( 1 c ) 3 , 3 c ( 1 c ) 2 , 3 c 2 ( 1 c ) , c 3 .
From Table 2 Equation (2) becomes,
P c ( p ) = C A A A P ( p ) + C B B B [ 1 P ( p ) ] ,
with
C A A A =   ( 1 c ) 3 + 2 3 3 c ( 1 c ) 2 + 1 3 3 c 2 ( 1 c ) =   ( 1 c ) ,
and
C B B B = 1 3 3 c ( 1 c ) 2 + 2 3 3 c 2 ( 1 c ) + c 3 = c .
Using Equations (4) and (5), the update Equation (3) writes as,
P c ( p ) = ( 1 c ) p 3 + 3 p 2 ( 1 p ) + c 3 p ( 1 p ) 2 + ( 1 p ) 3 ,
which shows that activation of contrarians could be also applied identically to individual agents out of the update groups. Equation (6) simplifies to,
P c ( p ) = ( 1 2 c ) 2 p t 3 + 3 p t 2 + c ,
which yields two attractors p A , p B and a tipping point p t defined as,
p A = ( 1 2 c ) + 1 8 c + 12 c 2 2 ( 1 2 c ) , p B = ( 1 2 c ) 1 8 c + 12 c 2 2 ( 1 2 c ) , p t = 1 2 ,
provided 0 c 1 6 0.1667 and 1 2 c 1 to ensure the existence of the square root.
However for 1 2 < c < 1 we have p A > 1 and p B < 0 making them non physical. At c= 1 6 and c = 1 2 , p A = p B = 1 2 . Variations of p A and p B as a function of c are shown in Figure 1. Indeed, for c 1 6 , p A and p B merge at the tipping point p t = 1 2 turning it to the unique attractor of the dynamics. The two regimes are shown in Figure 2 and Figure 3 with c = 0.08 p A = 0.11 ,   p B = 0.89 ,   p t = 0.50 and c = 0.20 p A = p B = p t = 0.50 . For the full range c > 1 6 any initial value p 0 ends up at 1 2 . For c = 1 2 we have p t + 1 = 1 2 for whatever value of p t .
The different regimes of opinion dynamics driven by a non zero proportion of contrarian behavior are exhibited in the series of Figure 2, Figure 3 and Figure 4.

3. Making Contrarians Asymmetric

Above symmetric contrarian behavior could be connected to the assumption that contrarianism is an individual feature produced by some social environment. It is individuals a psychological feature to oppose the surrounding leading opinion. It is, therefore, independent of opinions themselves, which justifies making contrarian behavior symmetrical with respect to A and B [18].
However, contrarianism could be also identified as arising from the cultural context in which a given opinion is rooted. Behind an opinion choice about a given issue would stand a global frame of thinking. In this case, contrarianism could become asymmetric, some social frames favoring more than others, an individual inclination to challenge the leading opinion adopting the minority view. Denoting x the probability that an A holding contrarian flips to B and y the probability that a B holding contrarian flips to A, we study the effect of such an asymmetry on the associated dynamics of opinion.
It is worth to notice that with symmetrical contrarians, the activation of contrarian shift was found identical when applied within the update groups or after the updates during the reshuffling prior to a new update. However, with asymmetric contrarians, while the equivalence is still valid for unanimity configurations A A A and B B B, it does not hold for configurations A A B and A B B.
Having x y makes simple analytical solving not possible to determine the values of the various attractors and tipping points as when x = y . However, it will be seen that the two regimes associated with symmetric contrarians, one with two attractors and a tipping point and another one with a single attractor, are still prevailing but with the loss of symmetry for the location of attractors and tipping points.
Accordingly, in the following of the paper, to keep the presentation clear, both regimes of the dynamics will be illustrated choosing a specific set of values for x and y, which exhibit the main shape of the corresponding phase diagram.

3.1. Post-Update Local Asymmetric Contrarians

When contrarians are activated after the update, Equation (6) is modified as,
P x , y 1 ( p ) = ( 1 x ) p 3 + 3 p 2 ( 1 p ) + y 3 p ( 1 p ) 2 + ( 1 p ) 3 ,
which shows that x y main effect is twofold removing the hung election scenario (single attractor at fifty–fifty) yielded by Equation (8) and turning the dynamics unbalanced. The shift of p A (p B ) towards p t becomes faster than the shift of p B (p A ) when x > y ( x < y ) since now p t is also moving departing from 1 2 to lower (higher) values.
For instance in the case x = 0.08 and y = 0.10 , Equation (9) yields p A = 0.15 ,   p t = 0.46 ,   p B = 0.89 . The coalescence phenomena observed with symmetrical contrarians is still valid but now it occurs at a value different from 1 2 as seen from the case x = 0.20 ,   y = 0.22 yielding p A = p t = p B = 0.57 . These two cases are exhibited in Figure 5.

3.2. Post-Update Global Asymmetric Contrarians

A variant of post-update contrarians has been studied in the symmetrical case considering the majority status of opinion at the global level instead of the local one [36]. Namely, a contrarian flip from A to B takes place with probability c only when p > 1 2 and from B to A with the same probability c only when p < 1 2 . In this case, a chaotic behavior is found around an attractor at fifty/fifty.
Within this global scheme, considering asymmetric contrarians makes a contrarian flip from A to B to occur with probability x only p > 1 2 . When p < 1 2 contrarian flips occur only from B to A with probability y. The associated update Equation writes as,
p > 1 2 P x 2 a ( p ) = ( 1 x ) p 3 + 3 p 2 ( 1 p ) , p < 1 2 P y 2 b ( ( p ) = p 3 + 3 p 2 ( 1 p ) + y ( 1 p ) 3 + 3 p ( 1 p ) 2 ,
which still produces a chaotic behavior but now at a value of p t 1 2 instead of p t = 1 2 .
Figure 6 illustrates the behavior from Equation (11) for the case x = 0.20 , y = 0.18 . The attractor p t = 1 2 has split in two limiting points located at p t a = 0.4 and p t b = 0.59 .

3.3. In-Group Asymmetric Contrarians

Activation of contrarians inside the update group opens two different possible schemes since A and B agents are not in the same situation as seen with the case A A B c (A c B B). The first scheme considers that B c (A c ) being already minority against A (B), keeps the same opinion B (A). On the other hand, the second scheme considers that B c (A c ) being minority, its contrarian feature is not activated and thus the agent follows the local majority turning A (B). However, for configurations A A A and B B B, there is no change with respect to the post-update scheme.

3.3.1. In-Group Majority/Minority Contrarians

The first scheme was investigated in [35] with for instance the update A A B c A A B. All configuration updates with independent contrarians are listed in Table 3 for groups of size 3 for the configuration A A B with 0–3 contrarians.
The update Equation (2) writes now,
P x , y 3 ( p ) = C A A A ( x ) p 3 + C A A B ( x , y ) 3 p 2 ( 1 p ) + C B B B ( y ) ( 1 p ) 3 + C A B B ( x , y ) 3 p ( 1 p ) 2 ,
with from Table 2,
C A A A ( x ) = ( 1 x ) 3 + 2 x ( 1 x ) 2 + x 2 ( 1 x ) + 0 x 3 = ( 1 x ) ,
and from Table 3,
C A A B ( x , y ) = ( 1 x ) 2 ( 1 y ) + 4 3 x ( 1 x ) ( 1 y ) + 2 3 ( 1 x ) 2 y + 1 3 x 2 ( 1 y ) + 2 3 x ( 1 x ) y = 1 2 3 x 1 3 y ,
By symmetry we get
C B B B ( y ) = y
and
C B B A ( x , y ) = 2 3 y + 1 3 x ,
to yield,
P x , y 3 ( p ) = ( 1 x ) p 3 + y ( 1 p ) 3 + 1 2 3 x 1 3 y ) 3 p 2 ( 1 p ) + 2 3 y + 1 3 x 3 p ( 1 p ) 2 ,
which writes,
P x , x 3 ( p ) = ( 1 x ) p 3 + 3 p 2 ( 1 p ) + x ( 1 p ) 3 + 3 p ( 1 p ) 2 ,
when x = y thus recovering Equation (6) with x = c .
For x = 0.08 and y = 0.10 , Equation (16) yields p A = 0.14 , p t = 0.48 , p B = 0.89 . The coalescence phenomena observed with symmetrical contrarians is still valid but now it occurs at a value different from 1 2 as seen from the case x = 0.20 , y = 0.22 yielding p A = p t = p B = 0.54 . These two cases are exhibited in Figure 7.

3.3.2. In-Group Majority Contrarians

While above, contrarians are always activated within the update group, we could consider that contrarianism is activated at a contrarian agent only when this agent finds itself in an initial local majority in a group. In this case, a contrarian B c (A c ) in a group A A B c (B B A c ) keeps dormant and acts as a regular floater following the local majority rule. According to this scheme A A B c (B B A c ) → A A A (B B B) instead of A A B (B B A) above. Table 3 is thus modified as shown in Table 4. Accordingly, Equation (13) becomes,
C A A B ( x , y ) = ( 1 x ) 2 ( 1 y ) + 4 3 x ( 1 x ) ( 1 y ) + ( 1 x ) 2 y + 1 3 x 2 ( 1 y ) + 4 3 x ( 1 x ) y + 1 3 x 2 y = 1 2 3 x ,
and by symmetry Equation (15) becomes,
C B B A ( x , y ) = 2 3 y ,
to yield,
P x , y 4 ( p ) = ( 1 x ) p 3 + y ( 1 p ) 3 + 1 2 3 x ) 3 p 2 ( 1 p ) + 2 3 y 3 p ( 1 p ) 2 .
For x = 0.08 and y = 0.10 , Equation (20) yields p A = 0.13 , p t = 0.47 , p B = 0.90 . The coalescence phenomena observed with symmetrical contrarians is still valid but now it occurs at a value different from 1 2 as seen from the case x = 0.20 , y = 0.22 yielding p A = p t = p B = 0.64 . These two cases are exhibited in Figure 8.

4. Discussion

Asymmetric contrarians were shown to generate four different schemes to account for the effect on the associated opinion dynamics. For update groups of size three, the various components to the update equations are listed in Table 5. Combined together, the respective update equations are,
P x , y 1 ( p ) = ( 1 x ) p 3 + 3 p 2 ( 1 p ) + y 3 p ( 1 p ) 2 + ( 1 p ) 3 , P x , 2 a ( p ) = ( 1 x ) p 3 + 3 p 2 ( 1 p ) p > 1 2 , P y 2 b ( p ) = p 3 + 3 p 2 ( 1 p ) + y 3 p ( 1 p ) 2 + ( 1 p ) 3 p < 1 2 , P x , y 3 ( p ) = ( 1 x ) p 3 + 1 2 3 x 1 3 y ) 3 p 2 ( 1 p ) + 2 3 y + 1 3 x 3 p ( 1 p ) 2 + y ( 1 p ) 3 , P x , y 4 ( p ) = ( 1 x ) p 3 + 1 2 3 x ) 3 p 2 ( 1 p ) + 2 3 y 3 p ( 1 p ) 2 + y ( 1 p ) 3 .
The common feature of four schemes is the breaking of the dynamics symmetry. In particular, the fifty-fifty attractor is shifted to a value different from fifty-fifty yielding a majority and a minority. Comparing schemes 1–4 shows some minor quantitative differences while scheme 4 is rather different as expected since it is connected to the global majority.
It is worth to notice that while P x , x 1 ( p ) = P x , x 3 ( p ) P x , x 4 ( p ) since indeed scheme 4 activates contrarian behavior only when an agent is at minority within an updating group before the update is completed.
To illustrate the differences between the various schemes 1–4, we show several cases of variations of the associated fixed points as a function of x for a fixed value of y. Figure 9 compares schemes 1 and 3 and 1 and 4 for y = 0.12 . Figure 10 compares schemes 1, 3 and 4 for y = 0 and y = 0.10 while Figure 11 does the same for y = 0.12 and y = 0.21 .
The left part of Figure 9 and right part of Figure 10 show that for scheme 3 the attractor p B increases slightly its value with increasing x, which looks paradoxical since contrarians are expected to decrease the overall support of their corresponding opinion. Indeed, scheme 3 maintains a contrarian being in a minority position to keep on its opinion. It implies that when support for A (B) is low, many configurations have A (B) in minority and thus preserves A (B) opinion. More A (B) contrarians for low p thus increases slightly the support for A (B).
Figure 12 provides an illustration of this paradox exhibiting the balance between the loss and gain from contrarian behavior as a function of p, from configurations A A B and A B B for schemes 3 (right) and 4 (left). The cases x = 0 , y = 0.10 and x = 0.12 , y = 0.10 are compared. It is seen that increasing x reduces the net balance at low p for scheme 4 as expected. However, the opposite occurs for scheme 3.
Indeed, looking at the 4 schemes associated with the update equations listed in Equation (21), it is seen that only P x , y 3 ( p ) has a positive x contribution from the minority configurations A B B, which explains the counterintuitive finding that at low values of p, increasing x increases a bit the A support. Figure 12 and Figure 13 shows both the amplitude of the negative contribution 2 3 x 3 p 2 ( 1 p ) from A contrarians (proportion x) and the amplitude of the positive contribution 1 3 x 3 p ( 1 p ) 2 for x = 0.20 . The positive contribution is seen to be bigger than the negative one for p < 1 3 .

5. Conclusions

We have introduced asymmetry in contrarian behavior with respect to two competing opinions A and B to account for inhomogeneity in the distribution of psychological traits within a given population in connexion to the ability to differ from the group. It hints at contrarian behavior being connected to the underlying cultural frame of social groups. Accordingly, symmetric contrarians are expected to prevail within homogeneous populations.
Asymmetry was shown to generate several different schemes for the contrarian activation. A rich variety of changes has been found about the values of the attractors and tipping point [37].
However, the differences by themselves are not of large amplitude establishing the robustness of contrarian behavior. Nevertheless, it is worth to emphasize that as soon as asymmetry is turned on, the salient feature of having a single attractor located precisely at fifty-fifty is waived off. However, now a slightly lower proportion of contrarian on its own side than on the other competing side is the key to eventually win a public debate which, in turn, can transform into an election victory. This result allows us to envision novel winning strategies focusing on the lowering of the proportion of contrarians on its own side.

Author Contributions

All authors have contributed equally to this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Noorazar, H.; Vixie, K.R.; Talebanpour, A.; Hu, Y. From classical to modern opinion dynamics. arXiv 2019, arXiv:1909.12089v2. [Google Scholar]
  2. Calvão, A.M.; Ramos, M.; Anteneodo, C. Role of the plurality rule in multiple choices. J. Stat. Mech. 2016, 2016, 023405. [Google Scholar] [CrossRef] [Green Version]
  3. Ajmone Marsan, G.; Bellomo, N.; Gibelli, L. Stochastic evolutionary differential games toward a systems theory of behavioral social dynamics. Math. Models Methods Appl. Sci. 2016, 26, 1051–1093. [Google Scholar] [CrossRef]
  4. Carbone, G.; Giannoccaro, I. Model of human collective decision-making in complex environments. Eur. Phys. J. 2015, 88, 339. [Google Scholar] [CrossRef] [Green Version]
  5. Sznajd-Weron, K.; Szwabiński, J.; Weron, R. Is the Person-Situation Debate Important for Agent-Based Modeling and Vice-Versa? PLoS ONE 2014, 9, e112203. [Google Scholar] [CrossRef] [PubMed]
  6. Javarone, M.A. Networks strategies in election campaigns. J. Stat. Mech. Theory Exp. 2014, P08013. [Google Scholar] [CrossRef] [Green Version]
  7. Lanchier, N.; Taylor, N. Galam’s bottom-up hierarchical system and public debate model revisited. Adv. Appl. Prob. 2013, 47, 668–692. [Google Scholar] [CrossRef]
  8. Goncalves, S.; Laguna, M.F.; Iglesias, J.R. Why, when, and how fast innovations are adopted. Eur. Phys. J. 2012, 85, 192. [Google Scholar] [CrossRef] [Green Version]
  9. Mobilia, M. Fixation and polarization in a three-species opinion dynamics model. Eur. Phys. Lett. 2011, 95, 50002. [Google Scholar] [CrossRef]
  10. Ellero, A.; Fasano, G.; Sorato, A. A modified Galam’s model for word-of-mouth information exchange. Physica A 2009, 388, 3901–3910. [Google Scholar] [CrossRef]
  11. Martins, A.C.R. Mobility and Social Network Effects on Extremist Opinions. Phys. Rev. 2008, 78, 036104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Diep, H.T.; Kaufman, M.; Kaufman, S. Dynamics of two-group conflicts: A statistical physics model. Physica A 2017, 469, 183–199. [Google Scholar] [CrossRef] [Green Version]
  13. Galam, S. Majority rule, hierarchical structures and democratic totalitarism: A statistical approach. J. Math. Psychol. 1986, 30, 426–434. [Google Scholar] [CrossRef]
  14. Galam, S.; Chopard, B.; Masselot, A.; Droz, M. Competing Species Dynamics. Eur. Phys. J. 1998, 4, 529–531. [Google Scholar] [CrossRef] [Green Version]
  15. Galam, S. Minority Opinion Spreading in Random Geometry. Eur. Phys. J. 2002, 25, 403–406. [Google Scholar] [CrossRef] [Green Version]
  16. Galam, S. The dynamics of minority opinion in democratic debate. Physica A 2004, 336, 56–62. [Google Scholar] [CrossRef]
  17. Galam, S. Heterogeneous beliefs, segregation, and extremism in the making of public opinions. Phys. Rev. E 2005, 71, 046123. [Google Scholar] [CrossRef] [Green Version]
  18. Galam, S. Contrarian deterministic effects on opinion dynamics, the hung elections scenario. Physica 2004, 333, 453. [Google Scholar] [CrossRef] [Green Version]
  19. Galam, S. From 2000 Bush-Gore to 2006 Italian elections: Voting at fifty-fifty and the contrarian effect. Qual. Quant. 2007, 41, 579–589. [Google Scholar] [CrossRef] [Green Version]
  20. Furioli, G.; Pulvirenti, A.; Terraneo, E.; Toscani, G. Fokker-Planck equations in the modelling of socio-economic phenomena. Math. Mod. Meth. Appl. Sci. 2017, 27, 115–158. [Google Scholar] [CrossRef] [Green Version]
  21. Albi, G.; Pareschi, L.; Toscani, G.; Zanella, M. Recent advances in opinion modeling: Control and social influence. In Active Particles, Volume 1: Theory, Models, Applications; Bellomo, N., Degond, P., Tadmor, E., Eds.; Springer: Basel, Switzerland, 2017; Chapter 2; pp. 49–98. [Google Scholar]
  22. Singh, P.; Sreenivasan, S.; Szymanski, B.K.; Korniss, G. Competing effects of social balance and influence. Phys. Rev. 2016, 93, 042306. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Cheon, T.; Morimoto, J. Balancer effects in opinion dynamics. Phys. Lett. 2016, 380, 429–434. [Google Scholar] [CrossRef] [Green Version]
  24. Siedlecki, P.; Szwabinski, J.; Weron, T. The interplay between conformity and anticonformity and its polarizing effect on society. JASSS 2016, 19, 9. [Google Scholar] [CrossRef] [Green Version]
  25. Lee, E.; Holme, P.; Lee, S.H. Modeling the dynamics of dissent. Physica A 2017, 486, 262. [Google Scholar] [CrossRef] [Green Version]
  26. Weisbuch, G. From Anti-Conformism to Extremism. J. Artif. Soc. Soc. Simul. 2015, 18, 1. [Google Scholar] [CrossRef] [Green Version]
  27. Crokidakis, N.; de Oliveira, P.M.C. Inflexibility and independence: Phase transitions in the majority-rule model. Phys. Rev. 2015, 92, 062122. [Google Scholar] [CrossRef] [Green Version]
  28. Bagnoli, F.; Rechtman, R. Bifurcations in models of a society of reasonable contrarians and conformists. Phys. Rev. 2015, 92, 042913. [Google Scholar] [CrossRef] [Green Version]
  29. Chacoma, A.; Zanette, D.H. Critical phenomena in the spreading of opinion consensus and disagreement. Papers Phys. 2014, 6, 060003. [Google Scholar] [CrossRef] [Green Version]
  30. Galam, S. Sociophysics: A Physicist’s Modeling of Psycho-Political Phenomena; Springer: New York, NY, USA, 2012. [Google Scholar]
  31. Castellano, C.; Fortunato, S.; Loreto, V. Statistical physics of social dynamics. Rev. Mod. Phys. 2009, 81, 591–646. [Google Scholar] [CrossRef] [Green Version]
  32. Galam, S. Sociophysics: A review of Galam models. Int. J. Mod. Phys. C 2008, 19, 409–440. [Google Scholar] [CrossRef]
  33. Chakrabarti, B.K.; Chakraborti, A.; Chatterjee, A. (Eds.) Econophysics and Sociophysics: Trends and Perspectives; Wiley-VCH Verlag: Weinheim, Germany, 2006. [Google Scholar]
  34. Sociophysics: A personal testimony. Physica A 2002, 336, 1375–1385.
  35. Jacobs, F.; Galam, S. Tow opinions dynamics generated by inlfexibles and non-contrarians and contrarians floaters. Adv. Complex Syst. 2019, 22, 1950008. [Google Scholar] [CrossRef] [Green Version]
  36. Borghesi, C.; Galam, S. Chaotic, staggered, and polarized dynamics in opinion forming: The contrarian effect. Phys. Rev. 2006, 73, 066118. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Pyko, N.S.; Pyko, S.A.; Markelov, O.A.; Karimov, A.I.; Butusov, D.N.; Zolotukhin, Y.V.; Uljanitski, Y.D.; Bogachev, M.I. Assessment of cooperativity in complex systems with non-periodical dynamics: Comparison of five mutual information metrics. Physica A 2018, 503, 1054–1072. [Google Scholar] [CrossRef]
Figure 1. Variations of p A (left) and p B (right) as a function of c are shown. The two attractors exist only in the range 0 c 1 .
Figure 1. Variations of p A (left) and p B (right) as a function of c are shown. The two attractors exist only in the range 0 c 1 .
Entropy 22 00025 g001
Figure 2. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0 (left) and c = 0.08 (right). The arrows show the dynamics from p 0 = 0.4 for several iterations. Associated attractors are p A = 0 ,   p B = 1 ,   p t = 0.50 for the first case and p A = 0.11 ,   p B = 0.89 ,   p t = 0.50 for the second case.
Figure 2. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0 (left) and c = 0.08 (right). The arrows show the dynamics from p 0 = 0.4 for several iterations. Associated attractors are p A = 0 ,   p B = 1 ,   p t = 0.50 for the first case and p A = 0.11 ,   p B = 0.89 ,   p t = 0.50 for the second case.
Entropy 22 00025 g002
Figure 3. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0.2 (left) and c = 0.5 (right). The arrows show the dynamics from p 0 = 0.1 (left) for several iterations. For the second case, any initial value p 0 leads to the single attractor within one iteration. For both cases one single attractor p A = p B = p t = 0.50 drives the dynamics.
Figure 3. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0.2 (left) and c = 0.5 (right). The arrows show the dynamics from p 0 = 0.1 (left) for several iterations. For the second case, any initial value p 0 leads to the single attractor within one iteration. For both cases one single attractor p A = p B = p t = 0.50 drives the dynamics.
Entropy 22 00025 g003
Figure 4. As soon as c > 0.5 the dynamics becomes alternating. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0.75 (left) and c = 0.9 (right). The arrows show the dynamics from p 0 = 0.1 (left) and p 0 = 0.35 (right) for several iterations. For the first case, any initial value p 0 leads to the single attractor p A = p B = p t = 0.50 but through jumps in respective supports for both opinions. For the second case, the iteration produces jumps around the tipping point p t = 0.50 but now alternating towards p A = 0.15 and p B = 0.85 .
Figure 4. As soon as c > 0.5 the dynamics becomes alternating. Using Equation (7) the expression p t + 1 is shown as a function of p t for c = 0.75 (left) and c = 0.9 (right). The arrows show the dynamics from p 0 = 0.1 (left) and p 0 = 0.35 (right) for several iterations. For the first case, any initial value p 0 leads to the single attractor p A = p B = p t = 0.50 but through jumps in respective supports for both opinions. For the second case, the iteration produces jumps around the tipping point p t = 0.50 but now alternating towards p A = 0.15 and p B = 0.85 .
Entropy 22 00025 g004
Figure 5. Using Equation (9) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.15 , p t = 0.46 , p B = 0.89 and p A = p t = p B = 0.57 .
Figure 5. Using Equation (9) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.15 , p t = 0.46 , p B = 0.89 and p A = p t = p B = 0.57 .
Entropy 22 00025 g005
Figure 6. Using Equation (11) the expression p t + 1 as a function of p t is shown for x = 0.20 , y = 0.18 (top). The associated evolution from p 0 = 0.26 is exhibited (bottom). The attractor p t = 1 2 has split in two limiting points located at p t a = 0.4 and p t b = 0.59 .
Figure 6. Using Equation (11) the expression p t + 1 as a function of p t is shown for x = 0.20 , y = 0.18 (top). The associated evolution from p 0 = 0.26 is exhibited (bottom). The attractor p t = 1 2 has split in two limiting points located at p t a = 0.4 and p t b = 0.59 .
Entropy 22 00025 g006
Figure 7. Using Equation (16) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.14 , p t = 0.48 , p B = 0.89 and p A = p t = p B = 0.54 .
Figure 7. Using Equation (16) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.14 , p t = 0.48 , p B = 0.89 and p A = p t = p B = 0.54 .
Entropy 22 00025 g007
Figure 8. Using Equation (20) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.13 , p t = 0.47 , p B = 0.90 and p A = p t = p B = 0.64 .
Figure 8. Using Equation (20) the expression p t + 1 as a function of p t is shown for x = 0.08 , y = 0.10 (top) and x = 0.20 , y = 0.22 (bottom). The associated fixed points are respectively p A = 0.13 , p t = 0.47 , p B = 0.90 and p A = p t = p B = 0.64 .
Entropy 22 00025 g008
Figure 9. Fixed points of scheme 1 and 3 (left) and schemes 1 and 4 (right) as a function of x for y = 0.12 . The two dotted vertical lines show the area with a tipping point between two attractors. Outside the lines, the dynamics is driven by one single attractor.
Figure 9. Fixed points of scheme 1 and 3 (left) and schemes 1 and 4 (right) as a function of x for y = 0.12 . The two dotted vertical lines show the area with a tipping point between two attractors. Outside the lines, the dynamics is driven by one single attractor.
Entropy 22 00025 g009
Figure 10. Fixed points of schemes 1–4 as a function of x for y = 0 (left) and y = 0.10 (right).
Figure 10. Fixed points of schemes 1–4 as a function of x for y = 0 (left) and y = 0.10 (right).
Entropy 22 00025 g010
Figure 11. Fixed points of schemes 1 and 3 and 4 as a function of x for y = 0.12 (left) and y = 0.21 (right).
Figure 11. Fixed points of schemes 1 and 3 and 4 as a function of x for y = 0.12 (left) and y = 0.21 (right).
Entropy 22 00025 g011
Figure 12. Net balance between gain minus loss from contrarian behavior for scheme 4 and 3 when x = 0 and y = 0.10 (left) and x = 0.12 and y = 0.10 (right). More A contrarians for low values of p increases the net gain for A.
Figure 12. Net balance between gain minus loss from contrarian behavior for scheme 4 and 3 when x = 0 and y = 0.10 (left) and x = 0.12 and y = 0.10 (right). More A contrarians for low values of p increases the net gain for A.
Entropy 22 00025 g012
Figure 13. Amplitude of the negative contribution 2 3 x 3 p 2 ( 1 p ) (Loss) from A contrarians (proportion x) and the amplitude of the positive contribution 1 3 x 3 p ( 1 p ) 2 (gain) for x = 0.20 . The positive contribution is seen to be bigger than the negative one for p < 1 3 .
Figure 13. Amplitude of the negative contribution 2 3 x 3 p 2 ( 1 p ) (Loss) from A contrarians (proportion x) and the amplitude of the positive contribution 1 3 x 3 p ( 1 p ) 2 (gain) for x = 0.20 . The positive contribution is seen to be bigger than the negative one for p < 1 3 .
Entropy 22 00025 g013
Table 1. Configuration lists the various distribution of A and B in a group of size r = 3 for a population with a proportion p in favor of A and ( 1 p ) in favor of B. Update shows the outcome of opinions once majority rule has been applied. Weight gives the corresponding contribution to opinion A. Probability gives the configuration probability.
Table 1. Configuration lists the various distribution of A and B in a group of size r = 3 for a population with a proportion p in favor of A and ( 1 p ) in favor of B. Update shows the outcome of opinions once majority rule has been applied. Weight gives the corresponding contribution to opinion A. Probability gives the configuration probability.
ConfigurationUpdateWeightProbability
A A A →A A A1 p 3
A A B, A B A, B A A →A A A1 3 p 2 ( 1 p )
A B B, B A B, B B A →B B B0 3 p ( 1 p ) 2
B B B →B B B0 ( 1 p ) 3
Table 2. Distribution of 0–3 contrarians in configuration A A A once the update has been completed. The associate weight contribution to A and the probability of occurrence are given.
Table 2. Distribution of 0–3 contrarians in configuration A A A once the update has been completed. The associate weight contribution to A and the probability of occurrence are given.
ContrariansConfigurationPost-UpdateWeightProbability
0A A A →A A A1 ( 1 c ) 3
1A c A A ( × 3 ) →B A A 2 3 3 c ( 1 c ) 2
2A c A c A ( × 3 ) →B B A 1 3 3 c 2 ( 1 c )
3A c A c A c B B B0 c 3
0B B B →B B B0 ( 1 c ) 3
1B c B B ( × 3 ) →A B B 1 3 3 c ( 1 c ) 2
2B c B c B ( × 3 ) →A A B 2 3 3 c 2 ( 1 c )
3B c B c B c A A A1 c 3
Table 3. Distribution of 0, 1, 2, 3 contrarians within a configuration A A B. Associated updates are shown with the weight contribution to A and the probability of occurrence.
Table 3. Distribution of 0, 1, 2, 3 contrarians within a configuration A A B. Associated updates are shown with the weight contribution to A and the probability of occurrence.
ContrariansConfigurationIn-Group Majority/MinorityWeightProba
0A A B →A A A1 ( 1 x ) 2 ( 1 y )
1A c A B →B A A 2 3 2 x ( 1 x ) ( 1 y )
1A A B c A A B 2 3 ( 1 x ) 2 y
2A c A c B →B B A 1 3 x 2 ( 1 y )
2A c A B c B A B 1 3 2 x ( 1 x ) y
3A c A c B c B B B0 x 2 y
Table 4. Distribution of 0, 1, 2 contrarians within a configuration A A B. Associated updates are shown with the weight contribution to A and the probability of occurrence.
Table 4. Distribution of 0, 1, 2 contrarians within a configuration A A B. Associated updates are shown with the weight contribution to A and the probability of occurrence.
ContrariansConfigurationIn-Group Majority/MinorityWeightProba
0A A B →A A A1 ( 1 x ) 2 ( 1 y )
1A c A B →B A A 2 3 2 x ( 1 x ) ( 1 y )
1A A B c A A A1 ( 1 x ) 2 y
2A c A c B →B B A 1 3 x 2 ( 1 y )
2A c A B c B A A 2 3 2 x ( 1 x ) y
3A c A c B c B B A 1 3 x 2 y
Table 5. Rescaling factors for the probability of occurrence of configurations A A A, A A B, A B B, B B B for the various contrarian schemes.
Table 5. Rescaling factors for the probability of occurrence of configurations A A A, A A B, A B B, B B B for the various contrarian schemes.
ConfigurationA A AA A BA B BB B B
Probability p 3 3 p 2 ( 1 p ) 3 p ( 1 p ) 2 ( 1 p ) 3
1post-update local majority ( 1 x ) ( 1 x ) yy
2apost-update global majority p > 1 2 ( 1 x ) ( 1 x ) 00
2bpost-update global majority p < 1 2 11yy
3In-groups local majority/minority ( 1 x ) ( 1 2 3 x 1 3 y ) ( 2 3 y + 1 3 x ) y
4In-groups local majority ( 1 x ) ( 1 2 3 x ) 2 3 y y

Share and Cite

MDPI and ACS Style

Galam, S.; Cheon, T. Asymmetric Contrarians in Opinion Dynamics. Entropy 2020, 22, 25. https://doi.org/10.3390/e22010025

AMA Style

Galam S, Cheon T. Asymmetric Contrarians in Opinion Dynamics. Entropy. 2020; 22(1):25. https://doi.org/10.3390/e22010025

Chicago/Turabian Style

Galam, Serge, and Taksu Cheon. 2020. "Asymmetric Contrarians in Opinion Dynamics" Entropy 22, no. 1: 25. https://doi.org/10.3390/e22010025

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop