I'm using the random forests algorithm in species distribution context. I have presence-absence as a dependent variable, so I decided to work with classification. At this point, I recovered the probability of presence this way :
predict(RF_MODEL, test_, type='prob')[,2]
If I want to compare the results of this model with other models "glm", "gam" and "brt", can we do this or we have to work with regression to have the estimation of occurrence probabilities as the other models?