Instruction on how to run the model supplied with the paper: Koshkina et al. "Integrated Species Distribution Models: Combining presence-only data and presence-absence data with imperfect detection"

Modified and used also in the paper:

Schank, C. J., M. V. Cove, M. J. Kelly, E. Mendoza, G. O’Farrill, R. Reyna-Hurtado, N. Meyer, C. A. Jordan, J. F. González-Maya, D. J. Lizcano, R. Moreno, M. T. Dobbins, V. Montalvo, C. Sáenz-Bolaños, E. C. Jimenez, N. Estrada, J. C. Cruz Díaz, J. Saenz, M. Spínola, A. Carver, J. Fort, C. K. Nielsen, F. Botello, G. Pozo Montuy, M. Rivero, J. A. de la Torre, E. Brenes-Mora, O. Godínez-Gómez, M. A. Wood, J. Gilbert, and J. A. Miller. 2017. Using a novel model approach to assess the distribution and conservation status of the endangered Baird’s tapir. Diversity and Distributions 23:1459–1471.

Steps:

1. Download all the covariate spatial layers, PB and SO data files and codes inside the R folder
2. Unarchive all files and save them in the same folder
3. Set this folder as working directory in R
4. Run file PBSO-yb_edited_Diego_v2
5. Include the name of your species of interest to get GBIF data
6. Chech the area of interest to download worldclim data