“…ENM were built in R environment (R Core Team, 2016), using the ‘biomod2’ modelling package, which permits to obtain the so‐called ‘Ensemble Models’ (EMs), the models obtained by merging single ENMs calculated by the different algorithms available within (Thuiller, Georges, & Engler, 2016). The ‘BIOMOD_EnsembleModeling’ function was used for this purpose, with single models obtained from Generalized Linear Models (set to type = ‘quadratic’ and interaction level = 3), Multiple Adaptive Regression Splines (set to type = ‘quadratic’ and interaction level = 3), Generalized Boosting Models (sometimes named BRTs, with number of trees set to 10,000, interaction depth = 3 and 10‐fold cross‐validation) and Maxent (Maxent.Phillips, maximum interactions = 5,000 and betamultiplier = 2), an approach often used to encompass and take advantage of different modelling techniques (Cerasoli et al., 2019; D'Alessandro, Iannella, Frasca, & Biondi, 2018; Iannella, D’Alessandro, & Biondi, 2018). Five sets of 1,000 pseudoabsences each were generated through the ‘sre’ (Surface Range Envelope, set to 0.05) algorithm, which calculates a linear envelope on the basis of selected predictors and selects pseudoabsences outside the set quantile, for all the reasons reported in Iannella, Cerasoli, D’Alessandro, Console, and Biondi (2018) and all the references within.…”