Função de aceleração que usa qual dentro de uma chamada solta em R

Eu tenho dois vetorese eg. Eu quero saber para cada elemento eme a porcentagem de elementos emg que são menores. Uma maneira de implementar isso em R é:

set.seed(21)
e <- rnorm(1e4)
g <- rnorm(1e4)
mf <- function(p,v) {100*length(which(v<=p))/length(v)}
mf.out <- sapply(X=e, FUN=mf, v=g)

Com grandee ougIsso leva muito tempo para ser executado. Como posso alterar ou adaptar este código para que isso seja mais rápido?

Note omf função acima é baseado no código domess função no pacote de dismo.

questionAnswers(1)

yourAnswerToTheQuestion