# Modeling Medical Services with Mobile Health Applications

Authors: Zhenfei Wang; Liying Zhang; Ling Ma; Bing Liu
Journal: Journal of Healthcare Engineering
 (2018)

---

## Abstract

The rapid development of mobile health technology (m-Health) provides unprecedented opportunities for improving health services. As the bridge between doctors and patients, mobile health applications enable patients to communicate with doctors through their smartphones, which is becoming more and more popular among people. To evaluate the influence of m-Health applications on the medical service market, we propose a medical service equilibrium model. The model can balance the supply of doctors and demand of patients and reflect possible options for both doctors and patients with or without m-Health applications in the medical service market. In the meantime, we analyze the behavior of patients and the activities of doctors to minimize patients’ full costs of healthcare and doctors’ futility. Then, we provide a resolution algorithm through mathematical reasoning. Lastly, based on artificially generated dataset, experiments are conducted to evaluate the medical services of m-Health applications.


---

## Body

## 1. Introduction

Healthcare performs an important role in people’s health and life quality. Healthcare spending in China is increasing rapidly. The country’s annual expenditure is projected to grow at an average rate of 11.8 percent a year from 2014 to 2018, reaching $892 billion by 2018 [1]. Despite these positive numbers, people still have to pay expensive costs when they go to hospitals for medical services. Besides, the people, especially for rural and remote area people who need medical services, have no alternative choice but travel to other cities to receive healthcare services and consultation. What is more, before consulting doctors, people may have to wait for a long time for diagnostic examinations. All of those will increase their commuting and waiting costs. In addition, many developing countries have a large population, while the medical staff are insufficient. For example, China only has 1.94 physicians and 0.83 health workers for every 1000 patients, far behind the Millennium Development Goal that every 1000 people need at least 2.5 medical staff for healthcare. [2] Because of the limited capacity and resources of the government and Information and Communication Technologies (ICT), the medical service system in China is unsound.Meanwhile, technologies have been advanced to improve the medical services, such as telemedicine [3] and mobile health [4]. Traditionally, patients always go to the doctors in hospital, which is not very convenient. This way can increase healthcare costs. Due to the advance of mobile and wireless communication technologies, people can receive medical services by their smartphones rather than having to go to hospitals. There are a total of 5 billion mobile phones in the world, of which more than 1.08 billion are smartphones. So, each of the 80% of the world population owns a smartphone [5]. The number of smartphones will soar to 16.8 billion globally by 2019 [6]. With the increase of mobile devices, various mobile health (m-Health) applications, such as MD, iTriage, UK’sNHS, Mayo Clinic’s Health App, and Spring Rain Doctor, have emerged in recent years. For example, Spring Rain Doctor, one main m-Health application in China, has reached 92 million users [7]. Spring Rain Doctor is the largest mobile doctor-patient communication platform in the world. In the medical service market, patients have the need for recovery. They can consult doctors and get services, expecting to spend as little as possible. As service providers, doctors will naturally want to get more pay. Taking Spring Rain Doctor as an example, for patients, they can get quick consultation services wherever and whenever they need. They could know more about the disease and treatment measures as well. It can help patients reduce the costs of time, space, and money and prevent excessive medical treatment. For doctors, Spring Rain Doctor can help doctors fully utilize their fragmented time. The doctors can increase their income by the convenient Internet communication and establish personal brand. In the meantime, doctors can decrease misdiagnosis rate by adding big data system auxiliary outside of doctor-patient multidirectional interaction. As the bridge between doctors and patients, m-Health applications will have a major effect on people’s health.The term m-Health was coined by Robert Istepanian as use of “emerging mobile communications and network technologies for healthcare” in 2005 [8]. m-Health is a subarea of eHealth. People can use mobile phones or mobile telecom equipment to receive health services and information. But m-Health cannot replace e-Health [9]. m-Health can help patients to receive health services using information and communication technologies. This convenient and efficient way can reduce treatment costs. The applications (apps) of m-Health can help people to change their unhealthy lifestyle and form good habits [10]. Steinhubl et al. [11] provided a comprehensive review of the emerging field of mobile health and the challenges we faced. Recently, many smartphone-based m-Health applications have been developed for medical purposes. Sclafani et al. [12] conducted an investigation of mobile tablet and applications by email. Nearly half of tablet users said that they would use mobile tablet and applications to get health services. Based on fundamental features of m-Health for physical activity, Ayubi et al. [13] developed an m-Health application successfully. According to relevant literatures of m-Health platform’s benefits and standards, Whiteside [14] examined some m-Health applications. In order to deal with alert messages, Kafeza et al. [15] developed a model for mobile devices to receive medical task alerts. In order to investigate how to interconnect biomedical sensors in a wireless body area network, Bao et al. [16] proposed a novel solution to tackle the problem of entity authentication in body area sensor network (BASN) for m-Health. To overcome the scarcity of follow-up health services for cancer, Tiancheng et al. [17] designed a smartphone app framework. The framework can make risk assessment for the patient’s health condition and supervise the follow-up health services. After analyzing the usage characteristics of a mobile health application, Garibay et al. [18] found that the app named Heartkeeper was useful to improve the quality of medical services for patients. After exploring the benefits and challenges of health mobile apps for patients with type 2 diabetes, Peng et al. [19] provided the acceptability, feasibility, and effectiveness of mobile apps for patients’ compliance behavior, self-management, and self-care ability. After describing the use and acceptability of m-Health by patients with HIV/tuberculosis and healthcare providers, Hirsch-Moverman et al. [20] found that the m-Health intervention was a low-tech, user-friendly intervention. To promote the development of m-Health, Jones et al. [21] proposed a model-driven design and development methodology. From the above studies, we can see that all of the studies have illuminated the roles of m-Health applications including medical service reminders, arranger, recorder, transmitter, mentor, and intervener. These m-Health applications can improve patients’ life quality. However, none of the studies have estimated the effect of communication of m-Health application and the impact on patients and doctors with or without m-Health applications.In fact, in many cities such as New York City and most cities in China, m-Health applications are widely adopted by doctors and patients. For example, a mobile health app is named My Medical™. Patients can store and send their medical records to doctors by pressing the button on it. And doctors also can give replies to patients by the app. Besides, doctors can communicate with patients by another app named Doximity, send secure HIPAA faxes, and follow information about their specialty [22]. In China, Spring Rain Doctor [23] allows patients to communicate with doctors about their health issues directly by their telephones.To estimate the impact of m-Health applications on these medical service markets, this paper proposes an equilibrium model. The model can balance the supply of doctors and demand of patients and reflect possible options of both doctors and patients with or without m-Health applications in the medical service market. Both doctors and patients have the right to use or not use the m-Health application and also have the right to accept or not accept each other. The earlier work by He and Shen [24] is adopted to describe the possible adoptions of both doctors and patients of the m-Health applications. The purpose of the patients is to minimize their own full costs of healthcare including the consultation fee and waiting and medical service time costs. Doctors can get as much pay as possible by minimizing their own futility.The main contributions of this paper are as follows:(i)
We redesign the futility function. The futility function in this paper mainly considers the waiting time of doctor, consulting income, and attractiveness of patient’s disease information to doctor.
(ii)
In this paper, both patient’s disease information and doctor registration information is considered in the matching process between doctors and patients, which is based on the need of patients and doctors.
(iii)
We construct a mobile health medical service market network based on a hypothetical distribution platform in the experiment. The results show that the model could balance the supply of doctors and demand of patients.The paper is structured as follows. The next section describes the problem and analyzes the behavior of patients and the activities of doctors to minimize patients’ full costs of healthcare and doctors’ futility. The following section defines the model and provides a resolution algorithm by mathematical reasoning. Numerical examples are presented in Section4 to demonstrate characteristics of the medical service system with the m-Health application. Lastly, Section 5 concludes the paper.


## 2. Characterization of Doctor-Patient Assignments

In this section, we first present the problem formulation and basic consideration of our problem. Then, we characterize our problem from two perspectives of both patients and doctors based on the earlier work of He and Shen [24].


### 2.1. Problem Formulation and Basic Consideration

#### 2.1.1. Problem Formulation

Given a mobile health medical service market network which consists of both doctors and patients, our research problem is how to assign doctors and patients, meet the supply of doctors and demand of patients within the network, and meanwhile minimize patients’ own full costs of healthcare and idle doctors’ own futility.We suppose that the m-Health application is widely available to doctors and patients. In the rest of this article, mode 1 represents the medical services that patients and doctors communicate with each other in the hospital. Mode 2 denotes the medical services that patients and doctors communicate with each other by utilizing the m-Health application. Both doctors and patients can choose each other through mode 1 and mode 2. In mode 2, a patient can send a request to a doctor through the m-Health application. A doctor can or cannot accept the request according to his/her spare time. If a doctor accepts the request, a consultation connecting process will be completed. After that, the doctor should communicate with the patient and give medical advice and then charge the patient for the consultation. The waiting time of doctors and patients is common in the consultation process. To simplify research, we only consider the common costs of the two models. In this paper, the patients’ full costs of healthcare consist of consultation fee and waiting and medical service time costs.


#### 2.1.2. Basic Definition and Setting

In the network, we need to make some assumptions and settings. We consider that medical service supply and patient demand are stationary in one unit period (i.e., one hour).(i)
I is the set of locations of doctors.
(ii)
J is the set of locations of patients.
(iii)
ϕ represents the set of D–P (doctor-patient) pairs.
(iv)
Pij is the number of patients’ demand between D–P pairs i,j∈φ in one unit time.
(v)
tij is the average medical service time between D–P pairs i,j∈φ in one unit time.
(vi)
Fij is the consultation fee between D–P pairs i,j∈φ of the two modes.The last three parameters are fixed and given.


### 2.2. Behavior of Patients

For the way to see a doctor, we give some notations which will be used in the mathematical model.(i)
Qijp1 is the number of patient demand between D–P pairs i,j using mode 1.
(ii)
Qijp2 is the number of patient demand between D–P pairs i,j utilizing mode 2, “p” stands for patient, and Qijp1+Qijp2=Pij.
(iii)
Wjp1 is the average waiting time of patients j using mode 1.
(iv)
Wijp2 is the average waiting time of patients consulting through mode 2.
(v)
Cijp1 is the full costs of healthcare for patients’ medical services using mode 1.
(vi)
Cijp2 is the full costs of healthcare for patients’ medical services through mode 2.Considering the above costs of two service modes, patients will choose the way of medical services to minimize their own healthcare costs. Therefore, we can get the patients’ healthcare costs of two modes. As shown in(1)Cijp1=Fij+μ1Wjp1+μ2tij,Cijp2=Fij+μ1Wijp2+μ2tij+μ1Z¯,where μ1 and μ2 are the cost values of patients’ waiting time and medical service time; Z¯ denotes the patients’ waiting time due to doctor’s departure during consultation process of mode 2.Therefore, the probabilities of a patient consulting a doctor through modes 1 and 2 are as follows:(1)
P1: the probability of a patient consulting a doctor through mode 1 is given by the logit model [24], as shown in
(2)P1=Qijp1Qijp1+Qijp2=exp−ω1Cijp1exp−ω1Cijp1+exp−ω1Cijp2∀i,j∈φ.
(2)
P2: the probability of a patient consulting a doctor through mode 2 is given by the logit model, as shown in
(3)P2=Qijp2Qijp1+Qijp2=exp−ω1Cijp2exp−ω1Cijp1+exp−ω1Cijp2∀i,j∈ϕ,where ω1 is a nonnegative parameter and reflects the degree of uncertainty about the patients’ demand.
From (1), (2), and (3), we get
(4)Qijp1=Pij⋅P1∀i,j∈φ,Qijp2=Pij⋅P2∀i,j∈φ.


### 2.3. Activities of Doctors

Accordingly, we also define some variables for doctors’ activities in the network.(i)
N is the doctors’ activities time, including occupied doctors’ time and idle doctors’ time,
(ii)
Tijod is the number of occupied doctor (number of persons/h) between D–P pairs i,j and “od” represents occupied doctors.
(iii)
∑ij∈φTijodtij is the total medical service time of all occupied doctors per unit time.
(iv)
Tijid2 is the number of idle doctors utilizing mode 2 to receive patients’ requests, and “id” stands for idle doctors.
(v)
Wijid2 is the average waiting time of doctors receiving patients utilizing mode 2.
(vi)
Tiid1 is the number of idle doctors utilizing mode 1 to receive patients’ requests.
(vii)
Wiid1 is the average waiting time of doctors receiving patients’ requests through mode 1.
(viii)
VT is the total unoccupied doctors’ service time in one unit time, VT=∑ij∈φTiid1Wiid1+∑ij∈φTijid2Wijid2+Z¯.
(ix)
Uijid1 is the futility of idle doctors accepting patients using mode 1.
(x)
Uijid2 is the futility of idle doctors accepting patients through mode 2.
(xi)
T~iid1 is the total number of idle doctors that receive the next patient using mode 1.
(xii)
T~iid2 is the number of idle doctors seeking the next patient through mode 2.Considering the total activity time of doctors, we can obtain(5)N=∑ij∈φTijodtij+∑ij∈φTiid1Wiid1+∑ij∈φTijid2Wijid2+Z¯.After completing one consulting activity, doctors will determine whether to receive the next patient to minimize their futility. The futility functions for idle doctors [25] are as follows:
(6)Uijid1=−Fi¯+πtWiid1+ti¯−yj¯∀i,j∈φ,Uijid2=−Fij+πtWijid2+tij+Z¯−yj∀i,j∈φ,where πt represents the doctor activity cost in one unit time; Fi¯ and ti¯ denote the average consulting income and occupied waiting time of mode 1, with Fi¯=∑i,j∈φFijQijp1/∑i,j∈φQijp1 and ti¯=∑i,j∈φtijQijp1/∑i,j∈φQijp1; yj is the attractiveness of patient j’s disease information to idle doctors that receive the next patient through mode 2; and yj¯ represents the average attractiveness of patient j’s disease information to idle doctors that receive the next patient through mode 1, with yj¯=∑i,j∈φyjQijp1/∑i,j∈φQijp1.After analyzing the doctors’ activity behavior and the futility functions, we can know that the idle doctors can accept the next patient in the hospital (mode 1) or through m-Health applications (mode 2). And we can get(7)T~iid1=∑i∈ITiid1,(8)T~iid2=∑i,j∈φTijid2,(9)∑i,j∈φTijod=T~iid1+T~iid2.Formula (9) dictates that for a given hour, the number of doctors completing medical services should equal the number of idle doctors.Then, we can obtain the probabilities of an idle doctor receiving the next patient through modes 1 and 2.(3)
P3: the probability of an idle doctor seeking the next patient through mode 1 is given by the logit model
(10)P3=T~iid1T~iid1+T~iid2=exp−ω2Liid1exp−ω2Liid1+exp−ω2Liid2.
(4)
P4: the probability of an idle doctor seeking the next patient through mode 2 is given by the logit model
(11)P4=T~iid2T~iid1+T~iid2=exp−ω2Liid2exp−ω2Liid1+exp−ω2Liid2,where ω2is the dispersion coefficient of the logit choice model.The probabilities of choosing patientj as the next patient by an idle doctor through modes 1 and 2 are as follows:
(5)
P5: the probability of choosing patient j as the next patient by an idle doctor through mode 1 is shown in
(12)P5=Tijid1T~iid1=exp−ω3Uijid1∑i,j∈φexp−ω3Uijid1.
(6)
P6: the probability of choosing patient j as the next patient by an idle doctor through mode 2 is shown in
(13)P6=Tijid2T~iid2=exp−ω4Uijid2∑i,j∈φexp−ω4Uijid2,where ω3 and ω4 are the dispersion coefficients of the logit choice model.From above formulas, we can get the log sums of the futility of idle doctors seeking their next patient through modes 1 and 2, that areLiid1=−1/ω3ln∑iexp−ω3Uijid1 and Liid2=−1/ω4ln∑i,j∈φexp−ω4Uijid2 and also ω3≫ω2 and ω4≫ω2 [26].In a stationary equilibrium state, no matter how a patient sees the doctor, the patient will receive medical services after waiting a certain time. So(14)∑i∈ITiid1=∑i,j∈φQijp1,(15)∑i∈ITijid2=∑i,j∈φQijp2.In the network, medical service process is a mutual selection process in which patients can search for doctors and doctors also can search for patients.Based on the Cobb–Douglas-type production function [25], we can get the waiting time function of mode 1.If0<α1,α2<1, we get
(16)Wjp1=Aj−1/α1∑i,j∈φQijp11−α1/α1Wiid1∑i∈ITiid1−α2/α1,where α1 is the meeting rate elasticity of the number of unserved patients, α2 is the meeting rate elasticity of the number of idle doctors, Aj is a zone-specific parameter, Wiid1∑i∈ITiid1 denotes the total waiting hours of the idle doctor using mode 1 of zone i, and ∑i,j∈ϕQijp1 is the unserved patient number of zone i.Ifα1=α2=1, we obtain Wjp1=1/AjWiid1∑i∈ITiid1.Similarly, we can also get the waiting time function of mode 2.Wijp2 is shown as
(17)Wijp2=A~−1/α1∑i,j∈ϕQijp21−α1/α1Wijid2∑i,j∈ϕTijid2−α2/α1,where A~ is a constant number to reflect the little-friction property.


## 2.1. Problem Formulation and Basic Consideration

### 2.1.1. Problem Formulation

Given a mobile health medical service market network which consists of both doctors and patients, our research problem is how to assign doctors and patients, meet the supply of doctors and demand of patients within the network, and meanwhile minimize patients’ own full costs of healthcare and idle doctors’ own futility.We suppose that the m-Health application is widely available to doctors and patients. In the rest of this article, mode 1 represents the medical services that patients and doctors communicate with each other in the hospital. Mode 2 denotes the medical services that patients and doctors communicate with each other by utilizing the m-Health application. Both doctors and patients can choose each other through mode 1 and mode 2. In mode 2, a patient can send a request to a doctor through the m-Health application. A doctor can or cannot accept the request according to his/her spare time. If a doctor accepts the request, a consultation connecting process will be completed. After that, the doctor should communicate with the patient and give medical advice and then charge the patient for the consultation. The waiting time of doctors and patients is common in the consultation process. To simplify research, we only consider the common costs of the two models. In this paper, the patients’ full costs of healthcare consist of consultation fee and waiting and medical service time costs.


### 2.1.2. Basic Definition and Setting

In the network, we need to make some assumptions and settings. We consider that medical service supply and patient demand are stationary in one unit period (i.e., one hour).(i)
I is the set of locations of doctors.
(ii)
J is the set of locations of patients.
(iii)
ϕ represents the set of D–P (doctor-patient) pairs.
(iv)
Pij is the number of patients’ demand between D–P pairs i,j∈φ in one unit time.
(v)
tij is the average medical service time between D–P pairs i,j∈φ in one unit time.
(vi)
Fij is the consultation fee between D–P pairs i,j∈φ of the two modes.The last three parameters are fixed and given.


## 2.1.1. Problem Formulation

Given a mobile health medical service market network which consists of both doctors and patients, our research problem is how to assign doctors and patients, meet the supply of doctors and demand of patients within the network, and meanwhile minimize patients’ own full costs of healthcare and idle doctors’ own futility.We suppose that the m-Health application is widely available to doctors and patients. In the rest of this article, mode 1 represents the medical services that patients and doctors communicate with each other in the hospital. Mode 2 denotes the medical services that patients and doctors communicate with each other by utilizing the m-Health application. Both doctors and patients can choose each other through mode 1 and mode 2. In mode 2, a patient can send a request to a doctor through the m-Health application. A doctor can or cannot accept the request according to his/her spare time. If a doctor accepts the request, a consultation connecting process will be completed. After that, the doctor should communicate with the patient and give medical advice and then charge the patient for the consultation. The waiting time of doctors and patients is common in the consultation process. To simplify research, we only consider the common costs of the two models. In this paper, the patients’ full costs of healthcare consist of consultation fee and waiting and medical service time costs.


## 2.1.2. Basic Definition and Setting

In the network, we need to make some assumptions and settings. We consider that medical service supply and patient demand are stationary in one unit period (i.e., one hour).(i)
I is the set of locations of doctors.
(ii)
J is the set of locations of patients.
(iii)
ϕ represents the set of D–P (doctor-patient) pairs.
(iv)
Pij is the number of patients’ demand between D–P pairs i,j∈φ in one unit time.
(v)
tij is the average medical service time between D–P pairs i,j∈φ in one unit time.
(vi)
Fij is the consultation fee between D–P pairs i,j∈φ of the two modes.The last three parameters are fixed and given.


## 2.2. Behavior of Patients

For the way to see a doctor, we give some notations which will be used in the mathematical model.(i)
Qijp1 is the number of patient demand between D–P pairs i,j using mode 1.
(ii)
Qijp2 is the number of patient demand between D–P pairs i,j utilizing mode 2, “p” stands for patient, and Qijp1+Qijp2=Pij.
(iii)
Wjp1 is the average waiting time of patients j using mode 1.
(iv)
Wijp2 is the average waiting time of patients consulting through mode 2.
(v)
Cijp1 is the full costs of healthcare for patients’ medical services using mode 1.
(vi)
Cijp2 is the full costs of healthcare for patients’ medical services through mode 2.Considering the above costs of two service modes, patients will choose the way of medical services to minimize their own healthcare costs. Therefore, we can get the patients’ healthcare costs of two modes. As shown in(1)Cijp1=Fij+μ1Wjp1+μ2tij,Cijp2=Fij+μ1Wijp2+μ2tij+μ1Z¯,where μ1 and μ2 are the cost values of patients’ waiting time and medical service time; Z¯ denotes the patients’ waiting time due to doctor’s departure during consultation process of mode 2.Therefore, the probabilities of a patient consulting a doctor through modes 1 and 2 are as follows:(1)
P1: the probability of a patient consulting a doctor through mode 1 is given by the logit model [24], as shown in
(2)P1=Qijp1Qijp1+Qijp2=exp−ω1Cijp1exp−ω1Cijp1+exp−ω1Cijp2∀i,j∈φ.
(2)
P2: the probability of a patient consulting a doctor through mode 2 is given by the logit model, as shown in
(3)P2=Qijp2Qijp1+Qijp2=exp−ω1Cijp2exp−ω1Cijp1+exp−ω1Cijp2∀i,j∈ϕ,where ω1 is a nonnegative parameter and reflects the degree of uncertainty about the patients’ demand.
From (1), (2), and (3), we get
(4)Qijp1=Pij⋅P1∀i,j∈φ,Qijp2=Pij⋅P2∀i,j∈φ.


## 2.3. Activities of Doctors

Accordingly, we also define some variables for doctors’ activities in the network.(i)
N is the doctors’ activities time, including occupied doctors’ time and idle doctors’ time,
(ii)
Tijod is the number of occupied doctor (number of persons/h) between D–P pairs i,j and “od” represents occupied doctors.
(iii)
∑ij∈φTijodtij is the total medical service time of all occupied doctors per unit time.
(iv)
Tijid2 is the number of idle doctors utilizing mode 2 to receive patients’ requests, and “id” stands for idle doctors.
(v)
Wijid2 is the average waiting time of doctors receiving patients utilizing mode 2.
(vi)
Tiid1 is the number of idle doctors utilizing mode 1 to receive patients’ requests.
(vii)
Wiid1 is the average waiting time of doctors receiving patients’ requests through mode 1.
(viii)
VT is the total unoccupied doctors’ service time in one unit time, VT=∑ij∈φTiid1Wiid1+∑ij∈φTijid2Wijid2+Z¯.
(ix)
Uijid1 is the futility of idle doctors accepting patients using mode 1.
(x)
Uijid2 is the futility of idle doctors accepting patients through mode 2.
(xi)
T~iid1 is the total number of idle doctors that receive the next patient using mode 1.
(xii)
T~iid2 is the number of idle doctors seeking the next patient through mode 2.Considering the total activity time of doctors, we can obtain(5)N=∑ij∈φTijodtij+∑ij∈φTiid1Wiid1+∑ij∈φTijid2Wijid2+Z¯.After completing one consulting activity, doctors will determine whether to receive the next patient to minimize their futility. The futility functions for idle doctors [25] are as follows:
(6)Uijid1=−Fi¯+πtWiid1+ti¯−yj¯∀i,j∈φ,Uijid2=−Fij+πtWijid2+tij+Z¯−yj∀i,j∈φ,where πt represents the doctor activity cost in one unit time; Fi¯ and ti¯ denote the average consulting income and occupied waiting time of mode 1, with Fi¯=∑i,j∈φFijQijp1/∑i,j∈φQijp1 and ti¯=∑i,j∈φtijQijp1/∑i,j∈φQijp1; yj is the attractiveness of patient j’s disease information to idle doctors that receive the next patient through mode 2; and yj¯ represents the average attractiveness of patient j’s disease information to idle doctors that receive the next patient through mode 1, with yj¯=∑i,j∈φyjQijp1/∑i,j∈φQijp1.After analyzing the doctors’ activity behavior and the futility functions, we can know that the idle doctors can accept the next patient in the hospital (mode 1) or through m-Health applications (mode 2). And we can get(7)T~iid1=∑i∈ITiid1,(8)T~iid2=∑i,j∈φTijid2,(9)∑i,j∈φTijod=T~iid1+T~iid2.Formula (9) dictates that for a given hour, the number of doctors completing medical services should equal the number of idle doctors.Then, we can obtain the probabilities of an idle doctor receiving the next patient through modes 1 and 2.(3)
P3: the probability of an idle doctor seeking the next patient through mode 1 is given by the logit model
(10)P3=T~iid1T~iid1+T~iid2=exp−ω2Liid1exp−ω2Liid1+exp−ω2Liid2.
(4)
P4: the probability of an idle doctor seeking the next patient through mode 2 is given by the logit model
(11)P4=T~iid2T~iid1+T~iid2=exp−ω2Liid2exp−ω2Liid1+exp−ω2Liid2,where ω2is the dispersion coefficient of the logit choice model.The probabilities of choosing patientj as the next patient by an idle doctor through modes 1 and 2 are as follows:
(5)
P5: the probability of choosing patient j as the next patient by an idle doctor through mode 1 is shown in
(12)P5=Tijid1T~iid1=exp−ω3Uijid1∑i,j∈φexp−ω3Uijid1.
(6)
P6: the probability of choosing patient j as the next patient by an idle doctor through mode 2 is shown in
(13)P6=Tijid2T~iid2=exp−ω4Uijid2∑i,j∈φexp−ω4Uijid2,where ω3 and ω4 are the dispersion coefficients of the logit choice model.From above formulas, we can get the log sums of the futility of idle doctors seeking their next patient through modes 1 and 2, that areLiid1=−1/ω3ln∑iexp−ω3Uijid1 and Liid2=−1/ω4ln∑i,j∈φexp−ω4Uijid2 and also ω3≫ω2 and ω4≫ω2 [26].In a stationary equilibrium state, no matter how a patient sees the doctor, the patient will receive medical services after waiting a certain time. So(14)∑i∈ITiid1=∑i,j∈φQijp1,(15)∑i∈ITijid2=∑i,j∈φQijp2.In the network, medical service process is a mutual selection process in which patients can search for doctors and doctors also can search for patients.Based on the Cobb–Douglas-type production function [25], we can get the waiting time function of mode 1.If0<α1,α2<1, we get
(16)Wjp1=Aj−1/α1∑i,j∈φQijp11−α1/α1Wiid1∑i∈ITiid1−α2/α1,where α1 is the meeting rate elasticity of the number of unserved patients, α2 is the meeting rate elasticity of the number of idle doctors, Aj is a zone-specific parameter, Wiid1∑i∈ITiid1 denotes the total waiting hours of the idle doctor using mode 1 of zone i, and ∑i,j∈ϕQijp1 is the unserved patient number of zone i.Ifα1=α2=1, we obtain Wjp1=1/AjWiid1∑i∈ITiid1.Similarly, we can also get the waiting time function of mode 2.Wijp2 is shown as
(17)Wijp2=A~−1/α1∑i,j∈ϕQijp21−α1/α1Wijid2∑i,j∈ϕTijid2−α2/α1,where A~ is a constant number to reflect the little-friction property.


## 3. Stationary Equilibrium of Medical Service

### 3.1. Equilibrium Definition

The rapid development of mobile health technology (m-Health) provides unprecedented opportunities for improving health services. In the medical service market, patients have the need for recovery. They can consult doctors and get services, expecting to spend as little as possible. As service providers, doctors will naturally want to get more pay. Therefore, we divide the medical service market into two modes with and without the m-Health application in this paper. We suppose that each patient and doctor have no market power to affect the two modes and the medical service market.In order to meet the supply of doctors and demand of patients within the network, we can get the equilibrium(18)Tijid=Qijp1+Qijp2,where (4) describes the choices of patients; (5), (6), (7), (8), (9), (10), (11), (12), (13), (14), and (15) characterize doctor activities in the network; and (16) and (17) calculate the patient waiting time of two modes.


### 3.2. Solution

If the medical service demand of patients, the average medical service time, and regulated consultation fee between each D–P pair are fixed and given, a stationary equilibrium in the medical service network can be achieved. The condition of the equilibrium is a series of nonlinear equations; so, we apply Brouwer’s fixed-point theorem to get the solution of the equilibrium.Theorem 1 (Brouwer’s fixed-point theorem) [27].
Every continuous function from a convex compact subset of a Euclidean space to itself has a fixed point.
It states that for any continuous functionψ:H→H mapping a compact convex set H into itself, there is a point θ∈H such that ψθ=θ.According to Brouwer’s fixed-point theorem, we need to do the following two steps:Step 1.
We need to define a convex compact set.According to the equilibrium, we can defineH=…,Qijp1,…,Qijp2,…∣Qijp1+Qijp2=Pij,∀i,j∈φ;Qijp1,Qijp2≥0 as the convex compact set.Step 2.
We need to construct a continuous function.We consider the artificial entropy-type minimization problem (AE) as the continuous function:(19)minT∑Α1T~iid2L2+∑Α2T~iid1L1+∑1ω4Tijid2L3+∑1ω3Tijid1L4+∑Tijid2Β1+∑Tijid1Β2,where A1=1/ω2−1/ω4, A2=1/ω2−1/ω3, B1=−Fij+πttij+Z¯−yj, B2=−Fi¯+πtti¯−yi¯, L2=lnT~iid2−1, L1=lnT~iid1−1, L3=lnTijid2−1, and L4=lnTijid1−1.From (9) and (18), we can get
(20)T~iid1+T~iid2=∑∀i,j∈φQijp1+Qijp2.From (7), (8), (14), and (15), we can get
(21)1ω2−1ω3lnT~iid1−bi+ci=0,1ω2−1ω4lnT~iid2−ai+ci=0,
(22)1ω4lnTijid2+−Fij+πttij+Z¯−yj+eij+ai=0,1ω3lnTijid1+−Fij¯+πtti¯−yi¯+fj+bi=0,where bi is the Lagrangian multiplier with constraint (7), ai with constraint (8), fj with constraint (14), eij with constraint (15), and ci with constraint (20).From (22) and B1=−Fij+πttij+Z¯−yj and B2=−Fi¯+πtti¯−yi¯, we can get
(23)exp−ω4ai=Tijid2exp−ω4B1+eij.From (8), we get
(24)exp−ω4ai=T~iid2∑i,j∈φexp−ω4B1+eij,ai=−1ω4ln∑i,j∈φexp−ω4B1+eij−1ω4lnT~iid2.Similarly, we can also obtainbi=−1/ω3ln∑i,j∈φexp−ω3B2+fj−1/ω3lnT~iid1.Then fromai, bi, and (21), we have
(25)1ω2lnT~iid1−1ω3ln∑i,j∈φexp−ω3B2+fj+ci=0,(26)1ω2lnT~iid2−1ω4ln∑i,j∈φexp−ω4B1+eij+ci=0.From (25), we can get
(27)T~iid1=expω2ω3ln∑i,j∈φexp−ω3B2+fj−ω2ci.From (26), we get T~iid2=expω2/ω4ln∑i,j∈φexp−ω4B1+eij−ω2ci.We interpretfj and eij as πtWiid1 and πtWijid2, respectively. It is straightforward to verify (25) and (26).Similarly, we can get(28)Tijid2=exp−ω4ai⋅exp−ω4Β1+eij,Tiid1=exp−ω3bi⋅exp−ω3B2+fj.In other words, given any…,Qijp1,…,Qijp2,…∈H, we can solve AE to obtain …,T~iid1,…,T~iid2,…,Tijid1,…,Tijid2,….If we addg to fj and eij and make Wiid1=fj/πt+g and Wijid2=eij/πt+g, we can obtain …,T~iid1,…,T~iid2,…,Tijid1,…,Tijid2,…Wiid1,…,Wijid2,…,Wip1,…,Wijp2,….The self-mapping is(29)ψQp:H→H,where
(30)Qp=…,Qijp1,…,Qijp2,…,ψijp1Qp=Pij⋅P1Qp∀i,j∈φ,ψijp2Qp=Pij⋅P2Qp∀i,j∈φ,P1Qp=exp−ω1Fij+μ1Wijp2Qp+μ2tijexp−ω1Fij+μ1Wjp1Qp+μ2tij+exp−ω1Fij+μ1Wijp2Qp+μ2tij+μZ¯∀i,j∈φ,P2Qp=exp−ω1Fij+μ1Wijp2Qp+μ1tij+μZ¯exp−ω1Fij+μ1Wjp1Qp+μ2tij+exp−ω1Fij+μ1Wijp2Qp+μ2tij+μZ¯∀i,j∈φ.For…Wip1,…,Wijp2,…,, mapping ψ is continuous. Wip1Qp and Wijp2Qp vary continuously with Qp; so, we can conclude that mapping ψ is a continuous function of Qp.Now, we provide a numerical algorithm to find the equilibrium.Algorithm 1: Find the equilibrium.
Data: H
Result: n
1. n←0;
2. forQ0p∈Hdo
3.  …Wip1Qnp,…,Wijp2Qnp,…,←AEQnp;
4.  Qij,n+1p1←ψijp1Qp
5.  Qij,n+1p2←ψijp2Qp
6. if∑i,j∈φQij,n+1p1−Qij,np12+Qij,n+1p2−Qij,np22<η
then
7.   n← n+1;
8.   go to step 3;
9.   end
10. end
11. return n;


## 3.1. Equilibrium Definition

The rapid development of mobile health technology (m-Health) provides unprecedented opportunities for improving health services. In the medical service market, patients have the need for recovery. They can consult doctors and get services, expecting to spend as little as possible. As service providers, doctors will naturally want to get more pay. Therefore, we divide the medical service market into two modes with and without the m-Health application in this paper. We suppose that each patient and doctor have no market power to affect the two modes and the medical service market.In order to meet the supply of doctors and demand of patients within the network, we can get the equilibrium(18)Tijid=Qijp1+Qijp2,where (4) describes the choices of patients; (5), (6), (7), (8), (9), (10), (11), (12), (13), (14), and (15) characterize doctor activities in the network; and (16) and (17) calculate the patient waiting time of two modes.


## 3.2. Solution

If the medical service demand of patients, the average medical service time, and regulated consultation fee between each D–P pair are fixed and given, a stationary equilibrium in the medical service network can be achieved. The condition of the equilibrium is a series of nonlinear equations; so, we apply Brouwer’s fixed-point theorem to get the solution of the equilibrium.Theorem 1 (Brouwer’s fixed-point theorem) [27].
Every continuous function from a convex compact subset of a Euclidean space to itself has a fixed point.
It states that for any continuous functionψ:H→H mapping a compact convex set H into itself, there is a point θ∈H such that ψθ=θ.According to Brouwer’s fixed-point theorem, we need to do the following two steps:Step 1.
We need to define a convex compact set.According to the equilibrium, we can defineH=…,Qijp1,…,Qijp2,…∣Qijp1+Qijp2=Pij,∀i,j∈φ;Qijp1,Qijp2≥0 as the convex compact set.Step 2.
We need to construct a continuous function.We consider the artificial entropy-type minimization problem (AE) as the continuous function:(19)minT∑Α1T~iid2L2+∑Α2T~iid1L1+∑1ω4Tijid2L3+∑1ω3Tijid1L4+∑Tijid2Β1+∑Tijid1Β2,where A1=1/ω2−1/ω4, A2=1/ω2−1/ω3, B1=−Fij+πttij+Z¯−yj, B2=−Fi¯+πtti¯−yi¯, L2=lnT~iid2−1, L1=lnT~iid1−1, L3=lnTijid2−1, and L4=lnTijid1−1.From (9) and (18), we can get
(20)T~iid1+T~iid2=∑∀i,j∈φQijp1+Qijp2.From (7), (8), (14), and (15), we can get
(21)1ω2−1ω3lnT~iid1−bi+ci=0,1ω2−1ω4lnT~iid2−ai+ci=0,
(22)1ω4lnTijid2+−Fij+πttij+Z¯−yj+eij+ai=0,1ω3lnTijid1+−Fij¯+πtti¯−yi¯+fj+bi=0,where bi is the Lagrangian multiplier with constraint (7), ai with constraint (8), fj with constraint (14), eij with constraint (15), and ci with constraint (20).From (22) and B1=−Fij+πttij+Z¯−yj and B2=−Fi¯+πtti¯−yi¯, we can get
(23)exp−ω4ai=Tijid2exp−ω4B1+eij.From (8), we get
(24)exp−ω4ai=T~iid2∑i,j∈φexp−ω4B1+eij,ai=−1ω4ln∑i,j∈φexp−ω4B1+eij−1ω4lnT~iid2.Similarly, we can also obtainbi=−1/ω3ln∑i,j∈φexp−ω3B2+fj−1/ω3lnT~iid1.Then fromai, bi, and (21), we have
(25)1ω2lnT~iid1−1ω3ln∑i,j∈φexp−ω3B2+fj+ci=0,(26)1ω2lnT~iid2−1ω4ln∑i,j∈φexp−ω4B1+eij+ci=0.From (25), we can get
(27)T~iid1=expω2ω3ln∑i,j∈φexp−ω3B2+fj−ω2ci.From (26), we get T~iid2=expω2/ω4ln∑i,j∈φexp−ω4B1+eij−ω2ci.We interpretfj and eij as πtWiid1 and πtWijid2, respectively. It is straightforward to verify (25) and (26).Similarly, we can get(28)Tijid2=exp−ω4ai⋅exp−ω4Β1+eij,Tiid1=exp−ω3bi⋅exp−ω3B2+fj.In other words, given any…,Qijp1,…,Qijp2,…∈H, we can solve AE to obtain …,T~iid1,…,T~iid2,…,Tijid1,…,Tijid2,….If we addg to fj and eij and make Wiid1=fj/πt+g and Wijid2=eij/πt+g, we can obtain …,T~iid1,…,T~iid2,…,Tijid1,…,Tijid2,…Wiid1,…,Wijid2,…,Wip1,…,Wijp2,….The self-mapping is(29)ψQp:H→H,where
(30)Qp=…,Qijp1,…,Qijp2,…,ψijp1Qp=Pij⋅P1Qp∀i,j∈φ,ψijp2Qp=Pij⋅P2Qp∀i,j∈φ,P1Qp=exp−ω1Fij+μ1Wijp2Qp+μ2tijexp−ω1Fij+μ1Wjp1Qp+μ2tij+exp−ω1Fij+μ1Wijp2Qp+μ2tij+μZ¯∀i,j∈φ,P2Qp=exp−ω1Fij+μ1Wijp2Qp+μ1tij+μZ¯exp−ω1Fij+μ1Wjp1Qp+μ2tij+exp−ω1Fij+μ1Wijp2Qp+μ2tij+μZ¯∀i,j∈φ.For…Wip1,…,Wijp2,…,, mapping ψ is continuous. Wip1Qp and Wijp2Qp vary continuously with Qp; so, we can conclude that mapping ψ is a continuous function of Qp.Now, we provide a numerical algorithm to find the equilibrium.Algorithm 1: Find the equilibrium.
Data: H
Result: n
1. n←0;
2. forQ0p∈Hdo
3.  …Wip1Qnp,…,Wijp2Qnp,…,←AEQnp;
4.  Qij,n+1p1←ψijp1Qp
5.  Qij,n+1p2←ψijp2Qp
6. if∑i,j∈φQij,n+1p1−Qij,np12+Qij,n+1p2−Qij,np22<η
then
7.   n← n+1;
8.   go to step 3;
9.   end
10. end
11. return n;


## 4. Numerical Example

A simple numerical example is given to demonstrate the equilibrium of medical services with or without the m-Health application. Since the medical service demand of patients, the average medical service time, and regulated consultation fee of each D–P pair are given, and the equilibrium always exists.In this section, we first describe the mobile health medical service market network that we need. Then, we analyze some parameters to investigate the impact on the equilibrium, such as the number of doctors, the dispersion coefficientω4, and profitability index λ. Last, we evaluate the influence of the m-Health application and compare two types of m-Health applications, one of which displays patient disease information and the other displays both patient disease information and doctor registration information.


### 4.1. Results

We adapt the paired network in [28] that was designed for children in rural Sichuan Province, China, for improving immunization based on a smart phone app, and get the health medical service market network that we need. The network that we need is a four-node network with 12 links connecting 12 D–P pairs.Table1 presents the medical service demand of patients for each D–P pair. Table 2 shows the average medical service time of each D–P pair. In the following analyses, medical service demand of patients, average medical service time, and regulated consultation fee of each D–P pair are fixed and given. The consultation fare between D–P pair i,j is equal to the medical service time multiplied by the time unit cost, that is, Fij=$60/h. We make dispersion coefficients ω1=0.2, ω2=0.2, ω3=0.5, and ω4=0.5. Other parameters include α1=1, α2=1,yj=0, z¯=0, A~=104, Ai=1, πt=$10/h, μ1=$20/h, and μ2=$10/h.Table 1
Medical service demand (services/h).
D
P
1
2
3
4
1
0
10
4
4
2
8
0
3
5
3
4
2
0
10
4
2
4
6
0Table 2
Average medical service time (h).
D
P
1
2
3
4
1
0.1
0.25
0.3
0.35
2
0.25
0.1
0.3
0.45
3
0.3
0.3
0.1
0.25
4
0.35
0.45
0.25
0.1We set five cases of the number of doctorsn equal to 20, 50, 80, 100, and 150, and set Qijp1=0.7Pij, Qijp2=0.3Pij, ∀i,j∈ϕ as the initial solution.For each of the five cases, the proposed algorithm shows a convergence with the error defined in Step 6 of the algorithm.Figure1 plots the sensitivity analysis of the number of doctors.Figure 1
Sensitivity analysis of the number of doctors.
(a)
(b)
(c)Figure1(a) shows the doctors’ waiting time and average waiting time of two modes as the number of doctors increases. From Figure 1(a), we can see that the doctors’ waiting time and average waiting time of two modes increase with the increase of the number of doctors. The increasing trends of doctors’ waiting time of two modes are almost the same when the number of doctors is from 20 to 50. From the value of 50, the doctors’ waiting time of mode 2 increases rapidly than the mode 1’s. Figure 1(b) shows the patients’ waiting time and the average waiting time of two modes as the number of doctors increases. As we can see from Figure 1(b), the patients’ waiting time and the average waiting time of both modes decrease as the number of doctors increases. Patient waiting time of mode 2 almost remains unchanged and nearly equals zero. Patient waiting time of mode 1 decreases rapidly. Figure 1(c) shows the doctor utilization rate and the average utilization rate of two modes as the number of doctors increases. The doctor utilization rate is equal to the doctors’ occupied medical service hours divided by the total activity time. From Figure 1(c), as the number of doctors increases, the utilization rates of doctors of two modes decrease. When the number of doctors is 50, the utilization rates of doctors of two modes are the same.Figure2 shows the comparison result of the doctor waiting time with different doctor groups in each pair. To further investigate the difference of the two modes, we set the number of doctors n=100 and n is fixed. The different doctor groups can be distinguished by the indexes above the columns. For example, index “1-2” denotes the group of doctors who utilize mode 2 to accept patients between D–P pair 1,2 and the index “non-app” denotes the group of doctors using mode 1. As you can see from Figure 2, the doctor’s waiting time is different among different groups at a particular pair and the doctors’ waiting time of mode 1 is almost equal to the average value. In addition, for mode 2, the more consulting time patients need, the longer waiting time the doctors will spend. For instance, the medical service time of D–P pairs 2–4 and 4–2 is the highest in the network. The doctor waiting time associated with these two D–P pairs is also the highest among all from Figure 2.Figure 2
Doctor waiting time withn=100.In order to reflect the net income of doctors, we define the profitability indexλ. Hence, we adjust the futility functions of idle doctors with consideration of λ, as shown below:
(31)Uijid1=λ−Fi¯+πtti¯+πtWiid1∀i,j∈φ,Uijid2=λ−Fij+πttij+πtWijid2∀i,j∈φ.As service providers, doctors will naturally want to get more pay and try their best to maximize their net income. If the value is big enough, doctors are most likely to wait for a long time only for a single consulting activity.To examine the variation of doctor and patient waiting time, we estimated the magnitude of change of dispersion parameterω4 and profitability parameter λ. The results are shown in Figure 3. We can see from Figure 3(a) that the waiting time of the doctor of mode 1 increases when ω4 increases. However, the waiting time of the doctor of mode 2 decreases and the extent of variation is small. From Figure 3(c), as the increase of profitability parameter λ, the waiting time of the doctor of mode 1 increases and the waiting time of the doctor of mode 2 decreases. When λ is about 1.1, the two values of the doctor and patient waiting time are equal. In the meantime, as we can see from Figure 3(b), the waiting time of the patient of mode 1 decreases with the adjust of ω4. But the waiting time of the patient of mode 2 also decreases and almost remains zero when ω4 is 2.5. The waiting time of the patient of mode 1 increases and the value of mode 2 decreases with the increase of λ.Figure 3
Sensitivity analysis ofω4 and λ.
(a)
(b)
(c)
(d)In many m-Health applications, patients need to provide disease information that will be sent to idle doctors. Meanwhile, doctors also need to be registered with expertise and availability details. With the patient disease information and doctor registration information revealed in the application, patients and doctors can choose each other according to their own information. Besides, the case of m-Health application displays both patient disease information and doctor registration information. We also provide two other cases: one is no m-Health application in the network and the other is the m-Health application only displaying patient disease information. As shown in Figure4, the indexes “traditional network,” “displaying patient,” and “displaying both patient and doctor” denote no m-Health applications available, the m-Health application available but only displaying patient disease information, and the application displaying both patient disease information and doctor registration information, respectively. We can see from Figure 4 that the average doctor waiting time of the three cases has little difference, as well as the average waiting time of patient. The doctor and patient waiting time in the case which displays both patient disease information and doctor registration information is longer than that of the case which only displays patient disease information.Figure 4
Doctor and patient waiting time.


### 4.2. Discussion

As shown in Figure1, the doctor waiting time of both modes and their average time increase while the utilization rates of doctors decrease with the increase of the number of doctors. The number of doctors has a significant impact on the waiting time and the utilization rates of doctors. Our analysis suggests that this is the result of the extra medical service time due to the increase of the number of doctors. And the extra medical service time is mostly devoted to waiting for patients. Therefore, the percentage of the occupied doctor hours decreases. In addition, we also can see that the doctor waiting time of mode 2 increases while the utilization rate of doctors of mode 2 decreases more rapidly than that of mode 1. From the above phenomenon, we may draw the conclusion that the doctor waiting time and utilization rate of mode 2 are more sensitive to the number of doctors in a small example. We may conclude that the service of mode 2 is more attractive to doctors than that of mode 1. With the increase of the number of doctors, the competition among the doctors of mode 2 is more intense. Besides, patient waiting time of mode 2 almost remains unchanged and nearly equals zero. This is most likely because the increase of the number of doctors makes it easier to see the doctors in hospital and communicating through the m-Health application exhibits little friction.For Figure2, our finding that the more the consulting time patients need, the longer the waiting time doctors will spend. This may be because idle doctors can make deterministic choices.We can analyze the futility value of mode 2.IfUijid2=Ui'j'id2, that is −Fij+πtWijid2+tij=−Fi'j'+πtWi'j'id2+ti'j'∀i,j,i',j'∈φ;i,j≠i',j'.If−Fij+πttij<−Fi'j'+πtti'j', then tij>ti'j' and Wijid2>Wi'j'id2. From the above formulas, if a D–P pair has a higher value of medical service time, then the doctor will also suffer longer waiting time. Because the doctors’ futility values of mode 1 are based on the medical service time and consulting fee, so, the doctor waiting time of mode 1 is nearly equal to the average time of each pair.The phenomenon is shown by the data of Figure3 and experimental studies indicate that idle doctors’ choices of choosing patient j as the next patient are more certain with the increase of dispersion parameter ω4. In order to get more income and reduce the cost of waiting time, more and more doctors choose mode 2 to receive patients. Therefore, the waiting time of patients decreases. Because the waiting time of patients is less than mode 1, mode 2 is more attractive to patients. Because the demand of patients is fixed and given in the network, idle doctors of mode 1 need to spend more time waiting for their next patient. Due to the little friction of mode 2, the waiting time of patient is almost zero. As idle doctors are inclined to maximize their net income and minimize idle doctors’ futility, more and more doctors will choose to use mode 2 with the increase of profitability index λ.As we can observe from Figure4, the medical service supply of patients can be cut down by introducing m-Health applications into the medical service market. Hiding the information of patient disease can increase the waiting time of patients of mode 2. However, providing disease information of patient can increase the attraction degree of mode 2. The average doctor waiting time of the three cases has little difference, as well as the average waiting time of patients. This could be because some patients may spend more time seeing the doctor using mode 2. On the contrary, patients spend less time seeing the doctor through mode 1.


## 4.1. Results

We adapt the paired network in [28] that was designed for children in rural Sichuan Province, China, for improving immunization based on a smart phone app, and get the health medical service market network that we need. The network that we need is a four-node network with 12 links connecting 12 D–P pairs.Table1 presents the medical service demand of patients for each D–P pair. Table 2 shows the average medical service time of each D–P pair. In the following analyses, medical service demand of patients, average medical service time, and regulated consultation fee of each D–P pair are fixed and given. The consultation fare between D–P pair i,j is equal to the medical service time multiplied by the time unit cost, that is, Fij=$60/h. We make dispersion coefficients ω1=0.2, ω2=0.2, ω3=0.5, and ω4=0.5. Other parameters include α1=1, α2=1,yj=0, z¯=0, A~=104, Ai=1, πt=$10/h, μ1=$20/h, and μ2=$10/h.Table 1
Medical service demand (services/h).
D
P
1
2
3
4
1
0
10
4
4
2
8
0
3
5
3
4
2
0
10
4
2
4
6
0Table 2
Average medical service time (h).
D
P
1
2
3
4
1
0.1
0.25
0.3
0.35
2
0.25
0.1
0.3
0.45
3
0.3
0.3
0.1
0.25
4
0.35
0.45
0.25
0.1We set five cases of the number of doctorsn equal to 20, 50, 80, 100, and 150, and set Qijp1=0.7Pij, Qijp2=0.3Pij, ∀i,j∈ϕ as the initial solution.For each of the five cases, the proposed algorithm shows a convergence with the error defined in Step 6 of the algorithm.Figure1 plots the sensitivity analysis of the number of doctors.Figure 1
Sensitivity analysis of the number of doctors.
(a)
(b)
(c)Figure1(a) shows the doctors’ waiting time and average waiting time of two modes as the number of doctors increases. From Figure 1(a), we can see that the doctors’ waiting time and average waiting time of two modes increase with the increase of the number of doctors. The increasing trends of doctors’ waiting time of two modes are almost the same when the number of doctors is from 20 to 50. From the value of 50, the doctors’ waiting time of mode 2 increases rapidly than the mode 1’s. Figure 1(b) shows the patients’ waiting time and the average waiting time of two modes as the number of doctors increases. As we can see from Figure 1(b), the patients’ waiting time and the average waiting time of both modes decrease as the number of doctors increases. Patient waiting time of mode 2 almost remains unchanged and nearly equals zero. Patient waiting time of mode 1 decreases rapidly. Figure 1(c) shows the doctor utilization rate and the average utilization rate of two modes as the number of doctors increases. The doctor utilization rate is equal to the doctors’ occupied medical service hours divided by the total activity time. From Figure 1(c), as the number of doctors increases, the utilization rates of doctors of two modes decrease. When the number of doctors is 50, the utilization rates of doctors of two modes are the same.Figure2 shows the comparison result of the doctor waiting time with different doctor groups in each pair. To further investigate the difference of the two modes, we set the number of doctors n=100 and n is fixed. The different doctor groups can be distinguished by the indexes above the columns. For example, index “1-2” denotes the group of doctors who utilize mode 2 to accept patients between D–P pair 1,2 and the index “non-app” denotes the group of doctors using mode 1. As you can see from Figure 2, the doctor’s waiting time is different among different groups at a particular pair and the doctors’ waiting time of mode 1 is almost equal to the average value. In addition, for mode 2, the more consulting time patients need, the longer waiting time the doctors will spend. For instance, the medical service time of D–P pairs 2–4 and 4–2 is the highest in the network. The doctor waiting time associated with these two D–P pairs is also the highest among all from Figure 2.Figure 2
Doctor waiting time withn=100.In order to reflect the net income of doctors, we define the profitability indexλ. Hence, we adjust the futility functions of idle doctors with consideration of λ, as shown below:
(31)Uijid1=λ−Fi¯+πtti¯+πtWiid1∀i,j∈φ,Uijid2=λ−Fij+πttij+πtWijid2∀i,j∈φ.As service providers, doctors will naturally want to get more pay and try their best to maximize their net income. If the value is big enough, doctors are most likely to wait for a long time only for a single consulting activity.To examine the variation of doctor and patient waiting time, we estimated the magnitude of change of dispersion parameterω4 and profitability parameter λ. The results are shown in Figure 3. We can see from Figure 3(a) that the waiting time of the doctor of mode 1 increases when ω4 increases. However, the waiting time of the doctor of mode 2 decreases and the extent of variation is small. From Figure 3(c), as the increase of profitability parameter λ, the waiting time of the doctor of mode 1 increases and the waiting time of the doctor of mode 2 decreases. When λ is about 1.1, the two values of the doctor and patient waiting time are equal. In the meantime, as we can see from Figure 3(b), the waiting time of the patient of mode 1 decreases with the adjust of ω4. But the waiting time of the patient of mode 2 also decreases and almost remains zero when ω4 is 2.5. The waiting time of the patient of mode 1 increases and the value of mode 2 decreases with the increase of λ.Figure 3
Sensitivity analysis ofω4 and λ.
(a)
(b)
(c)
(d)In many m-Health applications, patients need to provide disease information that will be sent to idle doctors. Meanwhile, doctors also need to be registered with expertise and availability details. With the patient disease information and doctor registration information revealed in the application, patients and doctors can choose each other according to their own information. Besides, the case of m-Health application displays both patient disease information and doctor registration information. We also provide two other cases: one is no m-Health application in the network and the other is the m-Health application only displaying patient disease information. As shown in Figure4, the indexes “traditional network,” “displaying patient,” and “displaying both patient and doctor” denote no m-Health applications available, the m-Health application available but only displaying patient disease information, and the application displaying both patient disease information and doctor registration information, respectively. We can see from Figure 4 that the average doctor waiting time of the three cases has little difference, as well as the average waiting time of patient. The doctor and patient waiting time in the case which displays both patient disease information and doctor registration information is longer than that of the case which only displays patient disease information.Figure 4
Doctor and patient waiting time.


## 4.2. Discussion

As shown in Figure1, the doctor waiting time of both modes and their average time increase while the utilization rates of doctors decrease with the increase of the number of doctors. The number of doctors has a significant impact on the waiting time and the utilization rates of doctors. Our analysis suggests that this is the result of the extra medical service time due to the increase of the number of doctors. And the extra medical service time is mostly devoted to waiting for patients. Therefore, the percentage of the occupied doctor hours decreases. In addition, we also can see that the doctor waiting time of mode 2 increases while the utilization rate of doctors of mode 2 decreases more rapidly than that of mode 1. From the above phenomenon, we may draw the conclusion that the doctor waiting time and utilization rate of mode 2 are more sensitive to the number of doctors in a small example. We may conclude that the service of mode 2 is more attractive to doctors than that of mode 1. With the increase of the number of doctors, the competition among the doctors of mode 2 is more intense. Besides, patient waiting time of mode 2 almost remains unchanged and nearly equals zero. This is most likely because the increase of the number of doctors makes it easier to see the doctors in hospital and communicating through the m-Health application exhibits little friction.For Figure2, our finding that the more the consulting time patients need, the longer the waiting time doctors will spend. This may be because idle doctors can make deterministic choices.We can analyze the futility value of mode 2.IfUijid2=Ui'j'id2, that is −Fij+πtWijid2+tij=−Fi'j'+πtWi'j'id2+ti'j'∀i,j,i',j'∈φ;i,j≠i',j'.If−Fij+πttij<−Fi'j'+πtti'j', then tij>ti'j' and Wijid2>Wi'j'id2. From the above formulas, if a D–P pair has a higher value of medical service time, then the doctor will also suffer longer waiting time. Because the doctors’ futility values of mode 1 are based on the medical service time and consulting fee, so, the doctor waiting time of mode 1 is nearly equal to the average time of each pair.The phenomenon is shown by the data of Figure3 and experimental studies indicate that idle doctors’ choices of choosing patient j as the next patient are more certain with the increase of dispersion parameter ω4. In order to get more income and reduce the cost of waiting time, more and more doctors choose mode 2 to receive patients. Therefore, the waiting time of patients decreases. Because the waiting time of patients is less than mode 1, mode 2 is more attractive to patients. Because the demand of patients is fixed and given in the network, idle doctors of mode 1 need to spend more time waiting for their next patient. Due to the little friction of mode 2, the waiting time of patient is almost zero. As idle doctors are inclined to maximize their net income and minimize idle doctors’ futility, more and more doctors will choose to use mode 2 with the increase of profitability index λ.As we can observe from Figure4, the medical service supply of patients can be cut down by introducing m-Health applications into the medical service market. Hiding the information of patient disease can increase the waiting time of patients of mode 2. However, providing disease information of patient can increase the attraction degree of mode 2. The average doctor waiting time of the three cases has little difference, as well as the average waiting time of patients. This could be because some patients may spend more time seeing the doctor using mode 2. On the contrary, patients spend less time seeing the doctor through mode 1.


## 5. Summary and Outlook

In the medical service market, patients have the need for recovery. They expect to get high-quality medical service and reduce healthcare spending. Their goals are typically to match doctor capacity with their demand and utilize medical personnel more efficiently without significantly increasing healthcare spending. As service providers, doctors are under pressure to improve healthcare quality and balance workloads. Meanwhile, they hope for a high quota reward. Thus, the supply and demand of medical services have become an increasingly essential component in a medical service market. In this study, we intend to propose an equilibrium model to evaluate the influence of m-Health applications on the medical service market. The model can balance the supply of doctors and demand of patients and reflect possible options for both doctors and patients with or without m-Health applications in the medical service market. In the meantime, we analyze the behavior of patients and the activities of doctors to minimize patients’ full costs of healthcare and doctors’ futility. After that, we provide a resolution algorithm by mathematical reasoning. Lastly, experiments are conducted on artificially generated dataset to evaluate the medical services of m-Health applications.The study does have some limitations. In order to provide a clear description of the medical services, we have assumed that the medical service time, the medical services demand of patients, and the consultation fee of each D–P pair are fixed and given. The m-Health applications are assumed to be widely used among patients. Our future research will relax those assumptions and so on.



---

*Source: 1385034-2018-04-10.xml*