A Fractional-Order SIR-C Cyber Rumor Propagation Prediction Model with a Clariﬁcation Mechanism

: As communication continues to develop, the high freedom and low cost of the communication network environment also make rumors spread more rapidly. If rumors are not clariﬁed and controlled in time, it is very easy to trigger mass panic and undermine social stability. Therefore, it is important to establish an efﬁcient model for rumor propagation. In this paper, the impact of rumor clariﬁers on the spread of rumors is considered and fractional order differentiation is introduced to solve the problem that traditional models do not take into account the "anomalous propagation" characteristics of information. A fractional-order Susceptible-Infected-Removal-Clarify (SIR-C) rumor propagation prediction model featuring the clariﬁcation mechanism is proposed. The existence and asymptotic stability conditions of the rumor-free equilibrium point (RFEP) E 0 ; the boundary equilibrium points (BEPs) E 1 and E 2 are also given. Finally, the stability conditions and practical cases are veriﬁed by numerical simulations. The experimental results conﬁrm the analysis of the theoretical study and the model ﬁts well with the real-world case data with just minor deviations. As a result, the model can play a positive and effective role in rumor propagation prediction.


Introduction
Cyber rumors [1] are rumors propagated through the Internet medium that are not based on facts and are offensive and purposeful. As communication technology has steadily developed, large social networking sites, represented by Twitter and Weibo, have become an important source of information reception for people. The China Internet Network Information Center's 50th Statistical Report [2] on the Development Status of the Internet in China states that as of June 2022, there were 1.051 billion Internet users in China, with 99.6% of those users accessing the Internet using mobile devices. This high degree of freedom, low cost of dissemination and the huge scale of users have also led to the "explosive" spread of cyber rumors. If cyber rumors are not clarified and controlled in a timely manner, but are instead spread excessively, they can easily trigger irrational behaviour such as panic and excessive anxiety among the public, disrupting public order, destabilising society and even affecting the credibility of the state. Hence, it is of significant theoretical and practical importance to analyse the mechanisms and laws of rumor propagation and to establish an effective model of cyber rumor propagation.
Research into rumor propagation models began in the 1960s with the DK model [3], which divided the audience into three categories according to rumor propagation: ignorants, spreaders and removers, and used stochastic process analysis to analyse rumor propagation. Subsequently, Maki and Thomson [4] proposed the MT model found on the DK model and analysed the model mathematically. These two models produce representative results and serve as the theoretical foundation for rumor propagation models. However, there are still differences in the transmission mechanisms and pathways between rumor propagation and disease propagation, and they are not fully applicable to the propagation of cyber rumors. Since then, in order to obtain a more precise description of the dynamic process and individual characteristics of rumor spread, a number of researchers have started to consider the use of network topology to improve the traditional epidemic model Susceptible-Infected-Removal (SIR), and have contributed significantly by using it to analyze the spread of cyber rumors. Among them, the influence of human behaviour on the rumor spreading process has received extensive attention and research from scholars. Examples include the forgetting mechanism [5][6][7], hesitation mechanism [8,9], and leader's opinion [10,11], etc. Some studies [12] have also found that the pattern of human activities can have an impact on the spread of information, making it slower.
In recent years, influenced by COVID- 19, there has been a renewed interest in infectious disease models and rumor propagation models, both domestically and internationally, and many scholars have extended the models for more in-depth study. The authors of [13] developed a model of information propagation with simultaneous censorship, sharing, collection, and suppression mechanisms. The authors of [14] provide an evolutionary game model that considers how individual choices affect the spread and management of rumors. The authors of [15] contend that the spread of rumors is influenced by the state of scientific knowledge. The authors of [16] provide a model with an influence mechanism, taking into account the likelihood of propagation in rumor transmission is not fixed, but is influenced by the number of current propagations. The authors of [17] address the current gap in research on the propagation power of rumors based on a function of two types of characteristics of content: false rumors and true rumors.
Although the above articles make an important contribution to rumor propagation, there are still problems with these models. Firstly, most of the models are built using the traditional integer-order. However, the integer-order model does not show the "anomalous propagation" of the actual information spread, that is, the rumor spreads explosively fast in the early stages and slows down in the later stages. The fractional-order model is a good solution to this problem. Fractional-order calculus provides a tool to describe the genetic and memory impacts of different materials. The memorability and heritability of fractional-order differentiation is theoretically proven and widely used [18][19][20][21]. For the rumor propagation model, memorability is an important feature of immunity. When a rumor spreads, people "remember" the clarifying information once they have received it in response to a successful immune response to the rumor, and this memorability plays an important role in stopping the spread of the rumor. Furthermore, the determination of future states throughout the spread of rumors depends on the historical process of propagation. The way in which a person reacts to external influences also depends on the experience he has accumulated in the past, and fractional-order differentiation is a very natural tool for modelling heritability, and is also considered to be the best method for modelling transmission [22]. Secondly, there are fewer existing rumor models that provide a mathematical analysis of their global asymptotic stable(GAS). Therefore, this study and gives local asymptotic stable(LAS) and GAS conditions for the rumor-free equilibrium point (RFEP) E 0 and the boundary equilibrium points (BEPs) E 1 , E 2 of the fractional-order model. Furthermore, most research make the assumption that rumors spread within a closed system and do not account for the inflow and outflow of the population. In fact, social networks are open platforms that need to take into account the impact of individual flows during the propagation process. Using the analysis just mentioned, this paper investigates a fractional-order Susceptible-Infected-Removal-Clarify (SIR-C) cyber rumor propagation prediction model with a clarification mechanism.
To sum up, the main contributions of this research are fourfold:

1.
We propose fractional-order differentiation to solve the problem of ''anomalous propagation" of cyber rumors;

2.
We propose a mechanism containing clarification considering the influence of clarifiers in the spread of cyber rumors; 3.
We propose LAS and GAS conditions for the RFEP (E 0 ) and the BEPs (E 1 , E 2 ) of the model; 4.
We consider the population inflow and outflow of rumor propagation in social networks.
The remainder of the paper is structured as follows. In Section 2, we describe a fractional-order SIR-C cyber rumor propagation prediction model featuring the clarification mechanism, and introduce the mathematical properties of the corresponding fractionalorder differential equations. In Section 3, the RFEP (E 0 ) and the BEPs (E 1 , E 2 ) are obtained and their stability is discussed. In Section 4, numerical simulations are used to verify the theoretical study and the usefulness of the model is verified by the actual event of ''Record low average maths score in 2022 National College Entrance Examnation(NEMT)". A brief conclusion of the paper is given in Section 5.

Formulation of Fractional-Order SIR-C Model with Basic Mathematical Properties
In this section, the model's fractional-order differential equations, the structure, and the basic mathematical properties are systematically introduced.

Conformable Fractional Derivative (CFD)
In the last few decades, fractional-order differential models have been popular in several domains such as engineering and science due to their unique memory effects and genetic properties [23], no longer applied only in pure mathematics. Typical definitions of fractional-order derivatives are the Riemann-Liouville (RL) [24], the Caputo [25] and the CFD [26]. The deficiency of RL and Caputo is that they lose the basic properties such as multiplication and chaining that general derivatives have [27]. In this paper, CFD is used, which satisfies the product, quotient and chain rules for non-linear derivatives and will make the model calculation easier. And it can be converted directly to an equation of integerorder, which facilitates subsequent model comparisons, which facilitates subsequent model comparisons [28]. The authors of [29,30] demonstrate the computational simplicity of CFD compared to other fractional order differentiations. CFD is a new well-behaved and simple definition that has the ability to successfully overcome some of the drawbacks of the traditional definition. So this paper takes advantage of CFD in applied problems.CFD fractional-order derivatives are defined as follows.
if f is α-differentiable in some (0, α), and lim t→0 + f α (t) exist, then define Sometimes, we will write f α (t) for D α f (t), to denote the CFD of f of order α.
We should be aware that D α (t p ) = pt p−α , Further, the definition of CFD is consistent with the classical definition of RL and Caputo on polynomials and has the multiplication and chain laws. Theorem 1. [26] Let α ∈ (0, 1], and f be α-differentiable at a point t > 0. If, in addition, f is differentiable. Then,

Formulation of the Model
It is well known that persuasive rumor-clearing messages have a significant effect in curbing the spread of rumors [14,31]. Moreover, social networking is an open-ended interactive platform, with users registering or logging out at any time, and a large number of users are attracted by rumors and flock to social networks to discuss and spread the topic. The rumor propagation model takes into account the inflow and outflow of people on the rumor topic, and combines the memory effect and genetic properties of fractional-order calculus to illustrate the ''anomalous propagation" of rumors. Finally, a fractional-order cyber rumor propagation prediction model with a clarification mechanism is formulated.
The model is based on the following assumptions: • In social networking sites, we assume that there is a fixed rate of population inflow and outflow, and the quantity of individuals entering each time unit remains constant, i.e. the inflow rate is Λ and the outflow rate is d; • The total number of transmission subjects N changes over time and is divided into three divisions: S (susceptible), which refers to nodes that not received information about the rumor; I (infected), which refers to nodes that receive the rumor and propagate it; C ( clarify), which refers to people who learn the truth and clear up the rumor; and R ( removal), which refers to nodes that immune to the cyber rumor and not propagate it; • Suppose that susceptible persons are converted to infected persons (transmitters, clarifiers) by a proportionality factor ( spread rate) β as a result of being affected by the spread of the rumor. The setting for the rumor transmission rate is β 1 and the rumor clarification rate is β 2 . The conversion rate of spreaders to clarifiers is set to µ; • Infected persons ( spreaders, clarifiers) are converted to immune persons by a scaling factor (removal rate) γ. The removal rate for spreaders is set to γ 1 and the removal rate for clarifiers is set to γ 2 ; • All parameters are constants with non-negative values between 0 and 1.
Based on the above assumptions, we can obtain the basic scheme of the fractional order SIR-C model, as seen in Figure 1. The specific details of the S, I, C and R nodes and the transformation relationships of the four nodes are described below: • S-nodes refer to those that have not received the cyber rumor and are not yet affected by it. Some S-nodes may receive rumor information and be influenced by it, listen to the rumor and transform into I-nodes with spread rate β 1 to spread cyber rumor. Some S-nodes may receive rumor clarification information, learn the truth, and transform into C-nodes at spread rate β 2 to spread the truth about the facts; • I-nodes refer to those that receive cyber rumors, listen to them and spread them. Some I-nodes are likely to learn the truth and transform to C-nodes with conversion rate µ, clarify the rumor and spread it. Some I-nodes may lose interest in the rumor topic as time passes and are transformed to an immune and indifferent state at a removal rate γ 1 ; • C-nodes refer to those that know the truth and clarify the rumor. When the C-nodes have finished clarifying the rumor, they will be immune to the rumor at a removal rate γ 2 and will no longer pay attention to it; • R-nodes refer to those that have received the information, lost interest in the rumor topic and will not have any further influence on the rumor topic.
Based on the above analysis of the model node transition rules, the differential equation for the fractional-order SIR-C model can be derived as Equation (4).
where, S(t), I(t), C(t) and R(t) are denoted as the quantities of S, I, R and C at moment t respectively. S(0), The number of individuals across all social networks is a in social networks is N(t), i.e., when, Equation (4) is substituted into Equation (5), we get Thus The positive variable set of the system (4) is Φ = (S, I, C, R) ∈ R + 4 : S + I + C + R Λ d . Obviously, R(t) does not affect the first three equations in system (4), and the R state node appears only in the fourth equation, which is determined by S, I and C. The dimensionality of system (4) can be reduced to eliminate the R state node in the differential equation. The simplified model is as follows.
Thus, the feasible region of the system (7) can be described as a closed positively invariant set, denoted by Ω = (S, I, C) ∈ R + 3 : A subsequent stability analysis of the fractional order SIR-C model will use system (7).

System Equilibrium Points and Stability Analysis
In this part, equilibrium points and basic reproduction number(BRN) of fractionalorder cyber rumor propagation models are computed, as well as the stability of cyber rumor propagation equilibrium states is discussed using Lyapunov's stability theorem and the Routh-Hurwitz criterion. This section focuses on the local and global stability of three equilibrium points, the RFEP (E 0 ) and the BEPs (E 1 , E 2 ).

RFEP(E 0 )
System (7) has a RFEP E 0 = (S 0 , I 0 , C 0 ), where there are neither rumor spreaders nor rumor clarifiers. Let System (7)'s rightmost component be zero. The new equation is expressed as follows: According to the definition of the RFEP, under the condition I 0 = 0, C 0 = 0. We obtain the unique RFEP.
Next, using the concept of the BRN R 0 [32] from the classical infectious disease model as an important parameter for determining whether a rumor can be quelled. R 0 is defined as the quantity of people a cyber rumor spreader can convert susceptible people into cyber rumor spreaders in the process of spreading. Subsequent analyses of stability were carried out based on R 0 .
The next generation matrix method [33] is applied to solve for R 0 , which is referred to as the next-generation matrix's spectral radius.
Let X(t) = (S(t), I(t), C(t)) T , system (7) is rewriteable as where, where, Therefore, the next generation FV −1 can be calculated as follows: The R 0 of system (7), which is the spectral radius of the matrix FV −1 , we obtain , where, .
Considering the upper triangular matrix's characteristics, the eigenvalues of the matrix J(E 0 ) can be obtained.

Theorem 3. For system
Proof. Let X(t) = (S(t), I(t), C(t)) T , According to the Lyapunov stability theorem, the Lyapunov function V(t) is constructed.
Next, by taking CFD in time of the function V(t) ,using the basic arithmetic properties of CFD and the define of BRN, we immediately obtain. where, from the basic inequality we know that.
The formula takes the equal sign when and only when dS(t) = Λ, i.e. S(t) = S 0 = Λ d . Therefore, we directly get that If R 01 < 1 and R 02 < 1, then, for all X(t) ∈ Ω, when X(t) = E 0 , we have D α V(t) < 0. According to the LaSalle invariant principle [35], E 0 is GAS.
According to the defintion of E 1 , under the condition I 0 = 0, brought into Equation (8). We obtain E 1 .
Futher, calculate the Jacobi matrix of F (x), Ψ(x) at E 1 , repectively, we have where, The BRN R 1 of system (7), which is the spectral radius of the matrix (F 1 V 1 ) −1 , we obtain Theorem 4. For system (7), if R 1 < 1. Then, the BEP E 1 is LAS.
Proof. The Jacobi matrix of system (7) at E 1 = (S 1 , I 1 , as follows. The characteristic equation of the matrix J(E 1 ) can be written as where, from the Equation (13), we can easily obtain the eigenvalues λ 1 .
Proof. Let X(t) = (S(t), I(t), C(t)) T , According to the Lyapunov stability theorem, the Lyapunov function V 1 (t) is constructed.
Next, by taking CFD in time of the function V 1 (t) ,using the basic arithmetic properties of CFD and the define of BRN, we immediately obtain. where, from the basic inequality we know that.
The formula takes the equal sign when and only when S 1 = S(t). Therefore, we directly get that where, If R 1 < 1, for all X(t) ∈ Ω, when X(t) = E 1 , we have D α V 1 (t) < 0. According to the LaSalle invariant principle, E 1 is GAS. From Equation (11), the E 1 exists if Overall, if R 02 > 1 and R 1 < 1. Then, the BEP E 1 exists and is GAS.
According to the defintion of E 2 , under the condition C 0 = 0, brought into Equation (8). We obtain the E 2 .
Futher, calculate the Jacobi matrix of F (x), Ψ(x) at E 2 , repectively, we have where, The BRN R 2 of system (7), which is the spectral radius of the matrix (F 2 V 2 ) −1 , we obtain Theorem 6. For system (7), if R 2 < 1. Then, the BEP E 2 is LAS.
Proof. The Jacobi matrix of system (7) at , 0 as follows.
The characteristic equation of the matrix J(E 2 ) can be obtained as Through simplification, we obtain where, from the Equation (17), we can easily obtain the eigenvalues λ 1 .
The remaining two eigenvalues λ 2 ans λ 3 satisfy the following equation.
Proof. Let X(t) = (S(t), I(t), C(t)) T , According to the Lyapunov stability theorem, the Lyapunov function V 2 (t) is constructed.
Next, by taking CFD in time of the function V 2 (t) ,using the basic arithmetic properties of CFD and the define of BRN, we immediately obtain. where, from the basic inequality we know that.
The formula takes the equal sign when and only when S 2 = S(t). Therefore, we directly get that where, If R 2 < 1, for all X(t) ∈ Ω, when X(t) = E 2 , we have D α V 2 (t) < 0. According to the LaSalle invariant principle, the BEP E 2 is GAS. From Equation (15), the E 2 exists if d(R 01 −1) β 1 > 0, i.e., R 01 > 1. Overall, if R 01 > 1 and R 2 < 1. Then, the BEP E 2 exists and is GAS.

Numerical Simulations and Discussions
In this part, Matlab is used as the simulation platform to numerically simulate the System (7) to verify the results of the RFEP (E 0 ) and the BEPs (E 1 , E 2 ) are in agreement with Theorems 3,5 and 7. The prediction results of the fractional-order SIR-C model and the traditional SIR model are compared to validate the feasibility of the model using the actual cyber rumor ''Record low average maths score in 2022 NEMT". The descriptions of the key parameters mentioned in this section are specified in Table 1. Using the data in Table 2, the global asymptotic stability of E 0 was simulated. where the initial value of the propagation density for each node in the system is set as: S(0) = 0.7, I(0) = 0.2, C(0) = 0.1, R(0) = 0. The different differential orders α ∈ {0.7, 1, 1.2} were set to verify the effect of different α on the fractional-order SIR-C rumor propagation model. The simulation results are shown in Figure 2.  According to the data in Table 2, it can be obtained that , R 01 = 0.833 < 1, and R 02 = 0.667 < 1. From Figure 2a, the system (7) will gradually converge to the E 0 = (1, 0, 0) as time goes on, i.e., E 0 is GAS. This is in agreement with Theorem 3.
Meanwhile, we observe that over time, S tends to 1 and I and C tend to 0. This means that the density of nodes spreading rumors will eventually shift to zero and cease to fluctuate, and the rumors will eventually be eliminated. That is, system (7) reaches a steady state where there are no rumor spreaders and clarifiers, only susceptible people. At the same time, it can be found that the fewer rumor spreaders there are, the fewer rumor clarifiers there will be, which is consistent with the actual rumor propagation process in social networking sites. The authors of [37] have analyzed the actual rumor spreading process of ''Does KFC sell rat?" in social networking sites. The analysis of the rumor messages and counter-rumor messages volumes obtained according to the conclusions of [37] is in agreement with the results of the theoretical variation of rumor spreaders and rumor clarifiers in this paper. The veracity of the results of this numerical simulation is further verified.
With Figure 2b, we observe that the fractional-order model is equivalent to that of the traditional integer-order system when the value of α is taken as 1. As the value of α decreases, the curve of rumor propagation converges more slowly. From the perspective of pubic opinion propagation, its interpretation suggests that such rumors take longer to eradicate.

Stability Simulation and Analysis of BEP E 1
Using the data in Table 3, the global asymptotic stability of E 1 was simulated. where the initial value of the propagation density for each node in the system is set as: S(0) = 0.8, I(0) = 0.1, C(0) = 0.1, R(0) = 0. The different differential orders α ∈ {0.6, 0.8, 1.1} were set to verify the effect of different α on the fractional-order SIR-C rumor propagation model. The simulation results are shown in Figure 3. Table 3. Parameter for the global asymptotic stability of E 1 .  According to the data in Table 3, it can be obtained that , R 02 = 1.90 > 1, and R 1 = 0.152 < 1. From Figure 3a, the system (7) will gradually converge to the E 1 = (0.525, 0, 0.45) as time goes on, i.e., E 1 is GAS. This is in agreement with Theorem 5.
Meanwhile, we observe that over time, S and C converge to a constant and I converges to 0. This means that the density of nodes spreading rumors will eventually shift to 0 and cease to fluctuate, and the rumors will eventually be eliminated. This means that system (7) reaches a steady state where there are no rumor spreaders, only susceptibles and clarifiers. It can also be found that the more rumor clarifiers there will be, the fewer rumor spreaders there will be, which is consistent with the actual rumor propagation process in social networking sites.
Through Figure 3b, we observe that the curve of rumor propagation converges more slowly as the value of α decreases, a situation that may affect the longer duration of rumor propagation.

Stability Simulation and Analysis of BEP E 2
Using the data in Table 4, the global asymptotic stability of E 2 was simulated. where the initial value of the propagation density for each node in the system is set as: S(0) = 0.8, I(0) = 0.1, C(0) = 0.1, R(0) = 0. The different differential orders α ∈ {0.9, 1, 1.3} were set to verify the effect of different α on the fractional-order SIR-C rumor propagation model. The simulation results are shown in Figure 4.    (7). (b) Consider the effect of different differential orders of α ∈ {0.9, 1, 1.3} on the spread of cyber rumors, for system (7).
According to the data in Table 4, it can be obtained that , R 01 = 1.667 > 1, and R 2 = 0.687 < 1. From Figure 4a, the system (7) will gradually converge to E 2 = (0.6, 0.38, 0) as time goes on, i.e. E 2 is GAS. This is in agreement with Theorem 7.
Meanwhile, we observe that over time, S and I converge to a constant and C converges to 0. This means that the density of nodes clarifying rumors will eventually shift to 0 and cease to fluctuate and rumors are continuously spread. That is, the system (7) reaches a steady state where there are no rumor clarifiers, only susceptibles and spreaders. Under this parameter condition, the rumor-clarifying mechanism cannot suppress rumors in a timely and effective manner.
Through Figure 4b, we observe that the curve of rumor propagation converges more slowly as the value of α decreases, a situation that may affect the longer duration of rumor propagation.

Experimental Simulation of Actual Cyber Rumor Propagation
Based on the fractional-order SIR-C model, the experimental simulation of the propagation process of "Record low average maths score in 2022 NEMT", an actual Internet rumor, was carried out, and the experiments were compared with the traditional integer-order SIR model to verify the effectiveness of the model.
The data for this experiment are the cyber rumors about ''Record low average maths score in 2022 NEMT" forwarded on social network platforms from 08:00 on June 22th to 20:00 on June 24th, 2022. Since this event involves keywords that are easily noticed by the society, such as education and examination, it has caused more controversies and concerns, so it was chosen as the validation experiment for this model.
Considering the obtained data, the initial value of the propagation density for each node in the system is assumed to be: S(0) = 0.98, I(0) = 0.01, C(0) = 0.01, R(0) = 0. Also, in order to obtain the best fit, the least squares method is used for data fitting in this paper. Fitting curve based on the actual data, the parameter values of each node can be calculated as shown in Table 5. The actual propagation process of rumors in social network platforms is usually ''anomalous propagation", but the integer-order model cannot predict ''anomalous propagation". Therefore, the fractional-order SIR-C model is proposed to split the prediction time, i.e., different differential orders are used to represent the propagation speed of rumors in the early, middle and late stages of propagation, so that the prediction curve is closer to the real data curve and the experimental results are more realistic.
According to the characteristics of each stage of the cyber rumor spreading process, the dichotomous search method is used to set different a values in different time stages. When t ∈ [0, 10], α = 1.2 is set to satisfy the accelerated spread of rumors in the early stage of rumor spread; when t ∈ (10, 20), α = 1 is set to satisfy the middle stage of rumor spread, where the speed of rumor propagation is reduced due to the presence of rumor clarifiers. When t ∈ (20, 60), α = 0.9 is set to satisfy the slowdown of rumor spread in the late stage of rumor spread as time passes and people lose interest in the rumor. The simulation diagram of the fractional-order SIR-C model results is shown in Figure 5a. In Figure 5b, for the traditional SIR model, rumor clarifiers do not exist in the whole rumor propagation process because there is no rumor clarification mechanism. However, in the model established in this paper, although there is a peak in the node density of cyber rumor propagation over time, it is significantly smaller than the peak in the traditional SIR model. Finally, with the action of rumor clarification mechanism and the passage of time, the rumor will be cleared eventually. We can find that the fractional-order SIR-C model proposed in this paper outperforms the existing SIR models. The model offers a more accurate approximation to actual data, and experimental results fit better with less error.

Conclusions
In this paper, based on the traditional SIR model, the rumor clarification mechanism is added and a fractional-order differential equation is introduced to establish a fractionalorder SIR-C rumor spread model that takes into account "anomalous propagation" of information. Next, the LAS and GAS of the RFEP (E 0 ) and the BEPs (E 1 , E 2 ) are analyzed by using the Routh-Hurwitz criterion and Lyapunov's stability theorem. Next, the validity of the above theoretical analysis is verified by numerical simulations, and a comparison with the existing SIR rumor propagation model is completed with actual rumor propagation examples to demonstrate that the model can play a positive and effective role in rumor propagation prediction. The model is optimized by constantly adjusting the parameters through the acquisition of the initial node densities of S, I, C and R in social networking sites at the time of rumor spread, as well as the real-time number of rumors spreading in the early stages. The model can eventually predict the trend of rumor propagation in the future.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data used to support the results of this study were the number of message retweets in social networking platforms.