2011-09-26 19 views

Odpowiedz

13

Wypróbuj cor(A,B, use = "pairwise.complete.obs"). To zignoruje NA w twoich obserwacjach.

Aby być statystycznie rygorystycznym, powinieneś także spojrzeć na # brakujących wpisów w swoich danych i sprawdzić, czy brakujące losowo założenie zawiera.

Edytuj 1: spójrz na ?cor, aby zobaczyć inne opcje dla parametru use.

4

Możesz rozważyć użycie funkcji rcorr w pakiecie Hmisc.

Jest bardzo szybki i zawiera tylko pełne sparowane obserwacje. Zwrócona obiekt zawiera matrycę

  1. punktów korelacji
  2. numerem obserwacji stosowane dla każdej wartości korelacji
  3. od wartości p dla każdego korelacji

przykładowy kod jest dostępny here:

Powiązane problemy