# Lomolino, M.V. (2000)  Ecology’s most general, yet protean pattern: the species-area relationship. Journal of Biogeography, 27, 17–26

CB<-read.csv("CarribTot.csv")
par(mfrow=c(2,2))
plot(CB)
plot(log(CB,10), xlab="log(Area)",ylab="log(Total species)") 
mod_CB<-lm(log(CB$Tot_butterflies,10) ~ log(CB$Area_square_km,10))
summary(mod_CB)
abline(mod_CB,col="blue",lwd=2)

####

a<-simulate(modelt,10)
plot(CB,pch=15)
for (i in 1:10) points(cbind(CB[1],10^a[i]),pch=3,col="red")