Using step selection functions to analyse human mobility using telemetry data in infectious disease epidemiology: a case study of leptospirosis
Figures
Map showing location of each study area in Salvador.
Each area includes symbology for stream (blue line), open sewer points (purple diamond), and domestic rubbish piles (orange triangle).
Schematic diagram showing what data sources are used in which model, and how models are linked with each other.
The blue sections represent phase 1, the individual-level models, whilst the orange section represents phase 2, the population-level model.
Descriptive diagram of step selection functions.
(A) Step lengths () and turning angles () are used to characterise an individual’s movements. (B) These parameters are used to create a set of available steps (grey dots) for every used step (black dots).
Graph showing results of final analyses.
(A) Results for distance-based rasters; values above 1 interpreted as increasing distance to points of reference. (B) Results for 20 m buffer-based rasters; values above 1 show movement within buffer zones. Each horizontal band represents a specific time period (right hand side y-axis label): all day (5 am to 9 pm, Tables 2 and 3), morning (5 am to 9 am), midday (9 am to 1 pm), afternoon (1 pm to 5 pm), and evening (5 pm to 9 pm). All data points include their corresponding 95% confidence intervals, some of which are too narrow to show up clearly.
Distribution of telemetry data provided by each individual across 24 hr periods (x-axis), separated into each of the four study areas (1: NVS, 2: ARE, 3: JSI, 4: CAL).
Overlapping areas represent multiple days. Vertical bars represent 5 am (left-hand bar) and 9 pm (right-hand bar), the period of analysis.
Distribution of Leptospirosis antibody status (serological status) by gender and age.
Distribution of nearest distance to each of the environmental factors being analysed (central stream, open sewer points, and domestic rubbish piles) by serological status (x-axis) and study area (1: NVS, 2: ARE, 3: JSI, 4: CAL).
NA represents the rest of households in the study area that did not take part in movement analysis.
Tables
Summary table comparing parent study participants and movement study participants.
| Parent study participants | Movement study participants | |||
|---|---|---|---|---|
| n=1086 | %/mean | n=128 | %/mean | |
| Study area | ||||
| 1—Nova Sussuarana | 297 | 27.3% | 32 | 25.0% |
| 2—Arenoso | 246 | 22.7% | 28 | 21.9% |
| 3—Jardim Santo Inacio | 278 | 25.6% | 35 | 27.3% |
| 4—Calabetao | 265 | 24.4% | 33 | 25.7% |
| Gender | ||||
| Female | 454 | 41.8% | 59 | 46.1% |
| Male | 632 | 58.2% | 69 | 53.9% |
| Age (mean±SD) | 32.2 | ±19.7 | 39.4 | ±15.4 |
| Leptospira antibody status | ||||
| Positive | 94 | 8.7% | 13 | 10.2% |
| Negative | 992 | 91.3% | 115 | 89.8% |
Proportion of tracked time (full-day period, 9 am to 5 pm) spent within each buffer.
Mean (standard deviation).
| River buffer | Open sewer buffer | Domestic rubbish buffer | ||
|---|---|---|---|---|
| Total | 0.54 (0.33) | 0.33 (0.31) | 0.07 (0.15) | |
| Area | 1—Nova Sussuarana | 0.54 (0.29) | 0.43 (0.29) | 0.15 (0.18) |
| 2—Arenoso | 0.65 (0.25) | 0.28 (0.29) | 0.03 (0.13) | |
| 3—Jardim Santo Inacio | 0.36 (0.35) | 0.41 (0.31) | 0.08 (0.18) | |
| 4—Calabetao | 0.67 (0.31) | 0.21 (0.28) | 0.01 (0.02) | |
| Gender | Female | 0.56 (0.32) | 0.37 (0.32) | 0.07 (0.15) |
| Male | 0.53 (0.34) | 0.30 (0.29) | 0.06 (0.16) | |
| Leptospirosis serological status | Negative | 0.55 (0.33) | 0.35 (0.31) | 0.07 (0.16) |
| Positive | 0.48 (0.29) | 0.25 (0.30) | 0.02 (0.05) |
Estimated differences (γ) in selection coefficients (β) for each environmental factor using distance-based rasters.
Values >1 represent increasing distance from points of reference.
| Community stream | Open sewers | Domestic rubbish piles | ||||
|---|---|---|---|---|---|---|
| Estimate | 95% CI | Estimate | 95% CI | Estimate | 95% CI | |
| Gender* | ||||||
| Male | (Ref) | – | (Ref) | – | (Ref) | – |
| Female | 0.98 | 0.97, 0.99 | 1.04 | 1.02, 1.06 | 0.99 | 0.98, 1.01 |
| Age† | ||||||
| Per year increase | 1.00 | 1.00, 1.00 | 1.00 | 1.00, 1.00 | 1.00 | 1.00, 1.00 |
| Leptospira serological status ‡ | ||||||
| Negative | (Ref) | – | (Ref) | – | (Ref) | – |
| Positive | 0.99 | 0.96, 1.01 | 1.03 | 1.00, 1.07 | 1.00 | 0.98, 1.02 |
-
(Ref) is the referrence group for the Odds Ratio.
-
*
Adjusted for age and study area.
-
†
Adjusted for gender and study area. Values represent increases by 1 year of age.
-
‡
Adjusted for gender, age, and study area.
Estimated differences (γ) in selection coefficients (β) for each environmental factor using 20 m buffers around each point of reference.
Values>1 represent movement within the buffer zone for each point of reference.
| Community stream | Open sewers | Domestic rubbish piles | ||||
|---|---|---|---|---|---|---|
| Estimate | 95% CI | Estimate | 95% CI | Estimate | 95% CI | |
| Gender* | ||||||
| Male | (Ref) | – | (Ref) | – | (Ref) | – |
| Female | 1.22 | 1.02, 1.46 | 0.95 | 0.80, 1.14 | 0.92 | 0.66, 1.27 |
| Age† | ||||||
| Per year increase | 1.00 | 1.00, 1.00 | 0.99 | 0.98, 1.00 | 1.00 | 0.99, 1.01 |
| Leptospira serological status ‡ | ||||||
| Negative | (Ref) | – | (Ref) | – | (Ref) | – |
| Positive | 0.89 | 0.67, 1.19 | 0.64 | 0.47, 0.87 | 0.85 | 0.48, 1.49 |
-
(Ref) is the referrence group for the Odds Ratio.
-
*
Adjusted for age and study area.
-
†
Adjusted for gender and study area. Values represent increases by 1 year of age.
-
‡
Adjusted for gender, age, and study area.
Demographic details of excluded individuals due to having less than 50 relocations.
| ID (anonymised) | Relocation below 50 | Period | Gender | Age group | Leptospirosis serological status |
|---|---|---|---|---|---|
| 60 | TRUE | 05–21 | Male | 50–54 | Neg |
| 91 | TRUE | 05–21 | Male | >55 | Neg |
| 15 | TRUE | 05–09 | Female | 45–49 | Neg |
| 60 | TRUE | 05–09 | Male | 50–54 | Neg |
| 81 | TRUE | 05–09 | Female | 50–54 | Neg |
| 91 | TRUE | 05–09 | Male | >55 | Neg |
| 108 | TRUE | 05–09 | Female | 50–54 | Neg |
| 109 | TRUE | 05–09 | Male | 20–24 | Neg |
| 128 | TRUE | 05–09 | Male | 25–29 | Neg |
| 129 | TRUE | 05–09 | Male | 40–44 | Pos |
| 15 | TRUE | 09–13 | Female | 45–49 | Neg |
| 24 | TRUE | 09–13 | Female | 50–54 | Neg |
| 60 | TRUE | 09–13 | Male | 50–54 | Neg |
| 70 | TRUE | 09–13 | Male | 35–39 | Neg |
| 71 | TRUE | 09–13 | Female | 35–39 | Neg |
| 76 | TRUE | 09–13 | Male | 35–39 | Neg |
| 91 | TRUE | 09–13 | Male | >55 | Neg |
| 108 | TRUE | 09–13 | Female | 50–54 | Neg |
| 109 | TRUE | 09–13 | Male | 20–24 | Neg |
| 128 | TRUE | 09–13 | Male | 25–29 | Neg |
| 129 | TRUE | 09–13 | Male | 40–44 | Pos |
| 24 | TRUE | 13–17 | Female | 50–54 | Neg |
| 60 | TRUE | 13–17 | Male | 50–54 | Neg |
| 71 | TRUE | 13–17 | Female | 35–39 | Neg |
| 76 | TRUE | 13–17 | Male | 35–39 | Neg |
| 91 | TRUE | 13–17 | Male | >55 | Neg |
| 5 | TRUE | 17–21 | Male | 35–39 | Neg |
| 7 | TRUE | 17–21 | Male | 50–54 | Neg |
| 18 | TRUE | 17–21 | Male | 30–34 | Neg |
| 22 | TRUE | 17–21 | Male | 45–49 | Neg |
| 24 | TRUE | 17–21 | Female | 50–54 | Neg |
| 27 | TRUE | 17–21 | Male | 30–34 | Neg |
| 30 | TRUE | 17–21 | Female | >55 | Neg |
| 60 | TRUE | 17–21 | Male | 50–54 | Neg |
| 71 | TRUE | 17–21 | Female | 35–39 | Neg |
| 91 | TRUE | 17–21 | Male | >55 | Neg |
| 114 | TRUE | 17–21 | Female | 40–44 | Neg |
Additional files
-
MDAR checklist
- https://cdn.elifesciences.org/articles/107153/elife-107153-mdarchecklist1-v1.docx
-
Source data 1
Data used in population-level models, showing step selection coefficients for each individual and each environmental factor.
- https://cdn.elifesciences.org/articles/107153/elife-107153-data1-v1.zip
-
Source code 1
R code used for individual-level step selection functions.
- https://cdn.elifesciences.org/articles/107153/elife-107153-code1-v1.zip
-
Source code 2
R code used for population-level models.
- https://cdn.elifesciences.org/articles/107153/elife-107153-code2-v1.zip