Multiple "spaces": using wildlife surveillance, climatic variables, and spatial statistics to identify and map a climatic niche for endemic plague in California, U.S.A.

Date repository last updated: 2024-11-05

Authors

• Ian D. Buller¹ - Corresponding Author - ORCID
• Gregory M. Hacker² - ORCID
• Mark G. Novak²
• James R. Tucker²
• A. Townsend Peterson³ - ORCID
• Lance A. Waller⁴ - Senior Author - ORCID

1. Environmental Health Sciences, James T. Laney School of Graduate Studies, Emory University, Atlanta, GA, 30322, USA
2. California Department of Public Health – Vector-Borne Disease Section, Sacramento, CA, 95814, USA
3. Biodiversity Institute, University of Kansas, Lawrence, KS, 66045, USA
4. Biostatistics and Bioinformatics, Rollins School of Public Health, Emory University, Atlanta, GA, 30322, USA

Project Details

We combine two analytic concepts (ecological niche modeling and spatial point processes) to identify plague-suitable climates and map their locations within California. Our approach uses serological samples from coyotes (Canis latrans), a sentinel species for plague, a zoonotic disease caused by the gram-negative bacterium, Yersinia pestis. The approach accounts for spatially heterogeneous sampling effort by considering locations of both seropositive (case) and seronegative (control) coyotes. Modification of Chapter 3 from Ian Buller's Doctoral Dissertation from Emory University (see Emory Theses and Dissertations Repository).

Project Timeframe

Time	Event
1981-2010	The Oregon State University Parameter-elevation Regression on Independent Slopes Model (PRISM) 30-year average climate normals at a 2.5 arcminute (~16 km2) resolution (see data availability section below).
1983-2015	The California Department of Public Health – Vector-Borne Disease Section (CDPH-VBDS) digitized coyote blood samples tested for Y. pestis antibodies (see data availability section below).
October 2016	Project Initiation
December 2020	The envi package in R published in the Comprehensive R Archive Network
November 2023	Initial manuscript submission to Spatial and Spatio-temporal Epidemiology for peer-review
October 2024	Manuscript accepted by Spatial and Spatio-temporal Epidemiology
November 2024	Manuscript published in Spatial and Spatio-temporal Epidemiology

R Scripts Included In This Repository

This repository includes R scripts used to calculate a spatial relative risk function in 'covariate space' and render the figures found in the following peer-reviewed manuscript:

Buller ID, Hacker GM, Novak MG, Tucker JR, Peterson AT, Waller LA. (2024) Multiple "spaces": using wildlife surveillance, climatic variables, and spatial statistics to identify and map a climatic niche for endemic plague in California, U.S.A. Spatial and Spatio-temporal Epidemiology, 51:100696 DOI: 10.1016/j.sste.2024.100696 PMID: XXXXXX

R Script	Description
Paths.R	Example paths for data. Must modify for your own system before beginning.
Preparation.R	Load settings, prepare data, run the log RR model, and process the results to generate the figures.
Figure1.R	Generate Figure 1
Figure2.R	Generate Figure 2
Figure3.R	Generate Figure 3
Supplemental1.R	Generate Supplemental Figure 1
Supplemental2.R	Generate Supplemental Figure 2
	Supplemental Figure 3 hand-generated, no code used
Supplemental4.R	Generate Supplemental Figure 4
	Supplemental Figure 5 hand-generated, no code used
Supplemental6.R	Generate Supplemental Figure 6
Supplemental7.R	Generate Supplemental Figure 7
Supplemental8.R	Generate Supplemental Figure 8

The repository also includes the code and resources to create the project hexagon sticker.

Getting Started

• Step 1: You must download the elevation BIL file at 4-km resolution from the PRISM data portal
• Step 2: Save the BIL file to the 'data' directory in this repository
• Step 3: Set your own file paths to the data in the 'Paths.R' file

Data Availability

Wildlife plague surveillance data from the California Department of Public Health – Vector-Borne Disease Section (CDPH-VBDS) available upon request to CDPH-VBDS Infectious Diseases Branch - Surveillance and Statistics Section. The Oregon State University Parameter-elevation Regression on Independent Slopes Model (PRISM) 30-year average climate normals are available through the prism package in R or directly from the PRISM data portal.

Questions?

For questions about the manuscript please e-mail the corresponding author Dr. Ian D. Buller.