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