Data and code from: Local exclusion and regional decline of an endemic Galápagos tree species (Psidium galapageium) by an invasive relative (P. guajava)
Data files
Jul 20, 2026 version files 28.51 MB
-
exclusion.zip
7.89 KB
-
niche_overlap.zip
466.54 KB
-
R_analyses.R
51.26 KB
-
README.md
6.38 KB
-
SDMs.zip
27.98 MB
Abstract
Invasive species can interact with native relatives in a variety of ways that may jeopardize their long-term coexistence. Here we show that interactions with an invasive species of guava (Psidium guajava) appear to be driving the local exclusion and regional decline of guayabillo (Psidium galapageium), a tree species endemic to the Galápagos archipelago. We find evidence consistent with recent historic exclusion of guayabillo from the highlands of San Cristóbal Island, signatures of ongoing demographic decline in sympatric populations at lower elevations, and evidence suggesting that the four coinhabited islands represent points along a time series of regional decline, with the extent of guayabillo decline depending on the date that guava was introduced to each island. Based on these results, we use the percentage of guava cover surrounding guayabillo populations to target populations that are at imminent risk of exclusion to aid in prioritizing management targets.
General Information
This data_availability directory contains all of the scripts and results necessary to recreate the analyses and plots used in the paper "Reatini et. al. 2026. Local exclusion and regional decline of an endemic Galápagos tree species (Psidium galapageium) by an invasive relative (P. guajava)."
Below is a description of the results files, the meaning of the column names within those files, and the corresponding analyses and plots for which they were used. A description of the scripts used to carry out the analyses follows.
Results Files
Below is a brief description of each of the results files and the figures they were used to produce:
The SDM subdirectory (SDMs.zip) contains the output of the final ensemble models for both species, including the species distribution models named gal_ensemble_cont.tif (continuous), gal_ensemble_bin.tif (binary), gua_ensemble_cont.tif (continuous),
gua_ensemble_bin.tif (binary). Presence data used to generate the SDMs can be found in galapageium.csv and guayava.csv; hyperparameters used to generate the models can be found in all_hyperparameters.xlsx, and model performance results files for both species can be found in gal_all_models_performance.csv and gua_all_models_performance.csv. All SDMs and associated results were generated via R_analyses.R
The niche_overlap subdirectory (niche_overlap.zip) contains the input data and results for all niche comparisons between species. Complete_distributions.csv includes all presence data for both species. gal_background_points.csv and
gua_background_points.csv contain all background points used in SDMs and were used to define the background environmental space on the islands. favorability.csv contains the difference in suitability scores (from the SDMs) for each presence location. All input data and associated niche overlap and comparison analyses were generated via R_analyses.R
The exclusion subdirectory (exclusion.zip) contains the input data for all exclusion analyses. favorability_ecos.csv contains the landcover classification type for each presence locality located with predicted sympatry. gal_CMCR_quadrats.csv and gal_GO_quadrats.csv contain the demographic structure data within the two sympatric and allopatric study sites. gal_fruit_set_perflower.csv and gua_fruit_set_perflower.csv
contain the input data for per-flower rate of fruit set for both species.
gal_seed_set.csv and gua_seed_set.csv contain the seed set input data used to assess reproductive fitness in sympatry. gal_decline.csv contains the input data on the relative position of P. galapageium populations within predicted sympatry, used to assess landscape-level patterns of exclusion. All of these input data and the associated results were processed and generated via R_analyses.R
Column Names
galapageium.csv and guayava.csv
- long - longitude
- lat - latitude
- pr_ab - presence classification (0 or 1)
all_hyperparameters.xlsx
- species - species for SDM model
- model - model type
- regmult - multiplier
- classes - classes of response modeled
- mtry - values modeled for "mtry" variable
- size - values modeled for "size" variable
- decay - values modeled for "decay" variable
gal_all_models_performance.csv and gal_all_models_performance.csv
- model_ID - model number
- model - model type
- threshold - threshold type for binary prediction of presence
- thr_value - threshold value for binary prediction of presence
- n_presences - number of presence points
- n_absences - number of absence points
- the following columns are mean values and standard deviations for the following statistics: TPR, TNR, SORENSEN, JACCARD, FPB, OR, TSS, AUC, BOYCE, IMAE
- regmult - value used for the "regmult" variable
- classes - value used for the "classes" variable
- mtry size - value used for the "mtry" variable
- decay - value used for the "decay" variable
complete_distributions.csv
- ID - individual ID
- source - source of datum
- species - species of datum
- island - island of occurrence
- Date - date collected
- latitude - latitude
- longitude - longitude
- altitude - elevation (m)
favorability.csv
- OBJECTID - point ID number
- species - species of occurrence
- island - island of occurrence
- dd_long - longitude
- dd_lat - latitude
- favorability_difference - difference in favorability score
gal_background_points.csv and gua_background_points.csv
- x - longitude
- y - latitude
- pr_ab - presence type
- .part - number of partition
- the following columns are the values extracted at the background point for the following variables: mean_diurnal_rn, temp_annual_rng, annual_precip, precip_coldest_qtr, dem, ecos
favorability_ecos.csv
- fav_species - species of occurrence
- FID_Galapagos_Agroecosystems_LandCover2018 - ID number of layer
- Level_1 - level 1 classification of land cover
- Level_2 - level 2 classification of land cover
- Level_3 - level 3 classification of land cover
- Level_4 - level 4 classification of land cover
- Level_5 - level 5 classification of land cover
- Island - island of occurrence
- Area_m2 - area in square meters
- Shape_Length - perimeter of polygon
- Shape_Area - area of polygon
gal_CRCM_quadrats.csv and gal_GO_quadrats
- height - height of plant in meters
- species - species of plant
- site - site ID
- quadrat - quadrat number
- type - type of site
gal_fruit_set_perflower.csv and gua_fruit_set_perflower.csv
- flower - flower number
- site - site ID
- type - type of site
- individual - individual ID
- branch - branch ID
- fruitset - fruitset condition (0 or 1)
gal_seed_set.csv and gua_seed_set.csv
- seeds - number of seeds
- species - species of fruit
- site - site ID
- type - type of site
- direction - side of island of site pair
gal_decline.csv
- OBJECTID - ID number of point
- species - species of point
- dd_long - longitude
- dd_lat - latitude
- Island - island of point
- centroid_to_edge - distance from centroid to edge of sympatry
- point_to_edge - distance from point to edge of sympatry
- normalized_centroid_to_point - relative position within sympatry
This dataset contains occurrence, population fitness, and geospatial data collected from field sites in the Galapagos Islands, Ecuador, for invasive guava (Psidium guajava) and endemic guayabillo (Psidium galapageium). The occurrence data includes collections by our team and known collections in public herbaria, as well as species abundance data collected from quadrats in field sites on San Cristobal Island. The population fitness data include measurements of fruit set and seed set from field sites on San Cristobal Island. The geospatial analyses - including species distribution models - were generated from the occurrence data described above. All statistical analyses were performed using scripts contained here.
