Как использовать Джулию для вычисления коэффициента корреляции Пирсона с p-значением? - PullRequest
Новая
       3

Как использовать Джулию для вычисления коэффициента корреляции Пирсона с p-значением?

0 голосов
/ 17 ноября 2018

Я ищу помощь для расчета коэффициента корреляции Пирсона с p-значением, используя язык Джулии. Аналогичная функция в Python: scipy.stats.pearson.

Функция Джулии, приведенная ниже, дает мне только корреляцию. Цените вашу помощь / подсказку о p-значении. 

using RDatasets, Statistics
iris = dataset("datasets", "iris");
Statistics.cor(iris.SepalLength, iris.SepalWidth)

1 Ответ

0 голосов
/ 17 ноября 2018

Я не знаю о существующей реализации, но вот двусторонний тест с H0, равным 0, с использованием преобразования Фишера : 

using Distributions

cortest(x,y) =
    if length(x) == length(y)
        2 * ccdf(Normal(), atanh(abs(cor(x, y))) * sqrt(length(x) - 3))
    else
        error("x and y have different lengths")
    end

или использования пакета HypothesisTests.jlНапример: 

using HypothesisTests

OneSampleZTest(atanh(cor(iris.SepalLength, iris.SepalWidth)),
               1, nrow(iris)-3)
...