Opinion Particles: Classical Physics and
arXiv:1307.3304v1 [physics.gen-ph] 12 Jul 2013
Opinion Dynamics
André C. R. Martins
NISC - EACH
Universidade de São Paulo, Brazil
July 15, 2013
Abstract
A model for Opinion Particles, based on Bayesian-inspired mod-
els of Opinion Dynamics such as the CODA model is presented. By
extending the discrete time characteristic of those models to continu-
ous time, a theory for the movement of opinion particles is obtained,
based only on inference ideas. This will allow inertia to be obtained
as a consequence of an extended CODA model. For the general case,
we will see that the likelihoods are associated with variables such as
velocity and acceleration of the particles. Newtonian forces are easily
defined and the relationship between a force and the equivalent like-
lihood provided. The case of the harmonic oscillator is solved as an
example, to illustrate clearly the relationship between Opinion Parti-
cles and Mechanics. Finally, possible paths to apply these results to
General Relativity are debated.
1 Introduction
The problem with integrating Quantum Mechanics and General Relativity
is so widely known that no introduction to the theme is actually required.
Both theories are very successful in the cases where we know they must be
applied. Yet, no agreement has been reached by the scientific community on
a good model for Quantum Gravity, despite the fact that efforts do exist, as,
per example, String Theory [1] or Loop Quantum Gravity [2]. A framework
1
that would allow both Quantum Mechanics and General Relativity to be
obtained from the same basic principles could, in principle, provide a impor-
tant new development in the efforts to understand how both theories could
be formulated in an unifying way. If that framework were to be based on
logical and informational principles, that would be a very nice bonus.
Efforts to understand Physics from an inference point of view are not new
and have provided some very interesting results [3, 4, 5, 6, 7]. They use en-
tropic principles based on the information one receives from experiments and
are compatible with Bayesian Statistics, in the cases where the information
comes as data. It is interesting to notice that Bayesian Statistics can actu-
ally be seen as an extension of Classical Logic where statements are not just
true or false, but can have different amounts of plausibility assigned to them,
according to what is currently known [8, 3]. That means that entropic and
informational approaches use very basic first principles of reasoning. While
these are interesting advances, they tend to apply inference principles only
to understand how an observer would interpret the results of an experiment.
In this paper, Bayesian Logic will be shown to possibly allow for a deeper
understanding of theoretical Physics than that allowed by just applying it to
the inferential process. The idea of the paper is to provide initial ideas in that
direction and to present far more questions than answers, while also illus-
trating how the general principles presented herein can be applied. This will
be accomplished by expanding on existing models from Opinion Dynamics.
In Opinion Dynamics problems [10, 11, 12, 13, 14, 15, 16], a framework
based on Bayesian methods has been proposed [17] as an extension of the
ideas presented in the Continuous Opinions and Discrete Actions (CODA)
model[15, 18]. It is interesting to notice that this framework is able to gen-
erate many of the other Opinion Dynamics models, both discrete [19] and
continuous [26], as particular or limit cases. Those are the ideas that will be
expanded here into a full model for opinion particles.
The paper is structured like this: First, a brief review of CODA model
will be presented. CODA is much simpler if we transform the probabilities p
to log-odds ν and, therefore, most of the analysis will be conducted using the
variables ν. In order to obtain a more realistic dynamics, I will proceed by
extending CODA to continuous time, thus defining opinion particles. These
particles move as a consequence of an inferential process based on the in-
formation they receive, like, per example, the position (and, possibly, other
variables) of the other particles. We will see that inertia can be obtained as
a consequence of this continuous time CODA, with no need to introduce any
2
kind of space structure. This means that Mach’s principle [21, 22, 23, 24, 25]
can be applied to Newtonian Mechanics, as inertia becomes a consequence
of interaction, with no need for introducing privileged inertial frames.
The relationship between likelihoods and the forces that act on the opin-
ion particles will be obtained and we will see how the harmonic oscillator
can be described in this framework. Through an inference where the full
information about the location of the other particles is used, I will show that
the results of the harmonic oscillator can be obtained from Beta function
likelihoods.
Finally, the problem of the geometric structure of the space of ν coordi-
nates will be briefly discussed. We will see it makes sense to identify ν with
space-time coordinates and obtain the equation that gives us the likelihood
equivalent to a parallel transport in a relativistic gravitational problem.
One warning is necessary before proceeding. While the name Opinion
Particles in a natural name, given they are inspired in Opinion Dynamics
models, it does not mean that the particles have any kind of opinion. The
name just reflects the fact that, as information about other particle arrives,
the first particle uses that information in some unknown process to determine
its trajectory. That this use can be compatible with principles of rationality
is an interesting fact, with no explanation or speculation attached to it.
2 CODA and Inertia
2.1 CODA
CODA model [15, 18] was developed to explore how human agents could
influence each other, changing their opinions when observing the choices of
other agents. It assumes that, in a situation where there are two possible
choices (or actions), each agent assigns a fixed probability α > 0.5 that each
one of its neighbors will have chosen the best alternative (repulsive forces
are also easy to introduce by allowing α < 0.5, when we have agents known
as contrarians [27, 28, 29]). Let the two choices be A and B, and pi (t) be
the probability agent i assigns at time t to the probability that A is the best
choice. CODA assumes a fixed likelihood α ≡ P (OAj |A), representing the
chance that, if A is indeed the best choice, when observing agent j, i will
observe j prefers A, indicated by OAj . If we don’t assume that the problem
is symmetrical in relation to both choices, that is, α ≡ P (OAj |A) 6= β ≡
3
P (OBj |B), a simple use of Bayes theorem will show how pi (t) is altered. Per
example, if agent j prefers A, we have
p(t) α
pi (t + 1|OAj ) = . (1)
1 − p(t) 1 − β
The model becomes much simpler if we use the log-odd function
p
ν ≡ ln( ) (2)
1−p
(where the agent index and time dependence were omitted). As p exists in
the interval 0 ≤ p ≤ 1, we have −∞ ≤ ν ≤ +∞. In the ν variable, a simple
additive model is obtained. That is, for the symmetric case where α = β,
ν(t + 1) = ν(t) ± C, (3)
with the plus sign corresponding to agent j preferring A and the minus sign
α
to the opposite choice and where C = ln 1−α (if α 6= β, the plus and minus
terms just have different sizes).
This symmetric case can be trivially normalized to ν ∗ so that when A is
observed, the agent adds +1 to ν ∗ , and when B is observed, −1 is added,
but we won’t normalize ν in this paper. The choice of the agent is defined
as the sign of ν ∗ , with positive signs indicating A is chosen. It is important
to notice that Equation 3 is actually general, even if we assume different
likelihoods. The expression for C can and will change, however, since C is a
simple constant due to the use of the simple likelihood of CODA model.
2.2 Continuous time
CODA model is defined in discrete time. Each observation of new data is
a discrete event and, from an inference point of view, it seems to make no
sense to talk about changes between the discrete observations. Of course,
this does not prevent us from extending the model to continuous time, from
a mathematical point of view. One might, if desired, think of it as the whole
information arriving in smaller bits.
Such an extension can be trivially obtained. If after ∆t = 1 we have
∆ν = C, we can assume a linear change in time so that ∆ν = C∆t. Actually,
any power of ∆t would give the correct discrete limit, and we will see bellow
that this simpler expression is not necessarily the best choice.
4
As we take ∆t → 0, we can write tentatively
dν
= C. (4)
dt
That is, C, while associated to the change of ν, would play the role of velocity
for the opinion particle. This role is basically correct if one looks for a
function C that describe the full movement. However, as we will see later, if
we divide the dynamics into an inertial component and another arising from
a force, Equation 4 will have to be replaced, as the force component will
assume a particle starting from rest. In that case, C will be associated with
the acceleration, with a change depending on ∆t2 , as we will see bellow.
α
In unidimensional CODA, supposing the choice was A, C = ln 1−β is
obtained from the ratio between the probability that A is chosen when A
is true (α) and the probability that A is chosen if B is true 1 − β. For an
opinion particle i that observes the position νj of another particle, we could
define an extension according to whether νj is negative or positive.
However, if we want to apply this theory to physical problems, the choice
of the zero for the ν-axis should be free. While a translation in ν does
correspond to a change of the probability p, this is a problem only if we hold
to an objective meaning of p. That objective view is not available to each
particle, though. From the point of view of an opinion particle, the central
position (p = 0.5 and ν = 0) is its location. That is, we can have a reference
frame ν(i) for each particle where the positions of the other particles, νj (i),
where j = 1, · · · , N represents each one of particles in the universe, are
measured.
What we gain from these choices is that each particle is influenced just
by the sum of of the influences to its right to the sum of the influences at its
left. In agreement with Sciama [22], these inertial forces obtained from this
CODA-like structure decay slower than the square of distance. Here, they
actually don’t decay at all. Of course, other functional forms can indeed
provide an inertial effect, but it is interesting that a simple extension of
CODA already does it.
3 Opinion Particles in d dimensions
The discussion at the previous Section was valid for the continuous extension
of the CODA model. But the class of possible dynamics based on proba-
5
bilistic reasoning is much larger. Basically, as long as we have information
arriving from some source and a likelihood about that information, this can
be used to update the probabilities associated with any parameter. In the
case of opinion particles, that simply means that the information will be
processed, causing the particle to move. If we retain the probabilistic in-
terpretation, that means a different evaluation about some parameter. If
we are only interested in obtaining equations of movement, we have nothing
more than information being used by a particle to alter its position. In any
case, the general case for Opinion Dynamics models on discrete time based
on Bayesian rules were presented before [17]. What we are interested here
are in models in continuous time that can be used in a d-dimensional space.
Therefore, at this point, we will not use some concepts like what how agent
j expresses its opinion. In this article, just the position, velocity and needed
parameters of j will be observed by particle i.
For that, we must return to the probabilistic representation and the Bayes
Theorem. Notice that while CODA models is simpler on the ν, it can just
as well be defined on the probability space of the p variable. This is a simple
choice of variables, distributions on ν have equivalents on p as well as any
other parameterization we choose. In order to compare with Mechanics, ν
just has the obvious advantage of existing in the range −∞ < µ < ∞, making
it easier to interpret and compare to Cartesian coordinates.
Just like regular space, opinions have several components, since they can
be about several issues [30]. If those issues are independent, they can be rep-
resented by the existence of a dimensional space, where the opinion particles
can move on. If we still assume the particles just decide to move whether to-
wards −∞ or towards ∞, we have something similar to the extended CODA
model in d dimensions. Assuming particle j is situates at ~x(j)
Therefore, we have a continuum of infinite possible values for ~x(j) as well
as, possible, other parameters. Representing the set of parameters by Θ,
we have, for the first coordinate p1 (t + 1) as a function of p1 (t), represented
simply as p1 for simplicity sake,
p1 f (Θ|A1)
p1 (t + 1) = , (5)
p1 f (Θ|A1) + (1 − p1 )f (Θ|B1)
where A1 corresponds to the possibility that p1 = 1 is the correct answer, and
B1 to p1 = 0. In a way, this might be considered equivalent to the Holographic
Principle [31], in the sense that each position is simply a representation of
6
the desired choices p = 0 or p = 1, that correspond to ν = −∞ and ν = +∞,
respectively.
Converting to ν1 , we have Equation 3 again, except now we have one such
equation for each coordinate, that is the vectorial equation
d~ν ~
=C (6)
dt
where
~ = ln f (~x|A1 )
C .
f (~x|B1 )
This is not complete, yet, however, as there the assumption that only
p = 0 or p = 1 are correct values is not necessarily correct. Even for real
opinions, it is perfectly reasonable to conclude that the best value for some
parameter is not in one of its extremes. In the case we have a continuum
of possible values for a best position (p or ν or any other), we will need a
probability distribution over the positions and it is this distribution that will
be updated, with the sum in the denominator of Equation 5 substituted by
an integration over all possible values. This case will not be addressed here,
since it is not necessary for the examples in this paper, but it is perfectly
reasonable that it might be needed for some specific forces and interactions.
4 Forces and opinions
Let’s assume for now that ν-space has the normal geometric structure of a
flat, Newtonian space-time and that there are three components of ν, equiv-
alent to the three Cartesian spatial directions. A similar supposition will be
made in the context of space-time coordinates.
We have obtained so far an equation for the velocities of the opinion
particles (Equation 4) and we have discussed how we can explain inertia in
a straight-forward way. To compare the movements of the opinion particles
with those of Newtonian particles we need to calculate the likelihoods that
match known forces. Notice that C is not actually a likelihood, but the
logarithm of the ratio of two likelihoods (or, if the likelihoods are defined
directly on ν, the ratio of two likelihods). The equivalence to a force can be
achieved by derivating Equation 4, to obtain the acceleration of the particle.
Since C, in general, depends on the relative position of the other particles,
it will change as the particle moves, the opinion particle will have a non-null
7
acceleration. And, given we define the particle mass in a consistent way, we
can derive Equation 4 in order to obtain a relation between a force F known
from Newtonian Mechanics and its associated likelihood. If C is associated
with the velocity of the particle, this is trivially given by, assuming a constant
mass m,
dC
F =m . (7)
dt
Integrating the expression provides us with the likelihood that is equivalent
to any force up to a constant term. We will see bellow that, in several cases,
C might be better identified with a force instead of the velocity. In those
cases, the Equation 7 is obviously not valid. Should C be associated with
∆t2 , it will play the role of an acceleration and its relation to the force is far
more obvious.
5 Velocities, inertia, and trust effects
One thing to notice is that, if C depended only on the relative position
of other particles, a given particle would have its trajectory from a fixed
point completely defined, with no possibility to have it, while at that point,
different velocities. This holds also for inertial forces and it would make
opinion particles behave in a way that is not compatible with the movement
of physical particles. While this is not a problem for Opinion Dynamics
problems, it would limit the applicability of the concept. Therefore, C must
have other types of dependencies, so that different velocities are allowed at
a point.
One possible way yo address the initial velocity problem we can introduce
trust between the particles [26, 20]. Basically, this might be seen as equiv-
alent to the problem of model choice in a Bayesian context [32], where the
likelihood is obtained from a sum of at least two different probability distribu-
tions. The term sum, of course, can mean, in the last sentence, an integration
over the values of a parameter. If one of the distributions is assumed to be
not related to a best position for the particle, as non-informative as possible,
that part of the distribution will not contribute to the movement. In a one-
dimensional model with two distributions [26], if the informative distribution
has a probability q of being the best model, the non-informative one will have
a probability 1 − q associated to it. This means that the movement of the
particle will be diminished by a factor q < 1 when compared with the model
8
with no non-informative part. Of course, as the particle moves, q may also
be updated and, as such, it does not have to be constant over the movement.
This means q can be seen as a measure of how much the information arriving
from the other particle is considered trustworthy or useless. This can cause
different velocities in any number of d dimensions, since q can be a function
of the angles in the d − 1 sphere.
In terms of implementation, however, the problem is much simpler, not re-
quiring a full model where distributions are specified and the hyper-parameter
for the probability of each distribution constantly updated. What happens is
that we already know the answer. The inertia component CI of any complete
particle description should be
~ I (t + dt) = ~v (t),
C (8)
where C~ I stands for the inertial likelihood and ~v(t) is the velocity of opinion
particle at time t.
6 Harmonic Oscillator
As an example of the application of the concept of using likelihoods and
opinion particles to describe Newtonian particles, let’s see how we can obtain
an harmonic oscillator. Equation 7 seems to provide a direct way to compute
a function C given a known force F , in this case, F = −kx. Basically, we
can use Newton laws to obtain how the force evolves in time and integrate
Equation
√ 7 in order to obtain a total function of the likelihoods given by
C ∝ 1 − x2 . While essentially correct, this solution misses some interesting
aspects of using opinion particles.
By integrating the force using the actual movement of the particle, both
the effects of inertia and the force get mixed in the function C. However, in
order to extend the ideas presented herein to General Relativity, it would be
useful to separate the inertial likelihood CI from the likelihood associated to
the force, CF .
In order to obtain the movement of a particle under an harmonic oscillator
force, let’s return to what happens in small time increments, from t1 to
t2 = t1 + ∆t. In this case, we can rewrite Equation 8 as
x(t1 ) − x(t0 )
CI (t2 ) = v(t1 ) = ,
∆t
9
where t0 = t1 − ∆t is the instant before t1 . Therefore, the change in position
from the inertial term is simply given by
x(t2 ) = x(t1 ) + v(t1 )∆t = x(t1 ) + [x(t1 ) − x(t0 )]. (9)
If CI were the only thing affecting the particle, its influence would be constant
and, therefore, the opinion particle would trivially follow a straight line.
It remains to introduce the effect of the force, by using Equation 7. We
have that
−kx(t1 )∆t = m∆C = m (C(t2 ) − C(t1 ))
Since C is basically the velocity of the particle, we would have, in principle,
x(t2 ) − x(t1 ) x(t1 ) − x(t0 )
− kx(t1 )∆t = m − . (10)
∆t ∆t
If we were to solve this, however, we would be back to using the full movement
of the particle and not just the effect of the force. If we mean to have just
the effect of CF from Equation 10, we must assume that the particle was
initially at rest, that is x(t1 )−x(t
∆t
0)
= 0. That way, we have
k
x(t2 ) = x(t1 ) − x(t1 )(∆t)2 . (11)
m
Combining the movements from inertia and the force, we have
k
x(t2 ) = x(t1 ) − x(t1 )(∆t)2 + [x(t1 ) − x(t0 )]. (12)
m
Figure 1 shows the time evolution obtained from Equation 12 as circles,
compared to the exact solution of the unidimensional harmonic oscillator,
represented as the line. We can see clearly that, despite its unconventional
form, Equation 12 does provide the correct dynamics for the problem.
7 Discrete Choices with Continuous Verbal-
ization
Finally, in order to better explore the formalism and also in order to make
some of its properties clearer, a variant of the CODA model will be discussed
10
Harmonic Oscillator
1.0
0.5
0.0
x
−1.0 −0.5
0 100 200 300 400
Time
Figure 1: Time evolution of the discrete harmonic oscillator, given by Equa-
tion 12 compared to the exact solution of the same problem. In the figure,
the discrete points are represented as circles and the exact solution is given
by the line. In the case shown, ∆t = 0.05 and ω = 1, while time in the figure
is given as multiples of ∆t.
11
where the communication between the agents is not a discrete spin value,
meaning the observed particle is at what side of the observing one. Instead,
the full probability pj particle j assigns to the possibility (measured in the
reference frame of i) that the right choice is x = +1 will be used in the in-
ference made by the particle i. Notice that the fact that the communication
is continuous does not imply that p should be. We still have a problem with
only two possible real choices p = −1 or p = +1. However, the continuous
probabilistic value is the communicated information. This distinction is a
very important albeit neglected one. In continuous opinion models, it is usu-
ally assumed that both the communication and the decision are continuous,
but that doesn’t have to be the case.
As the communication phase in the framework was changed, we need
now a new likelihood, that neighbor agent j will issue the value pj in the
case where p = −1 and in the case where p = +1, that is, functions f (pj |A)
and f (pj |B). Since all values for pi are limited to 0 ≤ pi ≤ 1, the simplest
choice is to take Beta distributions Be(pj |α, β) as priors.
1
Be(pj |α, β) = pα−1 (1 − pj )β−1
N(α, β) j
where N(α, β) is obtained from Gamma functions by
Γ(α)Γ(β)
N(α, β) = .
Γ(α + β))
As the Beta function is symmetric in α e β, if we want to keep the
symmetry between p = 0 and p = 1, we must have for the likelihoods that,
if f (pj |x = 1) = B(α, β), then f (pj |x = −1) = B(β, α). By applying the
Bayes Theorem to this problem, agent i, when observing pj , will update pi
to pi (t + 1),
pi pα−1
j (1 − pj )β−1
pi (t + 1) = . (13)
pi pjα−1 (1 − pj )β−1 + (1 − pi )pjβ−1 (1 − pj )α−1
If we adopt the same transformation of variables as in CODA model and
calculate νi we will see that the denominators cancel out and we have that
" α−β #
pi (t + 1) pi (t) pj (t)
ln = ln + ln (14)
1 − pi (t + 1) 1 − pi (t) 1 − pj (t)
12
Equation 14 can be rewritten in a more elegant fashion as
νi (t + 1) = νi (t) + (α − β)νj (15)
This is similar to the CODA dynamics, except that now, at each step, instead
of adding a term that is constant in size and only varies in sign, we add a
term proportional to the log-odds of the opinion of the neighbor.
However, we must also remember that, unlike the CODA case, that led us
to inertia, the objective here is to obtain an equivalent to a force component.
In Equation 11, we see that a force component enters as a term with ∆t2 ,
while, if we repeat the argument that associated C with velocity, we will just
have a ∆t. More importantly, very simple simulations of the limit ∆t → 0
show that a term ∆t2 is required in order for the dynamics to have a proper
limit. Per example, if, in Equation 12, the force term had just a linear
dependence on ∆t, the system would still oscillate in an apparently correct
way for a fixed ∆t, but with a frequency that would grow as ∆t → 0.
Of course, that might not seem reason enough. We certainly need to
understand why we can use a quadratic instead of a linear term, since C
seemed more easily associated with velocities. But we must remember that
we have divided the dynamics in two parts, a inertial part and one associated
with forces. The task of keeping the velocities the same as before is performed
by CI . And, as we have seen for the harmonic oscillator, the way CF works
is by providing the same change in position the particle would have if it were
at rest. That is, expanding the movement in a power series, the linear term,
giving the velocity due to CF must indeed be null. The first contribution,
therefore, is the acceleration, associated with ∆t2 .
That means we will have, from Equation 15 plus the inertial term, when
we start heading towards continuous time, the following equation for the
movement under the influence of Beta likelihoods:
νi (t + 1) = νi (t) + (α − β)νj ∆t2 + [νi (t) − νi (t − 1)], (16)
which is the same as the Equation 12 for the harmonic oscillator with ν in
the place of x, for a proper choice of α and β. The main difference is that
instead of νi the term in ∆t2 has νj . If one assumes, as mentioned before,
that a particle, at each instant, adjusts its coordinate system so that the
frame will return to the center, that is, we have an evolving referential frame
F where νiF (t) = 0, then νj is exactly the distance between i and j. That is,
13
in order to recover the results from Newtonian Mechanics, it seems we must
have translation invariance.
It is worth remembering that the models here assume the particles are
making an inference whether to move to −∞ or to +∞. This is not a
necessary characteristic of forces acting on opinion particles, they could very
well be trying to establish a value ν to move to. It is interesting, however,
that, for both inertial movement and the harmonic oscillator that was not
necessary.
8 Coordinate Changes and Metrics
We have seen the first details about how opinion particles relate to Newtonian
Mechanics. However, if we want the dynamics of the opinion particles to
be more general, it should also provide ways to describe central theories
of modern Physics, as General Relativity and Quantum Mechanics. Here,
however, I will just investigate if it makes sense for a likelihood to be used
in order in a way that corresponds to the geodesics in a space where points
are identified by the ν variable. The objective is simply to provide some
evidence that this framework could indeed be useful to describe gravity. The
already observed need for a local reference frame as well as the separation of
the movement in an inertial part and the part associated to other forces seem
to suggest that General Relativity should be a natural way to deal with ν
coordinates or any other coordinates we might decide to use to parameterize
the space-time.
As we have discussed, ν has no predefined geometric structure and, by
identifying a likelihood with movement on geodesics, we might be able to
get a new way to describe Gravity that could be compatible with General
Relativity. It should be noter that connections between distributions and
geometrical properties and curved spaces are not new in inferential problems
[33].
Since ν coordinates have no special geometric interpretation to them,
vectors in one point (per example, the velocity of the particle) have no natural
way to be defined as parallels at other point. The aim, if we want to find how
opinion particles can represent the movement of a particle in a gravitational
relativistic field, is to be able to define the likelihoods in terms of the metric
of the space-time.
14
8.1 Parallel Transport
So far, we have worked only with the spatial coordinates. We need to define
the model in the space-time and, for such, extending the likelihood to act
on a four-vector becomes necessary. In the Newtonian models, we had time
simply pass in equal infinitesimal steps. That means the quantity C α will
need a time component. From Equation 6, we can see that, in the discrete
case, the purpose of C α is to update the components of the position vector ν α .
That is, the time component for C α should be the time passed between the
two observations in the frame of reference of a particle moving on a geodesic
between the points. This is valid for any likelihoods in a space-time. In order
to differentiate the one that should be equivalent to gravity, I will denote it
by Gα from now on.
Since Gα is associated with how the coordinates of a particle change in
a gravitational field, it must correspond to the geodesics. That is, it must
have the same effect that parallel transport has on vectors, which means
that, despite its notation, it is not really a 4-vector field. This is easy to
understand, since Gα , as it was the case in Newtonian Mechanics, needs to
be not just a function of the position ν α but also of the velocity the particle
arrives at ν α and, therefore, might have different values at the same point
for different particles.
Following the same reasoning behind Equation 6, we must move from the
discretized probabilistic model to a continuum one. Here, however, we must
replace time for an invariant quantity, therefore, we need to differentiate the
position vector in relation to the distance element ds. In this case, Equation
6 becomes
dν α
= Gα . (17)
ds
A relation between the log-odds space and the coordinate space-time of
General Relativity is needed to proceed. Given the linear effect Gα already
has on ν coordinates, I will assume that both coordinate systems are the
same, that is that xα = ν α . This makes calculating the field Gα possible
for a given known gravitational problem, simply by observing that the term
dν α
ds
in Equation 17 can be obtained easily from the geodesic equation of
the movement of a particle (from any General Relativity textbook, as, per
example,[34])
duα
= −Γαβγ uβ uγ , (18)
ds
15
α
where uα = dνds = Gα .
Finally, it still needs to be investigated if the continuum limit would
work properly this way or would require a dependence in ds2 . This, however,
would mean a very simple change in Equation 18, since it is already the
second derivative of G that enters the equation.
9 Discussion
The work presented here is expected to provide a first step towards building
a complete theory for the movement of opinion particles. As an example of
its applicability in physical problems, we have shown how we can describe
the movement of a harmonic oscillator using this framework. In order for
more applications to Physics to become real, of course, other solutions that
match likelihoods and forces or metrics must be found. The exact functions
that do that translation will be studied in the future. I hope the functions
prove to be simple enough to allow analytical calculations, but there is no
reason for that to be the case.
In any case, the possibility of defining any forces (and, therefore, po-
tentials) and Gravitation in the same theoretical framework opens a new
possible path of exploration towards building an unified description of the
fundamental physical laws. This possibility alone makes the theory for the
movement of opinion particles quite worth exploring. One extra advantage
it already provides is in the area of Opinion Dynamics, with a way to trans-
form discrete opinion models that need simulation techniques into differential
equations that might have an analytic solution.
Of course, many problems are still open. Really describing gravity and
quantizing the theory are the main ones to be tackled in the near future.
Acknowledgments
The author would like to thank the Fundação de Amparo à Pesquisa do
Estado de São Paulo (FAPESP) for partial support to this research under
grant 2009/08186-0.
16
References
[1] Katrin Becker, Melanie Becker, and John H. Schwarz. String Theory
and M-Theory. Cambridge University Press, 2007.
[2] T. Thiemann. Modern Canonical Quantum General Relativity. Cam-
bridge University Press, 2007.
[3] E.T. Jaynes. Probability Theory: The Logic of Science. Cambridge,
Cambridge University Press, 2003.
[4] A. Caticha. Relative entropy and inductive inference. In G. Erickson and
Y. Zhai, editors, Bayesian Inference and Maximum Entropy Methods
in Science and Engineering, number 75 in AIP Conf. Proc. 707, 2004.
arXiv.org/abs/physics/0311093.
[5] A. Caticha and R. Preuss. Maximum entropy and bayesian data analysis:
Entropic prior distributions. Phys. Rev. E, 70:046127, 2004.
[6] Ariel Caticha and Adom Giffin. Updating probabilities. In
A. Mohammad-Djafari, editor, Bayesian Inference and Maximum En-
tropy Methods in Science and Engineering, volume 872 of AIP Conf.
Proc., page 31, 2007.
[7] Philip Goyal. Information physicstowards a new conception of physical
reality. Information, 3:567–594, 2012.
[8] R. T. Cox. The Algebra of Probable Inference. John Hopkins University
Press, 1961.
[9] C. Castellano, S. Fortunato, and V. Loreto. Statistical physics of social
dynamics. Reviews of Modern Physics, 81:591–646, 2009.
[10] S. Galam, Y. Gefen, and Y. Shapir. Sociophysics: A new approach of
sociological collective behavior: Mean-behavior description of a strike.
J. Math. Sociol., 9:1–13, 1982.
[11] S. Galam and S. Moscovici. Towards a theory of collective phenomena:
Consensus and attitude changes in groups. Eur. J. Soc. Psychol., 21:49–
74, 1991.
17
[12] K. Sznajd-Weron and J. Sznajd. Opinion evolution in a closed commu-
nity. Int. J. Mod. Phys. C, 11:1157, 2000.
[13] G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch. Mixing beliefs
among interacting agents. Adv. Compl. Sys., 3:87–98, 2000.
[14] R. Hegselmann and U. Krause. Opinion dynamics and bounded con-
fidence models, analysis and simulation. Journal of Artificial Societies
and Social Simulations, 5(3):3, 2002.
[15] André C. R. Martins. Continuous opinions and discrete actions in opin-
ion dynamics problems. Int. J. of Mod. Phys. C, 19(4):617–624, 2008.
[16] Serge Galam. Sociophysics: A Physicist’s Modeling of Psycho-political
Phenomena. Springer, 2012.
[17] André C. R. Martins. Bayesian updating as basis for opinion dynamics
models. AIP Conf. Proc., 1490:212–221, 2012.
[18] André C. R. Martins. Mobility and social network effects on extremist
opinions. Phys. Rev. E, 78:036104, 2008.
[19] André C. R. Martins. Discrete opinion models as a limit case of the
coda model. In Liu Yijun and Zhou Tao, editors, Social Physics Catena
(No.3), pages 146–157. Science Press, Beijing, 2012. arXiv:1201.4565v1.
[20] André C. R. Martins. Trust in the coda model: Opinion dynam-
ics and the reliability of other agents. Physics Letters A, in press
arXiv:1304.3518, 2013.
[21] Ernst Mach, T. J. McCormack, and Karl Menger. The Science of Me-
chanics: A Critical and Historical Account of Its Development. Open
Court Publishing Company, 1988.
[22] D. W. Sciama. On the origin of inertia. Royal Astronomical Society,
113:34–42, 1953.
[23] R. H. Dicke. Cosmology, mach’s principle and relativity. Am. J. Phys.,
31:500–509, 1963.
[24] Julian Barbour. The definition of mach’s principle. Found. Phys.,
40:1263–1284, 2010.
18
[25] Hanno Essén. Mechanics, cosmology and mach’s principle. Eur. J.
Phys., 34:139–145, 2013.
[26] André C. R. Martins. Bayesian updating rules in continuous opinion
dynamics models. Journal of Statistical Mechanics: Theory and Exper-
iment, 2009(02):P02017, 2009. arXiv:0807.4972v1.
[27] S. Galam. Contrarian deterministic effect: the hung elections scenario.
Physica A, 333:453–460, 2004.
[28] S. Galam. Heterogeneous beliefs, segregation, and extremism in the
making of public opinions. Physical Review E, 71:046123, 2005.
[29] André C. R. Martins and Cleber D. Kuba. The importance of disagree-
ing: Contrarians and extremism in the coda model. Adv. Compl. Sys.,
13:621–634, 2010.
[30] R. Vicente, André C. R. Martins, and N. Caticha. Opinion dynamics of
learning agents: Does seeking consensus lead to disagreement? Journal
of Statistical Mechanics: Theory and Experiment, 2009:P03015, 2009.
arXiv:0811.2099.
[31] Leonard Susskind. The world as a hologram. J. Math. Phys., 36:6377,
1995.
[32] Anthony O’Hagan. Kendall’s Advanced Theory of Statistics: Bayesian
Inference, volume 2B. Arnold, 1994.
[33] Shun ichi Amari and Hiroshi Nagaoka. Methods of Information Geom-
etry, volume 191 of Translations of Mathematical Monographs. Oxford
University Press, 2000.
[34] P. A. M. Dirac. General Theory of Relativity. Princeton University
Press, 1996.
19