Полихорическая корреляция с pvalue в R

Question

Я пытаюсь создать полихорическую корреляцию в R и получить ее значение p.Я использовал следующий код:

library(polycor)
polychor(x,y)

Хотя он предоставляет мне значение корреляции, мне все еще нужно значение, есть ли способ получить его с помощью этой функции?или я должен использовать другую функцию?

Спасибо!

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