Korelacja i odtwarzalność wyników

[pbiecek]

Krótka uwaga do Przewodnika:
Na str. 273 jest ciekawy przykład: test korelacji Pearsona vs. test korelacji Spearmana.
Ponieważ dobrą praktyką jest odtwarzalność wyników to warto podać ziarno generatora liczb losowych.

-------------------------------- R:
set.seed(3133); x <- rlnorm(50);
set.seed(3138); y <- x+rlnorm(50);
cor.test(x,y);

cor.test(x,y,method="spearman");

Może warto w tym miejscu zwrócić uwagę na:
[] gdy występują rangi wiązane test Spearmana nie może obliczyć wartości dokładnej p-wartości.
[] alternatywy - coin::spearman_test, pspearman::spearman.test.
[*] test - energy::dcor.ttest. Wartości mogą się powtarzać, wykrywa on nie tylko liniową korelacje ale także nieliniowe związki.