# Mathematical Modeling of Public Opinions: Parameter Estimation, Sensitivity Analysis, and Model Uncertainty Using an Agree-Disagree Opinion Model

Authors: Sara Bidah; Omar Zakary; Mostafa Rachik; Hanane Ferjouchia
Journal: Abstract and Applied Analysis
 (2020)

---

## Abstract

In this paper, we present a mathematical model that describes agree-disagree opinions during polls. We first present the different compartments of the model. Then, using the next-generation matrix method, we derive thresholds of the stability of equilibria. We consider two sets of data from the Reuters polling system regarding the approval rating of the U.S presidential in two terms. These two weekly polls data track the percentage of Americans who approve and disapprove of the way the President manages his work. To validate the reality of the underlying model, we use nonlinear least-squares regression to fit the model to actual data. In the first poll, we consider only 31 weeks to estimate the parameters of the model, and then, we compare the rest of the data with the outcome of the model over the remaining 21 weeks. We show that our model fits correctly the real data. The second poll data is collected for 115 weeks. We estimate again the parameters of the model, and we show that our model can predict the poll outcome in the next weeks, thus, whether the need for some control strategies or not. Finally, we also perform several computational and statistical experiments to validate the proposed model in this paper. To study the influence of various parameters on these thresholds and to identify the most influential parameters, sensitivity analysis is carried out to investigate the effect of the small perturbation near a parameter value on the value of the threshold. An uncertainty analysis is performed to evaluate the forecast inaccuracy in the outcome variable due to uncertainty in the estimation of the parameters.


---

## Body

## 1. Introduction

To attract the attention of citizens to participate in political life, candidate parties must provide certain benefit, because this task requires time, knowledge, motivation, and money. Citizen participation is then guaranteed if there are proportional benefits with political mobilizations, individual economic, psychological, and social costs [1].Public opinion is defined as the collective opinion of individuals within a society on a specific idea even if they have several perceptions. For example, in current American politics, many people tend to overestimate their opinions when they want to disagree with what they consider public opinion. This overestimation leads to false perceptions of public opinion, because all opinions, exaggerated or not, contribute to forming public opinion. Candidates and the media often use opinion polls before and during election campaigns to determine which candidates are first and which are likely to emerge victoriously. These polls are intention polls for a sample of potential voters, which reveal the expected voting quotas. Then, future results allow candidates to negotiate and discuss opinions based on a particular outcome. Public opinion polls are usually used to determine people’s position in political life, votes, and belonging and political orientations by answering questions about their opinions, personal characteristics, and their activities. These answers are collected, counted, mathematically analyzed, and interpreted. These analyses allow the candidate parties to improve their campaign strategies to change their positions in public opinion in the long term [2].In this article, we study how the preferences and expectations of voters are in competition and how they convince themselves in the face of the information of the candidates in opinion polls. We start by presenting the mathematical model that describes the evolution of opinions during polls to predict the probability of the result. Based on the next-generation matrix approach, we compute the thresholds of equilibrium stability. We consider two sets of data from the Reuters poll system regarding the approval rating of the US President in two terms. These two weekly polls track the percentage of Americans who agree and reject the way the president does business. To validate the underlying model, we use nonlinear least-squares regression to fit the model to the actual data. In the first survey, we only consider 31 weeks to estimate model parameters and then compare the rest of the data to the model result over the remaining 21 weeks. We show that our model fits real data correctly. The second survey data is collected for a period of 115 weeks. We again estimate model parameters and show that our model can predict the outcome of the survey in the coming weeks. Sensitivity analysis and model uncertainty are carried out to validate the underlying model.The development of more efficient models and the accuracy of the modeling process require reliable statistical methods, whereas sensitivity analysis is one of the most recently used methods [3–6], which is used in development decisions, recommendations, understanding and quantification of systems, verification of the validity and accuracy of the model and even in identification of the important parameter for other studies [7]. Sensitivity analysis is a tool to evaluate the effect of changes in the value of the input parameter on the output value of a simulation model. The method we use in this article uses the single-point derivative when an input parameter is changed while maintaining other fixed parameters. We calculate the sensitivity indices of the equilibrium stability thresholds to the model parameters. These indices tell us how important each parameter is to changing opinions. We perform also uncertainty analysis for the proposed model to illustrate the relationship between the model and the observations from the real data.The paper is organized as follows: Section2 introduces our new model, giving some details about interactions between the different compartments and parameters of the model. In Section 3, we derive basic reproduction numbers. Section 4 provides the results of the parameters’ estimation. In Section 5, sensitivity analysis is performed to identify the most important parameter in the proposed model for both polls. Section 6 provides the uncertainty analysis of the model parameters. and Section 7 concludes.


## 2. Presentation of the Model

Based on the fact that there are many scenarios involving binary decision, the poll model that we present here describes the opinions of agreement (or approval) and disagreement (or disapproval), regarding a candidate or idea, in polls. The model can also describe the political position of a specific candidate even on situations when there are more than two candidate parties, because we can always reduce the situation to two decisions. For instance, consider the four partsA, B, C, and D and we will investigate the political position of the party A, then we can examine the two subsets (A) and (B, C, and D) of the investigation, votes for A are considered as Agree opinions, and votes for B, C, or D are considered Disagree opinions. Or more simply, we consider the opinion poll of voters on the performance of the party studied.Without loss of generality, we devise a mathematical model describing the evolution of Agree and Disagree opinions during polls, and the types of surveys we consider are surveys that can be answered in agreement, disagreement, or otherwise. Thus, the targeted population by the poll is regrouped into three groups, and the model herewith has been formulated using compartments. Three compartments have been considered, and each of them has been described below.(1)
Indifferent (I) Individuals. Undecided or ambivalent people, or people who do not know about the poll yet or those who abstain from voting for personal reasons. This category has weak or nonexistent attitudes about the ideas, parties, or candidates and lacked any strong positive or negative associations.(2)
Agree/Approved (A) Individuals. People in agreement with the idea being studied.(3)
Disagree/Disapproved (D) Individuals. People in disagreement with the idea being studied.For the modeling processes, a set of assumptions has been used. These are as follows:(1)
The targeted population is well mixed; that is, the indifferent individuals are homogeneously spread throughout the entire population(2)
Recruitment and mortality are negligible under the temporal scale consideration; therefore, no individual recruited and no individual dies during the poll(3)
Individuals have the right to communicate with each other and can thus convince themselves(4)
People who are unsure, ambivalent or with mixed feelings about their opinion are indifferent(5)
People who abstain from voting are indifferentEveryone has their reasons for agreement or disagreement, and indifferent can be convinced by reasons of approved people at a rateβ1 or by reasons of disapproved individuals at a rate β2. Agreeing individuals can be convinced by disagreeing people at a rate α1, or disagreeing people convinced by agreeing individuals at a rate α2. People can abstain from voting or lose the interest without any direct contact with individuals from the opposite opinion group, then agreeing people become indifferent at the rate γ1, and disagreeing people become indifferent at the rate γ2. All contacts are modeled by the standard incidence rate. A flowchart describing the different interactions between the compartments of the model is presented in Figure 1.Figure 1
Flow chart for the model (1)-(3).All these assumptions and considerations are written as the following system of ordinary differential equations:(1)I′=−β1AIN−β2DIN+γ1A+γ2D,(2)A′=β1AIN+α1ADN−α2ADN−γ1A,(3)D′=β2DIN+α2ADN−α1ADN−γ2D,where I0≥0, A0≥0, and D0≥0, and N=I+A+D. Note that N′=I′+A′+D′=0; thus, the population size N is considered a constant in time. A summary of parameters description is given in Table 1.Table 1
Parameter descriptions.
ParameterDescriptionValuesβ1Indifferent to Agree transmission rate0.4842β2Indifferent to Disagree transmission rate0.5776α1Disagree to Agree transmission rate0.7439α2Agree to Disagree transmission rate0.2990γ1Interest lose factor of Agree individuals0.4498γ2Interest lose factor of Disagree individuals0.1699I0Initial Indifferent population6.9A0Initial Agree population45.8D0Initial Disagree population47.3Sensitivity analysis allows us to measure the relative change in a state variable when a parameter changes while the other parameters are retained fixed at constant values. Later, we carry out a local sensitivity analysis of the thresholds to the model parameters. Next, we determine the covariance matrix which provides uncertainties on the estimation of parameters and captures relations in the measurement uncertainties based on two sets of data. To do this, we start by computing the model thresholds in the next section.


## 3. Thresholds

In epidemiology, the basic reproductive numberR0 (or epidemic threshold) is defined as the average number of secondary cases of an infection produced by a “typical” infected individual during his entire life as infectious when introduced in a population of susceptibles [8–12]. The threshold R0 has been shown to be mathematically characterized in terms of infection transmission as a “demographic process”, but offspring production is not seen as giving birth in a demographic sense, but it causes new infections through transmission [13]. Thus, the infection process can be considered a “consecutive generation of infected individuals”. The following growing generations indicate a growing population (i.e., an epidemic), and the growth factor for each generation indicates the potential for growth. So the mathematical characterization of R0 is this growth factor [13]. Generally, if R0>1, an epidemic occurs whereas if R0<1, there will probably be no epidemic.Following this definition, we define the thresholds as follows:RD0 is the average number of new disagreements produced by an individual disagreeing introduced in a population of indifferent people, during the period in which he or she was in this opinion. And, RA0 is the average number of new agreements produced by an individual in agreement and that was introduced into a population of indifferent people during the period in which he or she was in that opinion.For the analysis of epidemic models, the first step is to calculate the disease-free equilibrium (DFE). This equilibrium point is then used to calculate the basic reproductive number using the next-generation matrix. In this contribution, letRX0 be the threshold of growth of the opinion X (either X=A “Agree” or X=D “Disagree”). Then, RD0 is the threshold associated with the Disagree-free equilibrium while RA0 is the one associated with the Agree-free equilibrium.To use the next-generation matrix method, first, we need to determine the Disagree-free equilibrium and the Agree-free equilibrium, and second, we can derive associated thresholds. Recall that ([8]) the points of equilibrium of the systems (1)–(3) are the solutions of
(4)I′=A′=D′=0.For the Disagree-free equilibrium, when there is no negative opinion, that is whenD=0, we have
(5)I∗=γ1β1N,A∗=Nβ1−γ1β1,where I∗ and A∗ represent the numbers of indifferent and agreeing individuals, respectively, in the absence of disagreeing people. Hence, for the system governed by (1)–(3), the disagree-free equilibrium is
(6)e1=γ1β1N,Nβ1−γ1β1,0.Following the second-generation approach [13], we compute the threshold RD0 associated to the disagree-free equilibrium, which is
(7)RD0=α2β1−α2γ1+β2γ1α1β1−α1γ1+β1γ2.Regarding the Agree-free equilibrium when there is no positive opinion, i.e., if we putA=0, we have
(8)I∗=γ2β2N,D∗=Nβ2−γ2β2,where I∗ and D∗ represent the numbers of indifferent and disagreeing individuals, respectively, in the absence of agreeing people.Therefore, for the system governed by (1)–(3), the agree-free equilibrium is
(9)e2=γ2β2N,0,Nβ2−γ2β2By following also the second-generation approach [13], we compute the threshold RA0 associated to the Agree-free equilibrium, which is
(10)RA0=α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1.


## 4. Parameter Estimation

### 4.1. Methods

The approval rating is the percentage determined by the survey process that indicates the percentage of respondents who agree with a particular individual or program. Habitually, a level of approval is given to a politician based on answers to a public opinion poll in which persons’ samples are invited whether they approve or disapprove this particular political person. Due to the availability of data in the form of proportions, we need to work in terms of the proportions of Indifferent, Approved, and Disapproved individuals; thus, leti=I/N, a=A/N, and d=D/N denote the fractions of the classes I, A, and D in the population, respectively. By replacing I by i, A by a, and D by d, Equations (1)-(3) can be written as
(11)i′=−β1ai−β2di+γ1a+γ2d,a′=β1ai+α1ad−α2ad−γ1a,d′=β2di+α2ad−α1ad−γ2d,where i+a+d=1.To verify the reality of this model, we use nonlinear least-squares regression to fit the model to actual observations by following the method used in [15], where the authors estimated the parameters of their influenza epidemic models. The least-squares method, which involves minimizing the sum of the squared differences between the measurements and the model predictions, is used to estimate the unknown parameters for the two influenza epidemic models presented in [15].Therefore, the following process has been followed for parameter estimation: using MATLAB ode45 routine, the system of ordinary differential equations is solved numerically, with initial chosen values for parameters and state variables from Table1. Model outcomes are compared with the field data, and the Levenberg–Marquardt optimization algorithm determines a new set of parameters’ values with the model outcomes in a better fit to the field data [16, 17]. After new parameters’ values are determined by this optimizer, the system of ordinary differential equations is solved numerically using these new parameters’ value and the model outcomes are compared again with the field data. This iteration process between parameter updating and numerical solutions of the system of ordinary differential equations using the Runge–Kutta method [18, 19] continue till convergence criteria for the parameters are met. In this process of estimating parameters, about one thousand values are chosen using a random process for each of the parameters to be estimated.


### 4.2. First Poll

To estimate the model parameters, we use the data from a poll carried out by the Reuters polling system entitled “Approval of Obama’s handling of job as president” from [14]. The data in Table 2 are collected for approximately 52 weeks, from January 17, 2016, to January 1, 2019. The poll is answered with the following three questions: approve, disapprove, and mixed feelings. Here, we consider the third class of mixed feeling as indifferents, see Figure 2.Table 2
Data of the first Poll from [14], where the first week is 17/01/2016 with I=3.7, A=45.3, and D=51.1.
DateIndifferentAgreeDisagreeDateIndifferentAgreeDisagreeDateIndifferentAgreeDisagree24/01/20164.746.948.422/05/20165.848.545.818/09/20165.246.947.831/01/20164.446.948.729/05/20164.151.544.425/09/20165.448.945.707/02/20163.445.651.105/06/2016449.346.702/10/20164.550.54514/02/20164.544.950.612/06/20165.448.546.109/10/20164.550.844.721/02/20164.445.849.819/06/20164.949.14616/10/20164.349.74628/02/20165.645.249.226/06/20165.448.246.423/10/20165.249.84506/03/2016644.849.203/07/20165.449.145.530/10/20164.250.745.113/03/20164.645.549.910/07/20166.247.846.106/11/20164.351.344.420/03/20164.845.349.917/07/20165.945.348.913/11/20165.250.344.427/03/20164.946.448.824/07/2016644.449.720/11/20165.952.141.903/04/20164.848.346.931/07/20164.747.447.927/11/20164.953.24210/04/20164.648.147.307/08/20164.55045.404/12/20164.452.74317/04/20165.548.945.614/08/20165.147.747.211/12/20165.153.441.524/04/20164.446.748.921/08/20166.244.84918/12/20163.355.541.201/05/20164.349.646.228/08/20166.545.947.625/12/20164.452.443.208/05/20163.746.649.704/09/20165.34945.701/01/20173.751.444.915/05/20164.847.747.511/09/20164.248.447.408/01/20173.154.742.2Figure 2
Data of the first poll from [14], entitled: Overall, do you approve or disapprove about the way Obama is handling his job as President? see Table 2 for more details.To validate our model, we chose as a period of parameter estimation from 17/01/2016 to 21/08/2016 (31 weeks), see Figure3. Thus, the estimated parameters’ values are used to validate our model, by plotting the three states of the model with the rest of the real data from 28/08/2016 to 08/01/2017 (20 weeks), see Figure 4. It can be seen from this figure that the model fits correctly the real data, and it can predict the position of people opinions at every next point of time and then the outcome of the survey.Figure 3
The Agree-Disagree model fits for the poll results data within 30 weeks. (a) The solid line depicts the model output forI, while the crosses show the indifferent individuals from data. (b) The solid line depicts the model output for A, while the crosses show the Agree individuals from data. (c) The solid line depicts the model output for D, while the crosses show the Disagree individuals from data.Figure 4
Validation of the Agree-Disagree model with the poll results data until the 50th week.Based on the estimated parameters in the first period, and by results of [8], we can see which equilibrium point is stable; then, we can predict the outcome of the long-term survey. By calculating the thresholds corresponding to this poll, we use the parameters’ values from Table 3 and we get
(12)RA0=0.9943<1,RD0=0.9365<1.Table 3
Estimated parameters value for models (1)-(3) for the first poll.
Parametersβ1β2α1α2γ1γ2Values0.24970.89990.52340.54330.00010.0571And we can see thatβ1>γ1 and β2>γ2, which means that the model presents two equilibria in stability at the same time, that is, the Disagree-free equilibrium e1=0.0004,0.9996,0 and the Agree-free equilibrium e2=0.0635,0,0.9365, see [8] for more details. However, only the knowledge of the stability of the equilibria does not guarantee the knowledge of the result of the survey, because the convergence of the states of the model towards this equilibrium can take a long time. Voting could end before this convergence. In the situation of the bistability of equilibria, the knowledge of the poll outcome could be more difficult.Despite that, these results provide an insight into the position of people’s opinions and the course of events, to help the authorities take the necessary measures to bring the situation under control. For example, we can see that by having the data for 30 weeks, we can mathematically estimate the result of the survey before 20 weeks from the end of the poll, see Figure4. Thus, our results are more effective in predicting the outcome of the poll. We can see from Figure 4 that the Disagree-free equilibrium is the more attractive, because D started to tend towards 0. This means that after this decrease, there is no more need for interventions.


### 4.3. Second Poll

A poll is carried out also by the Reuters polling system entitled “Overall, do you approve or disapprove about the way Donald Trump is handling his job as President?” from [20]. The data in Table 4 are collected for approximately 115 weeks, from January 29, 2017, to April 7, 2019. The survey is answered also with the following three questions: approve, disapprove, and mixed feelings. These data are plotted against time in Figure 5. It is clear that Figure 5 depicts a positive equilibrium e∗, see [8] for more details about equilibria.Table 4
Data of the second poll from [20].
DateIndifferentAgreeDisagreeDateIndifferentAgreeDisagreeDateIndifferentAgreeDisagree29/01/20176.945.847.322/10/20174.636.159.322/07/20186415305/02/20175.347.547.229/10/20174.33560.829/07/2018542.352.712/02/20173.94551.105/11/20174.436.159.505/08/20184.941.753.419/02/20176.946.646.612/11/20175.137.757.212/08/2018540.454.626/02/20175.74648.219/11/20174.637.45819/08/20184.742.652.605/03/20175.548.645.926/11/20175.237.157.726/08/20185.938.355.912/03/20176.443.849.803/12/20174.335.460.302/09/20185.639.554.919/03/20175.84747.210/12/20174.134.861.109/09/20186.140.353.626/03/20176.846.24717/12/20175.536.158.416/09/20184.24253.702/04/20174.645.649.824/12/20175.238.85623/09/20185.641.652.809/04/20174.845.849.407/01/20184.541.853.730/09/20185.840.753.516/04/20176.140.753.214/01/20184.640.954.507/10/20186.438.455.223/04/20175.542.152.421/01/20184.736.658.714/10/20186.142.95130/04/20175.642.951.528/01/20184.337.85821/10/20186.741.152.207/05/20175.241.553.304/02/20185.540.753.928/10/20185.242.652.214/05/20176.237.356.511/02/20185.340.454.304/11/20185.741.952.421/05/20175.238.756.118/02/20186.140.253.711/11/20184.442.153.528/05/20174.539.456.225/02/2018540.654.418/11/20185.443.750.904/06/20174.837.25804/03/20185.440.254.425/11/20184.644.750.711/06/20173.838.357.911/03/2018639.654.502/12/20186.141.95218/06/20176.137.356.718/03/20184.739.555.909/12/20185.941.552.625/06/20176.837.655.725/03/20184.638.45716/12/20183.845.850.402/07/20174.139.156.801/04/20184.539.955.723/12/20184.941.753.409/07/20173.939.556.708/04/20185.637.35730/12/20185.440.154.516/07/20177.136.556.315/04/20183.640.356.106/01/20194.94055.123/07/20174.436.858.822/04/20185.138.156.813/01/20194.1415530/07/20174.536.459.129/04/20184.241.454.320/01/20194.540.55506/08/20175.437.756.906/05/20182.745.451.927/01/20194.139.356.613/08/20176.135.558.413/05/20182.543.354.203/02/20194.239.156.620/08/2017535.459.520/05/20182.343.554.210/02/20194385827/08/2017437.958.127/05/20184.243.452.517/02/20194.141.554.403/09/20175.139.255.803/06/20185.339.755.124/02/20194.740.654.610/09/20176.137.955.910/06/20185.64054.403/03/20194.73956.417/09/20175.737.556.817/06/20184.342.852.910/03/20195.142.152.924/09/20174.738.356.924/06/2018540.954.117/03/20194.340.75501/10/20175.435.858.701/07/2018540.554.424/03/20195.641.652.908/10/20174.73758.308/07/20183.940.655.531/03/20194.94253.115/10/20174.835.959.415/07/20185.639.155.307/04/2019539.355.7Figure 5
Data of the second poll from [20] entitled “Overall, do you approve or disapprove about the way Donald Trump is handling his job as President?,” see Table 4 for more details.When theit remains below 10% and dt is greater than at from the beginning of March 2017 until the end of the poll on 07 April 2019, we will use the data plotted in Figure 5, to estimate the parameters of models (11)-(13) for this poll.Estimated parameter values of models (11)-(13) for this poll are given in Table 5. By calculating the thresholds corresponding to this poll, we use the parameters’ values from this Table, we get
(13)RA0=1.1404>1,RD0=6.1346>1.Table 5
Estimated parameters value for models (1)-(3) for the second poll.
Parametersβ1β2α1α2γ1γ2Values0.12000.89000.48580.33620.08980.0001And we can see thatβ1γ2−β2γ1=−0.0799<0 and α=α1−α2>0, which means that the model presents the equilibrium e∗=0.0651,0.3868,0.5480 in stability, see the proof of Proposition 5 and Remark 2 in [8]. In this case, it may be noted that some control strategies are necessary to change the course of events prior to Election Day.Figure6 depicts the effectiveness of the estimation of the parameters, while in (a) we can see the indifferent individual curve of the model plotted with the solid line while the data of Indifferent individuals is plotted with small crosses. In (b), a comparison of the Agree people function of the model plotted with the solid line and the data of Agree individuals is plotted with small crosses. Also in (c), the estimated function of Disapproved people is plotted with the solid line and the data of disagree individuals is plotted with crosses. In Figure 7, the three functions of the model and the three states of the real data are plotted with different patterns to show that the model fits the data correctly.Figure 6
The Agree-Disagree model fits for the second poll data. (a) The solid line depicts the model output forI, while the small crosses show the indifferent data. (b) The solid line depicts the model output for A, while the small crosses show the Agree data. (c) The solid line depicts the model output for D, while the small crosses show the Disagree data.Figure 7
Three states of models (1)-(3) simulated with estimated parameter values from Table 5 and the real data of the second poll.


## 4.1. Methods

The approval rating is the percentage determined by the survey process that indicates the percentage of respondents who agree with a particular individual or program. Habitually, a level of approval is given to a politician based on answers to a public opinion poll in which persons’ samples are invited whether they approve or disapprove this particular political person. Due to the availability of data in the form of proportions, we need to work in terms of the proportions of Indifferent, Approved, and Disapproved individuals; thus, leti=I/N, a=A/N, and d=D/N denote the fractions of the classes I, A, and D in the population, respectively. By replacing I by i, A by a, and D by d, Equations (1)-(3) can be written as
(11)i′=−β1ai−β2di+γ1a+γ2d,a′=β1ai+α1ad−α2ad−γ1a,d′=β2di+α2ad−α1ad−γ2d,where i+a+d=1.To verify the reality of this model, we use nonlinear least-squares regression to fit the model to actual observations by following the method used in [15], where the authors estimated the parameters of their influenza epidemic models. The least-squares method, which involves minimizing the sum of the squared differences between the measurements and the model predictions, is used to estimate the unknown parameters for the two influenza epidemic models presented in [15].Therefore, the following process has been followed for parameter estimation: using MATLAB ode45 routine, the system of ordinary differential equations is solved numerically, with initial chosen values for parameters and state variables from Table1. Model outcomes are compared with the field data, and the Levenberg–Marquardt optimization algorithm determines a new set of parameters’ values with the model outcomes in a better fit to the field data [16, 17]. After new parameters’ values are determined by this optimizer, the system of ordinary differential equations is solved numerically using these new parameters’ value and the model outcomes are compared again with the field data. This iteration process between parameter updating and numerical solutions of the system of ordinary differential equations using the Runge–Kutta method [18, 19] continue till convergence criteria for the parameters are met. In this process of estimating parameters, about one thousand values are chosen using a random process for each of the parameters to be estimated.


## 4.2. First Poll

To estimate the model parameters, we use the data from a poll carried out by the Reuters polling system entitled “Approval of Obama’s handling of job as president” from [14]. The data in Table 2 are collected for approximately 52 weeks, from January 17, 2016, to January 1, 2019. The poll is answered with the following three questions: approve, disapprove, and mixed feelings. Here, we consider the third class of mixed feeling as indifferents, see Figure 2.Table 2
Data of the first Poll from [14], where the first week is 17/01/2016 with I=3.7, A=45.3, and D=51.1.
DateIndifferentAgreeDisagreeDateIndifferentAgreeDisagreeDateIndifferentAgreeDisagree24/01/20164.746.948.422/05/20165.848.545.818/09/20165.246.947.831/01/20164.446.948.729/05/20164.151.544.425/09/20165.448.945.707/02/20163.445.651.105/06/2016449.346.702/10/20164.550.54514/02/20164.544.950.612/06/20165.448.546.109/10/20164.550.844.721/02/20164.445.849.819/06/20164.949.14616/10/20164.349.74628/02/20165.645.249.226/06/20165.448.246.423/10/20165.249.84506/03/2016644.849.203/07/20165.449.145.530/10/20164.250.745.113/03/20164.645.549.910/07/20166.247.846.106/11/20164.351.344.420/03/20164.845.349.917/07/20165.945.348.913/11/20165.250.344.427/03/20164.946.448.824/07/2016644.449.720/11/20165.952.141.903/04/20164.848.346.931/07/20164.747.447.927/11/20164.953.24210/04/20164.648.147.307/08/20164.55045.404/12/20164.452.74317/04/20165.548.945.614/08/20165.147.747.211/12/20165.153.441.524/04/20164.446.748.921/08/20166.244.84918/12/20163.355.541.201/05/20164.349.646.228/08/20166.545.947.625/12/20164.452.443.208/05/20163.746.649.704/09/20165.34945.701/01/20173.751.444.915/05/20164.847.747.511/09/20164.248.447.408/01/20173.154.742.2Figure 2
Data of the first poll from [14], entitled: Overall, do you approve or disapprove about the way Obama is handling his job as President? see Table 2 for more details.To validate our model, we chose as a period of parameter estimation from 17/01/2016 to 21/08/2016 (31 weeks), see Figure3. Thus, the estimated parameters’ values are used to validate our model, by plotting the three states of the model with the rest of the real data from 28/08/2016 to 08/01/2017 (20 weeks), see Figure 4. It can be seen from this figure that the model fits correctly the real data, and it can predict the position of people opinions at every next point of time and then the outcome of the survey.Figure 3
The Agree-Disagree model fits for the poll results data within 30 weeks. (a) The solid line depicts the model output forI, while the crosses show the indifferent individuals from data. (b) The solid line depicts the model output for A, while the crosses show the Agree individuals from data. (c) The solid line depicts the model output for D, while the crosses show the Disagree individuals from data.Figure 4
Validation of the Agree-Disagree model with the poll results data until the 50th week.Based on the estimated parameters in the first period, and by results of [8], we can see which equilibrium point is stable; then, we can predict the outcome of the long-term survey. By calculating the thresholds corresponding to this poll, we use the parameters’ values from Table 3 and we get
(12)RA0=0.9943<1,RD0=0.9365<1.Table 3
Estimated parameters value for models (1)-(3) for the first poll.
Parametersβ1β2α1α2γ1γ2Values0.24970.89990.52340.54330.00010.0571And we can see thatβ1>γ1 and β2>γ2, which means that the model presents two equilibria in stability at the same time, that is, the Disagree-free equilibrium e1=0.0004,0.9996,0 and the Agree-free equilibrium e2=0.0635,0,0.9365, see [8] for more details. However, only the knowledge of the stability of the equilibria does not guarantee the knowledge of the result of the survey, because the convergence of the states of the model towards this equilibrium can take a long time. Voting could end before this convergence. In the situation of the bistability of equilibria, the knowledge of the poll outcome could be more difficult.Despite that, these results provide an insight into the position of people’s opinions and the course of events, to help the authorities take the necessary measures to bring the situation under control. For example, we can see that by having the data for 30 weeks, we can mathematically estimate the result of the survey before 20 weeks from the end of the poll, see Figure4. Thus, our results are more effective in predicting the outcome of the poll. We can see from Figure 4 that the Disagree-free equilibrium is the more attractive, because D started to tend towards 0. This means that after this decrease, there is no more need for interventions.


## 4.3. Second Poll

A poll is carried out also by the Reuters polling system entitled “Overall, do you approve or disapprove about the way Donald Trump is handling his job as President?” from [20]. The data in Table 4 are collected for approximately 115 weeks, from January 29, 2017, to April 7, 2019. The survey is answered also with the following three questions: approve, disapprove, and mixed feelings. These data are plotted against time in Figure 5. It is clear that Figure 5 depicts a positive equilibrium e∗, see [8] for more details about equilibria.Table 4
Data of the second poll from [20].
DateIndifferentAgreeDisagreeDateIndifferentAgreeDisagreeDateIndifferentAgreeDisagree29/01/20176.945.847.322/10/20174.636.159.322/07/20186415305/02/20175.347.547.229/10/20174.33560.829/07/2018542.352.712/02/20173.94551.105/11/20174.436.159.505/08/20184.941.753.419/02/20176.946.646.612/11/20175.137.757.212/08/2018540.454.626/02/20175.74648.219/11/20174.637.45819/08/20184.742.652.605/03/20175.548.645.926/11/20175.237.157.726/08/20185.938.355.912/03/20176.443.849.803/12/20174.335.460.302/09/20185.639.554.919/03/20175.84747.210/12/20174.134.861.109/09/20186.140.353.626/03/20176.846.24717/12/20175.536.158.416/09/20184.24253.702/04/20174.645.649.824/12/20175.238.85623/09/20185.641.652.809/04/20174.845.849.407/01/20184.541.853.730/09/20185.840.753.516/04/20176.140.753.214/01/20184.640.954.507/10/20186.438.455.223/04/20175.542.152.421/01/20184.736.658.714/10/20186.142.95130/04/20175.642.951.528/01/20184.337.85821/10/20186.741.152.207/05/20175.241.553.304/02/20185.540.753.928/10/20185.242.652.214/05/20176.237.356.511/02/20185.340.454.304/11/20185.741.952.421/05/20175.238.756.118/02/20186.140.253.711/11/20184.442.153.528/05/20174.539.456.225/02/2018540.654.418/11/20185.443.750.904/06/20174.837.25804/03/20185.440.254.425/11/20184.644.750.711/06/20173.838.357.911/03/2018639.654.502/12/20186.141.95218/06/20176.137.356.718/03/20184.739.555.909/12/20185.941.552.625/06/20176.837.655.725/03/20184.638.45716/12/20183.845.850.402/07/20174.139.156.801/04/20184.539.955.723/12/20184.941.753.409/07/20173.939.556.708/04/20185.637.35730/12/20185.440.154.516/07/20177.136.556.315/04/20183.640.356.106/01/20194.94055.123/07/20174.436.858.822/04/20185.138.156.813/01/20194.1415530/07/20174.536.459.129/04/20184.241.454.320/01/20194.540.55506/08/20175.437.756.906/05/20182.745.451.927/01/20194.139.356.613/08/20176.135.558.413/05/20182.543.354.203/02/20194.239.156.620/08/2017535.459.520/05/20182.343.554.210/02/20194385827/08/2017437.958.127/05/20184.243.452.517/02/20194.141.554.403/09/20175.139.255.803/06/20185.339.755.124/02/20194.740.654.610/09/20176.137.955.910/06/20185.64054.403/03/20194.73956.417/09/20175.737.556.817/06/20184.342.852.910/03/20195.142.152.924/09/20174.738.356.924/06/2018540.954.117/03/20194.340.75501/10/20175.435.858.701/07/2018540.554.424/03/20195.641.652.908/10/20174.73758.308/07/20183.940.655.531/03/20194.94253.115/10/20174.835.959.415/07/20185.639.155.307/04/2019539.355.7Figure 5
Data of the second poll from [20] entitled “Overall, do you approve or disapprove about the way Donald Trump is handling his job as President?,” see Table 4 for more details.When theit remains below 10% and dt is greater than at from the beginning of March 2017 until the end of the poll on 07 April 2019, we will use the data plotted in Figure 5, to estimate the parameters of models (11)-(13) for this poll.Estimated parameter values of models (11)-(13) for this poll are given in Table 5. By calculating the thresholds corresponding to this poll, we use the parameters’ values from this Table, we get
(13)RA0=1.1404>1,RD0=6.1346>1.Table 5
Estimated parameters value for models (1)-(3) for the second poll.
Parametersβ1β2α1α2γ1γ2Values0.12000.89000.48580.33620.08980.0001And we can see thatβ1γ2−β2γ1=−0.0799<0 and α=α1−α2>0, which means that the model presents the equilibrium e∗=0.0651,0.3868,0.5480 in stability, see the proof of Proposition 5 and Remark 2 in [8]. In this case, it may be noted that some control strategies are necessary to change the course of events prior to Election Day.Figure6 depicts the effectiveness of the estimation of the parameters, while in (a) we can see the indifferent individual curve of the model plotted with the solid line while the data of Indifferent individuals is plotted with small crosses. In (b), a comparison of the Agree people function of the model plotted with the solid line and the data of Agree individuals is plotted with small crosses. Also in (c), the estimated function of Disapproved people is plotted with the solid line and the data of disagree individuals is plotted with crosses. In Figure 7, the three functions of the model and the three states of the real data are plotted with different patterns to show that the model fits the data correctly.Figure 6
The Agree-Disagree model fits for the second poll data. (a) The solid line depicts the model output forI, while the small crosses show the indifferent data. (b) The solid line depicts the model output for A, while the small crosses show the Agree data. (c) The solid line depicts the model output for D, while the small crosses show the Disagree data.Figure 7
Three states of models (1)-(3) simulated with estimated parameter values from Table 5 and the real data of the second poll.


## 5. Sensitivity Analysis

Sensitivity analysis is measured using the sensitivity index. Sensitivity indices allow us to measure the relative change in a state variable when a parameter changes while the other parameters are retained fixed at constant values. The normalized forward sensitivity index of a variable to a parameter is the ratio of the relative change in the variable to the relative change in the parameter. When the variable is a differentiable function of the parameter, the sensitivity index may be alternatively defined using partial derivatives.We carry out local sensitivity indices ofRD0 and RA0 to the model parameters. A sensitivity index of 1 indicates that a variation of 10% of the parameter causes a variation of 10% of the final value of RD0 (or RA0). Indices less than 0 show a negative fluctuation, while those greater than 0 show a positive variation. These results help us to identify the most influential parameter on each threshold.


### 5.1. Sensitivity Analysis ofRD0

The thresholdRD0 is a function of the six parameters β1, β2, α1, α2, γ1, and γ2 where
(14)RD0=α2β1−α2γ1+β2γ1α1β1−α1γ1+β1γ2.A small perturbationδx to a parameter x and the corresponding change in RD0 as δRD0 is given by
(15)δRD0=RD0x+δx−RD0x=RD0x+δx−RD0xδxδx≅δx∂RD0∂x.Thus, the normalized sensitivity indexϕx is defined as
(16)ϕx=∂RD0/RD0∂x/x=xRD0.∂RD0∂x.


#### 5.1.1. The First Poll

The normalized sensitivity indices for the 6 parameters are obtained as(17)ϕβ1=−β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=−0.00062401,ϕβ2=β2γ1α2β1−α2γ1+β2γ1=0.00066317,ϕα1=−α1β1−γ1α1β1−α1γ1+β1γ2=−0.9016,ϕα2=α2β1−γ1α2β1−α2γ1+β2γ1=0.99934,ϕγ1=β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=0.00062401,ϕγ2=−β1γ2α1β1−α1γ1+β1γ2=−0.98399.Sinceϕβ1=−0.00062401, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0062401%. As ϕγ1=+0.00062401, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.0062401%. Similarly, as ϕα1=−0.9016, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.016%. Then again, as ϕβ2=+0.00066317, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.0066317%. And since ϕα2=+0.99934, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) RD0 by 9.9934%. Finally, as ϕγ2=−0.98399, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) RD0 by 9.8399%.


#### 5.1.2. The Second Poll

The normalized sensitivity indices for the 6 parameters using the values of the second poll are obtained as(18)ϕβ1=−3.5232,ϕβ2=0.8873,ϕα1=−0.9992,ϕα2=0.1127,ϕγ1=3.5232,ϕγ2=−8.1726×10−4.Sinceϕβ1=−3.5232, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 35.232%. As ϕγ1=+3.5232, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 35.232%. Similarly, as ϕα1=−0.9992, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.992%. Then again, as ϕβ2=+0.8873, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 8.873%. And since ϕα2=+0.1127, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 1.127%. Finally, as ϕγ2=−8.1726×10−4, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 8.1726×10−3%.


### 5.2. Sensitivity Analysis ofRA0

Also the thresholdRA0 is a function of the six parameters β1, β2, α1, α2, γ1, and γ2 where
(19)RA0=α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1.


#### 5.2.1. The First Poll

Using the same procedure as before, sensitivity indices for the 6 parameters are obtained as(20)ϕβ1=β1γ2α1β2−α1γ2+β1γ2=0.03131,ϕβ2=−β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=−0.033418,ϕα1=α1β2−γ2α1β2−α1γ2+β1γ2=0.96869,ϕα2=−α2β2−γ2α2β2−α2γ2+β2γ1=−0.9998,ϕγ1=−β2γ1α2β2−α2γ2+β2γ1=−0.00019649,ϕγ2=β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=0.033418.Asϕα1=+0.96869, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 9.6869%. And since ϕα2=−0.9998, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.998%. Similarly, As ϕγ1=−0.00019649, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0019649%. Since ϕβ1=+0.03131, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.3131%. Then again, as ϕβ2=−0.33418, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 0.33418%. Finally, as ϕγ2=+0.033418, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.33418%.


#### 5.2.2. The Second Poll

Using the same procedure as before and using the values corresponding to the second poll, sensitivity indices for the 6 parameters are obtained as(21)ϕβ1=2.7757×10−5,ϕβ2=−4.0700×10−6,ϕα1=1,ϕα2=−0.7892,ϕγ1=−0.2108,ϕγ2=4.0700×10−6.Asϕα1=+1, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 10%. And since ϕα2=−0.7892, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 7.892%. Similarly, as ϕγ1=−0.2108, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 2.108%. Since ϕβ1=+2.7757×10−5, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 2.7757×10−4%. Then again, as ϕβ2=−4.0700×10−6, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 4.0700×10−5%. Finally, as ϕγ2=+4.0700×10−6, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 4.0700×10−5%.A summary of sensitivity indices forRD0 and RA0 corresponding to the first poll is given in Figure 8, and to the second poll is given in Figure 9.Figure 8
Sensitivity indices for the first poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.Figure 9
Sensitivity indices for the second poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.


## 5.1. Sensitivity Analysis ofRD0

The thresholdRD0 is a function of the six parameters β1, β2, α1, α2, γ1, and γ2 where
(14)RD0=α2β1−α2γ1+β2γ1α1β1−α1γ1+β1γ2.A small perturbationδx to a parameter x and the corresponding change in RD0 as δRD0 is given by
(15)δRD0=RD0x+δx−RD0x=RD0x+δx−RD0xδxδx≅δx∂RD0∂x.Thus, the normalized sensitivity indexϕx is defined as
(16)ϕx=∂RD0/RD0∂x/x=xRD0.∂RD0∂x.


### 5.1.1. The First Poll

The normalized sensitivity indices for the 6 parameters are obtained as(17)ϕβ1=−β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=−0.00062401,ϕβ2=β2γ1α2β1−α2γ1+β2γ1=0.00066317,ϕα1=−α1β1−γ1α1β1−α1γ1+β1γ2=−0.9016,ϕα2=α2β1−γ1α2β1−α2γ1+β2γ1=0.99934,ϕγ1=β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=0.00062401,ϕγ2=−β1γ2α1β1−α1γ1+β1γ2=−0.98399.Sinceϕβ1=−0.00062401, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0062401%. As ϕγ1=+0.00062401, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.0062401%. Similarly, as ϕα1=−0.9016, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.016%. Then again, as ϕβ2=+0.00066317, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.0066317%. And since ϕα2=+0.99934, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) RD0 by 9.9934%. Finally, as ϕγ2=−0.98399, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) RD0 by 9.8399%.


### 5.1.2. The Second Poll

The normalized sensitivity indices for the 6 parameters using the values of the second poll are obtained as(18)ϕβ1=−3.5232,ϕβ2=0.8873,ϕα1=−0.9992,ϕα2=0.1127,ϕγ1=3.5232,ϕγ2=−8.1726×10−4.Sinceϕβ1=−3.5232, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 35.232%. As ϕγ1=+3.5232, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 35.232%. Similarly, as ϕα1=−0.9992, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.992%. Then again, as ϕβ2=+0.8873, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 8.873%. And since ϕα2=+0.1127, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 1.127%. Finally, as ϕγ2=−8.1726×10−4, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 8.1726×10−3%.


## 5.1.1. The First Poll

The normalized sensitivity indices for the 6 parameters are obtained as(17)ϕβ1=−β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=−0.00062401,ϕβ2=β2γ1α2β1−α2γ1+β2γ1=0.00066317,ϕα1=−α1β1−γ1α1β1−α1γ1+β1γ2=−0.9016,ϕα2=α2β1−γ1α2β1−α2γ1+β2γ1=0.99934,ϕγ1=β1γ1α1β2−α2γ2+β2γ2α1β1−α1γ1+β1γ2α2β1−α2γ1+β2γ1=0.00062401,ϕγ2=−β1γ2α1β1−α1γ1+β1γ2=−0.98399.Sinceϕβ1=−0.00062401, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0062401%. As ϕγ1=+0.00062401, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.0062401%. Similarly, as ϕα1=−0.9016, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.016%. Then again, as ϕβ2=+0.00066317, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.0066317%. And since ϕα2=+0.99934, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) RD0 by 9.9934%. Finally, as ϕγ2=−0.98399, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) RD0 by 9.8399%.


## 5.1.2. The Second Poll

The normalized sensitivity indices for the 6 parameters using the values of the second poll are obtained as(18)ϕβ1=−3.5232,ϕβ2=0.8873,ϕα1=−0.9992,ϕα2=0.1127,ϕγ1=3.5232,ϕγ2=−8.1726×10−4.Sinceϕβ1=−3.5232, thus, β1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 35.232%. As ϕγ1=+3.5232, thus, γ1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 35.232%. Similarly, as ϕα1=−0.9992, then, α1 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.992%. Then again, as ϕβ2=+0.8873, thus, β2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 8.873%. And since ϕα2=+0.1127, thus, α2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 1.127%. Finally, as ϕγ2=−8.1726×10−4, then, γ2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 8.1726×10−3%.


## 5.2. Sensitivity Analysis ofRA0

Also the thresholdRA0 is a function of the six parameters β1, β2, α1, α2, γ1, and γ2 where
(19)RA0=α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1.


### 5.2.1. The First Poll

Using the same procedure as before, sensitivity indices for the 6 parameters are obtained as(20)ϕβ1=β1γ2α1β2−α1γ2+β1γ2=0.03131,ϕβ2=−β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=−0.033418,ϕα1=α1β2−γ2α1β2−α1γ2+β1γ2=0.96869,ϕα2=−α2β2−γ2α2β2−α2γ2+β2γ1=−0.9998,ϕγ1=−β2γ1α2β2−α2γ2+β2γ1=−0.00019649,ϕγ2=β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=0.033418.Asϕα1=+0.96869, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 9.6869%. And since ϕα2=−0.9998, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.998%. Similarly, As ϕγ1=−0.00019649, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0019649%. Since ϕβ1=+0.03131, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.3131%. Then again, as ϕβ2=−0.33418, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 0.33418%. Finally, as ϕγ2=+0.033418, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.33418%.


### 5.2.2. The Second Poll

Using the same procedure as before and using the values corresponding to the second poll, sensitivity indices for the 6 parameters are obtained as(21)ϕβ1=2.7757×10−5,ϕβ2=−4.0700×10−6,ϕα1=1,ϕα2=−0.7892,ϕγ1=−0.2108,ϕγ2=4.0700×10−6.Asϕα1=+1, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 10%. And since ϕα2=−0.7892, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 7.892%. Similarly, as ϕγ1=−0.2108, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 2.108%. Since ϕβ1=+2.7757×10−5, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 2.7757×10−4%. Then again, as ϕβ2=−4.0700×10−6, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 4.0700×10−5%. Finally, as ϕγ2=+4.0700×10−6, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 4.0700×10−5%.A summary of sensitivity indices forRD0 and RA0 corresponding to the first poll is given in Figure 8, and to the second poll is given in Figure 9.Figure 8
Sensitivity indices for the first poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.Figure 9
Sensitivity indices for the second poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.


## 5.2.1. The First Poll

Using the same procedure as before, sensitivity indices for the 6 parameters are obtained as(20)ϕβ1=β1γ2α1β2−α1γ2+β1γ2=0.03131,ϕβ2=−β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=−0.033418,ϕα1=α1β2−γ2α1β2−α1γ2+β1γ2=0.96869,ϕα2=−α2β2−γ2α2β2−α2γ2+β2γ1=−0.9998,ϕγ1=−β2γ1α2β2−α2γ2+β2γ1=−0.00019649,ϕγ2=β2γ2α2β1−α1γ1+β1γ1α1β2−α1γ2+β1γ2α2β2−α2γ2+β2γ1=0.033418.Asϕα1=+0.96869, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 9.6869%. And since ϕα2=−0.9998, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 9.998%. Similarly, As ϕγ1=−0.00019649, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 0.0019649%. Since ϕβ1=+0.03131, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 0.3131%. Then again, as ϕβ2=−0.33418, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 0.33418%. Finally, as ϕγ2=+0.033418, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 0.33418%.


## 5.2.2. The Second Poll

Using the same procedure as before and using the values corresponding to the second poll, sensitivity indices for the 6 parameters are obtained as(21)ϕβ1=2.7757×10−5,ϕβ2=−4.0700×10−6,ϕα1=1,ϕα2=−0.7892,ϕγ1=−0.2108,ϕγ2=4.0700×10−6.Asϕα1=+1, then, α1 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 10%. And since ϕα2=−0.7892, thus, α2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 7.892%. Similarly, as ϕγ1=−0.2108, thus, γ1 is decreasing (or increasing) by a 10% increase (or decreases) in RD0 by 2.108%. Since ϕβ1=+2.7757×10−5, thus, β1 is decreasing (or increasing) by a 10% decrease (or increases) in RD0 by 2.7757×10−4%. Then again, as ϕβ2=−4.0700×10−6, thus, β2 is increasing (or decreasing) by a 10% decrease (or increases) in RD0 by 4.0700×10−5%. Finally, as ϕγ2=+4.0700×10−6, then, γ2 is increasing (or decreasing) by a 10% increase (or decreases) in RD0 by 4.0700×10−5%.A summary of sensitivity indices forRD0 and RA0 corresponding to the first poll is given in Figure 8, and to the second poll is given in Figure 9.Figure 8
Sensitivity indices for the first poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.Figure 9
Sensitivity indices for the second poll. (a) Sensitivity indices forRD0. (b) Sensitivity indices for RA0.


## 6. Model Uncertainty

A covariance matrix estimates the variability of the model parameters and hence random disturbances in the output. Thus, covariance matrices contain information about the uncertainties in the model output. All covariance matrices are symmetrical. The absolute value giving accuracy information is contained in the diagonal elements of the covariance matrix. All other elements in covariance matrices can be used to study the interrelationships. If the covariance between any two coefficients is positive, then the values of the coefficients tend to vary in a positive way. On the other hand, if the covariance between any two coefficients is negative, then the coefficients’ values tend to move in opposite directions.


### 6.1. The First Poll

The variance of each parameters is shown in the leading diagonal elements of the covariance matrix∑1 given in the Appendix. The parameter α2 has the largest variance (1019328818) followed by the parameter α1 (1019328472), the parameter β2 (5.9071), the parameter β1 (1.2486), the parameter γ2 (0.0152), and the parameter γ1 (0.0034).Covariance matrices contain negative and positive elements; negative covariance means that the parameters’ values on those parameters tend to move in opposite directions which means that if one parameter’s value increases, then other parameter’s value starts decreasing, while positive covariances mean that these parameters’ values tend to vary in a positive way which means that if one parameter’s value increases, then other parameter’s value also increases. The covariance matrix∑1 contains only positive elements. This means that these parameters’ values vary in a positive way.As shown in the covariance matrix∑1, the covariance between transmission coefficients β1 and β2 is 0.9856, the covariance between transmission coefficients β1 and α1 is 2703.3959, the covariance between transmission coefficients β1 and α2 is 2703.3971, the covariance between the transmission coefficient β1 and the parameter γ1 is 0.0609, and the covariance between the transmission coefficient β1 and the parameter γ2 is 0.0494.The covariance between the transmission coefficientβ2 and the parameter α1 is 27274.0064, the covariance between the transmission coefficient β2 and the parameter α2 is 27274.0112, the covariance between the transmission coefficient β2 and the parameter γ1 is 0.0467, and the covariance between the transmission coefficient β2 and the parameter γ2 is 0.2942.The covariance between transmission coefficientsα1 and α2 is 1019328645, the covariance between transmission coefficients α1 and γ1 is 72.0217, and the covariance between transmission coefficients α1 and γ2 is 1483.2081.The covariance between transmission coefficientsα2 and γ1 is 72.0206, and the covariance between transmission coefficients α2 and γ2 is 1483.2095. Finally, the covariance between transmission coefficients γ1 and γ2 is 0.0018.


### 6.2. The Second Poll

The parameterα1 has the largest variance (10927264.36) followed by the parameter α2 (10927020.6), the parameter β1 (0.002112), the parameter β2 (0.001307), the parameter γ2 (0.001046), and the parameter γ1 (0.000787).The covariance matrix∑2 given in the Appendix contains negative and positive elements.As shown in the covariance matrix∑2, the covariance between transmission coefficients β1 and β2 is positive (0.000104), the covariance between transmission coefficients β1 and α1 is positive (64.542655), the covariance between transmission coefficients β1 and α2 is positive (64.540946), the covariance between the transmission coefficient β1 and the parameter γ1 is positive (0.001141), the covariance between the transmission coefficient β2 and the parameter γ2 is positive (0.000968), the covariance between transmission coefficients α1 and α2 is positive (10927142.47), the covariance between the transmission coefficient α1 and the parameter γ1 is positive (55.366504), and the covariance between the transmission coefficient α2 and the parameter γ1 is positive (55.3651). This means that these parameters’ values vary in a positive way.The covariance between the transmission coefficientβ1 and the parameter γ1 is negative (−0.000676), the covariance between transmission coefficients β2 and α1 is negative (−67.701944), the covariance between transmission coefficients β2 and α2 is negative (−67.700578), the covariance between the transmission coefficient β2 and the parameter γ1 is negative (−0.000362), the covariance between the transmission coefficient α1 and the parameter γ2 is negative (−70.976139), the covariance between the transmission coefficient α2 and the parameter γ2 is negative (−70.974401), and the covariance between parameters γ1 and the parameter γ2 is negative (−0.00067). This reflect that these parameters’ values tend to vary in a negative way. Table 6 shows the summary of covariance of the Agree-Disagree model parameters for the poll data.Table 6
Summary of covariance.
CovarianceDirectionCovarianceDirectionCovarianceDirectionCβ1,β2+veCβ2,α1−veCα1,γ1+veCβ1,α1+veCβ2,α2−veCα1,γ2−veCβ1,α2+veCβ2,γ1−veCα2,γ1+veCβ1,γ1+veCβ2,γ2+veCα2,γ2−veCβ1,γ2−veCα1,α2+veCγ1,γ2+ve


## 6.1. The First Poll

The variance of each parameters is shown in the leading diagonal elements of the covariance matrix∑1 given in the Appendix. The parameter α2 has the largest variance (1019328818) followed by the parameter α1 (1019328472), the parameter β2 (5.9071), the parameter β1 (1.2486), the parameter γ2 (0.0152), and the parameter γ1 (0.0034).Covariance matrices contain negative and positive elements; negative covariance means that the parameters’ values on those parameters tend to move in opposite directions which means that if one parameter’s value increases, then other parameter’s value starts decreasing, while positive covariances mean that these parameters’ values tend to vary in a positive way which means that if one parameter’s value increases, then other parameter’s value also increases. The covariance matrix∑1 contains only positive elements. This means that these parameters’ values vary in a positive way.As shown in the covariance matrix∑1, the covariance between transmission coefficients β1 and β2 is 0.9856, the covariance between transmission coefficients β1 and α1 is 2703.3959, the covariance between transmission coefficients β1 and α2 is 2703.3971, the covariance between the transmission coefficient β1 and the parameter γ1 is 0.0609, and the covariance between the transmission coefficient β1 and the parameter γ2 is 0.0494.The covariance between the transmission coefficientβ2 and the parameter α1 is 27274.0064, the covariance between the transmission coefficient β2 and the parameter α2 is 27274.0112, the covariance between the transmission coefficient β2 and the parameter γ1 is 0.0467, and the covariance between the transmission coefficient β2 and the parameter γ2 is 0.2942.The covariance between transmission coefficientsα1 and α2 is 1019328645, the covariance between transmission coefficients α1 and γ1 is 72.0217, and the covariance between transmission coefficients α1 and γ2 is 1483.2081.The covariance between transmission coefficientsα2 and γ1 is 72.0206, and the covariance between transmission coefficients α2 and γ2 is 1483.2095. Finally, the covariance between transmission coefficients γ1 and γ2 is 0.0018.


## 6.2. The Second Poll

The parameterα1 has the largest variance (10927264.36) followed by the parameter α2 (10927020.6), the parameter β1 (0.002112), the parameter β2 (0.001307), the parameter γ2 (0.001046), and the parameter γ1 (0.000787).The covariance matrix∑2 given in the Appendix contains negative and positive elements.As shown in the covariance matrix∑2, the covariance between transmission coefficients β1 and β2 is positive (0.000104), the covariance between transmission coefficients β1 and α1 is positive (64.542655), the covariance between transmission coefficients β1 and α2 is positive (64.540946), the covariance between the transmission coefficient β1 and the parameter γ1 is positive (0.001141), the covariance between the transmission coefficient β2 and the parameter γ2 is positive (0.000968), the covariance between transmission coefficients α1 and α2 is positive (10927142.47), the covariance between the transmission coefficient α1 and the parameter γ1 is positive (55.366504), and the covariance between the transmission coefficient α2 and the parameter γ1 is positive (55.3651). This means that these parameters’ values vary in a positive way.The covariance between the transmission coefficientβ1 and the parameter γ1 is negative (−0.000676), the covariance between transmission coefficients β2 and α1 is negative (−67.701944), the covariance between transmission coefficients β2 and α2 is negative (−67.700578), the covariance between the transmission coefficient β2 and the parameter γ1 is negative (−0.000362), the covariance between the transmission coefficient α1 and the parameter γ2 is negative (−70.976139), the covariance between the transmission coefficient α2 and the parameter γ2 is negative (−70.974401), and the covariance between parameters γ1 and the parameter γ2 is negative (−0.00067). This reflect that these parameters’ values tend to vary in a negative way. Table 6 shows the summary of covariance of the Agree-Disagree model parameters for the poll data.Table 6
Summary of covariance.
CovarianceDirectionCovarianceDirectionCovarianceDirectionCβ1,β2+veCβ2,α1−veCα1,γ1+veCβ1,α1+veCβ2,α2−veCα1,γ2−veCβ1,α2+veCβ2,γ1−veCα2,γ1+veCβ1,γ1+veCβ2,γ2+veCα2,γ2−veCβ1,γ2−veCα1,α2+veCγ1,γ2+ve


## 7. Conclusion

In this paper, an IAD model-type compartmental model has been considered to explore agree-disagree opinions during polls. The next-generation matrix method is used to derive thresholds of equilibrium stabilityRA0 and RD0. Two sets of data from the Reuters polling system regarding the approval rating of the U.S presidential in two terms are used to estimate the parameters and to validate the underlying model.To identify the most influential parameter in the proposed model, local sensitivity analysis is carried out. We calculated sensitivity indices based on the estimated parameters to identify the most influential parameters for both polls. We found in the first poll that the most important parameters in determining the thresholdsRD0 are α1, disagree to agree transmission rate, α2, agree to disagree transmission rate, and the parameter γ2, of the loss of interest of the disagreeing people, while the most influential parameters on the threshold RA0 are α1, disagree to agree transmission rate, and α2, agree to disagree transmission rate.Regarding the second poll, we found that the most important parameters in determining the thresholdsRD0 are β1, the polarization parameter of indifferent people by approved people, the parameter, γ1, of the loss of interest of the approved people, and α1, disagree to agree transmission rate, while the most influential parameters on the threshold RA0 are α1, disagree to agree transmission rate, α2, agree to disagree transmission rate, and the parameter γ1, of the loss of interest of the agree people. Uncertainty analysis is carried out to determine the relationship between the different parameters of the model in both polls.



---

*Source: 1837364-2020-07-02.xml*