Una forma muy simple es usar el argumento use = "pairwise.complete.obs"
en la función cor()
. Con ese argumento se calcula la correlación solamente en aquellas filas en las que las dos columnas tienen valores válidos. El valor defecto es use = "everything"
y en ese caso se conservan los NA
y las operaciones aritméticas de cor
los reciben y propagan. El resultado es que el coeficiente de correlación es NA
siempre que haya al menos un valor NA
.
La ventaja que le veo es que no necesitas crear (y rastrear) estructuras de datos adicionales.
En tu caso:
cor(qog$bci_bci, qog$undp_hdi, use = "pairwise.complete.obs")
En la ayuda de cor
bajo el apartado Details puedes encontrar todas las opciones (bastante rudimentarias) para manejar casos perdidos con esta función. Lamentablemente no está traducida al español.
PD: suscribo plenamente el comentario de @Patricio Moracho sobre la importancia de diagnosticar los NA
presentes en los datos. Eliminarlos es una solución técnicamente simple y en algunos casos analíticamente correcta. Pero conlleva el riesgo de introducir sesgos a los datos u omitir cuestiones importantes como que los NA
provienen de un procesos generador de datos diferente a los datos válidos.