Я использую приведенный ниже код для расчета карты корреляции между двумя наборами данных. Этот код работал нормально, и я получил следующие результаты:! [Введите описание изображения здесь]! [Введите описание изображения здесь] [1].
Я также хотел бы получить другую карту, показывающую, сколько пар использовалось при вычислении каждого пикселя, поэтому я получаю карту N a long с картой корреляции. согласно Полю Хиемстра эта функция дала cor и N:
cor_withN = function(...) {
cor_obj = cor.test(...)
print(sprintf("N = %s", cor_obj$parameter + 2))
return(data.frame(cor = cor_obj$estimate, N = cor_obj$parameter + 2))
}
cor_withN(runif(100), runif(100))
[1] "N = 100"
cor N
cor 0.1718225 100
когда я просто заменил cor на cor_withN, я получил эту ошибку:
Error in cor.test.default(...) : not enough finite observations
Как я могу включить эту функцию в свой код, чтобы получить две карты корреляции и N значений?