R中带有p值的多元相关

Question

我正在尝试在R中创建多色相关并获取其p值。我一直在使用以下代码：

library(polycor)
polychor(x,y)

尽管它为我提供了相关的值，但我仍然需要pvalue，有什么方法可以使用此函数获取它？还是应该使用其他功能？

谢谢！

Answer 1

目前尚不清楚您正在寻找什么测试。

但是，如果将std.err设置为TRUE，则会得到二元正态性检验。那是你所追求的吗？

library(polycor)
x <- rnorm(10)
y <- runif(10)
ans <- polychor(x, y, std.err=TRUE)
print(ans)
#
#Polychoric Correlation, 2-step est. = 0.2633 (0.3001)
#Test of bivariate normality: Chisquare = 45.39, df = 80, p = 0.9994

运行?polychor了解更多信息。

您可以使用以下方法提取p值：

pchisq(ans$chisq, ans$df, lower.tail = FALSE)
#0.9994