This is an ugly R script to performing a parametric bootstrap approach to simulate no direct effect of the geographical character on the estimated rates as implemented in Pinto-Ledezma et al. (2017).
Pinto‐Ledezma, JN., L. Simon, JAF. Diniz‐Filho & F. Villalobos (2017) The geographical diversification of Furnariides: the role of forest versus open habitats in driving species richness gradients. Journal of Biogeography, 44, 1683-1693.