Coordination of a Dual-Channel Pharmaceutical Supply Chain Based on the Susceptible-Infected-Susceptible Epidemic Model

Abstract

With the continuous development of Internet, online pharmaceutical channels in many countries have seen rapid expansion. As a result, pharmaceutical supply chain participants can adopt dual channels, namely, both online channels and offline channels. As online channels compete with traditional offline channels, it is of great relevance to study the potential conflicts and coordination between them, which is the focus of this paper. Specifically, this article develops a susceptible-infected-susceptible epidemic model of the dual channels for a pharmaceutical supply chain. Our main findings are that in a competitive situation, there is a positive stable equilibrium. Furthermore, increasing the rate of influence of offline transmission, online transmission, and cross transmission will improve sales. Moreover, improving the transmission influence rate will turn more potential customers into purchasers, increase channel sales, and achieve dual channel coordination. We then conduct numerical analysis to illustrate and complement the findings from the model. Finally, we provide managerial insights for implementing successful dual-channel pharmaceutical supply chains.

1. Introduction

With the continuous development of the Internet, online pharmaceutical channels in many countries have expanded rapidly. The United States was one of the first countries to initiate pharmaceutical sales on the Internet. As early as 1991, the United States launched the online pharmacy website certification program and became the pioneering country in online pharmaceuticals. This has led to the current prevalent coexistence of business to business (B2B), business to consumer (B2C), and third-party pharmaceutical Internet trading platforms in the United States. Among them, B2B mode is the most popular one, where suppliers, hospitals, pharmacies, and clinics are integrated. B2C services such as online diagnosis and treatment and online pharmacies also has a large market share. At present, there are more than 1000 online pharmacies and more than 20,000 health-related websites in the United States, and online channels of pharmaceutical sales have become the main growth point of e-commerce in the United States. Furthermore, the UK, Germany, Italy, Switzerland, Japan, and many other countries have a large proportion of online trading. In this paper, we define online channels and offline channels together as dual channels, which means that pharmaceutical supply chain participants can adopt both online channels and offline channels.

Pharmaceutical supply chains are socio-technical systems designed to align firms in order to enable improved health [1]. When a pharmaceutical supply chain involves both online and offline channels, its participants face conflicts in choosing the channel. Channel conflicts can be divided into internal conflicts and external conflicts. Internal conflicts are the contradiction between different departments of the same enterprise. External conflicts refer to the contradiction between two or more independent enterprises. External conflicts can be divided into horizontal conflicts, vertical conflicts, and dual-channel conflicts. Horizontal conflicts are the conflicts between different enterprises at the same level, such as channel conflicts between distributors in different pharmaceutical supply chains. Vertical conflicts are the conflicts between different levels of enterprises in a supply chain, which are more common than horizontal conflicts. For example, the conflicts between distributors and manufacturers, and the conflicts between distributors and retailers. Dual-channel conflicts are the conflicts between different channels, such as online and offline channels. Channel conflicts have both positive and negative effects. Multiple channels may lead to cross-selling, competition for customers and resources, and channel conflicts. As online channels compete with the traditional offline channels, it is therefore of great significance to study the conflicts between online and offline channels.

A number of existing studies focus on dual-channel price coordination problems [2,3,4]. Some studies use supply chain contracts to coordinate the dual-channel [5,6,7], assuming two-stage supply chain distribution systems. For example, reference [8] discussed a manufacturer with two channels, based on the theory of consumer utility. The authors found that the demand of each channel was only affected by prices. Reference [9] considered coordination from three potential perspectives: manufacturer Stackelberg, retailer Stackelberg and Nash equilibrium. The authors showed that price discount could improve the profits of the supply chain members. Reference [10] also employed game theory to study a manufacturer’ dual-channel pricing strategy and showed that a single wholesale price contract could not coordinate channels. However, changing the retail price within a certain range could achieve coordination. References [11] and [12] reviewed the research on pharmaceutical supply chain management and argued that related research on dual channel is rare. Reference [13] pointed out that drug users increasingly purchase drugs online, and purity and adulteration do not vary considerably between drugs purchased online and offline for most substances, while online prices are mostly higher than offline prices. Reference [14] presented a model of multichannel distribution between two separate companies, particularly in the area of retail organizations.

Moreover, scholars have studied the stability of markets and channels borrowing the knowledge from the field of Biology to study dual-channel coordination problem. Reference [15] explored the origins of different channel and investigated the efficiency and evolution of the channel change under the generalized Lotka–Volterra model. Reference [16] used a logistic equation accurately describing the growth law of populations of the channel. The authors also studies population sustainability of a channel [17,18]. Reference [19] studied the self-governing susceptible-infected-susceptible (SIS) epidemic model of two competing channels. In [20], a dual- channel competition model was built and the evolution and stability of consumers’ purchasing behavior in different channels was discussed. Those studies established the foundation for the SIS model for pharmaceutical supply chains and for stability analysis and how to derive the population at different channels.

The rest of this paper is organized as follows. In Section 2, a dual-channel pharmaceutical supply chain structure consisting of an online channel and an offline channel is discussed, and the main mode of each channel is introduced. Section 3 describes the construction of the SIS model of the dual channel for the pharmaceutical supply chain. Section 4 discusses the two-dimensional and four-dimensional equilibrium points and the stability of equilibrium points in the SIS model. Finally, the paper concludes with simulation analyses and discussions.

2. Dual-Channel Pharmaceutical Supply Chain Structure

In the case of pharmaceuticals, the dual-channel system has both positive and negative effects. The offline channel of the pharmaceutical supply chain includes pharmaceutical manufacturers, distributors, retailers, doctors, and patients. The online channel also comprises trading platforms as channel members, including self-operated platforms and third-party platforms. Different countries have different regulations on pharmaceutical distribution. For example, online B2B pharmaceutical transactions is very common among Italian pharmaceutical enterprises, but online B2C transactions are prohibited by regulations. Pharmacies are the main sales channels in the United States, France, Germany, and some other countries, while China sells pharmaceuticals to patients mainly through hospitals.

Table 1. Offline pharmaceutical supply chain structure.

	Offline Channel	Structure
B2B	Agency Model	Pharmaceutical manufacturers sell pharmaceuticals through agents or distributors
	Direct Model	Pharmaceutical manufacturers directly sell pharmaceuticals to medical institutions and pharmacies without agents or distributors
B2C	Hospital Channel	Patients purchase medicine from hospital pharmacy after a doctor prescribes prescription
	Pharmacy Channel	Consumers directly purchase OTC pharmaceuticals, medical devices, health care products, etc., or purchase prescription pharmaceuticals according to prescriptions in the pharmacy

Source: Summarized according to references.

and

Table 2. Online pharmaceutical supply chain structure.

	Online Channel	Structure
B2B	Self-operated	Enterprises build online platforms for pharmaceutical sales to other enterprises
	Third-party platform	E-commerce platform build by third-party for business to business transation
B2C	Self-operated	Enterprises build online platforms for pharmaceutical sales to customers
	Third-party Platform	E-commerce platform build by third-party for business to consumer transation
O2O	Self-operated	Pharmaceutical retailer B2C + Online information release + offline pharmaceutical delivery and payment
	Third-party platform	Online ordering and offline delivery established by third-party platforms

Source: Summarized according to references.

summarize various offline and online pharmaceutical supply chain structures.

In this context, although offline and online channels together expand the scope of consumers and increase market share, they present various advantages and disadvantages to each member of the supply chain. Online channels help pharmaceutical manufacturers obtain more profits, but simultaneously, the profits of distributors decrease, because of direct online transactions between the manufacturers and downstream retailers. As far as pharmaceutical retailers are concerned, on the one hand, online channels increase the profits of enterprises, but on the other, offline channels sales transfer to online channels, thereby damaging the profits of offline channels.

Although online channels cannot provide a face-to-face professional service, with the development of electronic payment and logistics, they facilitate one-stop purchasing for buyers. As professional consultants are also available online, consumers can purchase directly through online channels. In terms of marketing, online channels pay more attention to the maintenance of the website, the choice of online products, the improvement of logistics services, etc. Online channels thus reduce on-site marketing costs. For B2B, the online channel can not only make available information on demand, supply, transfer, and bidding related to pharmaceuticals, but consumers can also compare pharmaceuticals prices, place orders, and so on. Logistics information can be shared among manufacturers, distributors, retailers, and consumers. Manufacturers, distributors, and retailers can also rationally plan their manufacturing and inventory accordingly.

In the traditional offline channels, consumers obtain pharmaceuticals information mainly by the introduction of pharmaceuticals packaging and through pharmacists. The purchase choice is also based on the brand effect or the recommendation of the pharmacist. Contrastingly, for consumers on online B2C channels, a simple click on the web page can provide information on pharmaceuticals, such as brands, prices, and detailed assessments. Consumers on online channels can thus make the most economical and effective purchase choices based on this information. As a result, competition between offline and online channels may lead to price chaos and cost increases, and it is therefore essential to understand the competition between the two channels.

Therefore, in this study, we use a dual-population infectious disease model to simulate the sales process of identical or similar products sold by a pharmaceutical supply chain using offline and online channels simultaneously.

3. The Susceptible-Infected-Susceptible Epidemic Model for Pharmaceutical Dual-Channel Supply Chain

The susceptible-infected-susceptible (SIS) epidemic model, sometimes also called the infectious disease model or the SIS model, is a mathematical model to study the process of transmission of infectious diseases. The infectious disease model has been widely accepted after more than a century of development. It was initially used to study the impacts of diseases [21,22], and was then gradually applied to various other fields such as risk transfer [23], rumor spreading [24], and information spreading [25].

The decision-making process of consumers in the dual channels of the pharmaceutical supply chain is very similar to the transmission of infectious diseases. The susceptible people are those who are likely to be infected: “infected people”, refers to the customers who have already purchased the pharmaceutical’s products through the dual channels; the cured people are those who have recovered through treatment and will be removed from the infected people group. They may acquire immune capacity and no longer buy pharmaceutical products. In this paper, the infectious disease model is used to study consumers’ decision-making behavior and to explore the interaction between the two channels, so as to provide suggestions for the dual-channel development of the pharmaceutical supply chain.

3.1. The Infectious Disease Model in a Dual-Channel Pharmaceutical Supply Chain

3.1.1. The Consumer Purchase Process

(1) Demand formation. There are two kinds of consumers: individual and corporate. Individual consumers purchase pharmaceutical products when they suffer from a disease or have some health-care need; corporate consumers may purchase pharmaceuticals when it is out of stock or inventory falls to the reorder point.

(2) Understanding and comparing pharmaceuticals. To make the most reasonable and effective purchase decision, consumers need to know the function, price, and other information in advance to choose and purchase pharmaceutical products. The purchase criteria can be a measure of the effectiveness of each potential product for the consumer. It can also represent the utility for the customer from purchasing a pharmaceutical candidate from each potential channel.

(3) Deciding on a purchase choice and consumers form purchases. Through pharmaceuticals comparison, consumers choose a pharmaceutical product and buy it from one of the channels. If there is no urgent need for the product or if it does not meet the requirements, the consumers choose not to buy it.

(4) Purchase or repeat purchase. Individual consumers terminate purchasing behavior because their disease is cured or they repeat the purchase because the disease is not cured; corporate consumers stop purchasing because the pharmaceutical sales are lower than expected, or they may re-purchase because of good sales.

For clarity, the customers who form purchase demand, compare pharmaceuticals, or decide on a purchase choice are referred to as potential customers. The consumers who form a purchase behavior are referred to as purchased customers. Finally, the consumers who no longer need to purchase are referred to as cured customers.

3.1.2. Comparison between the Infectious Disease Model and the Dual-Channel Model of Pharmaceutical Supply Chain

As shown in Figure 1, in the two-population infectious disease model, the infectious diseases exist in both populations. There are three different groups in each population: the susceptible, the infected, and the cured. There is a cross-infection phenomenon, which means an infected person in the population 1 can infect a susceptible person in the population 2, or vice versa. In the dual-channel system of pharmaceuticals enterprises, as shown in Figure 2, the online and offline channels are similar to the two populations in the transmission process of infectious diseases. Purchasing the same pharmaceutical is similar to spreading infectious diseases in the two populations. There are three relationships among populations:

• Competition, which refers to the phenomenon that two populations living in the same environment compete for ecological resources;
• Mutually beneficial symbiosis, which refers to the mutual benefit of two populations and when they are separated, neither of them can live independently;
• Parasitism, which refers to the fact that two populations benefit from and suffer due to each other. They become parasitic and host populations, respectively.

In the infectious disease model, a susceptible person is someone who has not been infected but is more likely to be infected. This corresponds to a potential consumer who has not yet purchased pharmaceuticals from the dual channels. An infected person is someone who has been infected with the disease, corresponding to a customer who has already purchased the pharmaceutical in the pharmaceutical supply chain. A cured person is someone who has recovered through treatment. S/he may stop buying because of acquiring immunity, or s/he may buy again because of reinfection.

3.2. Models and Assumptions

When people are infected with a disease, such as meningitis, gonorrhea, etc., they can recover and do not acquire immunity after rehabilitation. The single-channel SIS infectious disease model is shown in Figure 3.

In the single channel SIS model, the population is a closed loop, regardless of birth or death. Hence, $N=S(t)+I(t)$ is a constant. $\gamma I$ indicates the number of recoveries and $\gamma$ indicates recovery rate of infected persons. In the SIS model, the following differential equation is satisfied:

$$\{\begin{array}{l}\frac{dS}{dt}=-\beta SI+\gamma I\\ \frac{dI}{dt}=\beta SI-\gamma I\end{array}$$

(1)

To simplify the model, we make the following assumptions:

• The pharmaceutical products sold in the online and offline channels are similar or identical;
• At a given time, a consumer either has purchased a product or has demand for it;
• At a given time, a consumer only buys pharmaceuticals from one channel;
• At a given time, consumers’ preferences for one buying channel can change to the other; that is, online channel consumers can transform to the offline channel, or vice versa.

As shown in Figure 4, a dual-channel infectious disease SIS model for pharmaceutical supply chain is established. Unlike a single-channel SIS infectious disease model, dual-channel models have mutual influence and transformation in offline and online channels.

3.3. Results

In the dual channels, for rational consumers, the sales of one channel will affect the sales volume of the other channel. Therefore, the two channels are generally in a competitive and substitutable relationship. If an enterprise wants to maximize the total profit or revenue from the two channels, they must coordinate. In this paper, the SIS model is adopted to study the channel stability, and the coordination between the online and offline channels. We have the following two propositions:

Proposition 1.

Only in the competitive situation, there is a stable equilibrium point.

Proposition 2.

The positive impact of the purchased customers on the potential customers increases the sales volume.

3.4. The Dual-Channel Infectious Disease Sis Model

Based on the SIS epidemic model proposed by [19], we construct the competitive dual-channel infectious disease SIS model for the pharmaceutical compound channel, as shown in the Equation (A1).

4. Model Solution

4.1. Equilibrium Point

We solve (Equation (A1)) and predict the purchase situation of customers in the dual channels of pharmaceutical companies based on the equilibrium point. In this way, we can judge the trend and stability of the two channels. Since N_i, S_I, I_i are linearly related, as well as I and I_i, and the main parameter is the amount of customer purchase, the system can be simplified into a four-dimensional system of (N₁, I₁, N₂, I₂), show in the Equation (A2).

4.1.1. Two-Dimensional Equilibrium Point

N₁ and N₂ are relatively independent. Thus, we can form a two-dimensional system for N₁ and N₂, shown in Equation (A3). We can then obtain the equilibrium points $O{\left(0,0\right)}^T$, $P{\left(K_1,0\right)}^T$ and $Q{\left(0,K_2\right)}^T$, which are in area G: $G=\left\{{\left(N_1,N_2\right)}^T|K_i\ge N_i\ge 0,i=1,2\right\}$. Positive equilibrium points $M{\left(N_{1E},N_{2E}\right)}^T$, ${\left(N_{1E},N_{2E}\right)}^T$ are the positive solution of Equation (A4), shown in Appendix A.

We have $N_{1E}=\frac{\left(m+r_1\right)r_2{K}_1}{r_1{r}_2-mn}$, $N_{2E}=\frac{\left(n+r_2\right)r_1{K}_2}{r_1{r}_2-mn}$, and $r_1{r}_2-mn\ne 0$, $\left(m+r_1\right)r_2{K}_1\ne 0$, $\left(n+r_2\right)r_1{K}_2\ne 0$. Since $K_i\ge N_i\ge 0,i=1,2$, we have $K_1\ge N_1\ge 0$, $K_2\ge N_2\ge 0$, then $K_1\ge \frac{\left(m+r_1\right)r_2{K}_1}{r_1{r}_2-mn}\ge 0$, $K_2\ge \frac{\left(n+r_2\right)r_1{K}_2}{r_1{r}_2-mn}\ge 0$, that is $1\ge \frac{{\mathrm{r}}_1{r}_2+m{r}_2}{r_1{r}_2-mn}\ge 0$, $1\ge \frac{{\mathrm{r}}_1{r}_2+n{r}_1}{r_1{r}_2-mn}\ge 0$.

The existence conditions of the positive equilibrium point M are:

• When $r_1{r}_2-mn>0$, that is $r_1{r}_2>mn$, then $m{r}_2\le -mn,n{r}_1\le -mn,\left(m+r_1\right)r_2{K}_1>0$, $\left(n+r_2\right)r_1{K}_2>0$. At this point, $r_1>-m,r_2>-n$; Therefore, m, n are negative, and $r_1>-m,r_2>-n$.
• When $r_1{r}_2-mn<0$, that is $r_1{r}_2<mn$, then $r_1<-m,r_2<-n$, where m and n are negative.

The influence rate between the two channels is negative, indicating that there will be a positive equilibrium point when the channels are in competition. Based on the above analysis, we further obtained the following results:

• When $\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}$ or $\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}$, $m<0$ and $n<0$. The two-dimensional system has four equilibrium points in Area G: $O{\left(0,0\right)}^T$, $P{\left(K_1,0\right)}^T$, $Q{\left(0,K_2\right)}^T$ and $M{\left(N_{1E},N_{2E}\right)}^T$.
• When $\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1<-m\mathrm{or}{r}_2<-n\end{array}$ or $\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1>-m\mathrm{or}{r}_2>-n\end{array}$, the two-dimensional system has three equilibrium points in Area G: $O{\left(0,0\right)}^T$, $P{\left(K_1,0\right)}^T$ and $Q{\left(0,K_2\right)}^T$.

Figure 5 shows the phase diagram of the two-dimensional system of N₁ and N₂:

When $\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}$ or $\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}$, M is the saddle point. $\stackrel{⌢}{OM}$ and $\stackrel{⌢}{MS}$ are the dividing lines and both tend to M, dividing the area into X and Y.

4.1.2. Four-Dimensional Equilibrium Point

The equation shown is Equation (A5). The equilibrium points $P_0{\left(0,0,0,0\right)}^T$, $P_1{\left(K_1,0,0,0\right)}^T$ and $P_3{\left(0,0,K_2,0\right)}^T$ are persistent in area G^’: $G^{’}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T|K_i\ge N_i\ge I_i\ge 0,i=1,2\right\}$. Since the two-dimensional system is included in the four-dimensional one, the equilibrium solution of the two-dimensional system is used to solve the four-dimensional system’s. According to the definition of the threshold ($R_0=\frac{\beta N}{\gamma }$) in the infectious disease model, we let $R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}\mathrm{and}{R}_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}$.

• When N₁ = K₁, N₂ = 0, since $0\le I_2\le N_2$, I₂ = 0. According to Equation (A6), ${\beta }_{11}{I}_1\left(K_1-I_1\right)-{\gamma }_1{I}_1=0$. When I₁ = 0, it coincides with point $P_1{\left(K_1,0,0,0\right)}^T$, then $I_1=\frac{{\beta }_{11}{K}_1-{\gamma }_1}{{\beta }_{11}}$; Since $I_1=\frac{{\beta }_{11}{K}_1-{\gamma }_1}{{\beta }_{11}}>0$, ${\beta }_{11}{K}_1-{\gamma }_1>0$. That is, when $R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}>1$, there is an equilibrium point $P_2{\left(K_1,I_{1E},0,0\right)}^T$, in which $I_{1E}=\frac{{\beta }_{11}{K}_1-{\gamma }_1}{{\beta }_{11}}=K_1\left(1-\frac{{\gamma }_1}{{\beta }_{11}{K}_1}\right)=K_1\left(1-\frac{1}{R_1}\right)$
• When N₁ = 0, N₂ = K₂ and $R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}>0$, there is equilibrium point $P_4{\left(0,0,K_2,I_{2E}\right)}^T$, in which $I_{2E}=\frac{{\beta }_{22}{K}_2-{\gamma }_2}{{\beta }_{22}}=K_1\left(1-\frac{1}{R_2}\right)$
• When $N_1=N_{1E}=\frac{\left(m+r_1\right)r_2{K}_1}{r_1{r}_2-mn}$, $N_2=N_{2E}=\frac{\left(n+r_2\right)r_1{K}_2}{r_1{r}_2-mn}$. According to Equation (A5), both m and n are smaller than 0, and Equation (A6) in the Appendix A is obtained. When I₁ = 0 and I₂ = 0, there is an equilibrium point $P_5{\left(N_{1E},0,N_{2E},0\right)}^T$, which satisfies $\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}$ or $\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}$
• When $I_1{I}_2\ne 0$, there is an equilibrium point $P_6{\left(N_{1E},N_{1EE},N_{2E},N_{2EE}\right)}^T$, which satisfies $\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}$ or $\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}$. We have $A-{\beta }_{11}{N}_{1E}>0$, $B-{\beta }_{22}{N}_{2E}>0$ and ${\prod }_1<{\prod }_2$, where ${\prod }_1=\frac{A-{\beta }_{11}{N}_{1E}}{{\beta }_{12}{N}_{1E}};{\prod }_2=\frac{{\beta }_{21}{N}_{2E}}{B-{\beta }_{22}{N}_{2E}}$; $A={\gamma }_1-r_1\left(1-\frac{N_{1E}}{K_1}\right)-m\frac{N_{2E}}{K_2};B={\gamma }_2-r_2\left(1-\frac{N_{2E}}{K_2}\right)-n\frac{N_{1E}}{K_1}$.

Thus, the following equilibrium points are obtained in different regions: $P_0{\left(0,0,0,0\right)}^T$, $P_1{\left(K_1,0,0,0\right)}^T$, $P_2{\left(K_1,I_{1E},0,0\right)}^T$, $P_3{\left(0,0,K_2,0\right)}^T$, $P_4{\left(0,0,K_2,I_{2E}\right)}^T$, $P_5{\left(N_{1E},0,N_{2E},0\right)}^T$ and $P_6{\left(N_{1E},N_{1EE},N_{2E},N_{2EE}\right)}^T$. The stability of each equilibrium point in Region $G^{’}$ is discussed below.

4.2. Stability Analysis

4.2.1. Stability of Equilibrium Points

Following [20], we use the Lyapunov Method and the Jacobi Matrix to explore the existence conditions and stability of the equilibrium point of the model.

Table 3. Stability of equilibrium point of the four-dimensional system.

Existence Condition			Stability of Equilibrium Points
			$P_0$ is unstable in ${\mathsf{\Omega }}_1$
$\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}$	$R_1<1$		$P_1$ is asymptotically stable in ${\mathsf{\Omega }}_2$
	$R_1>1$		$P_2$ exists, asymptotically stable in ${\mathsf{\Omega }}_3$
	$R_2<1$		$P_3$ is asymptotically stable in ${\mathsf{\Omega }}_6$
	$R_2>1$		$P_4$ exists, asymptotically stable in ${\mathsf{\Omega }}_7$
	$A-{\beta }_{11}{N}_{1E}>0$$B-{\beta }_{22}{N}_{2E}>0$	${\prod }_1>{\prod }_2$	$P_5$ exists, asymptotically stable in ${\mathsf{\Omega }}_{12}$
		${\prod }_1<{\prod }_2$	$P_6$ exists, asymptotically stable in ${\mathsf{\Omega }}_{13}$
$\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}$	$A-{\beta }_{11}{N}_{1E}>0$$B-{\beta }_{22}{N}_{2E}>0$	${\prod }_1>{\prod }_2$	$P_5$ exists, asymptotically stable in ${\mathsf{\Omega }}_{10}$
		${\prod }_1<{\prod }_2$	$P_6$ exists, asymptotically stable in ${\mathsf{\Omega }}_{11}$
$\{\begin{array}{l}{r}_1{r}_2>mn\\ r_2<-n\end{array}\mathrm{or}\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1>-m\end{array}$	$R_1<1$		$P_1$ is asymptotically stable in ${\mathsf{\Omega }}_4$
	$R_1>1$		$P_2$ exists, asymptotically stable in ${\mathsf{\Omega }}_5$
$\{\begin{array}{l}{r}_1{r}_2<mn\\ r_2>-n\end{array}\mathrm{or}\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1<-m\end{array}$	$R_2<1$		$P_3$ is asymptotically stable in ${\mathsf{\Omega }}_8$
	$R_2>1$		$P_4$ exists, asymptotically stable in ${\mathsf{\Omega }}_9$

lists the stability of equilibrium points under different conditions.

The meaning of the symbols in the table:

${\mathsf{\Omega }}_1=G^{’}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T|K_i\ge N_i\ge I_i\ge 0,i=1,2\right\}$	${\mathsf{\Omega }}_2=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in Y\right\}$
${\mathsf{\Omega }}_3=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in Y,I_1\ne 0\right\}$	${\mathsf{\Omega }}_4=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|N_1\ne 0\right\}$
${\mathsf{\Omega }}_5=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|N_1\ne 0,I_1\ne 0\right\}$	${\mathsf{\Omega }}_6=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in X\right\}$
${\mathsf{\Omega }}_7=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in X,I_2\ne 0\right\}$	${\mathsf{\Omega }}_8=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|N_2\ne 0\right\}$
${\mathsf{\Omega }}_9=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|N_2\ne 0,I_2\ne 0\right\}$	${\mathsf{\Omega }}_{10}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|N_1\ne 0,N_2\ne 0\right\}$
${\mathsf{\Omega }}_{11}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|\begin{array}{l}{N}_1\ne 0,I_1\ne 0,\\ N_2\ne 0,I_2\ne 0\end{array}\right\}$	${\mathsf{\Omega }}_{12}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in \stackrel{⌢}{OM}\cup \stackrel{⌢}{MS}\right\}$
${\mathsf{\Omega }}_{13}=\left\{{\left(N_1,I_1,N_2,I_2\right)}^T\in G^{’}|{\left(N_1,N_2\right)}^T\in \stackrel{⌢}{OM}\cup \stackrel{⌢}{MS},I_1\ne 0,I_2\ne 0\right\}.$

X and Y are regions composed of N₁ and N₂ in the two-dimensional system.

4.2.2. Stability of the Dual-Channel Pharmaceutical Supply Chain

Table 4. Stability of the dual-channel pharmaceutical supply chain.

Existence Condition		Offline	Online		Total Sales
		Demand	Sales	Demand	Sales
$\mathbf{No}\mathbf{Constraint}\mathbf{in}\mathbf{Area}{\mathit{G}}^{\prime }$		Does Not Approach 0	Do Not Approach 0	Does Not Approach 0	Do Not Approach 0	Do Not Approach 0
$\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1<-m\\ r_2<-n\end{array}m,n<0$	$R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}<1$; ${\beta }_{11}$ smaller, ${\gamma }_1$ larger; in region ${\mathsf{\Omega }}_2$	K1a	0	0	0	0
	$R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}>1$; ${\beta }_{11}$ larger, ${\gamma }_1$ smaller; in region ${\mathsf{\Omega }}_3$	K1	I1E	0	0	I1E
	$R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}<1$; ${\beta }_{22}$ smaller, ${\gamma }_2$ larger; in region ${\mathsf{\Omega }}_6$	0	0	K2	0	0
	$R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}>1$; ${\beta }_{22}$ larger, ${\gamma }_2$ smaller; in region ${\mathsf{\Omega }}_7$	0	0	K2	I2E	I2E
	${\beta }_{11}$, ${\beta }_{12}$, ${\beta }_{21}$, ${\beta }_{22}$ smaller; in region ${\mathsf{\Omega }}_{12}$	N1E	0	N2E	0	0
	${\beta }_{11}$, ${\beta }_{12}$, ${\beta }_{21}$, ${\beta }_{22}$ larger; in region ${\mathsf{\Omega }}_{13}$	N1E	I1EE	N2E	I2EE	I1EE + I2EE
$\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1>-m\\ r_2>-n\end{array}m,n<0$	${\beta }_{11}$, ${\beta }_{12}$, ${\beta }_{21}$, ${\beta }_{22}$ smaller; in region ${\mathsf{\Omega }}_{10}$	N1E	0	N2E	0	0
	${\beta }_{11}$, ${\beta }_{12}$, ${\beta }_{21}$, ${\beta }_{22}$ larger; in region ${\mathsf{\Omega }}_{11}$	N1E	I1EE	N2E	I2EE	I1EE + I2EE
$\{\begin{array}{l}{r}_1{r}_2>mn\\ r_2<-n\end{array}\mathrm{or}\{\begin{array}{l}{r}_1{r}_2<mn\\ r_1>-m\end{array}$(Unable to judge the positive and negative of m and n)	$R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}<1$; ${\beta }_{11}$ smaller, ${\gamma }_1\mathrm{larger};\mathrm{in}\mathrm{region}$	K1	0	0	0	0
	$R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}>1$; ${\beta }_{11}$ smaller, ${\gamma }_1$ larger; in region ${\mathsf{\Omega }}_5$	K1	I1E	0	0	I1E
$\{\begin{array}{l}{r}_1{r}_2<mn\\ r_2>-n\end{array}\mathrm{or}\{\begin{array}{l}{r}_1{r}_2>mn\\ r_1<-m\end{array}$(Unable to judge the positive and negative of m and n)	$R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}<1$; ${\beta }_{22}$ smaller, ${\gamma }_2$ larger; in region ${\mathsf{\Omega }}_8$	0	0	K2	0	0
	$R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}>1$; ${\beta }_{22}$ smaller, ${\gamma }_2$ larger; in region ${\mathsf{\Omega }}_9$	0	0	K2	I2E	I2E

shows that when $R_1>1$ (${\beta }_{11}$ is larger, ${\gamma }_1$ is smaller), offline sales remain at a stable level, while online sales are poor, and the total sales of the company tended to be offline sales. When $R_2>1$ (${\beta }_{22}$ is larger, ${\gamma }_2$ is smaller), then online sales remain stable and offline sales are poor, and the total sales tend to be online sales. When $B-{\beta }_{22}{N}_{2E}>0$, ${\prod }_1=\frac{A-{\beta }_{11}{N}_{1E}}{{\beta }_{12}{N}_{1E}}<{\prod }_2=\frac{{\beta }_{21}{N}_{2E}}{B-{\beta }_{22}{N}_{2E}}$, ${\beta }_{11}$, ${\beta }_{12}$, ${\beta }_{21}$, ${\beta }_{22}$ are larger, both online and offline sales are at a stable level, and total dual-channel sales tend to be stable.

5. Simulation Analysis

This paper mainly discusses the change of the total sales volume of pharmaceuticals in the dual channels, and explores the possible coordination. When $P_0$, $P_1$, $P_3$, $P_5$ are in an unstable condition, the sales volume tends to zero, a situation that companies would wish to avoid. Simulation analysis is used to analyze the coordination of the dual channel pharmaceutical supply chain.. Due to the instability of $P_0$, it will not be discussed further. A positive balance only occurs when there is competition between the two channels.

5.1. Analyses under a Stable $P_1{\left(K_1,0,0,0\right)}^T$

$R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}<1$ exists under the two circumstances of $P_1{\left(K_1,0,0,0\right)}^T$. When it meets the requirement that $R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}>1$, $P_1$ can be converted to $P_2$. This means if the ${\beta }_{11}$ increases and ${\gamma }_1$ decreases, the offline sales volume improve and the total sales volume of the enterprise improves as well. Assume that ${\gamma }_1={\gamma }_2=1$, $K_1=K_2=1000$, $N_1\left(0\right)=N_2\left(0\right)=100,m=-1.2,n=-1.5,t=300,$ ${\gamma }_1={\gamma }_2=0.01$, ${\beta }_{12}=0.05,{\beta }_{22}=0.05,{\beta }_{21}=0.05,$ ${\beta }_{11}=0.000005;0.00003;0.00006;0.00009,$ it meets the requirement that $R_1=\frac{{\beta }_{11}{K}_1}{{\gamma }_1}=0.5<1$.

As shown in Figure 6, it presents a change in the total sales volume of dual-channel pharmaceutical companies at time t. The total sales volume of the dual channel increases with the increase of the offline transmission influence rate ${\beta }_{11}$, achieving the coordination of the two channels.

5.2. Analyses under a Stabile of $P_3{\left(0,0,K_2,0\right)}^T$

When it meets the requirement that $R_2=\frac{{\beta }_{22}{K}_2}{{\gamma }_2}>1$, $P_3{\left(0,0,K_2,0\right)}^T$ can be converted to $P_4{\left(0,0,K_2,I_{2E}\right)}^T$. This means that if the ${\beta }_{22}$ increases and ${\gamma }_2$ decreases, online sales volumes improve, as do the total sales volumes of the enterprise. Assume that ${\gamma }_1={\gamma }_2=1$, $K_1=K_2=1000$, $N_1\left(0\right)=N_2\left(0\right)=100,m=-1.2,n=-0.5,t=10,$ ${\gamma }_1={\gamma }_2=0.01$, ${\beta }_{11}=0.00005$, ${\beta }_{12}=0.00005$ ${\beta }_{22}=0.000005;0.00003;0.00006;0.00009$ and ${\beta }_{21}=0.00005,$

As shown in Figure 7, the total sales volume of the dual-channel increases with the increase of the online transmission influence rate ${\beta }_{22}$, achieving the coordination of the two channels.

5.3. Analyses under a Stabile $P_5{\left(N_{1E},0,N_{2E},0\right)}^T$

When it meets the requirement ${\prod }_1=\frac{A-{\beta }_{11}{N}_{1E}}{{\beta }_{12}{N}_{1E}}<{\prod }_2=\frac{{\beta }_{21}{N}_{2E}}{B-{\beta }_{22}{N}_{2E}}$, $P_5{\left(N_{1E},0,N_{2E},0\right)}^T$ can be converted to $P_6{\left(N_{1E},N_{1EE},N_{2E},N_{2EE}\right)}^T$. This means that if the ${\beta }_{12}$ and ${\beta }_{21}$ increase, both offline and online sales volumes improve and the total sales volume of the enterprise improves as well. Assume that ${\gamma }_1={\gamma }_2=1$, $K_1=K_2=1000$, $N_1\left(0\right)=N_2\left(0\right)=100,m=-0.8,n=-0.5,t=10,$ ${\gamma }_1={\gamma }_2=0.01$, ${\beta }_{11}=0.00002$, ${\beta }_{22}=0.00001$, ${\beta }_{12}=0.000005;0.00001;0.00002;0.00005$, and ${\beta }_{22}=0.000003;0.00001;0.00002;0.00005$,

As shown in Figure 8, the total sales volume of the dual channels increases with the increase of the offline and online transmission influence rate ${\beta }_{12}$ and ${\beta }_{21}$. Thus, dual channel coordination is achieved.

6. Conclusions

This paper first discussed the offline and online channels of a pharmaceutical supply chain. It then constructed the susceptible-infected-susceptible epidemic model for the dual-channel supply chain. Finally, it derived the equilibrium points and conducted simulation analyses. Our results show that the increase in offline transmission influence rate ${\beta }_{11}$, online transmission influence rate ${\beta }_{22}$, offline transmission influence rate ${\beta }_{12}$, and offline transmission influence rate ${\beta }_{21}$ can improve the total sales in the pharmaceutical supply chain. Based on our results, the following managerial insights are obtained:

• Offline channels can reduce prices, enhance service quality, and attract consumers who prefer offline channels to become buyers. Enterprises can use revenue-sharing supply chain contracts and other methods to obtain lower wholesale prices and lower pharmaceutical prices, thereby achieving sales growth. They can also provide comprehensive pre-sales and after-sales services.
• Online channels can provide online consultation and diversified distribution services. Professionals can provide consulting services for online consumers. For customers with time constraints, pharmaceutical enterprises can cooperate with express logistics enterprises. For more price-sensitive customers, as long as the reliability of shipping is guaranteed, they would choose cheaper modes of transportation.
• Online and offline channels should differentiate between their target customers. The online channel should target young customers, who have a higher acceptance of new things. In contrast, offline channels should target older people as customers as they are more willing to visit an offline store. Online and offline channels should also differentiate brands and packaging, and control the price gap within a reasonable range.
• Dual-channel pharmaceutical supply chains need to achieve the mutual coordination of online and offline flow and information sharing. Through remote diagnosis and treatment and O2O in certain regions, we can improve the information sharing of offline and online channels, including marketing strategies, market information, information distribution, and so on. In this way, dual-channel pharmaceutical supply chains can realize the coordination of marketing and logistics.

There are several limitations in this research. Especially, the potential factors in the dual-channel SIS model are not considered comprehensively. There are still many potential problems to be further explored. First, future research can conduct more detailed quantitative analysis on the factor of the influence rate. Second, actual data can be collected based on which our models and results can be tested empirically. Third, issues like network design, channel selection, and fulfillment service contracts can be addressed for dual-channel pharmaceutical supply chains [26,27].