2011-09-26 17 views

Respuesta

13

Trate cor(A,B, use = "pairwise.complete.obs"). Eso ignorará las NA en tus observaciones.

Para ser estadísticamente riguroso, también debe mirar el n. ° de entradas faltantes en sus datos y observar si se cumple la aserción faltante al azar.

Edición 1: Eche un vistazo a ?cor para ver otras opciones para el parámetro use.

4

Usted podría considerar el uso de la función Rcorr en el paquete Hmisc.

Es muy rápido, y solo incluye observaciones completas por pares. El objeto devuelto contiene una matriz

  1. de las puntuaciones de correlación
  2. con el número de observación utilizada para cada valor de correlación
  3. de un valor de p para cada correlación

un código de ejemplo está disponible here:

Cuestiones relacionadas