Question

我正在寻求使用Julia语言来计算p值的Pearson相关系数的帮助。 Python中的类似函数是scipy.stats.pearson。

下面的Julia函数仅给我相关性。感谢您对p值部分的帮助/提示。

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

Answer 1

我不知道现有的实现，但这是使用Fisher transformation进行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)

如何使用Julia来计算p值的皮尔逊相关系数？

1 个答案: