
==== Front
Heliyon
Heliyon
Heliyon
2405-8440
Elsevier

S2405-8440(24)11927-5
10.1016/j.heliyon.2024.e35896
e35896
Research Article
Modelling the role of tourism in the spread of HIV: A case study from Malaysia
Apenteng Ofosuhene O. ofap@sund.ku.dk
⁎
Rasmussen Philip
Conrady Beate
Department of Veterinary and Animal Sciences, University of Copenhagen, Denmark
⁎ Corresponding author. ofap@sund.ku.dk
08 8 2024
30 8 2024
08 8 2024
10 16 e358961 8 2024
6 8 2024
© 2024 The Authors
2024
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
This study aimed to assess the role of tourism in the spread of human immunodeficiency virus (HIV) using Malaysian epidemiological data on HIV and acquired immunodeficiency syndrome (AIDS) incidence from 1986 to 2004. A population-level mathematical model was developed with the following compartments: the population susceptible to HIV infection, the clinically confirmed HIV-positive population, the population diagnosed with AIDS, and the tourist population. Additionally, newborns infected with HIV were considered. Sensitivity analyses and variations in fixed parameter values were used to explore the effect of changes to various parameter values on HIV incidence in the model. It was determined that variations in the rate of HIV-positive inbound tourist entries and the rate of foreign tourist exits (i.e., the duration of time tourists spent in Malaysia) significantly impacted the predicted incidence of HIV and AIDS in Malaysia. The model's fit to observed HIV and AIDS incidence was evaluated, resulting in adjusted R2 values of 53.3% and 53.2% for HIV and AIDS, respectively. Furthermore, the reproduction number (R0) was also calculated to quantify the stability of the HIV endemicity in Malaysia. The findings suggest that a steady-state level of HIV in Malaysia is achievable based on the low value of R0 = 0.0136, and the disease-free equilibrium was stable from the negative eigenvalues obtained, which is encouraging from a public health perspective. The Partial Rank Correlation Coefficient (PRCC) values between the proportion of newborns born HIV-positive, the rate of Malaysian tourist entries returning home after contracting HIV, and the rate of foreign tourist exits have a significant impact on the R0. The methods provide a framework for epidemiological modelling of HIV spread through transient population groups. The model results suggest that the role of tourism should not be overlooked within the set of available measures to mitigate the spread of HIV.

Keywords

AIDS
Compartment model
HIV
Malaysia
Reproduction number
Tourism
==== Body
pmc1 Introduction

HIV is a dangerous infection that attacks the immune system [1,2]. The virus mostly targets immune system CD4+ T cells [3], which are central to achieving a regulated effective immune response to pathogens [3]. HIV infection begins with little to no symptoms, and in the early months, it is characterized by slight changes to the immune system, until seroconversion occurs to the point where HIV-specific antibodies can be detected [4]. While the outcome of infection varies widely across individuals, the progression of the disease is typically slow [5], often taking several years to move from the primary infection stage to the symptoms of advanced HIV disease, immunosuppression [4], and increased vulnerability to co-infections [6]. Sexual transmission is the main pathway for the spread of HIV [7], with unprotected anal or vaginal sex greatly increasing the risk of infection [[8], [9], [10]]. Breastfeeding and pregnancy can also pass the virus from mother to child [11], with contaminated organ donations, infusions with contaminated blood or blood products, and injections with contaminated needles being additional transmission pathways [[12], [13], [14], [15]]. HIV-positive people are also at risk for persistent infections from other sexually transmitted infections (STIs), such as gonorrhea, chlamydia, herpes, and syphilis, while STIs increase HIV transmission efficiency through both an increased likelihood of transmission and an increased susceptibility to HIV infection [16,17]. Modelling of co-infection between HIV and Hepatitis C virus (HCV) has been investigated to understand how fractional calculus can be used to model infectious diseases [18].

Globally, tourism is among the fastest-growing industries [19] and revenue from this sector is an essential contributor to many economies, including Malaysia [20]. The Malaysian government began to take the expansion of the tourism industry seriously in the 1980s, with the establishment of The Ministry of Culture, Arts and Tourism in 1987, which later evolved into the Ministry of Tourism [21]. The role of tourism in Malaysia's economic well-being was further emphasized with the adoption of the nation's Third Industrial Master Plan (IMP3), which identified tourism services as a non-governmental sub-sector to be targeted for more significant development and promotion [21,22]. Tourism revenue increased from 4% of Gross domestic product (GDP) in the early 1990s to 16% in 2014. The contribution of tourism to the Malaysian economy has steadily increased, with almost 15% of the nation's total employment attributable to tourism-related activities [21,23]. The country saw more than a 140% increase in tourist arrivals during the same period. Since the year 2000, Malaysia's global tourism ranking has improved, rising from 17th place in 2000 to 10th place in 2012 [21,24]. In addition, over 27 million foreign tourists visited Malaysia in 2016, bringing in USD23 billion from tourism-related activities [23]. However, during the COVID period, tourist entries to Malaysia dropped to 4.3 million in 2020 in response to globally imposed restrictions on international travel [25], and the country's earnings from foreign tourism fell to USD3.0 billion that year [26]. However, in 2022, when travel restrictions eased, tourist entries rebounded to over 10 million [27].

Despite its critical economic role and associated benefits, it is vital to consider the potential for infectious disease transmission between tourists and the populations in their destination countries [28]. Several studies have suggested that touristic regions are associated with increased HIV risk among destination populations. For example, Ketshabile et al. (2007) investigated the relationship between tourism and HIV and the final stage of HIV, which is acquired immunodeficiency syndrome (AIDS), in southern Africa [29], and the role of tour guides/operators has been investigated [30]. Padilla et al. (2010) argued that the direction of future studies will be influenced by the change in the environmental perspective on sexual health in Caribbean tourist destinations and the incidence of HIV/AIDS cases [31]. By conducting a qualitative study of policymakers in the Dominican Republic, Guilamo-Ramos et al. (2013) explored the policy environment of prevention efforts regarding tourism-related HIV in the Caribbean [32]. A study by De Matos et al. (2013) suggested that, in Central Brazil, most sexually transmitted diseases have increased due to the sex trade through tourist-aimed prostitution [33]. In Malaysia, up to 70,000 new STIs are reported annually [34,35], and therefore, the sexual behavior of tourists can play a role in the transmission and spread of STIs and HIV, especially when unprotected sexual activity takes place [36,37]. A better understanding of the role of tourists within the transmission dynamics of HIV, particularly in a nation like Malaysia where tourism is a focal point of economic strategy, would benefit both human health and the economy.

It is known that mathematical models can help public health authorities plan for disease preparedness and mitigation measures [[38], [39], [40], [41]]. However, few mathematical models focus on the impact of tourism on HIV transmission. Some studies have investigated the role of tourists in spreading the disease [42,43], while others have used simulation models to determine the effectiveness of different approaches to controlling HIV transmission caused by tourism [[43], [44], [45]]. However, existing mathematical simulations have some limitations. They do not consider that children born to infected parents may be carriers of the disease, the birth rate within the susceptible population is independent of vertical transmission, and the study periods are not long enough. Therefore, this study aims to investigate the role of tourism in the spread of HIV with conderation for these omissions.

2 Materials and methods

2.1 Data

The epidemiological data from Malaysia between 1986 and 2004 was used to develop and validate a compartmental epidemiological model [46] (section 2.2). The data included the incidence of clinically confirmed HIV cases, the incidence of AIDS cases, and the number of tourist entries into Malaysia. The model was developed based on the reported data that three individuals acquired HIV in 1986, which made up the original HIV cases compartment I(0). Additionally, one person developed AIDS due to their HIV infection, making up the AIDS compartment A(0). The duration of data was between 1986 and 2004.The susceptible compartment S(0) had 16,329,396 members, which was the entire domestic population of Malaysia. The HIV compartment I(0) had three cases, and the AIDS compartment A(0) had only one case. In 1986, there were 3,200,000 tourist entries into Malaysia [47].

2.2 Description of the model

The study used a compartment model based on the SIR (susceptible, infectious, recovered/removed) framework (Fig. 1), commonly used in epidemiological modelling studies. The proposed HIV spread model in this study consists of four components: the susceptible population S(t), the clinically confirmed HIV-positive population I(t), and the population diagnosed with AIDS A(t). Assuming that both I(t) and A(t) were sexually active, newborns who were HIV-positive joined the HIV compartment at a rate of b. The proportion of newborns born HIV-positive is represented by ϕ [48], while b is the natural birth rate. Hence, (1−ϕ)b(I+A) denotes the proportion of newborns susceptible to HIV infection. The rate of susceptibility to HIV transmission through sexual contact is denoted by β. It is assumed that the natural death rate μ remains constant across all compartments. The parameter d represents the rate of AIDS-related death, while the parameter α represents the rate of transfer from the HIV-positive compartment to the AIDS-positive compartment.Fig. 1 Schematic model for assessing tourism impact on HIV/AIDS incidence. S = population susceptible to HIV infection;I = population with clinically-confirmed HIV; A = population diagnosed with AIDS; ϕ = the proportion newborns born HIV-positive; b = the natural birthrate; (1−ϕ)b(I+A) = the fraction of newborns susceptible to HIV infection; β = the rate of susceptibility to HIV infection; μ = the natural death rate; d = AIDS-related death rate; α = the rate of population transfer from HIV-positive to AIDS-positive; δ = the rate of inbound tourist entries with susceptibility to HIV infection; δ1 = the rate of HIV-positive inbound tourist entries; δ2 = the rate of AIDS-positive inbound tourist entries; σ = the rate of foreign tourist exits; ρ = the rate of Malaysians returning home after contracting HIV abroad with unconfirmed infection status.

Fig. 1

In a study by Apenteng (2017) [43], it was observed that tourists entered three different compartments at different rates. The compartments were susceptible, HIV-positive, and AIDS-positive. The rates of entry were defined as δ, δ1, and δ2, respectively, where δ is the rate of susceptibility of tourists, δ1 is the rate of HIV of the tourists, and finally, and δ2 is the rate of the AIDS of tourist. Similarly, tourists exited each compartment at a rate of σ, representing the rate of foreign tourist exits. The study also considered the proportion of Malaysian tourists who traveled abroad and contracted HIV before returning home, denoted as ρ. To estimate the population in each compartment, differential equation (1) through (3) below were used to calculate the rates of change.(1) dSdt=δT+ρS+(1−ϕ)b(I+A)−β(I+A)NS−(σ+μ)S

(2) dIdt=δ1T+ϕb(I+A)+β(I+A)NS−(α+σ+μ)I

(3) dAdt=δ2T+αI−(σ+d+μ)A

2.3 Analysis of the model

In this section, we analyze the model and the equilibrium points with regard to the stability, instability of the model and the basic reproduction number.

2.3.1 Equilibrium and stability analysis

We analyze the equilibrium points to find out whether there are changes or constants. To determine this disases free equilibrium, (1), (2), (3) becomes dSdt=dIdt=dAdt=0. For information, see [[49], [50], [51]](4) δT+ρS*+(1−ϕ)b(I*+A*)−β(I*+A*)NS*−(σ+μ)S*=0

(5) δ1T+ϕb(I*+A*)+β(I*+A*)NS−(α+σ+μ)I*=0

(6) δ2T+αI*−(σ+d+μ)A*=0

From (6), we obtain,(7) A*=δ2T+αI*(σ+d+μ)

by substituting (7) into (4) and (5), respectively, and by making solving for S* and I* we obtainS*=(σ+d+μ)(α+σ+μ)N−(σ+d+μ)δ1TN−ϕbδ2TN−(σ+d+μ)ϕbTN−ϕbαTN(σ+d+μ)β+βαN+βδ2TN

I*=(σ+d+μ)δTN+ρS*N(σ+d+μ)−(1−ϕ)b(σ+d+μ)δ1TN−βSδ2T−(σ+d+μ)(σ+μ)N(σ+d+μ)β+βαN+βδ2TN

The inflow of tourists suggested by (1), (2), (3) does not show a disease-free equilibrium. The endemic equilibrium E*=(S*,I*,A*) occurs whenever there is HIV and AIDS found in the population. Now, from (2), the presence of E* is shown below(8) F(S)=δ1T+ϕb(I+A)+β(I+A)NS−(α+σ+μ)I

Now from (8), S≤δTμ+σ−ρ has a positive root range between 0 and δTμ+σ−ρ. By substituting when S=0, (5) becomes:(9) F(S)=δ1T+ϕb(I+A)‐(α+σ+μ)I

and when S=δTμ+σ−ρ, (8) gives:(10) F(δTμ+σ−ρ)=dIdt=δ1T+ϕb(I+A)+β(I+A)N(δTμ+σ−ρ)−(α+σ+μ)I

Now, when we differentiate (8) with respect to S, the final solution becomes:(11) dF(S)dt=β(I+A)N

2.3.2 Local stability of the equilibrium points

To determine the local stability of the endemicity of the disease, we use the eigenvalues of the characteristic equation of the Jacobian matrix for (1), (2), (3), J(S*,I*,A*)=J(E*) This is given by:(12) J(E*)=[ρ−β(I+A)N−(σ+μ)(1−ϕ)b−βSN(1−ϕ)b−βSNβ(I+A)Nϕb+βSN−(α+σ+μ)ϕb+βSN0α−(α+d+μ)]

The J(E*) corresponding equation of the characteristic polynomial equation is given by(13) f(λ)=(λ1−ρ+β(I+A)N+(σ+μ))(λ2−ϕb−βSN+(α+σ+μ))(λ3+(α+d+μ))=0

whereλ1=ρ−β(I+A)N−(σ+μ),

λ2=ϕb+βSN−(α+σ+μ),

λ3=−(α+d+μ).

2.4 Calculation of the reproduction number

The basic reproduction number R0 (i.e., epidemiologic metric used to describe the contagiousness or transmissibility of infectious agents), is the anticipated number of new HIV-positive cases that a typical infected person introduces into a group with a specific percentage of protection [52,53]. The model's parameters and variables were assumed to be non-negative. From equation (1) through (3), when the disease is absent at equilibrium equation (1) becomes S(t)=δTμ+σ−ρ. However, under the dynamic restrictions of (1) through (3), the region Σ={(x=S,I,A)∈R+3:S≤δTμ+σ−ρ} is positively invariant. So, for the initial starting point x∈R+3, the trajectory lies in Σ. The analysis in the current study is limited to this region Σ. The existence, uniqueness, and continuation results for equation (1) through (3) in this region are described in other studies [[54], [55], [56], [57], [58]]. From this, we now write the following two matrices:(14) F=(βSNβSN00)

(15) V=(−ϕb+(α+σ+μ)−ϕb−α(α+d+μ))

The inverse multiplication of (14) and (15) will give the spectral radius of the matrix of the basic reproduction number as βSN(α+d+μ)(α+d+μ)(α+σ+μ−ϕb)−ϕαb at equilibrium S(t)=δTμ+σ−ρ.

Therefore, the final basic reproduction number isR0=βδT(α+d+μ)N(μ+σ−ρ)(α+d+μ)(α+σ+μ−ϕb)−ϕαb

2.5 Model calibration

The model's parameters were calibrated to reflect Malaysian epidemiological data on HIV and AIDS incidence. Assume that θ represents the epidemiological parameters estimated and that θ is the independent and additive Gaussian prior. Let ϕ be the error, such that ϕ∼N(0,σ2). The Delayed-Rejection Adaptive Metropolis algorithm [59] was used to solve for θ, the posterior distribution of the defined model. The standard errors of the fitted parameters were approximated by taking the square root of the diagonal of the generalized inverse of the Hessian matrix of the fitted parameters [59]. Lastly, the model accuracy was evaluated by calculating the adjusted R2 [60,61] using the following equations:(16) R2=1−∑i=1n(yi−yˆi)2∑i=1n(yi−y‾i)2

and(17) Radj2=1−(1−R2)(n−1n−p)

where, yˆi, yˆi, and y‾i represent the observed data values, the predicted data values, and the mean of the observed data values, respectively, and n and p represent the number of data points and the number of independent variables.

2.6 Sensitivity analysis

Using the package 'Flexible Modelling Environment' (FME) [59] in R programming language [62] within RStudio, sensitivity analyses were conducted to identify the parameters whose values significantly impacted the model's predictions. Specifically, the sensFun function within the 'FME' package estimates the sensitivity of the model's outputs to parameter values within a set of sensitivity functions [59]. These functions estimate the instantaneous rate of change of the model's output relative changes to the value of its parameters. As described in another study [59], the function generated a matrix of normalized, dimensionless sensitivities whose (i,j)th element Sij contains the following:(18) ∂yi∂Θj*wΘjwyi

where, yi is an output variable and Θj is a parameter, and wyi and wΘj are the scaling weights of yi and Θj, respectively [58]. The L1 and L2 norms, where the L1 norm is defined as the sum of the absolute values of the vector (i.e., L1=∑|Sij|/n), and the L2 norm is defined as the square root of the sum of the squared vector values (i.e., L2=(Sij2)/n), can then be compared across parameters to determine the relative sensitivity of the model's outputs to variations in the parameter values. Specifically, the larger the magnitude of these norms, the more significant the impact of variations in the parameter's value on the model's output. Additionally, the effects of changes in β and δ1 on the spread of HIV were explored by comparing the number of diagnosed cases of HIV from the period 1986 to 2004 across repeated simulations using a range of values for these parameters.

To determine which of these parameters affects the R0, we subsequently conducted a sensitive analysis based on the R0 for each estimated parameter. This was done by performing the Partial Rank Correlation Coefficient (PRCC) by calculating the PRCC values for each estimated parameter based on the R0 [63,64].

3 Results

The model was fitted to the Malaysian data, resulting in, among other estimates, an estimated HIV transmission rate of 4.767x10−1 per susceptible individual, a rate of compartmental transfer from HIV-positive to AIDS-positive of 1.885x10−1 per HIV-positive individual, and an exit rate of foreign tourists of 5.208x10−1 (Table 1). The epidemiological data indicates that HIV incidence in Malaysia peaked around 2001, with a subsequent decline, followed by a peak in AIDS incidence in 2004. The model's fitted parameters reflect these epidemiological data (Fig. 2) with adjusted R2 values of 53.3% and 53.2% for HIV and AIDS, respectively.Table 1 Estimated parameters for the compartmental model based on Malaysian human immunodeficiency virus (HIV), acquired immunodeficiency syndrome (AIDS), and tourism data from 1986 to 2004.

Table 1Parameters	Definition	Estimate (95%CI)	
β	Transmission rate	4.767x10−1 (4.735x10−1 - 4.798x10−1)	
α	Rate of population transfer from HIV-positive to AIDS-positive	1.086x10−1 (1.081x10−1 - 1.092x10−1)	
μ	Natural death rate	4.56x10−3 (Source: [65])	
σ	Rate of foreign tourist exits	5.208x10−1 (5.176x10−1 - 5.239x10−1)	
δ	Rate of inbound tourist entries with susceptibility to HIV infection	9.964x10−1 (Source: [46])	
δ1	Rate of HIV-positive inbound tourist entries	1.992x10−3 (Source: [46])	
δ2	Rate of AIDS-positive inbound tourist entries	1.992x10−8 (Source: [46])	
ρ	Rate of Malaysian tourist entries returning home after contracting HIV abroad	9.837x10−9 (9.425x10−9 - 9.946x10−9)	
d	AIDS-related death rate	1.23x10−4 (Source: [66])	
b	Natural birth rate	2.91x10−3 (Source: [65])	
ϕ	The proportion of newborns born HIV-positive	4.429x10−6(3.257x10−5 - 4.143x10−5)	

Fig. 2 Reported cases of human immunodeficiency virus (HIV), corresponding to the Infected compartment (left), and acquired immunodeficiency syndrome (AIDS), corresponding to the AIDS compartment (right) in Malaysia (1986–2004) compared to the results of the fitted compartment model. Observations are shown in red while fitted model is shown in a solid black line.

Fig. 2

Table 2 presents the results of sensitivity analyses using L1 and L2 norms. Based on both norms, the predicted number of HIV cases was found to be most sensitive to changes in α, the rate of population transfer from HIV-positive to AIDS-positive, and σ, the rate of foreign tourist exits (i.e., the duration of time foreign tourists spends in Malaysia). The parameter with the next highest impact was β, followed by ρ, the rate of Malaysian tourist entries returning home after contracting HIV abroad, which had the least impact.Table 2 L1 and L2 norms were calculated to determine the sensitivity of the compartmental model's predicted number of human immunodeficiency virus (HIV) cases to variations in parameter values.

Table 2Parameter	L1 norm	L2 norm	
β	1.00x10−3	1.50x10−3	
α	3.30x10−1	5.10x10−1	
σ	2.20x10−1	2.60x10−1	
ρ	5.10x10−9	1.00x10−8	
ϕ	7.10x10−3	1.00x10−2	

According to Fig. 4 and Table 3, if the baseline foreign tourist exit rate (σ) is increased by 50%, there will be an estimated 218,377 cumulative HIV cases in the Malaysian population. This is equivalent to a ten-fold increase from the current estimated number of HIV cases, which is 82,255.Table 3 The effect of variations in δ1, the rate at which HIV-positive tourists entered the Malaysian population, and σ, the rate of foreign tourists exits, on the cumulative predicted number of HIV cases in Malaysia between 1986 and 2004 based on the compartmental model.

Table 3Changes inδ1	HIV cases	
No change (baseline)	82,255	
50% increase	122,234	
40% decrease	66,059	
60% decrease	49,743	
80% decrease	33,303	
Changes inσ	HIV cases	
No change (baseline)	82,254	
50% increase	218,377	
40% decrease	273,012	
60% decrease	123,452	
80% decrease	48,475	

Conversely, a reduction in σ translates to a significant decrease in the predicted number of HIV cases, with, for example, an 80% decrease from the baseline σ resulting in 48,475 predicted cases of HIV, a decrease of over 50% in the number of HIV cases. Similarly, the effects of δ1 on the number of predicted HIV cases is also shown in Fig. 4. A 50% increase in δ1 will result in 122,234 new cases of HIV in the Malaysian population compared to the baseline prediction.

Table 4 and Fig. 5 can be utilized to determine the most significant estimated parameters for the basic reproduction number output, considering the PRCC values. The bigger the absolute PRCC value, the greater the impact of the parameter on R0.Table 4 Partial Rank Correlation Coefficients (PRCC) of the estimated parameters with R0.

Table 4Parameters	Definition	PRCC values	
β	Transmission rate	0.304	
α	Rate of population transfer from HIV-positive to AIDS-positive	0.023	
σ	Rate of foreign tourist exits	−0.355	
ρ	Rate of Malaysian tourist entries returning home after contracting HIV abroad	−0.445	
ϕ	The proportion of newborns born HIV-positive	0.391	

4 Discussion and conclusions

A population-level mathematical model was developed to investigate the role of tourism in the spread of HIV in Malaysia and fitted to epidemiological data on Malaysian HIV and AIDS with explicit consideration for newborns born with HIV. It was determined that variations in the model's tourism-related parameters, such as the rate of foreign tourist exits, which reflects the duration of touristic visits, and the rate of HIV-positive inbound tourist entries, which reflects the overall rate of touristic entries, significantly impacted the predicted cumulative incidence of HIV in Malaysia (Table 2, Table 3, and Fig. 3). Understanding the role of tourism in the spread of infectious diseases is crucial for national efforts to control or prevent the spread of HIV, particularly for countries that rely heavily on tourism, such as Malaysia. Additionally, this research suggests that a steady-state level of HIV in Malaysia is achievable, from a disease transmission perspective, with an estimated basic reproduction number R0<1.Fig. 3 Reported cases of human immunodeficiency virus (HIV), corresponding to the Infected compartment (top) and acquired immunodeficiency syndrome (AIDS), corresponding to the AIDS compartment (bottom) in Malaysia (1986–2004), where the light gray shade represents the upper confidence interval and dark gray shade represents the lower confidence interval (Min-Max).

Fig. 3

Fig. 4 Graphical comparison of the effect of variations in δ1, the rate at which tourists moved to the HIV-positive compartment (top) and σ, the rate of foreign tourist exits (bottom) on the predicted number of HIV cases in Malaysia between 1986 and 2004 based on the compartmental model.

Fig. 4

Fig. 5 PRCC for the basic reproduction number with each of the estimated parameters.

Fig. 5

Mathematical modelling is a valuable tool in researching different aspects of HIV [[67], [68], [69]], such as transmission between sexual partners [[70], [71], [72]] and intravenous drug users [73], as well as cell-to-cell spread [74]. It can also help prevent drug-resistant HIV and aid in the study the infection's pathogenesis and treatment [75]. Compartmental models are particularly relevant to this study's results, as they extend the work of Apenteng and Ismail (2014) and Apenteng (2017) by incorporating the proportion of newborns infected with HIV (as shown in Table 1) [42,43].

Specifically, by expanding on the Malaysian compartmental HIV model [42,43] to include the proportion of newborns ϕb(I+A) and, thus, the proportion of newborns born with HIV (1-ϕ)b(I+A), the current model more accurately reflects HIV transmission dynamics and yields different parameter estimations. For example, in the study by Apenteng (2017) [43], the estimated parameters for the rate of inbound tourists entering the susceptible, HIV, and AIDS compartments were 9.98x10−2, 3.19x10−5, and 5.86x10−2, respectively (Table 1).

This study proposes a new method to evaluate the impact of tourism on HIV transmission. The framework used in this method also estimates epidemiological parameters that can help develop surveillance and control efforts by identifying the factors that strongly influence HIV transmission. Based on the parameter vector's L1 and L2 norms, the sensitivity analyses suggested that the rate of foreign tourist exits (σ) was the most sensitive factor, meaning that the duration of time that foreign tourists spend in Malaysia has a significant impact on the predicted number of HIV cases. Despite identifying the potential role of touristic flows in HIV transmission, the formulation of an appropriate policy response is complex and may remain difficult to achieve. This complexity is apparent in compartmental models of HIV transmission in other regions and countries.

Positive PRCC values signify the positive correlation between an estimated parameter and R0, as seen with ∅,α, and β. In contrast, negative PRCC values signify the negative correlation between an estimated parameter and R0 , as seen with ρ and σ (for more details, see Table 4 and Fig. 5). These PRCC values give information on potentially effective mitigation plans to prevent the spread of HIV and AIDS. To determine the eigenvalues, we substituted the values of the estimated parameters and the initial values of the compartments; we have λ1=−0.5253, λ2=−0.2372 , and λ3=−0.5254 from (13). This shows that the disease-free equilibrium is stable since all the eigenvalues have negative parts or values [76].

As an example, Nyabadza and Mukandavire (2011) developed a compartmental model in 2011 to predict the trend of HIV transmission in South Africa. They used simulations to fit their model [38], similar to the approach described in this paper. Their findings suggested that an individual's knowledge about HIV/AIDS, which is typically obtained through treatment provided during HIV testing and screening, can help reduce the likelihood of transmission. However, based on a further exploration of the impact of screening and treatment on the dynamics of HIV transmission in South Africa, the same authors suggested that, despite 1.6 million people undergoing HIV screening annually, screening and surveillance efforts alone will not stop the spread of HIV. Further, Nyabadza and Mukandavire (2011) estimated that the HIV epidemic in South Africa had an R0 value of less than one, which means that the disease could not spread rapidly. However, they also noted that HIV could continue to spread without bounds unless there is a change in sexual behaviour in the country. To prevent new HIV infections, the authors suggested that people should reduce the number of sexual partners, use condoms, educate individuals about HIV, and change their sexual conduct overall [46]. This current study suggests that a steady-state level of HIV in Malaysia is achievable, from a disease transmission perspective, with an estimated reproduction number R0 of 0.0136 (Table 1). In epidemiological terms, it is expected that transmission will wane and eventually stop when R0≤1 and conversely, disease transmission will persist or increase indefinitely when R0>1 [54,77]. It has been observed that HIV and AIDS are prevalent diseases in Malaysia. However, taking preventive measures such as reducing the number of sexual partners and using condoms can help in phasing out the disease. Nevertheless, it is essential to consider tourism as an ongoing risk for introducing diseases, as demonstrated in this study.

Bacaër et al. (2006) created a model to study the HIV/AIDS epidemic among injectable drug users and sex workers in Kunming, Yunnan, the provincial capital of China [40]. They developed a compartmental model in which individuals went through three stages, namely susceptible S, HIV-positive-AIDS-negative I, and AIDS-positive A compartments. Interestingly, their estimated basic reproduction number was 32, which means that, despite their assumption that individuals in the AIDS compartment did not have contact with others, the public health system was far more unstable than the one modeled in the current study. Bozkurt and Peker (2014) proposed that educating individuals about HIV and the risks of engaging in unsafe sexual behavior or any other risky behavior that increases the risk of HIV transmission would have the most significant impact on lowering the infection rate and prevalence [68]. Their compartmental model suggested that immigration is positively related to the spread of the disease. At the same time, a combination of alerting individuals of their HIV-positive status through screening and contact tracing can result in increased abstinence and, therefore, a reduction in spread. However, if contact tracing is omitted, the spread will increase regardless of status notification through testing and screening.

Naresh et al. (2008) conducted a study to understand the impact of HIV-positive immigration on different population sizes. Their research primarily focused on how the disease spreads through sexual intercourse and assumed that all infected individuals eventually progress to the AIDS compartment. Their findings suggested that, by limiting the number of sexual partners, the transmission of the disease can be reduced over time. In other words, the study highlights the importance of taking preventive measures to control the spread of HIV, especially in communities with high-risk factors. Benotsch et al. (2006) studied visitors to a popular tourist destination in the US. They found that these travelers were at a higher risk of contracting HIV and could potentially spread the virus to their home populations [78]. Zöldi et al. (2017) performed a non-compartmental study from 1995 to 2015 in Finland and examined all cases of STIs related to foreign travel. The authors discovered that out of approximately 9,000 STI cases, 2,300 (26%) were linked to travel. The study estimated the total risks and analyzed trends across different geographical areas [79].

Tourism is an essential source of income for many countries around the world. Therefore, when developing public health policies to prevent the spread of diseases, it is crucial to take into account the role of tourism. Some countries, such as the United States, Canada, Australia, New Zealand, Spain, Croatia, Denmark, and others, require visitors to undergo a medical examination before entering their territory, the requirement being based on the length of stay. This policy is valid as this study suggests that the incidence of HIV and other transmissible diseases is related to the duration of stay (σ or the rate of foreign tourist exit). Recently, Selvanathan et al. (2021) found that tourism significantly contributed to the spread of COVID-19 within and across countries due to the movement of individuals [80]. The COVID-19 pandemic resulted in the widespread adoption of travel restrictions [81], reflecting the vital role of transient populations such as international travelers and tourists [82,83]. Travel restrictions may be effective in reducing the transmission of HIV, however, imposing travel restrictions based on an individual's HIV status may be unethical. Therefore, public health policy should focus on developing and promoting education initiatives aimed at reducing the transmission of HIV among at-risk populations. In addition, comprehensive surveillance and contact tracing capacities should be developed to control the transmission of HIV further.

Finally, it is important to recognize the limitations of this study. We had to estimate the rate of Malaysians returning home after contracting unconfirmed HIV abroad, as well as the rate of foreign tourist exits since the actual values of these parameters were not available in the literature or the epidemiological data used in this study. The findings of this study may be negatively impacted by the lack of comprehensive epidemiological data from the Malaysian HIV outbreak up to this point. This is because the estimations of the parameters were affected by the amount of noise in the data [84,85]. Furthermore, while the average number of days spent in Malaysia was incorporated into the proposed model as a fixed parameter, the actual period a tourist stays in Malaysia may vary significantly in real-world conditions. Nevertheless, the parameters δ, δ1, δ2, and σ indicate that tourism had a substantial impact on HIV transmission in Malaysia from 1986 to 2004 (Table 3), based on the available data, and it is possible that by collecting more observations, we can improve the model's accuracy.

The study highlights the need for the country to take significant steps in quantifying the trend and magnitude of tourism. Developing platforms that will enable public health officials to implement effective control strategies is also crucial. The study highlights the need for the country to take significant steps in quantifying the trend and magnitude of tourism. It is also crucial to develop platforms that will enable public health officials to implement effective control strategies. The findings of this study suggest that visitors to Malaysia have played a role in the spread of HIV and AIDS infection among the Malaysian population from 1986 to 2004. The findings suggest that a steady-state level of HIV in Malaysia is achievable based on the low value of R0 = 0.0136, which is encouraging from a public health perspective. Moreover, the modelling framework can be used to explore the roles of other population groups and transmission parameters. For instance, future modelling studies should consider factors such as engaging in sexual activity with multiple partners and not protecting oneself against HIV.

Data availability statement

Data supporting the results presented in this manuscript is available as stated above.

CRediT authorship contribution statement

Ofosuhene O. Apenteng: Conceptualization, Data curation, Formal analysis, Methodology, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. Philip Rasmussen: Conceptualization, Validation, Visualization, Writing – review & editing. Beate Conrady: Validation, Visualization, Writing – review & editing.

Declaration of competing interest

All authors declared that no competing interests exist.
==== Refs
References

1 Gray M. Nieburg P. Clinics R.D. Pediatric human immunodeficiency virus continuum of care: a concise review of evidence-based practice Pediatr. Clin. 64 4 2017 879 891
2 Anderson R.M. Medley G.F. May R.M. Johnson A.M. A preliminary study of the transmission dynamics of the human immunodeficiency virus (HIV), the causative agent of AIDS http://imammb.oxfordjournals.org/ 1986
3 Bitnun A. Samson L. Chun T.W. Kakkar F. Brophy J. Murray D. Justement S. Soudeyns H. Ostrowski M. Mujib S. Harrigan P.R. Kim J. Sandstrom P. Read S.E. Early initiation of combination antiretroviral therapy in HIV-1 Infected newborns can achieve sustained virologic suppression with low frequency of CD4+ T cells carrying HIV in peripheral blood Clin. Infect. Dis. 59 2014 10.1093/cid/ciu432
4 Naif H.M. Pathogenesis of HIV infection Infect. Dis. Rep. 5 2013 10.4081/idr.2013.s1.e6
5 Cunningham A.L. Li S. Juarez J. Lynch G. Alali M. Naif H. The level of HIV infection of macrophages is determined by interaction of viral and host cell genotypes J. Leukoc. Biol. 68 2000 10.1189/jlb.68.3.311
6 Palefsky J.M. Holly E.A. Chapter 6: immunosuppression and co-infection with HIV J. Natl. Cancer Inst. Monogr. 2003 10.1093/oxfordjournals.jncimonographs.a003481
7 Galvin S.R. Cohen M.S. The role of sexually transmitted diseases in HIV transmission Nat. Rev. Microbiol. 2 2004 10.1038/nrmicro794
8 Valverde E.E. Painter T. Heffelfinger J.D. Schulden J.D. Chavez P. DiNenno E.A. Migration patterns and characteristics of sexual partners associated with unprotected sexual intercourse among hispanic immigrant and migrant women in the United States J Immigr Minor Health 17 2014 10.1007/s10903-014-0132-6
9 Choi K.H. Hudes E.S. Steward W.T. Social discrimination, concurrent sexual partnerships, and hiv risk among men who have sex with men in Shanghai, China AIDS Behav. 12 2008 10.1007/s10461-008-9394-0
10 Bond L. Wheeler D.P. Millett G.A. LaPollo A.B. Carson L.F. Liau A. Black men who have sex with men and the association of down-low identity with HIV risk behavior Am. J. Publ. Health 99 Suppl 1 2009 10.2105/AJPH.2007.127217
11 Greene S. Ion A. Elston D. Kwaramba G. Smith S. Carvalhal A. Loutfy M. “Why aren't you breastfeeding?”: how mothers living with HIV talk about infant feeding in a “breast is best” world Health Care Women Int. 36 2015 883 901 10.1080/07399332.2014.888720 24527767
12 Hadland S.E. Werb D. Kerr T. Fu E. Wang H. Montaner J.S. Wood E. Childhood sexual abuse and risk for initiating injection drug use: a prospective cohort study Prev. Med. 55 2012 10.1016/j.ypmed.2012.08.015
13 Valjent E. Bertran-Gonzalez J. Aubier B. Greengard P. Hervé D. Girault J.A. Mechanisms of locomotor sensitization to drugs of abuse in a two-injection protocol Neuropsychopharmacology 35 2010 10.1038/npp.2009.143
14 Ulibarri M.D. Strathdee S.A. Ulloa E.C. Lozada R. Fraga M.A. Magis-Rodríguez C. De La Torre A. Amaro H. O'Campo P. Patterson T.L. Injection drug use as a mediator between client-perpetrated abuse and HIV status among female sex workers in two Mexico-US border cities AIDS Behav. 15 2011 10.1007/s10461-009-9595-1
15 Prangnell A. Imtiaz S. Karamouzian M. Hayashi K. Childhood abuse as a risk factor for injection drug use: a systematic review of observational studies Drug Alcohol Rev. 39 2020 10.1111/dar.13001
16 McIntosh E.D.G. Development of vaccines against the sexually transmitted infections gonorrhoea, syphilis, Chlamydia, herpes simplex virus, human immunodeficiency virus and Zika virus Ther Adv Vaccines Immunother 8 2020 10.1177/2515135520923887
17 Du M. Yan W. Jing W. Qin C. Liu Q. Liu M. Liu J. Increasing incidence rates of sexually transmitted infections from 2010 to 2019: an analysis of temporal trends by geographical regions and age groups from the 2019 Global Burden of Disease Study BMC Infect. Dis. 22 2022 10.1186/s12879-022-07544-7
18 Naik P.A. Yeolekar B.M. Qureshi S. Yeolekar M. Madzvamuse A. Modeling and analysis of the fractional-order epidemic model to investigate mutual influence in HIV/HCV co-infection Nonlinear Dynam. 2024 10.1007/s11071-024-09653-1
19 Bauer I. The health impact of tourism on local and indigenous populations in resource-poor countries Trav. Med. Infect. Dis. 6 2008 10.1016/j.tmaid.2008.05.005
20 Jalil A. Mahmood T. Idrees M. Tourism-growth nexus in Pakistan: evidence from ARDL bounds tests Econ. Modell. 35 2013 10.1016/j.econmod.2013.06.034
21 Habibi F. The determinants of inbound tourism to Malaysia: a panel data analysis Curr. Issues Tourism 20 2017 10.1080/13683500.2016.1145630
22 T. and I.-M. Ministry of Investment Third Industrial Master Plan (IMP3) 2006-2020 2006
23 Shahbaz M. Benkraiem R. Miloudi A. Tiwari A.K. Tourism-induced financial development in Malaysia: new evidence from the tourism development index Tourism Econ. 25 2019 10.1177/1354816618806123
24 Kibona I. Mahera W. Makinde D. Mango J. A deterministic model of HIV/AIDS with vertical transmission in the presence of infected immigrants Int. J. Phys. Sci. 6 2011 5383 5398 10.5897/IJPS11.452
25 Malaysian Immigration Department Tourist Arrivals to Malaysia by Country of Nationality December 2020 2020
26 Jong M.C. Puah C.H. Arip M.A. Assessing the impacts of tourism climate index and holiday climate index on tourism demand in Malaysia Asian J. Bus. Res. 13 2023 10.14707/ajbr.220142
27 Ab Rahman S.A. Bachok S. Mahamod L.H. Effects of travel characteristics on tourism expenditure: a case study of Malaysia young outbound tourists Planning Malaysia 21 2023 10.21837/PM.V21I25.1233
28 Bauer I. Understanding sexual relationships between tourists and locals in Cuzco/Peru Trav. Med. Infect. Dis. 5 2007 10.1016/j.tmaid.2007.06.004
29 Ketshabile L.Simeon HIV and AIDS as a Threat to Southern African Tourism 2007 Cape Peninsula University of Technology
30 Wright E.R. Travel, tourism, and HIV risk among older adults J. Acquir. Immune Defic. Syndr. 1988 33 10.1097/00126334-200306012-00023 2003
31 Padilla M.B. Guilamo-Ramos V. Bouris A. Reyes A.M. HIV/AIDS and tourism in the Caribbean: an ecological systems perspective Am. J. Publ. Health 100 2010 10.2105/AJPH.2009.161968
32 Guilamo-Ramos V. Padilla M. Cedar A.L. Lee J. Robles G. HIV sexual risk behavior and family dynamics in a Dominican tourism town Arch. Sex. Behav. 42 2013 10.1007/s10508-012-0064-y
33 de Matos M.A. Caetano K.A.A. França D.D. da S. Pinheiro R.S. de Moraes L.C. Teles S.A. Vulnerability to sexually transmitted infections in women who sell sex on the route of prostitution and sex tourism in central Brazil Rev Lat Am Enfermagem 21 2013 10.1590/S0104-11692013000400011
34 MOH The Global AIDS Response Progress Report 2014 Malaysia, (2015
35 Chuah F.L.H. Tan S.T. Yeo J. Legido-Quigley H. The health needs and access barriers among refugees and asylum-seekers in Malaysia: a qualitative study Int. J. Equity Health 17 2018 10.1186/s12939-018-0833-x
36 Richens J. Sexually transmitted infections and HIV among travellers: a review Trav. Med. Infect. Dis. 4 2006 184 195 10.1016/j.tmaid.2005.06.010
37 Uusküla A. McMahon J.M. Raag M. Silm S. Rüütel K. Talu A. Abel-Ollo K. Ahas R. Des Jarlais D.C. Emergent properties of HIV risk among injection drug users in Tallinn, Estonia: synthesis of individual and neighbourhood-level factors Sex. Transm. Infect. 86 Suppl 3 2010 10.1136/sti.2009.040212
38 Nyabadza F. Mukandavire Z. Modelling the HIV/AIDS epidemic trends in South Africa: insights from a simple mathematical model Nonlinear Anal. Real World Appl. 12 4 2011 2019 2104
39 Nyabadza F. Mukandavire Z. Modelling HIV/AIDS in the presence of an HIV testing and screening campaign J. Theor. Biol. 280 1 2011 167 179 21536051
40 Bacaër N. Abdurahman X. Ye J. Modeling the HIV/AIDS epidemic among injecting drug users and sex workers in Kunming, China Bull. Math. Biol. 68 2006 10.1007/s11538-005-9051-y
41 Hyman J.M. Li J. Ann Stanley E. The differential infectivity and staged progression models for the transmission of HIV Math. Biosci. 155 1999 10.1016/S0025-5564(98)10057-3
42 Apenteng O.O. Ismail N.A. Modelling the impact of international travellers on the trend of HIV epidemic Lecture Notes in Engineering and Computer Science 2014
43 Apenteng O.O. Assessing the effects of tourism on the spread of HIV and AIDS in Malaysia using susceptible infected removed models http://studentsrepo.um.edu.my/7980/9/ofosuhene.pdf 2017
44 Jamal T. Budke C. Tourism in a world with pandemics: local-global responsibility and action J. Tourism Futur. 6 2020 10.1108/JTF-02-2020-0014
45 Ullah S. Altanji M. Khan M.A. Alshaheri A. Sumelka W. The dynamics of HIV/AIDS model with fractal-fractional Caputo derivative Fractals 31 2023 10.1142/S0218348X23400157
46 ari Ngadiman, Malaysia The Global AIDS Response Progress Report 2012 2012
47 Nicholas B. The Impact of the Existence of Ulu Skrang Roads upon Tourism Activities and Rural Livelihood 2005
48 Apenteng O.O. Ismail N.A. A Markov chain Monte Carlo approach to estimate AIDS after HIV infection PLoS One 10 2015 10.1371/journal.pone.0131950
49 Oduro B. Apenteng O.O. Nkansah H. Assessing the effect of fungicide treatment on cocoa black pod disease in Ghana: insight from mathematical modeling, Statistics Optimization and Information Computing 8 2020 10.19139/soic-2310-5070-686
50 Apenteng O.O. Ismail N.A. Modelling the impact of migration on HIV persistency in Ghana, statistics, optimization and information computing 7 10.19139/soic.v7i1.552 2019
51 Naik P.A. Global dynamics of a fractional-order SIR epidemic model with memory Int. J. Biomath. (IJB) 13 2020 10.1142/S1793524520500710
52 Van Den Driessche P. Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission www.elsevier.com/locate/mbs 2002
53 P. Van Den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, n.d. www.elsevier.com/locate/mbs (accessed July 9, 2020).
54 Van Den Driessche P. Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Math. Biosci. 180 2002 29 48 10.1016/S0025-5564(02)00108-6 12387915
55 Diekmann O. Heesterbeek J.A.P. Roberts M.G. The construction of next-generation matrices for compartmental epidemic models 10.1098/rsif.2009.0386 2009
56 Roberts M.G. Heesterbeek J.A.P. Model-consistent estimation of the basic reproduction number from the incidence of an emerging infection J. Math. Biol. 55 2007 803 816 10.1007/s00285-007-0112-8 17684743
57 Zhai X. Li W. Wei F. Mao X. Dynamics of an HIV/AIDS transmission model with protection awareness and fluctuations Chaos, Solit. Fractals 169 2023 10.1016/j.chaos.2023.113224
58 Li W. Cai S. Zhai X. Ou J. Zheng K. Wei F. Mao X. Transmission dynamics of symptom-dependent HIV/AIDS models Math. Biosci. Eng. 21 2024 10.3934/mbe.2024079
59 Soetaert K. Petzoldt T. Inverse modelling, sensitivity and Monte Carlo analysis in R using package FME J. Stat. Softw. 33 3 2010
60 Bhandari A. Difference between R-squared and adjusted R-squared Analytics Vidhya 2020
61 Miles J. Squared R. Adjusted R ‐squared Encyclopedia of Statistics in Behavioral Science 2005 10.1002/0470013192.bsa526
62 Team R.C. R: A Language and Environment for Statistical Computing 2021 R Foundation for Statistical Computing
63 Mangongo Y.T. Bukweli J.D.K. Kampempe J.D.B. Mabela R.M. Munganga J.M.W. Stability and global sensitivity analysis of the transmission dynamics of malaria with relapse and ignorant infected humans Phys. Scripta 97 2022 10.1088/1402-4896/ac4862
64 Bürli C. Harbrecht H. Odermatt P. Sayasone S. Chitnis N. Mathematical analysis of the transmission dynamics of the liver fluke, Opisthorchis viverrini J. Theor. Biol. 439 2018 10.1016/j.jtbi.2017.11.020
65 Lan W.S. Ismail H.I.M. A Study on under Five Deaths in Malaysia in the Year 2006 2008 Ministry of Health of Malaysia 37 38
66 Darani W.N.S.W.M. Chen X.W. Samsudin E.Z. Nor F.M. Ismail I. Determinants of successful human immunodeficiency virus treatment outcomes: a linkage of national data sources in Malaysia Malays. J. Med. Sci. 30 2023 10.21315/mjms2023.30.1.15
67 Cassels S. Clark S.J. Morris M. Mathematical models for HIV transmission dynamics: tools for social and behavioral science research J Acquir Immune Defic Syndr 1988 10.1097/QAI.0b013e3181605da3 2008
68 Bozkurt F. Peker F. Mathematical modelling of HIV epidemic and stability analysis Adv. Differ. Equ. 2014 2014 10.1186/1687-1847-2014-95
69 Anderson R.M. The role of mathematical models in the study of HIV transmission and the epidemiology of AIDS J. Acquir. Immune Defic. Syndr. 1988 1 1988
70 Rao D.W. Wheatley M.M. Goodreau S.M. Enns E.A. Partnership dynamics in mathematical models and implications for representation of sexually transmitted infections: a review Ann. Epidemiol. 59 2021 10.1016/j.annepidem.2021.04.012
71 Johnson L.F. Dorrington R.E. Bradshaw D. Pillay-Van Wyk V. Rehle T.M. Sexual behaviour patterns in South Africa and their association with the spread of HIV: insights from a mathematical model Demogr. Res. 21 2009 10.4054/DemRes.2009.21.11
72 Punyacharoensin N. Edmunds W.J. De Angelis D. White R.G. Mathematical models for the study of HIV spread and control amongst men who have sex with men Eur. J. Epidemiol. 26 2011 10.1007/s10654-011-9614-1
73 Greenhalgh D. Mathematical modelling of the spread of HIV/AIDS amongst injecting drug users Math. Med. Biol. 14 1997 10.1093/imammb/14.1.11
74 Culshaw R.V. Ruan S. Webb G. A mathematical model of cell-to-cell spread of HIV-1 that includes a time delay J. Math. Biol. 46 2003 10.1007/s00285-002-0191-5
75 Wodarz D. Nowak M.A. Mathematical models of HIV pathogenesis and treatment Bioessays 24 2002 10.1002/bies.10196
76 Fudolig M. Howard R. The local stability of a modified multi-strain SIR model for emerging viral strains PLoS One 15 2020 10.1371/journal.pone.0243408
77 Jones J.H. Notes on R 0 1 the Basic Reproduction Number in a Nutshell 2013
78 Benotsch E.G. Mikytuck J.J. Ragsdale K. Pinkerton S.D. Sexual risk and HIV acquisition among men who have sex with men travelers to Key West, Florida: a mathematical modeling analysis AIDS Patient Care STDS 20 2006 10.1089/apc.2006.20.549
79 Zöldi V. Sane J. Kantele A. Rimhanen-Finne R. Salmenlinna S. Lyytikäinen O. Destination specific risks of acquisition of notifiable food- and waterborne infections or sexually transmitted infections among Finnish international travellers, 1995–2015 Trav. Med. Infect. Dis. 25 2018 10.1016/j.tmaid.2017.10.006
80 Selvanathan E.A. Jayasinghe M. Selvanathan S. International tourism and infectious disease transmission nexus: a cross-country and regional study J. Trav. Res. 61 2022 10.1177/00472875211048932
81 Naik P.A. Farman M. Zehra A. Nisar K.S. Hincal E. Analysis and modeling with fractal-fractional operator for an epidemic model with reference to COVID-19 modeling Partial Differential Equations in Applied Mathematics 10 2024 10.1016/j.padiff.2024.100663
82 Fotiadis A. Polyzos S. Huan T.C.T.C. The good, the bad and the ugly on COVID-19 tourism recovery Ann. Tourism Res. 87 2021 10.1016/j.annals.2020.103117
83 Bou-Karroum L. Khabsa J. Jabbour M. Hilal N. Haidar Z. Abi Khalil P. Khalek R.A. Assaf J. Honein-AbouHaidar G. Samra C.A. Hneiny L. Al-Awlaqi S. Hanefeld J. El-Jardali F. Akl E.A. El Bcheraoui C. Public health effects of travel-related policies on the COVID-19 pandemic: a mixed-methods systematic review J. Infect. 83 2021 10.1016/j.jinf.2021.07.017
84 Gallagher M.T. Wain R.A.J. Dari S. Whitty J.P. Smith D.J. Non-identifiability of parameters for a class of shear-thinning rheological models, with implications for haematological fluid dynamics J. Biomech. 85 2019 10.1016/j.jbiomech.2019.01.036
85 Chen S. McMullan G. Faruqi A.R. Murshudov G.N. Short J.M. Scheres S.H.W. Henderson R. High-resolution noise substitution to measure overfitting and validate resolution in 3D structure determination by single particle electron cryomicroscopy Ultramicroscopy 135 2013 10.1016/j.ultramic.2013.06.004
