medelev <- median(elev(dkmodel$cusp))
basins <- dkcrop(dkmodel$cusp, which(elev(dkmodel$cusp) < medelev))
cusps <- dkcrop(dkmodel$cusp, which(elev(dkmodel$cusp) > medelev))
rrate(dkmodel$cusp, basins)
rrate(dkmodel$cusp, cusps)
Run the code above in your browser using DataLab