2011-09-26 23 views

risposta

13

Prova cor(A,B, use = "pairwise.complete.obs"). Questo ignorerà le NA nelle tue osservazioni.

Per essere statisticamente rigorosi, è necessario anche esaminare il numero di voci mancanti nei dati e verificare se il presupposto mancante a caso è valido.

Modifica 1: Dai un'occhiata a ?cor per visualizzare altre opzioni per il parametro use.

4

È possibile utilizzare la funzione rcorr nel pacchetto Hmisc.

È molto veloce e include solo osservazioni complete a coppie. L'oggetto restituito contiene una matrice

  1. di punteggi di correlazione
  2. con il numero di osservazione utilizzato per ogni valore di correlazione
  3. di un p-valore per ogni correlazione

qualche esempio di codice è disponibile here:

Problemi correlati