[R-br] sugestao de melhorias em codigo?

Prezados, existe algo obvio (em termos de R, visto que gostaria de evitar o uso de codigo compilado) que eu esteja esquecendo em termos de otimizacao para a seguinte operacao? Ela se resume a: R[i,j] = mean( X[i,j] < Y[,j]) O CMR segue: set.seed(1) X <- matrix(rnorm(100), nc=5) Y <- matrix(rnorm(50), nc=5) f = function(i, x, y){ n = ncol(x) sapply(1:n, function(j, v, ref) mean(ref[,j] > v[j]), x[i,], y) } R = t(sapply(1:nrow(X), f, X, Y) O problema que surge para mim refere-se ao tempo de execucao quando tenho q lidar com X com dimensoes 5e6 x 1e3.... pois, naturalmente, a solucao acima levara' um bom tempo para finalizar... obrigado, benilton