Prehospital emergency care patient satisfaction scale [PECPSS] for care provided by emergency medical teams: Scale development and validation

Junpei Haruna; Nobuyasu Hayasaka; Yukiko Taguchi; Saori Muranaka; Sachi Niiyama; Hirotoshi Inamura; Shuji Uemura; Keigo Sawamoto; Hirotoshi Mizuno; Nobuaki Himuro; Eichi Narimatsu; Junpei Haruna; Nobuyasu Hayasaka; Yukiko Taguchi; Saori Muranaka; Sachi Niiyama; Hirotoshi Inamura; Shuji Uemura; Keigo Sawamoto; Hirotoshi Mizuno; Nobuaki Himuro; Eichi Narimatsu

doi:10.3934/publichealth.2023011

AIMS Public Health

2023, Volume 10, Issue 1: 129-144. doi: 10.3934/publichealth.2023011

Previous Article Next Article

Research article

Prehospital emergency care patient satisfaction scale [PECPSS] for care provided by emergency medical teams: Scale development and validation

1.
Department of Intensive Care Medicine, School of Medicine, Sapporo Medical University, South-1, West-16, Chuo-ku, Sapporo, Hokkaido, 060–8543, Japan
2.
Emergency Sector of Ebetsu City Fire Department, Yoyogicho, 80–8, Nopporo, Hokkaido, 069–0817, Japan
3.
Department of Nursing, School of Health Sciences, Sapporo Medical University, South-1, West-17, Chuo-ku, Sapporo, Hokkaido, 060–8556, Japan
4.
Department of Advanced Critical Care and Emergency Center, Sapporo Medical University Hospital, South-1, West-16, Chuo-ku, Sapporo, Hokkaido, 060–8543, Japan
5.
Department of Pharmacy, Sapporo Medical University Hospital, South-1, West-16, Chuo-ku, Sapporo, Hokkaido, 060–8543, Japan
6.
Department of Emergency Medicine, Sapporo Medical University School of Medicine, South 1, West 16, Chuo-ku, Sapporo, Hokkaido, 060–8543, Japan
7.
Department of Public Health, School of Medicine, Sapporo Medical University, Sapporo, Japan
8.
Department of Emergency Medical Services, Life Flight and Disaster medicine, Sapporo Medical University, South 1, West 16, Chuo-ku, Sapporo, Hokkaido, 060–8543, Japan

Received: 01 November 2022 Revised: 28 January 2023 Accepted: 21 February 2023 Published: 01 March 2023

The purpose of this study was to develop and validate an emergency medical technician (EMT) care patient satisfaction scale to measure patient satisfaction with prehospital emergency care. To date, patient satisfaction surveys of EMTs have been performed subjectively, e using each facility's questionnaire, without the use of a validated patient satisfaction scale. However, no specific scale has been devised to assess patient satisfaction with EMTs. The study population comprised patients who used an ambulance between November 2020 and May 2021 (N = 202). A survey instrument was administered to participants who provided informed consent. In the process of validating the patient satisfaction scale, an exploratory factor analysis (EFA) of construct validity was performed. The results of the EFA showed a factor structure consisting of five factors: “teamwork”, “explanation and communication”, “physical treatment and psychological support”, “quickness of transport”, and “environment in the ambulance”. In addition, domain and summary scores showed good internal reliability (Cronbach's range = 0.82–0.94). The patient satisfaction scale developed in this study was designed and validated considering the role of EMTs and patients' needs for prehospital care. This scale may be useful in the development of assessments and interventions to improve patient satisfaction with EMTs.

Keywords:

Citation: Junpei Haruna, Nobuyasu Hayasaka, Yukiko Taguchi, Saori Muranaka, Sachi Niiyama, Hirotoshi Inamura, Shuji Uemura, Keigo Sawamoto, Hirotoshi Mizuno, Nobuaki Himuro, Eichi Narimatsu. Prehospital emergency care patient satisfaction scale [PECPSS] for care provided by emergency medical teams: Scale development and validation[J]. AIMS Public Health, 2023, 10(1): 129-144. doi: 10.3934/publichealth.2023011

Related Papers:

[1]	Deniz Sevinç . Volatility spillovers among MIST stock markets. Data Science in Finance and Economics, 2022, 2(2): 80-95. doi: 10.3934/DSFE.2022004
[2]	Nurun Najwa Bahari, Hafizah Bahaludin, Munira Ismail, Fatimah Abdul Razak . Network, correlation, and community structure of the financial sector of Bursa Malaysia before, during, and after COVID-19. Data Science in Finance and Economics, 2024, 4(3): 362-387. doi: 10.3934/DSFE.2024016
[3]	Moses Khumalo, Hopolang Mashele, Modisane Seitshiro . Quantification of the stock market value at risk by using FIAPARCH, HYGARCH and FIGARCH models. Data Science in Finance and Economics, 2023, 3(4): 380-400. doi: 10.3934/DSFE.2023022
[4]	Jiamin Cheng, Yuanying Jiang . How can carbon markets drive the development of renewable energy sector? Empirical evidence from China. Data Science in Finance and Economics, 2024, 4(2): 249-269. doi: 10.3934/DSFE.2024010
[5]	Kuo-Shing Chen . Interlinkages between Bitcoin, green financial assets, oil, and emerging stock markets. Data Science in Finance and Economics, 2024, 4(1): 160-187. doi: 10.3934/DSFE.2024006
[6]	Dachen Sheng, Opale Guyot . Market power, internal and external monitoring, and firm distress in the Chinese market. Data Science in Finance and Economics, 2024, 4(2): 285-308. doi: 10.3934/DSFE.2024012
[7]	Zimei Huang, Zhenghui Li . What reflects investor sentiment? Empirical evidence from China. Data Science in Finance and Economics, 2021, 1(3): 235-252. doi: 10.3934/DSFE.2021013
[8]	Katleho Makatjane, Ntebogang Moroke . Examining stylized facts and trends of FTSE/JSE TOP40: a parametric and Non-Parametric approach. Data Science in Finance and Economics, 2022, 2(3): 294-320. doi: 10.3934/DSFE.2022015
[9]	Samuel Asante Gyamerah, Collins Abaitey . Modelling and forecasting the volatility of bitcoin futures: the role of distributional assumption in GARCH models. Data Science in Finance and Economics, 2022, 2(3): 321-334. doi: 10.3934/DSFE.2022016
[10]	Nitesha Dwarika . The risk-return relationship in South Africa: tail optimization of the GARCH-M approach. Data Science in Finance and Economics, 2022, 2(4): 391-415. doi: 10.3934/DSFE.2022020

Abstract

1. Introduction

When searching the literature, we noticed that several propositions have been made by several researchers to find kernels that can be used to obtain fractional differential operators. The main reason for this is that, real world problems show signs of processes resembling behaviors of some mathematical functions. Riemann, Liouville, Cauchy and Abel works have led to a fractional calculus with a power law kernel. Their work was later modified by Caputo; this version has been used in many fields of science because of its ability to allow classical initial conditions ^[1]. Prabhakar suggested a different kernel as a product of power-law and the generalized Mittag-Leffler function with three parameters. This version has also attracted the attention of many researchers, studies have been done on theory as well as on applications. Indeed, the two kernels have their specific values, for example, power-law helps only to replicate processes exhibiting power-law behavior, while the product of the power-law and the generalized three-parameters Mittag-Leffler has also its field of application ^[2]. As nature is complex, a new kernel was suggested by Caputo and Fabrizio, a special exponential kernel with Delta Dirac properties. A differential operator that is in fashion nowadays due to its ability to replicate processes following fading memory. Indeed, this kernel brought a new shift into fractional calculus, as the concept of the fractional derivative with a non-singular kernel was introduced ^[3]. A point made by some researcher about the non-fractionality of the kernel led to a new kernel, the generalized Mittag-Leffler function with one parameter. This version was suggested by Atangana and Baleanu, constituting another new step forward in the field of fractional calculus. The operators have been used in several fields of studies with great success ^[4]. In ^[5], the basic theory of fractional differentiation, existence-uniqueness theorems and analytical numerical methods of solution of fractional differential equations are presented. In ^[6], the authors examined Noether's theorems of fractional generalized Birkhoffian systems in terms of classical and combined Caputo derivatives. The development of the time-fractional damage model for hyperelastic body is considered in ^[7]. While looking at nature and its complexities, one can with no doubt conclude that these suggested kernels are not enough to predict all complex behaviors of our universe. On this note, one will proceed to search for a different kernel or modified kernel, or a class of functions that will be used to introduce new differential operators. Sabatier recently presented some variants of kernels that will also open new doors for investigation ^[8]. In addition to these outstanding contributions, several other ideas have been suggested, for example the concept of short memory was suggested and a fractional derivative in the Caputo sense has been defined for different values of fractional orders. The idea is to have a different type of variable-order derivative unlike the well-known version that considers a fractional orders. The idea is initiated to have a different type of variable order derivative unlike the well-known version that considers a fractional-order to be a function of time. This case was suggested by Wu et al. ^[9] and applied in chaos. On the other hand, researchers have noticed that several real-world problems exhibit processes with different behaviors as a function of time and space. A particular case is a passage from deterministic to stochastic, or from power law to exponential decay. It was noted that existing differential operators may not be able to account for these behaviors, thus piecewise differential and integral operators were introduced to deal with problems exhibiting crossover behaviors ^[10]. The main aim of this note is to provide a critical analysis of the possible applications, advantages, and disadvantages of these two concepts.

2. Motivation for piecewise derivatives

We will illustrate the motivation with some examples.

Death body decay in different temperatures. Consider a corpse found in a snowy place, and assume that such body has been found after 20 days. The corpse is taken and brought to the house and kept in a normal temperature for a few days and later put in a mortuary and buried. The main aim here is to replicate the process of decay. The first part will provide a very low decay. The second part will provide a fast decay, the third part will again be a fast process and finally, a slow process. Indeed, slow and fast processes can be characterized by some mathematical functions. The power-law function is

$\begin{equation} t^{-\alpha }. \end{equation}$

(2.1)

The exponential decay function with a Dirac delta property is

$\begin{equation} \frac{1}{1-\alpha }\exp \left[ -\frac{\alpha }{1-\alpha }t\right] . \end{equation}$

(2.2)

However, a crossover from fast to slow decay processes can be modeled using the following Mittag-Leffler function with Delta-Dirac property

$\begin{equation} \frac{1}{1-\alpha }E_{\alpha }\left[ -\frac{\alpha }{1-\alpha }t^{\alpha } \right] . \end{equation}$

(2.3)

The process can be divided into several intervals to capture each behavior. In the first part, one can have

$\begin{equation} _{0}^{C}D_{t}^{\alpha }y\left( t\right) = -\lambda _{1}y\left( t\right) \;{\rm{ if }}\;0 < t\leq T_{1}. \end{equation}$

(2.4)

The second processes can be

$\begin{equation} _{T_{1}}^{CF}D_{t}^{\alpha }y\left( t\right) = -\lambda _{2}y\left( t\right) \;{\rm{ if }}\;T_{1} < t\leq T_{2}. \end{equation}$

(2.5)

The two last parts will be characterized by

$\begin{equation} _{T_{2}}^{ABC}D_{t}^{\alpha }y\left( t\right) = -\lambda _{3}y\left( t\right) \; {\rm{ if }}\;T_{2} < t\leq T. \end{equation}$

(2.6)

Thus the whole process will be a system with the following crossover behaviors

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }y\left( t\right) = -\lambda _{1}y\left( t\right) \;{\rm{ if }}\;0 < t\leq T_{1}, \\ _{T_{1}}^{CF}D_{t}^{\alpha }y\left( t\right) = -\lambda _{2}y\left( t\right) \;{\rm{ if }}\;T_{1} < t\leq T_{2}, \\ _{T_{2}}^{ABC}D_{t}^{\alpha }y\left( t\right) = -\lambda _{3}y\left( t\right) \;{\rm{ if }}\;T_{2} < t\leq T. \end{array} \right. \end{equation}$

(2.7)

Therefore in general the concept of piecewise derivative has been introduced.

3. Short memory concept

In this section, we present the definitions of both short memory and piecewise differentiation. For the short memory case, the idea has already been discussed; for example a paper published by Deng in 2007 has already discussed the short memory principle ^[11], which has also been called the fixed memory principle and logarithmic memory principle. Wu et al. ^[9] presented the concept of a fractional variable-order derivative where the fractional order changes within an interval. The suggested definition is given below as:

$\begin{equation} \left\{ \begin{array}{c} _{t_{0}}^{C}D_{t}^{\alpha _{0}}y\left( t\right)=f\left( y,t\right) \text{ for }t\in \left[ t_{0},t_{1}\right] , \\ _{t_{1}}^{C}D_{t}^{\alpha _{1}}y\left( t\right)=f\left( y,t\right) \text{ for }t\in \left[ t_{1},t_{2}\right] , \\ _{t_{2}}^{C}D_{t}^{\alpha _{2}}y\left( t\right)=f\left( y,t\right) \text{ for }t\in \left[ t_{2},t_{3}\right] . \end{array} \right. \end{equation}$

(3.1)

Here Caputo power law derivative is used, which is known to have a singularity at the origin for a class of functions. The author did not give explicitly space of the functions as the derivative was defined in general, therefore, if we assume a class of the functions for which $\frac{dy}{dt}$ is continuous then at each $t=t_{i},$ $f\left( y,t\right)$ should be zero.

On the other hand, the piecewise concept is defined as follows.

$1)$ A piecewise derivative of $f$ with classical, power-law and stochastic processes is given by

$\begin{equation} \left\{ \begin{array}{c} Dy\left( t\right) = f\left( t,y\left( t\right) \right) \;{\rm{ if }}\;0 < t \leq t_{1}, \\ _{t_{1}}^{C}D_{t}^{\alpha }y\left( t\right) = f\left( t,y\left( t\right) \right) \;{\rm{ if }}\;t_{1} < t\leq t_{2}, \\ dy\left( t\right) = f\left( t,y\left( t\right) \right) dt+\sigma y\left( t\right) dB\left( t\right)\; {\rm{ if }}\;t_{2} < t \leq t_{3}. \end{array} \right. \end{equation}$

(3.2)

$2)$ A piecewise derivative of $f$ with power, exponential and Mittag-Leffler laws is given as

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }y\left( t\right) = f\left( t,y\left( t\right) \right) \;{\rm{ if }}\;0 < t \leq t_{1}, \\ _{t_{1}}^{CF}D_{t}^{\alpha }y\left( t\right) = f\left( t,y\left( t\right) \right) \;{\rm{ if }}\;t_{1} < t \leq t_{2}, \\ _{t_{1}}^{ABC}D_{t}^{\alpha }y\left( t\right) = f\left( t,y\left( t\right) \right) \;{\rm{ if }}\;t_{2} < t \leq t_{3}. \end{array} \right. \end{equation}$

(3.3)

Several definitions can be found in ^[10] $.$

In the next section, we will provide a deep analysis of these concepts.The aim of the section is to present advantages and disadvantages of these two independent concepts. We note that

$\begin{eqnarray} _{0}^{C}D_{t}^{\alpha }{\rm{ }}f\left( t\right) & = &\frac{1}{\Gamma \left( 1-\alpha \right) }\int\limits_{0}^{t}\frac{d}{d\tau }f\left( \tau \right) \left( t-\tau \right) ^{-\alpha }d\tau ,{\rm{ \ }} \end{eqnarray}$

(3.4)

$\begin{eqnarray} _{0}^{CF}D_{t}^{\alpha }f\left( t\right) & = &\frac{1}{1-\alpha } \int\limits_{0}^{t}\frac{d}{d\tau }f\left( \tau \right) \exp \left[ -\frac{ \alpha }{1-\alpha }\left( t-\tau \right) \right] d\tau , \end{eqnarray}$

(3.5)

$\begin{eqnarray} _{0}^{ABC}D_{t}^{\alpha }f\left( t\right) & = &\frac{1}{1-\alpha } \int\limits_{0}^{t}\frac{d}{d\tau }f\left( \tau \right) E_{\alpha }\left[ - \frac{\alpha }{1-\alpha }\left( t-\tau \right) ^{\alpha }\right] d\tau . \end{eqnarray}$

(3.6)

4. Analysis and difference

Here, we present a deep analysis that will help the readers to see the difference between the two concepts. Let us start with short memory. A deep look inside the short memory concept defined by gives the following:

$1)$ The short memory principle considers a change in constant variable order; however, it uses a single kernel. Indeed, this can be viewed as variable order where the order is changed within shorter intervals. Nonetheless, it describes some process and in the case of Caputo derivation, this process will only describe a power-law process that is, this process will not have crossover behavior.

$2)$ The concept is unable to capture classical processes because the fractional order should change in each interval.

On a serious note, however, if we have $f\left(t, y\left(t\right) \right) \neq 0, \forall t\in \left[ 0, T\right].$ So at each $t_{j}$

$\begin{equation} \lim\limits_{t\rightarrow t_{j}}{\rm{ }}_{t_{j}}^{C}D_{t}^{\alpha }y\left( t\right) = f\left( t,y\left( t\right) \right) \mid _{t = t_{j}}. \end{equation}$

(4.1)

Since $\frac{dy}{dt}$ is continuous,

$\begin{equation} \lim\limits_{t\rightarrow t_{j}}{\rm{ }}_{t_{j}}^{C}D_{t}^{\alpha }y\left( t\right) = f\left( t_{j},y\left( t_{j}\right) \right) \neq 0. \end{equation}$

(4.2)

However the right hand side produces $f\left( t,y\left( t\right) \right) \neq 0$ , which is a contradiction. It infers that $f\left( t_{j},y\left( t_{j}\right) \right) =0.$ This may be used to explain the process with a renewal force that follows the power-law process. For example, the trajectory of a flea with a constant jump. There are fewest real world problems that present these behaviors. However, due to the power-law singularity the definition suggested by Wu may have some problems at the boundaries. Nonetheless, as a result of using non-singular kernels, the concept of short memory will be well-posed since the conditions at the boundaries will be well controlled and the renewal processes would be well-posed

On the other hand, however, the concept of piecewise was introduced for different purposes. The following example will give light to the situation. We consider evaluating the velocity of water within a geological formation with fractures. We record the velocity as a function of time. The velocity will obviously be low in the rock's matrix; however as the water reaches the fracture, there will be a crossover behavior because the velocity will suddenly increase. The first part of this process follows behaviors resembling a declining process and, later, an almost constant high velocity. Now if one was to withdraw water from such a geological formation and record the water level, one will observe a fading groundwater level in earlier times; then when water is being taken from the fracture, we observe a steady groundwater level. Thus, there is a crossover behavior from the fading process to long range, which can be captured using an exponential function, then the power law. Therefore, the differential operator able to replicate this process is

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }h\left( r,t\right) = f\left( t,r,h\left( r,t\right) \right) \;{\rm{ if }}\;0 < t\leq T_{1}, \\ _{T_{1}}^{CF}D_{t}^{\alpha }h\left( r,t\right) = f_{1}\left( t,r,h\left( r,t\right) \right) \;{\rm{ if }}\;T_{1} < t\leq T. \end{array} \right. \end{equation}$

(4.3)

Therefore, one can see a clear difference and objectives of both concepts. Piecewise differential operators are thus conceived to capture processes exhibiting different patterns, or crossover behaviors. Here the order does not change within an interval as changing order does not change the process but only the memory with the same pattern, rather, the kernels change to bring into the mathematical formulation the crossover behaviors by each kernel.

For example, using the piecewise derivative the decay equation

$\begin{equation} \left\{ \begin{array}{c} y\prime \left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;0 < t\leq T_{1}, \\ _{T_{1}}^{C}D_{t}^{\alpha }y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;T_{1} < t\leq T_{2}, \\ dy\left( t\right) = -\lambda y\left( t\right) dt+\sigma y\left( t\right) dB\left( t\right) \;{\rm{ if }}\;T_{2} < t \leq T \end{array} \right. \end{equation}$

(4.4)

which will lead to decay with crossover behaviors. From $0 < t < T_{1}$ , we shall observe a normal decay process. From $T_{1} < t < T_{2}$ , we shall observe power-law decay or the Mittag-Leffler decay. From $T_{2} < t < T$ , we shall observe decay with randomness. On the other hand, with the short memory

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha _{1}}y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;0 < t \leq T_{1}, \\ _{T_{1}}^{C}D_{t}^{\alpha _{2}}y\left( t\right) = -\lambda y\left( t\right) \; {\rm{ if }}\;T_{1} < t \leq T_{2}, \\ _{T_{2}}^{C}D_{t}^{\alpha _{3}}y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;T_{2} < t \leq T \end{array} \right. \end{equation}$

(4.5)

which leads only to a power-law process or the Mittag-Leffler process.

5. Application to real world problems

In this section, we shall present clearly the difference between the two concepts by applying them to some problems. We shall start with a simple problem about the decay problem with two intervals.

Example 1. We consider a simple decay model within $\left[ 0, T_{1} \right]$ and $\left[ T_{1}, T_{2}\right].$ In the case of short memory, we have

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha _{1}}y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ _{T_{1}}^{C}D_{t}^{\alpha _{2}}y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right]. \end{array} \right. \end{equation}$

(5.1)

In the case of piecewise, we can consider the following

$\begin{equation} \left\{ \begin{array}{c} \frac{dy\left( t\right) }{dt} = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ _{T_{1}}^{C}D_{t}^{\alpha }y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. \end{equation}$

(5.2)

$\begin{equation} \left\{ \begin{array}{c} _{0}^{ABC}D_{t}^{\alpha }y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ \frac{dy\left( t\right) }{dt} = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. \end{equation}$

(5.3)

$\begin{equation} \left\{ \begin{array}{c} _{0}^{CF}D_{t}^{\alpha }y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ _{T_{1}}^{C}D_{t}^{\alpha }y\left( t\right) = -\lambda y\left( t\right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right]. \end{array} \right. \end{equation}$

(5.4)

Using the Laplace transform on the above system yields in the case of short memory

$\begin{equation} \left\{ \begin{array}{c} y\left( t\right) = y\left( 0\right) E_{\alpha _{1}}\left[ -\lambda t^{\alpha _{1}}\right] \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ y\left( t\right) = y\left( T_{1}\right) E_{\alpha _{2}}\left[ -\lambda \left( t-T_{1}\right) ^{\alpha _{2}}\right] \;{\rm{ if }}\;t\in \left[ T_{1},T_{2} \right]. \end{array} \right. \end{equation}$

(5.5)

In the case of piecewise, we have the following solution

$\begin{equation} \left\{ \begin{array}{c} y\left( t\right) = y\left( 0\right) \exp \left[ -\lambda t\right] \;{\rm{ if }}\; t\in \left[ 0,T_{1}\right] , \\ y\left( t\right) = y\left( T_{1}\right) E_{\alpha }\left[ -\lambda \left( t-T_{1}\right) ^{\alpha }\right] \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. \end{equation}$

(5.6)

$\begin{equation} \left\{ \begin{array}{c} y\left( t\right) = y\left( 0\right) \;{\rm{ if }}\;t = 0, \\ y\left( t\right) = \frac{y\left( 0\right) }{1+\lambda \left( 1-\alpha \right) }E_{\alpha }\left[ -\frac{\lambda \alpha }{1+\lambda \left( 1-\alpha \right) }t^{\alpha }\right] b\;{\rm{ if }}\;t\in \left] 0,T_{1}\right] , \\ y\left( t\right) = y\left( T_{1}\right) \exp \left[ -\lambda \left( t-T_{1}\right) \right] \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. \end{equation}$

(5.7)

$\begin{equation} \left\{ \begin{array}{c} y\left( t\right) = y\left( 0\right) \;{\rm{ if }}\;t = 0, \\ y\left( t\right) = \frac{y\left( 0\right) }{1+\lambda \left( 1-\alpha \right) }\exp \left[ -\frac{\lambda \alpha }{1+\lambda \left( 1-\alpha \right) } t^{\alpha }\right] b\;{\rm{ if }}\;t\in \left] 0,T_{1}\right] , \\ y\left( t\right) = y\left( T_{1}\right) E_{\alpha }\left[ -\lambda \left( t-T_{1}\right) ^{\alpha }\right] \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right]. \end{array} \right. \end{equation}$

(5.8)

$T_{1} = 1$ where the crossover occurs and $T_{2} = 2.4$ is the final time. The numerical simulations are presented in Figures 1–4 below.

Figure 1. Decay model with short memory for

$\lambda = 3, \alpha _{1} = 0.9$ and

$\alpha _{2} = 0.6$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. Decay model with piecewise derivative for

$\lambda = 3$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. Decay model with piecewise derivative for

$\lambda = 3$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. Decay model with piecewise derivative for

$\lambda = 3$ .

DownLoad: Full-Size Img PowerPoint

Example 2. We consider the spread of an infectious disease. We consider a Susceptible, Exposed, Infected and Recovered (SEIR) model. It was proven by Atangana and Seda that such a model may not be able to predict waves ^[12]. Thus to introduce waves, in particular two waves, we consider the spread to be in the periods $\left[ 0, T_{1} \right]$ and $\left[ T_{1}, T_{2}\right].$ In the case of short memory, we have

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }S = \mu N-\mu S-\frac{\beta SI}{N} \\ _{0}^{C}D_{t}^{\alpha }E = \frac{\beta SI}{N}-\left( \mu +a\right) E \\ _{0}^{C}D_{t}^{\alpha }I = aE-\left( \gamma +\mu \right) I \\ _{0}^{C}D_{t}^{\alpha }R = \gamma I-\mu R \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{T_{1}}^{C}D_{t}^{\beta }S = \mu N-\mu S-\frac{\beta SI}{N} \\ _{T_{1}}^{C}D_{t}^{\beta }E = \frac{\beta SI}{N}-\left( \mu +a\right) E \\ _{T_{1}}^{C}D_{t}^{\beta }I = aE-\left( \gamma +\mu \right) I \\ _{T_{1}}^{C}D_{t}^{\beta }R = \gamma I-\mu R \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.9)

In the case of piecewise several scenarios can be obtained. So we can have

$\begin{eqnarray} \left\{ \begin{array}{c} \overset{\cdot }{S} = \mu N-\mu S-\frac{\beta SI}{N} \\ \overset{\cdot }{E} = \frac{\beta SI}{N}-\left( \mu +a\right) E \\ \overset{\cdot }{I} = aE-\left( \gamma +\mu \right) I \\ \overset{\cdot }{R} = \gamma I-\mu R \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} dS = \left( \mu N-\mu S-\frac{\beta SI}{N}\right) dt+\sigma _{1}SdB_{1}\left( t\right) \\ dE = \left( \frac{\beta SI}{N}-\left( \mu +a\right) E\right) dt+\sigma _{2}EdB_{2}\left( t\right) \\ dI = \left( aE-\left( \gamma +\mu \right) I\right) dt+\sigma _{3}IdB_{3}\left( t\right) \\ dR = \left( \gamma I-\mu R\right) dt+\sigma _{4}RdB_{4}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] \end{eqnarray}$

(5.10)

$\begin{eqnarray} \left\{ \begin{array}{c} dS = \left( \mu N-\mu S-\frac{\beta SI}{N}\right) dt+\sigma _{1}SdB_{1}\left( t\right) \\ dE = \left( \frac{\beta SI}{N}-\left( \mu +a\right) E\right) dt+\sigma _{2}EdB_{2}\left( t\right) \\ dI = \left( aE-\left( \gamma +\mu \right) I\right) dt+\sigma _{3}IdB_{3}\left( t\right) \\ dR = \left( \gamma I-\mu R\right) dt+\sigma _{4}RdB_{4}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{T_{1}}^{C}D_{t}^{\alpha }S = \mu N-\mu S-\frac{\beta SI}{N} \\ _{T_{1}}^{C}D_{t}^{\alpha }E = \frac{\beta SI}{N}-\left( \mu +a\right) E \\ _{T_{1}}^{C}D_{t}^{\alpha }I = aE-\left( \gamma +\mu \right) I \\ _{T_{1}}^{C}D_{t}^{\alpha }R = \gamma I-\mu R \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.11)

Noting that several more scenarios can be considered, however, we will only consider these two in the case of piecewise. Before proceeding with the analysis of these models, we shall first provide an interpolation of each case in terms of wave behaviors. In the case of short memory, the first and second waves show the non-Gaussian distribution associated to the power-law variables $t^{-\alpha _{1}}$ and $t^{-\alpha _{2}}$ . Therefore, one will observe different long tails for $t^{-\alpha _{1}}$ and $t^{-\alpha _{2}}$ . On the other hand, however, the first considered model shows that the first wave has a normal distribution, while the second wave has random behavior. The second model shows that the first wave has a long tail spread while the second has random behavior.

We now present the numerical solutions of the considered models. We shall consider the numerical scheme suggested by Ghanbari et al. ^[13] for fractional cases. The short memory case can be simplified to

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }y_{i}\left( t\right) = F_{i}\left( t,y_{i}\left( t\right) \right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ _{T_{1}}^{C}D_{t}^{\beta }y_{i}\left( t\right) = F_{i}\left( t,y_{i}\left( t\right) \right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right]. \end{array} \right. \end{equation}$

(5.12)

Applying the Ghanbari method yields

$\begin{equation} \left\{ \begin{array}{c} y_{i}^{n_{1}} = y_{i}\left( 0\right) +h^{\alpha }\left( \begin{array}{c} \widetilde{a}_{n_{1}}^{\left( \alpha \right) }F_{i}\left( t_{0},y_{i}\left( t_{0}\right) \right) \\ +\sum\limits_{k = 2}^{n_{1}}\widetilde{a}_{n_{1}-k}^{\left( \alpha \right) }F_{i}\left( t_{k},y_{i}\left( t_{k}\right) \right) \end{array} \right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ y_{i}^{n_{2}} = y_{i}\left( T_{1}\right) +h^{\beta }\left( \begin{array}{c} \widetilde{a}_{n_{2}}^{\left( \beta \right) }F_{i}\left( T_{1},y_{i}\left( T_{1}\right) \right) \\ +\sum\limits_{j = n_{1}+2}^{n_{2}}a_{n_{2}-j}^{\left( \beta \right) }F_{i}\left( t_{j},y_{i}\left( t_{j}\right) \right) \end{array} \right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. \end{equation}$

(5.13)

where

$\begin{eqnarray} \widetilde{a}_{n_{1}}^{\left( \alpha \right) } & = &\frac{\left( n_{1}-1\right) ^{\alpha +1}-n_{1}^{\alpha }\left( n_{1}-\alpha -1\right) }{ \Gamma \left( \alpha +2\right) }, \\ a_{n_{2}}^{\left( \beta \right) } & = &\left\{ \begin{array}{c} \frac{1}{\Gamma \left( \beta +2\right) }\;{\rm{ if }}\;n = 0, \\ \frac{\left( n_{2}-1\right) ^{\beta +1}-2n_{2}^{\beta +1}+\left( n_{2}+1\right) ^{\beta +1}}{\Gamma \left( \beta +2\right) }\;{\rm{ if }}\;n\geq 1. \end{array} \right. \end{eqnarray}$

(5.14)

and

$\begin{eqnarray} \widetilde{a}_{n_{2}}^{\left( \beta \right) } &=&\frac{\left( n_{2}-1\right) ^{\beta +1}-n_{2}^{\beta }\left( n_{2}-\beta -1\right) }{\Gamma \left( \beta +2\right) }, \\ \theta _{k}^{\left( \beta \right) } &=&\frac{\left( k-1\right) ^{\beta +1}-2k^{\beta +1}+\left( k+1\right) ^{\beta +1}}{\Gamma \left( \beta +2\right) }\text{ if }k=n_{1}+1,,..,n_{2}-1. \notag \end{eqnarray}$

(5.15)

In the case of the piecewise first model, we have

$\left\{ \begin{matrix} y_{i}^{{{n}_{1}}}={{y}_{i}}\left( 0 \right)+\sum_{j=2}^{{{n}_{1}}}{\left[ \begin{matrix} \frac{23}{12}{{F}_{i}}\left( {{t}_{j}},{{y}_{i}}\left( {{t}_{j}} \right) \right)\Delta t \\ -\frac{4}{3}{{F}_{i}}\left( {{t}_{j-1}},{{y}_{i}}\left( {{t}_{j-1}} \right) \right)\Delta t \\ +\frac{5}{12}{{F}_{i}}\left( {{t}_{j-2}},{{y}_{i}}\left( {{t}_{j-2}} \right) \right)\Delta t \\ \end{matrix} \right]}\ \text{if}\ t\in \left[ 0,{{T}_{1}} \right], \\ \left\{ \begin{matrix} y_{i}^{{{n}_{2}}}=y\left( {{T}_{1}} \right)+\sum_{j={{n}_{1}}+3}^{{{n}_{2}}}{\left[ \begin{matrix} \frac{23}{12}{{F}_{i}}\left( {{t}_{j}},{{y}_{i}}\left( {{t}_{j}} \right) \right)\Delta t \\ -\frac{4}{3}{{F}_{i}}\left( {{t}_{j-1}},{{y}_{i}}\left( {{t}_{j-1}} \right) \right)\Delta t \\ +\frac{5}{12}{{F}_{i}}\left( {{t}_{j-2}},{{y}_{i}}\left( {{t}_{j-2}} \right) \right)\Delta t \\ \end{matrix} \right]}\ \text{if}\ t\in \left[ {{T}_{1}},{{T}_{2}} \right]. \\ +{{\sigma }_{i}}{{y}_{i}}\left( B_{i}^{{{n}_{2}}+1}-B_{i}^{{{n}_{2}}} \right) \\ \end{matrix} \right. \\ \end{matrix} \right.$

(5.16)

For Case 2, we obtain

$\begin{equation} \left\{ \begin{array}{c} \left\{ \begin{array}{c} y_{i}^{n_{1}} = y_{i}\left( 0\right) +\sum_{j = 2}^{n_{1}}\left[ \begin{array}{c} \frac{23}{12}F_{i}\left( t_{j},y_{i}\left( t_{j}\right) \right) \Delta t \\ -\frac{4}{3}F_{i}\left( t_{j-1},y_{i}\left( t_{j-1}\right) \right) \Delta t \\ +\frac{5}{12}F_{i}\left( t_{j-2},y_{i}\left( t_{j-2}\right) \right) \Delta t \end{array} \right] {\rm{ }} \\ +\sigma _{i}\sum_{j = 1}^{n_{1}}y_{j}\left( B_{i}^{j}-B_{i}^{j-1}\right) \end{array} \right. \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ y_{i}^{n_{2}} = y_{i}\left( T_{1}\right) +h^{\alpha }\left( \begin{array}{c} \widetilde{a}_{n_{1}}^{\left( \alpha \right) }F_{i}\left( t_{n_{1}},y_{i}\left( t_{n_{1}}\right) \right) \\ +\sum_{k = 0}^{n_{2}}\widetilde{a}_{n_{1}+k}^{\left( \alpha \right) }F_{i}\left( t_{k},y_{i}\left( t_{k}\right) \right) \end{array} \right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right]. \end{array} \right. \end{equation}$

(5.17)

The numerical solutions for short memory are depicted in Figure 5 with the following parameters

$\begin{equation} N = 1000,\mu = 0.01,\beta = 0.6,a = 0.2,\gamma = 0.03, \end{equation}$

(5.18)

Figure 5. Numerical solutions of Covid-19 model with short memory for

$\alpha _{1} = 0.95, \alpha _{2} = 0.75$ .

DownLoad: Full-Size Img PowerPoint

and the following initial conditions

$\begin{eqnarray} S\left( 0\right) & = &1000,E\left( 0\right) = 25,I\left( 50\right) = 22,R\left( 50\right) = 5,S\left( 50\right) = 39 \\ E\left( 50\right) & = &67,I\left( 50\right) = 505,R\left( 50\right) = 390 \end{eqnarray}$

(5.19)

The numerical solutions for Cases 1 and 2 are respectively depicted in Figures 6 and 7 with the following parameters

$\begin{equation} N = 1000,\mu = 0.01,\beta = 0.6,a = 0.2,\gamma = 0.03, \end{equation}$

(5.20)

Figure 6. Numerical solutions of Covid-19 model under the conditions of Case 1

$\sigma _{1} = 0.000005, \sigma _{2} = 0.008, \sigma _{3} = 0.0021, \sigma _{4} = 0.009$ and

$\alpha = 0.7$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. Numerical solutions of Covid-19 model under the conditions of Case 2

$\sigma _{1} = 0.17, \sigma _{2} = 0.25, \sigma _{3} = 0.21, \sigma _{4} = 0.19$ and

$\alpha = 0.7$ .

DownLoad: Full-Size Img PowerPoint

and the following initial conditions

$\begin{eqnarray} S\left( 0\right) & = &1000,E\left( 0\right) = 25,I\left( 50\right) = 22,R\left( 50\right) = 5,S\left( 50\right) = 35 \\ E\left( 50\right) & = &44,I\left( 50\right) = 424,R\left( 50\right) = 527,\alpha = 0.7 \end{eqnarray}$

(5.21)

demonstrates the short memory effect, where the first and second parts present different long tails behaviors characterized by $\alpha _{1}$ and $\alpha _{2}$ ; however, both clearly depict the same process which is the power law process. On the other hand, Figures 6 and 7 show the clear effects of the piecewise differential operators; the first part in Figure 6 shows normal distribution, which is characteristic of the classical differential operator, while the second part clearly shows the effect of randomness. The complete model therefore shows crossover behavior from deterministic to stochastic with no steady state. In Figure 6, the process is different. We start with randomness then end up with long tail behaviors as a result of the power law kernel.

Example 3. We now consider a chaotic problem in particular the 8-wing attractor ^[14]. This model has been studied by several researchers ^[15], the model is given as

$\begin{eqnarray} \overset{\cdot }{x} & = &a\left( y-x\right) +f\left( t\right) yz \\ \overset{\cdot }{y} & = &cx+dy-xz \\ \overset{\cdot }{z} & = &-bz+xy \end{eqnarray}$

(5.22)

where the function $f\left(t\right) = M$ $sgn\left(\sin \left(wt\right) \right) +K.$ The constant $w$ is known to be a switch frequency and $M$ and $K$ are constant parameters.

In the case of short memory, we have

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha _{1}}x = a\left( y-x\right) +f\left( t\right) yz \\ _{0}^{C}D_{t}^{\alpha _{1}}y = cx+dy-xz \\ _{0}^{C}D_{t}^{\alpha _{1}}z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{T_{1}}^{C}D_{t}^{\alpha _{2}}x = a\left( y-x\right) +f\left( t\right) yz \\ _{T_{1}}^{C}D_{t}^{\alpha _{2}}y = cx+dy-xz \\ _{T_{1}}^{C}D_{t}^{\alpha _{2}}z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.23)

In the case of piecewise, we can consider three cases. For Case 1, we can write

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{ABC}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{0}^{ABC}D_{t}^{\alpha }y = cx+dy-xz \\ _{0}^{ABC}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} dx = \left( a\left( y-x\right) +f\left( t\right) yz\right) dt+\sigma _{1}xdB_{1}\left( t\right) \\ dy = \left( cx+dy-xz\right) dt+\sigma _{2}ydB_{2}\left( t\right) \\ dz = \left( -bz+xy\right) dt+\sigma _{3}zdB_{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.24)

For Case 2, we have

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{0}^{C}D_{t}^{\alpha }y = cx+dy-xz \\ _{0}^{C}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} dx = \left( a\left( y-x\right) +f\left( t\right) yz\right) dt+\sigma _{1}xdB_{1}\left( t\right) \\ dy = \left( cx+dy-xz\right) dt+\sigma _{2}ydB_{2}\left( t\right) \\ dz = \left( -bz+xy\right) dt+\sigma _{3}zdB_{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.25)

For Case 3, we can get

$\begin{eqnarray} \left\{ \begin{array}{c} \overset{\cdot }{x} = a\left( y-x\right) +f\left( t\right) yz \\ \overset{\cdot }{y} = cx+dy-xz \\ \overset{\cdot }{z} = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} dx = \left( a\left( y-x\right) +f\left( t\right) yz\right) dt+\sigma _{1}xdB_{1}\left( t\right) \\ dy = \left( cx+dy-xz\right) dt+\sigma _{2}ydB_{2}\left( t\right) \\ dz = \left( -bz+xy\right) dt+\sigma _{3}zdB_{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] . \end{eqnarray}$

(5.26)

We shall now present the numerical solution of each case. In the case of short memory, one can use the numerical method based on Lagrange interpolation ^[16]. We can first simplify the model as follows.

$\begin{equation} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha _{1}}M_{i}\left( t\right) = \Pi \left( t,M_{i}\left( t\right) \right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ _{T_{1}}^{C}D_{t}^{\alpha _{2}}M_{i}\left( t\right) = \Pi \left( t,M_{i}\left( t\right) \right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. {\rm{ }}. \end{equation}$

(5.27)

Thus, applying such method to the above system yields

$\begin{equation} \left\{ \begin{array}{c} M_{i}^{n_{1}} = M_{i}\left( 0\right) +\left[ \begin{array}{c} \frac{\left( \Delta t\right) ^{\alpha _{1}}}{\Gamma \left( \alpha _{1}+2\right) }\sum\limits_{q = 0}^{n_{1}}\Pi \left( t_{q},M_{i}^{q}\right) \\ \times \left[ \begin{array}{c} \left( n_{1}-q+1\right) ^{\alpha _{1}}\left( n_{1}-q+2+\alpha _{1}\right) \\ -\left( n_{1}-q\right) ^{\alpha _{1}}\left( n_{1}-q+2+2\alpha _{1}\right) \end{array} \right] \\ -\frac{\left( \Delta t\right) ^{\alpha _{1}}}{\Gamma \left( \alpha _{1}+2\right) }\sum\limits_{q = 1}^{n_{1}}\Pi \left( t_{q-1},M_{i}^{q-1}\right) \\ \times \left[ \begin{array}{c} \left( n_{1}-q+1\right) ^{\alpha _{1}+1} \\ -\left( n_{1}-q\right) ^{\alpha _{1}}\left( n_{1}-q+1+\alpha _{1}\right) \end{array} \right] \end{array} \right] \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ M_{i}^{n_{2}} = M_{i}\left( T_{1}\right) +\left[ \begin{array}{c} \frac{\left( \Delta t\right) ^{\alpha _{2}}}{\Gamma \left( \alpha _{2}+2\right) }\sum\limits_{k = n_{1}+1}^{n_{2}}\Pi \left( t_{k},M_{i}^{k}\right) \\ \times \left[ \begin{array}{c} \left( n_{2}-q+1\right) ^{\alpha _{2}}\left( n_{2}-q+2+\alpha _{2}\right) \\ -\left( n_{2}-q\right) ^{\alpha _{2}}\left( n_{2}-q+2+2\alpha _{2}\right) \end{array} \right] \\ -\frac{\left( \Delta t\right) ^{\alpha _{2}}}{\Gamma \left( \alpha _{2}+2\right) }\sum\limits_{k = n_{1}+2}^{n_{2}}\Pi \left( t_{k-1},M_{i}^{k-1}\right) \\ \times \left[ \begin{array}{c} \left( n_{2}-q+1\right) ^{\alpha _{2}+1} \\ -\left( n_{2}-q\right) ^{\alpha _{2}}\left( n_{2}-q+1+\alpha _{2}\right) \end{array} \right] \end{array} \right] \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. {\rm{ }}. \end{equation}$

(5.28)

With the piecewise case, we obtain the following numerical solution by using a scheme based on Newton polynomial interpolation ^[17]. Also, we simplify the model as follows.

$\begin{equation} \left\{ \begin{array}{c} _{0}^{ABC}D_{t}^{\alpha _{1}}M_{i}\left( t\right) = \Pi \left( t,M_{i}\left( t\right) \right) \;{\rm{ if }}\;t\in \left[ 0,T_{1}\right] , \\ dM_{i}\left( t\right) = \Pi \left( t,M_{i}\left( t\right) \right) dt+\sigma _{i}M_{i}dB_{i}\left( t\right) \;{\rm{ if }}\;t\in \left[ T_{1},T_{2}\right] \end{array} \right. {\rm{ }}. \end{equation}$

(5.29)

Then applying this method yields

$\begin{equation} \begin{array}{c} M_{i}^{n_{1}} = M_{i}\left( 0\right) +\frac{1-\alpha }{AB\left( \alpha \right) }\Pi \left( t_{n_{1}},M_{i}^{n_{1}}\right) +\frac{\alpha \left( \Delta t\right) ^{\alpha }}{AB\left( \alpha \right) \Gamma \left( \alpha +1\right) } \sum\limits_{k = 2}^{n_{1}}\Pi \left( t_{k-2},M_{i}^{k-2}\right) \\ \times \left\{ \left( n_{1}-k+1\right) ^{\alpha }-\left( n_{1}-k\right) ^{\alpha }\right\} \\ +\frac{\alpha \left( \Delta t\right) ^{\alpha }}{AB\left( \alpha \right) \Gamma \left( \alpha +2\right) }\sum\limits_{k = 2}^{n_{1}}\left[ \Pi \left( t_{k-1},M_{i}^{k-1}\right) -\Pi \left( t_{k-2},M_{i}^{k-2}\right) \right] \\ \times \left\{ \begin{array}{c} \left( n_{1}-k+1\right) ^{\alpha }\left( n_{1}-k+3+2\alpha \right) \\ -\left( n_{1}-k\right) ^{\alpha }\left( n_{1}-k+3+3\alpha \right) \end{array} \right\} \\ +\frac{\alpha \left( \Delta t\right) ^{\alpha }}{2AB\left( \alpha \right) \Gamma \left( \alpha +3\right) }\sum\limits_{k = 2}^{n_{1}}\left[ \Pi \left( t_{k},M_{i}^{k}\right) -2\Pi \left( t_{k-1},M_{i}^{k-1}\right) +\Pi \left( t_{k-2},M_{i}^{k-2}\right) \right] \\ \times \left\{ \begin{array}{c} \left( n_{1}-k+1\right) ^{\alpha }\left[ \begin{array}{c} 2\left( n_{1}-k\right) ^{2}+\left( 3\alpha +10\right) \left( n_{1}-k\right) \\ +2\alpha ^{2}+9\alpha +12 \end{array} \right] \\ -\left( n_{1}-k\right) ^{\alpha }\left[ \begin{array}{c} 2\left( n_{1}-k\right) ^{2}+\left( 5\alpha +10\right) \left( n_{1}-k\right) \\ +6\alpha ^{2}+18\alpha +12 \end{array} \right] \end{array} \right\} ,0\leq t\leq T_{1} \\ \\ M_{i}^{n_{2}+4} = M_{i}^{n_{2}+3}+\sum\limits_{j_{3} = n_{1}+3}^{n_{3}}\left[ \begin{array}{c} \frac{23}{12}\Pi \left( t_{n_{2}+3},M_{i}^{n_{2}+3}\right) -\frac{4}{3}\Pi \left( t_{n_{2}+2},M_{i}^{n_{2}+2}\right) \\ +\frac{5}{12}\Pi \left( t_{n_{2}+1},M_{i}^{n_{2}+1}\right) \end{array} \right] \Delta t \\ +\sigma _{i}M_{i}^{n_{2}}\left( B_{i}^{n_{2}+1}-B_{i}^{n_{2}+1}\right) ,T_{1}\leq t\leq T_{2} \end{array} \end{equation}$

(5.30)

For Case 2, we again apply the scheme with the Newton polynomial to obtain

$\begin{equation} \left\{ \begin{array}{c} M_{i}^{n_{1}} = M_{i}\left( 0\right) +\frac{\left( \Delta t\right) ^{\alpha }}{ \Gamma \left( \alpha +1\right) }\sum\limits_{k = 2}^{n_{1}}\Pi \left( t_{k-2},M_{i}^{k-2}\right) \\ \times \left\{ \left( n_{1}-k+1\right) ^{\alpha }-\left( n_{1}-k\right) ^{\alpha }\right\} \\ +\frac{\left( \Delta t\right) ^{\alpha }}{\Gamma \left( \alpha +2\right) } \sum\limits_{k = 2}^{n_{1}}\left[ \Pi \left( t_{k-1},M_{i}^{k-1}\right) -\Pi \left( t_{k-2},M_{i}^{k-2}\right) \right] , \\ \times \left\{ \begin{array}{c} \left( n_{1}-k+1\right) ^{\alpha }\left( n_{1}-k+3+2\alpha \right) \\ -\left( n_{1}-k\right) ^{\alpha }\left( n_{1}-k+3+3\alpha \right) \end{array} \right\} \\ +\frac{\left( \Delta t\right) ^{\alpha }}{2\Gamma \left( \alpha +3\right) } \sum\limits_{k = 2}^{n_{1}}\left[ \Pi \left( t_{k},M_{i}^{k}\right) -2\Pi \left( t_{k-1},M_{i}^{k-1}\right) +\Pi \left( t_{k-2},M_{i}^{k-2}\right) \right] \\ \times \left\{ \begin{array}{c} \left( n_{1}-k+1\right) ^{\alpha }\left[ \begin{array}{c} 2\left( n_{1}-k\right) ^{2}+\left( 3\alpha +10\right) \left( n_{1}-k\right) \\ +2\alpha ^{2}+9\alpha +12 \end{array} \right] \\ -\left( n_{1}-k\right) ^{\alpha }\left[ \begin{array}{c} 2\left( n_{1}-k\right) ^{2}+\left( 5\alpha +10\right) \left( n_{1}-k\right) \\ +6\alpha ^{2}+18\alpha +12 \end{array} \right] \end{array} \right\} ,0\leq t\leq T_{1} \\ \\ M_{i}^{n_{2}+4} = M_{i}^{n_{2}+3}+\sum\limits_{j_{3} = n_{2}+3}^{n_{3}}\left[ \begin{array}{c} \frac{23}{12}\Pi \left( t_{n_{2}+3},M_{i}^{n_{2}+3}\right) -\frac{4}{3}\Pi \left( t_{n_{2}+2},M_{i}^{n_{2}+2}\right) \\ +\frac{5}{12}\Pi \left( t_{n_{2}+1},M_{i}^{n_{2}+1}\right) \end{array} \right] \Delta t \\ +\sigma _{i}M_{i}\left( C_{n_{2}}\right) \left( B_{i}^{n_{2}+1}-B_{i}^{n_{2}+1}\right) ,T_{1}\leq t\leq T_{2} \end{array} \right. \end{equation}$

(5.31)

where $C_{n_{2}}\in \left[ t_{n_{2}+4}, t_{n_{2}+3}\right]$ . For Case 3, we use the simple Adams-Bashforth methods as follows

$\begin{equation} \left\{ \begin{array}{c} M_{i}^{n_{2}+3} = M_{i}^{n_{2}+2}+\left[ \begin{array}{c} \frac{23}{12}\Pi \left( t_{n_{2}+2},M_{i}^{n_{2}+2}\right) -\frac{4}{3}\Pi \left( t_{n_{2}+1},M_{i}^{n_{2}+1}\right) \\ +\frac{5}{12}\Pi \left( t_{n_{2}},M_{i}^{n_{2}}\right) \end{array} \right] \Delta t,0\leq t\leq T_{1} \\ \\ M_{i}^{n_{2}+4} = M_{i}^{n_{2}+3}+\left[ \begin{array}{c} \frac{23}{12}\Pi \left( t_{n_{2}+3},M_{i}^{n_{2}+3}\right) -\frac{4}{3}\Pi \left( t_{n_{2}+2},M_{i}^{n_{2}+2}\right) \\ +\frac{5}{12}\Pi \left( t_{n_{2}+1},M_{i}^{n_{2}+1}\right) \end{array} \right] \Delta t \\ +\sigma _{i}M_{i}\left( C_{n_{2}}\right) \left( B_{i}^{n_{2}+1}-B_{i}^{n_{2}}\right) ,T_{1}\leq t\leq T_{2} \end{array} \right. \end{equation}$

(5.32)

where $C_{n_{2}}\in \left[ t_{n_{2}+4}, t_{n_{2}+3}\right].$

Numerical solutions are obtained for the following parameters

$\begin{equation} a = 14,b = 43,c = -1,d = 16,M = 9,K = 10,w = \frac{2\pi }{50}. \end{equation}$

(5.33)

The numerical solutions for short memory are depicted in for the initial data $x\left(0\right) = 22.3, y\left(0\right) = 4.9$ and $z\left(0\right) = 2.6.$

Figure 8. Numerical solutions for 8-wing attractor with short memory for

$\alpha _{1} = 0.96$ and

$\alpha _{2} = 0.98$ .

DownLoad: Full-Size Img PowerPoint

For piecewise, we consider the following stochastic-deterministic problem with power-law and Mittag-Leffler kernels:

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz+\sigma _{1}xB\prime _{1}\left( t\right) \\ _{0}^{C}D_{t}^{\alpha }y = cx+dy-xz+\sigma _{2}yB\prime _{2}\left( t\right) \\ _{0}^{C}D_{t}^{\alpha }z = -bz+xy+\sigma _{3}zB\prime _{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }y = cx+dy-xz \\ _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] \end{eqnarray}$

(5.34)

with the stochastic constants

$\begin{equation} \sigma _{1} = 0.031;\sigma _{2} = 0.035;\sigma _{3} = 0.032. \end{equation}$

(5.35)

The numerical simulation for such a problem can be performed as illustrated in Figure 9 by using the same initial data and parameters.

Figure 9. Stochastic-deterministic 8-wing attractor with piecewise derivative for

$\alpha = 0.99$ .

DownLoad: Full-Size Img PowerPoint

For piecewise, we consider the following stochastic-deterministic problem with power-law and Mittag-Leffler kernels:

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{ABC}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz+\sigma _{1}xB\prime _{1}\left( t\right) \\ _{0}^{ABC}D_{t}^{\alpha }y = cx+dy-xz+\sigma _{2}yB\prime _{2}\left( t\right) \\ _{0}^{ABC}D_{t}^{\alpha }z = -bz+xy+\sigma _{3}zB\prime _{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{{{T}_{1}}}^{C}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{{{T}_{1}}}^{C}D_{t}^{\alpha }y = cx+dy-xz \\ _{{{T}_{1}}}^{C}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] \end{eqnarray}$

(5.36)

with the following stochastic constants

$\begin{equation} \sigma _{1} = 0.031;\sigma _{2} = 0.035;\sigma _{3} = 0.032. \end{equation}$

(5.37)

The numerical simulation for such a problem can be performed as shown in Figure 10 by using the same initial data and parameters.

Figure 10. Stochastic-deterministic 8-wing attractor with piecewise derivative for

$\alpha = 0.94$ .

DownLoad: Full-Size Img PowerPoint

Also, we consider the following deterministic-stochastic problem with power-law and Mittag-Leffler kernels:

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{C}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{0}^{C}D_{t}^{\alpha }y = cx+dy-xz \\ _{0}^{C}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz+\sigma _{1}xB\prime _{1}\left( t\right) \\ _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }y = cx+dy-xz+\sigma _{2}yB\prime _{2}\left( t\right) \\ _{{{T}_{1}}}^{ABC}D_{t}^{\alpha }z = -bz+xy+\sigma _{3}zB\prime _{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] \end{eqnarray}$

(5.38)

with the following stochastic constants:

$\begin{equation} \sigma _{1} = 0.02;\sigma _{2} = 0.012;\sigma _{3} = 0.021. \end{equation}$

(5.39)

The numerical solutions for the above problem can be depicted as shown in Figure 11 using same initial data and parameters.

Figure 11. Deterministic-stochastic 8-wing attractor with piecewise derivative for

$\alpha = 0.99$ .

DownLoad: Full-Size Img PowerPoint

Figure 7 shows the results obtained by applying the power-law short memory concept. To show the difference within the interval, we have opted to consider the first part of the interval to be in blue and the second to be in red. However, a quick look at the results shows clearly the effect of a power law, no great change from one interval to another is observed since both depict long-tail behaviors as results of the power-law kernel. In this case, there is no crossover behavior, just a repetition of different long tail behaviors. On the other hand, in Figures 8 and 9, there is a clear change in patterns where the first patterns show deterministic behaviors in particular, long-tails due to power law in Figure 8, and the second part shows randomness, there is, therefore, a clear crossover from power law to randomness. While, in Figure 9, we have two crossovers, the first is due to the Mittag-Leffler kernel that shows a change from stretched exponential to power law and the second is stochastic.

Also, we consider deterministic-stochastic problem with power-law and Mittag-Leffler kernel

$\begin{eqnarray} \left\{ \begin{array}{c} _{0}^{ABC}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz \\ _{0}^{ABC}D_{t}^{\alpha }y = cx+dy-xz \\ _{0}^{ABC}D_{t}^{\alpha }z = -bz+xy \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ 0,T_{1}\right] , \\ \left\{ \begin{array}{c} _{{{T}_{1}}}^{C}D_{t}^{\alpha }x = a\left( y-x\right) +f\left( t\right) yz+\sigma _{1}xB\prime _{1}\left( t\right) \\ _{{{T}_{1}}}^{C}D_{t}^{\alpha }y = cx+dy-xz+\sigma _{2}yB\prime _{2}\left( t\right) \\ _{{{T}_{1}}}^{C}D_{t}^{\alpha }z = -bz+xy+\sigma _{3}zB\prime _{3}\left( t\right) \end{array} \right. \;{\rm{ if }}\;t &\in &\left[ T_{1},T_{2}\right] \end{eqnarray}$

(5.40)

with the following stochastic constants:

$\begin{equation} \sigma _{1} = 0.02;\sigma _{2} = 0.012;\sigma _{3} = 0.021. \end{equation}$

(5.41)

The numerical solutions for the above problem can be depicted as shown in Figure 12 by using the same initial data and parameters.

Figure 12. Deterministic-stochastic 8-Wing attractor with piecewise derivative for

$\alpha = 0.9$ .

DownLoad: Full-Size Img PowerPoint

6. Conclusions

In the last decades, researchers have devoted their attention to better understanding complex realworld problems even on a small scale. They have therefore developed several methods, different differential, and integral operators. To understand the process by which different long tails occur in different intervals, the concept of short memory was introduced. This concept defines power-law derivatives in different intervals and each interval has its own order. This order accounts for the long tail associated with that interval. However, the process is scale-invariant in terms of patterns. On the other hand, because there are many real-world problems that exhibit passages from one process to another, for example, a passage from power law to randomness, the concept of piecewise was introduced. In this work, a detailed analysis was given to show their different fields of applications and show also, when short memory and piecewise concepts can be applied. This paper thus helps researchers to identify what problem is suitable for short memory and piecewise.

Conflict of interest

The authors declare there is no conflict of interest.

Acknowledgments

We thank the Ebetsu City Fire Department for their cooperation in the survey. This project received the Emergency Medical Services Foundation grant. The funders had no role in study design, data collection and analysis, the decision to publish, or preparation of the manuscript.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	Kamimura A, Ashby J, Myers K, et al. (2015) Satisfaction with healthcare services among free clinic patients. J Community Health 40: 62-72. https://doi.org/10.1007/s10900-014-9897-8
[2]	Marshall GN, Hays RD (1994) The patient satisfaction questionnaire short-form (PSQ-18). Santa Monica, CA: Rand.
[3]	Fufa BD, Negao EB (2019) Satisfaction of Outpatient Service Consumers and Associated Factors Towards the Health Service Given at Jimma Medical Center, South West Ethiopia. Patient Relat Outcome Meas 10: 347-354. https://doi.org/10.2147/PROM.S220404
[4]	Safran DG, Taira DA, Rogers WH, et al. (1998) Linking primary care performance to outcomes of care. J Fam Pract 47: 213-220.
[5]	Fitzpatrick R, Hopkins A (1993) Measurement of Patients' Care Satisfaction with Their Care.Royal college of physician of London.
[6]	Guldvog B (1999) Can patient satisfaction improve health among patients with angina pectoris?. Int J Qual Health Care 11: 233-240. https://doi.org/10.1093/intqhc/11.3.233
[7]	Rahmqvist M, Bara AC (2010) Patient characteristics and quality dimensions related to patient satisfaction. Int J Qual Health Care 22: 86-92. https://doi.org/10.1093/intqhc/mzq009
[8]	Shah MN (2006) The formation of the emergency medical services system. Am J Public Health 96: 414-423. https://doi.org/10.2105/AJPH.2004.048793
[9]	Tanigawa K, Tanaka K (2006) Emergency medical service systems in Japan: past, present, and future. Resuscitation 69: 365-370. https://doi.org/10.1016/j.resuscitation.2006.04.001
[10]	Rivard MK, Cash RE, Mercer CB, et al. (2021) Demography of the National Emergency Medical Services Workforce: A Description of Those Providing Patient Care in the Prehospital Setting. Prehosp Emerg Care 25: 213-220. https://doi.org/10.1080/10903127.2020.1737282
[11]	Bernard AW, Lindsell CJ, Handel DA, et al. (2007) Postal survey methodology to assess patient satisfaction in a suburban emergency medical services system: an observational study. BMC Emerg Med 7: 5. https://doi.org/10.1186/1471-227X-7-5
[12]	Doering GT (1998) Customer care. Patient satisfaction in the prehospital setting. Emerg Med Serv 27: 71-74.
[13]	Hudak PL, Wright JG (2000) The characteristics of patient satisfaction measures. Spine 25: 3167-3177. https://doi.org/10.1097/00007632-200012150-00012
[14]	Sitzia J (1999) How valid and reliable are patient satisfaction data? An analysis of 195 studies. Int J Qual Health Care 11: 319-328. https://doi.org/10.1093/intqhc/11.4.319
[15]	Davis BA, Bush HA (2003) Patient satisfaction of emergency nursing care in the United States, Slovenia, and Australia. J Nurs Care Qual 18: 267-274. https://doi.org/10.1097/00001786-200310000-00004
[16]	Romero-García M, de la Cueva-Ariza L, Jover-Sancho C, et al. (2013) La percepción del paciente crítico sobre los cuidados enfermeros: una aproximación al concepto de satisfacción. Enferm Intensiva 24: 51-62. https://doi.org/10.1016/j.enfi.2012.09.003
[17]	Anthoine E, Moret L, Regnault A, et al. (2014) Sample size used to validate a scale: a review of publications on newly-developed patient reported outcomes measures. Health Qual Life Outcomes 12: 176. https://doi.org/10.1186/s12955-014-0176-2
[18]	Bilsbury CD, Richman A (2002) A staging approach to measuring patient-centred subjective outcomes. Acta Psychiatr Scand 106: 5-40. https://doi.org/10.1046/j.0001-690x.2002.00052.x
[19]	Davis LL (1992) Instrument review: Getting the most from a panel of experts. Appl Nurs Res 5: 194-197. https://doi.org/10.1016/S0897-1897(05)80008-4
[20]	Polit DF, Beck CT (2006) The content validity index: are you sure you know what's being reported? Critique and recommendations. Res Nurs Health 29: 489-497. https://doi.org/10.1002/nur.20147
[21]	Mokkink LB, Terwee CB, Patrick DL, et al. (2012) COSMIN checklist manual. Amsterdam: University Medical Center.
[22]	Herdman M, Gudex C, Lloyd A, et al. (2011) Development and preliminary testing of the new five-level version of EQ-5D (EQ-5D-5L). Qual Life Res 20: 1727-1736. https://doi.org/10.1007/s11136-011-9903-x
[23]	Janssen MF, Pickard AS, Golicki D, et al. (2013) Measurement properties of the EQ-5D-5L compared to the EQ-5D-3L across eight patient groups: a multi-country study. Qual Life Res 22: 1717-1727. https://doi.org/10.1007/s11136-012-0322-4
[24]	Shiroiwa T, Fukuda T, Ikeda S, et al. (2016) Japanese population norms for preference-based measures: EQ-5D-3L, EQ-5D-5L, and SF-6D. Qual Life Res 25: 707-719. https://doi.org/10.1007/s11136-015-1108-2
[25]	Revicki DA (2004) Patient assessment of treatment satisfaction: methods and practical issues. Gut 53: iv40-iv44. https://doi.org/10.1136/gut.2003.034322
[26]	Vardanjani AE, Moghaddasi J, Rafiei H, et al. (2014) Satisfaction from pre-hospital emergency medical services: A cross sectional study in Shahrekord. J Clin Rheumatol .
[27]	Bahadori M, Ghardashi F, Izadi AR, et al. (2016) Pre-Hospital Emergency in Iran: A Systematic Review. Trauma Mon 21: e31382. https://doi.org/10.5812/traumamon.31382
[28]	Jansson J, Larsson M, Nilsson J (2021) Advanced paramedics and nurses can deliver safe and effective pre-hospital and in-hospital emergency care: An integrative review. Nurs Open 8: 2385-2405. https://doi.org/10.1002/nop2.866
[29]	García-Alfranca F, Puig A, Galup C, et al. (2018) Patient Satisfaction with Pre-Hospital Emergency Services. A Qualitative Study Comparing Professionals' and Patients' Views. Int J Environ Res Public Health 15: 233. https://doi.org/10.3390/ijerph15020233
[30]	Ferketich S (1991) Focus on psychometrics. Aspects of item analysis. Res Nurs Health 14: 165-168. https://doi.org/10.1002/nur.4770140211
[31]	Bland JM, Altman DG (1997) Cronbach's alpha. BMJ 314: 572. https://doi.org/10.1136/bmj.314.7080.572
[32]	Nunnally JC (1994) Psychometric theory 3E. Tata McGraw-hill education .
[33]	Kline P (2014) An easy guide to factor analysis. Routledge: . https://doi.org/10.4324/9781315788135
[34]	Romero-García M, de la Cueva-Ariza L, Benito-Aracil L, et al. (2018) Nursing Intensive-Care Satisfaction Scale [NICSS]: Development and validation of a patient-centred instrument. J Adv Nurs 74: 1423-1435. https://doi.org/10.1111/jan.13546
[35]	Atkinson MJ, Sinha A, Hass SL, et al. (2004) Validation of a general measure of treatment satisfaction, the Treatment Satisfaction Questionnaire for Medication (TSQM), using a national panel study of chronic disease. Health Qual Life Outcomes 2: 12. https://doi.org/10.1186/1477-7525-2-12
[36]	Wiebe S, Wahby S, Lawal OA, et al. (2020) Development and validation of the Epilepsy Surgery Satisfaction Questionnaire (ESSQ-19). Epilepsia 61: 2729-2738. https://doi.org/10.1111/epi.16709
[37]	Will KK, Johnson ML, Lamb G (2019) Team-Based Care and Patient Satisfaction in the Hospital Setting: A Systematic Review. J Patient Cent Res Rev 6: 158-171. https://doi.org/10.17294/2330-0698.1695
[38]	Ng JHY, Luk BHK (2019) Patient satisfaction: Concept analysis in the healthcare context. Patient Educ Couns 102: 790-796. https://doi.org/10.1016/j.pec.2018.11.013
[39]	Moore L (1999) Measuring quality and effectiveness of prehospital EMS. Prehosp Emerg Care 3: 325-331. https://doi.org/10.1080/10903129908958963
[40]	Reu RC, Dirks B, Richter D, et al. (1997) Emergency telephone calls and response actions at an German emergency base. Resuscitation 2: 199. https://doi.org/10.1016/S0300-9572(97)84273-6
[41]	Patterson PD, Weaver MD, Weaver SJ, et al. (2012) Measuring teamwork and conflict among emergency medical technician personnel. Prehosp Emerg Care 16: 98-108. https://doi.org/10.3109/10903127.2011.616260
[42]	Herzberg S, Hansen M, Schoonover A, et al. (2019) Association between measured teamwork and medical errors: an observational study of prehospital care in the USA. BMJ Open 9: e025314. https://doi.org/10.1136/bmjopen-2018-025314
[43]	Patterson PD, Weaver MD, Hostler D (2017) Teams and teamwork in emergency Medical Services. Human factors and ergonomics of prehospital emergency care.CRC Press 95-108. https://doi.org/10.1201/9781315280172-7
[44]	Leiter MP, Harvie P, Frizzell C (1998) The correspondence of patient satisfaction and nurse burnout. Soc Sci Med 47: 1611-1617. https://doi.org/10.1016/s0277-9536(98)00207-x
[45]	Curka PA, Pepe PE, Zachariah BS, et al. (1995) Incidence, source, and nature of complaints received in a large, urban emergency medical services system. Acad Emerg Med 2: 508-512. https://doi.org/10.1111/j.1553-2712.1995.tb03249.x
[46]	Colwell CB, Pons PT, Pi R (2003) Complaints against an EMS system. J Emerg Med 25: 403-408. https://doi.org/10.1016/j.jemermed.2003.02.004
[47]	Kuisma M, Määttä T, Hakala T, et al. (2003) Customer satisfaction measurement in emergency medical services. Acad Emerg Med 10: 812-815. https://doi.org/10.1111/j.1553-2712.2003.tb00080.x
[48]	Boudreaux ED, O'Hea EL (2004) Patient satisfaction in the Emergency Department: a review of the literature and implications for practice. J Emerg Med 26: 13-26. https://doi.org/10.1016/j.jemermed.2003.04.003
[49]	Søvsø MB, Kløjgaard TA, Hansen PA, et al. (2019) Repeated ambulance use is associated with chronic diseases-a population-based historic cohort study of patients' symptoms and diagnoses. Scand J Trauma Resusc Emerg Med 27: 46. https://doi.org/10.1186/s13049-019-0624-4
[50]	Pedersen CK, Stengaard C, Friesgaard K, et al. (2019) Chest pain in the ambulance; prevalence, causes and outcome - a retrospective cohort study. Scand J Trauma Resusc Emerg Med 27: 84. https://doi.org/10.1186/s13049-019-0659-6
[51]	Elmqvist C, Fridlund B, Ekebergh M (2008) More than medical treatment: the patient's first encounter with prehospital emergency care. Int Emerg Nurs 16: 185-192. https://doi.org/10.1016/j.ienj.2008.04.003
[52]	Holmberg M, Forslund K, Wahlberg AC, et al. (2014) To surrender in dependence of another: the relationship with the ambulance clinicians as experienced by patients. Scand J Caring Sci 28: 544-551. https://doi.org/10.1111/scs.12079
[53]	Ahl C, Nyström M (2012) To handle the unexpected–The meaning of caring in pre-hospital emergency care. Int Emerg Nurs 20: 33-41. https://doi.org/10.1016/j.ienj.2011.03.001
[54]	Rantala A, Behm L, Rosén H (2019) Quality Is in the Eye of the Beholder—A Focus Group Study from the Perspective of Ambulance Clinicians, Physicians, and Managers. Healthc Pap 7: 41. https://doi.org/10.3390/healthcare7010041
[55]	Bowman GS (2001) Emotions and illness. J Adv Nurs 34: 256-263. https://doi.org/10.1046/j.1365-2648.2001.01752.x
[56]	Togher FJ, O'Cathain A, Phung VH, et al. (2015) Reassurance as a key outcome valued by emergency ambulance service users: a qualitative interview study. Health Expect 18: 2951-2961. https://doi.org/10.1111/hex.12279
[57]	Persse DE, Jarvis JL, Corpening J, et al. (2004) Customer satisfaction in a large urban fire department emergency medical services system. Acad Emerg Med 11: 106-110. https://doi.org/10.1197/j.aem.2003.07.020
[58]	Moshinskie JF (1997) Reader attitudes toward journal-based continuing education for EMS personnel. Prehosp Disaster Med 12: 258-263. https://doi.org/10.1017/S1049023X00037766
[59]	Becker J, Hugelius K (2021) Driving the ambulance: an essential component of emergency medical services: an integrative review. BMC Emerg Med 21: 160. https://doi.org/10.1186/s12873-021-00554-9
[60]	Batbaatar E, Dorjdagva J, Luvsannyam A, et al. (2017) Determinants of patient satisfaction: a systematic review. Perspect Public Health 137: 89-101. https://doi.org/10.1177/1757913916634136
[61]	Naidu A (2009) Factors affecting patient satisfaction and healthcare quality. Int J Health Care Qual Assur 22: 366-381. https://doi.org/10.1108/09526860910964834
[62]	Aléx J, Karlsson S, Saveman BI (2013) Patients' experiences of cold exposure during ambulance care. Scand J Trauma Resusc Emerg Med 21: 44. https://doi.org/10.1186/1757-7241-21-44
[63]	Ligtenberg JJM, Arnold LG, Stienstra Y, et al. (2005) Quality of interhospital transport of critically ill patients: a prospective audit. Crit Care 9: R446-451. https://doi.org/10.1186/cc3749
[64]	Rantala A, Ekwall A, Forsberg A (2020) Significant others' perceptions of being taken seriously by the Swedish Ambulance Service when the patient is assessed as nonurgent. Scand J Caring Sci 34: 1028-1037. https://doi.org/10.1111/scs.12811
[65]	Doohan I, Saveman BI (2015) Need for compassion in prehospital and emergency care: a qualitative study on bus crash survivors' experiences. Int Emerg Nurs 23: 115-119. https://doi.org/10.1016/j.ienj.2014.08.008
[66]	Wireklint Sundström B, Bremer A, Lindström V, et al. (2019) Caring science research in the ambulance services: an integrative systematic review. Scand J Caring Sci 33: 3-33. https://doi.org/10.1111/scs.12607
[67]	O'Cathain A, Coleman P, Nicholl J (2008) Characteristics of the emergency and urgent care system important to patients: a qualitative study. J Health Serv Res Policy 13: 19-25. https://doi.org/10.1258/jhsrp.2007.007097
[68]	Buller MK, Buller DB (1987) Physicians' communication style and patient satisfaction. J Health Soc Behav 28: 375-388. https://doi.org/10.2307/2136791
[69]	Tyser AR, Gaffney CJ, Zhang C, et al. (2018) The Association of Patient Satisfaction with Pain, Anxiety, and Self-Reported Physical Function. J Bone Joint Surg Am 100: 1811-1818. https://doi.org/10.2106/JBJS.17.00372
[70]	Lozano Calderón SA, Paiva A, Ring D (2008) Patient satisfaction after open carpal tunnel release correlates with depression. J Hand Surg Am 33: 303-307. https://doi.org/10.1016/j.jhsa.2007.11.025
[71]	Kelly AM (2000) Patient satisfaction with pain management does not correlate with initial or discharge VAS pain score, verbal pain rating at discharge, or change in VAS score in the Emergency Department. J Emerg Med 19: 113-116. https://doi.org/10.1016/s0736-4679(00)00219-5
[72]	Shirai K, Matsumoto R, Horii J (2022) Questionnaire survey of the occupants of high-rise residential buildings in Sapporo after the 2018 Hokkaido Eastern Iburi Earthquake. J Asian Archit Build 21: 884-899. https://doi.org/10.1080/13467581.2021.1900856
[73]	Ikeda H, Kubo T, Sasaki T, et al. (2018) Cross-sectional Internet-based survey of Japanese permanent daytime workers' sleep and daily rest periods. J Occup Health 60: 229-235. https://doi.org/10.1539/joh.17-0165-OA
[74]	Takeuchi M, Showa S, Kitazawa K, et al. (2018) Living alone is associated with an increased risk of institutionalization in older men: A follow-up study in Hamanaka Town of Hokkaido, Japan. Geriatr Gerontol Int 18: 867-872. https://doi.org/10.1111/ggi.13267
[75]	Leppink J, Pérez-Fuster P (2017) We need more replication research-A case for test-retest reliability. Perspect Med Educ 6: 158-164. https://doi.org/10.1007/s40037-017-0347-z

This article has been cited by:

1.	Rodrigo Moreira, Larissa Ferreira Rodrigues Moreira, Flávio de Oliveira Silva, Brazilian Selic Rate Forecasting with Deep Neural Networks, 2024, 0927-7099, 10.1007/s10614-024-10597-2
2.	Michele Bufalo, Claudia Ceci, Giuseppe Orlando, Addressing the financial impact of natural disasters in the era of climate change, 2024, 73, 10629408, 102152, 10.1016/j.najef.2024.102152
3.	Guo-Hui Yang, Si-Qi Ma, Xiao-Dong Bian, Jiang-Cheng Li, Pierluigi Vellucci, The roles of liquidity and delay in financial markets based on an optimal forecasting model, 2023, 18, 1932-6203, e0290869, 10.1371/journal.pone.0290869
4.	Giacomo Ascione, Michele Bufalo, Giuseppe Orlando, Modeling volatility of disaster-affected populations: A non-homogeneous geometric-skew Brownian motion approach, 2024, 130, 10075704, 107761, 10.1016/j.cnsns.2023.107761

Reader Comments

Your name:*

Email:*
© 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Public Health

3.1 4.8

Metrics

Article views(2505) PDF downloads(125) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Public Health

Prehospital emergency care patient satisfaction scale [PECPSS] for care provided by emergency medical teams: Scale development and validation

Related Papers:

Abstract

1. Introduction

2. Motivation for piecewise derivatives

3. Short memory concept

4. Analysis and difference

5. Application to real world problems

6. Conclusions

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. Motivation for piecewise derivatives

3. Short memory concept

4. Analysis and difference

5. Application to real world problems

6. Conclusions

Conflict of interest

References

AIMS Public Health

Prehospital emergency care patient satisfaction scale [PECPSS] for care provided by emergency medical teams: Scale development and validation

Related Papers:

Abstract

1. Introduction

2. Motivation for piecewise derivatives

3. Short memory concept

4. Analysis and difference

5. Application to real world problems

6. Conclusions

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. Motivation for piecewise derivatives

3. Short memory concept

4. Analysis and difference

5. Application to real world problems

6. Conclusions

Conflict of interest

References