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.
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.