Estoy replicando la correlación de spearman
en R
, en conjunto con el Rcpp
para usar las funciones de C++
, pero tengo un error consistente:
Cada vez que tomo matrices con filas mayores a las (alrededor) de 45 - 50 mil, mi correlación (delimitada teóricamente entre -1 y 1) se distorsiona terriblemente y puede ir tanto para cifras menores a -1 como mayores a 1. He corroborado que cada vez que uso matrices con filas menores al umbral que les comento, todo va bien. Supongo debe haber algún límite de trabajo en mi computadora (16 GB de RAM), ya que no encuentro otra explicación.
Solo para comentarles, una matriz de 60 000 filas, la recorto a 45 000 y el resultado de mi función Rcpp
cae dentro del rango normal de una correlación. Le permito trabajar con todas sus files y termina siendo -18.
Si alguien tuviera algún comentario le agradecería muchísimo.