Вычисление корреляции между матрицей — R и Python возвращают разные результаты

Предположим, что матрицы X и Y имеют размер 2x3 и 2x2 соответственно. Функция cor в R возвращает матрицу 3x2, а функция numpy.corrcoef в Python возвращает матрицу 5x5. Примеры ниже:

Р:

X<-matrix(c(0.2,0.5,0.1,0.7,0.5,0.3), nrow=2, ncol=3)
Y<-matrix(c(0.2,0.3,0.6,0.7), nrow=2)
cor(X,Y)
     [,1] [,2]
[1,]    1    1
[2,]    1    1
[3,]   -1   -1

Питон:

X = np.array([[0.2,0.5], [0.1, 0.7], [0.5,0.3]], ndmin=2).T
Y = np.array([[0.2,0.3],[0.6,0.7]], ndmin=2).T
corr = np.corrcoef(X, Y, rowvar=False)

array([[ 1.,  1., -1.,  1.,  1.],
       [ 1.,  1., -1.,  1.,  1.],
       [-1., -1.,  1., -1., -1.],
       [ 1.,  1., -1.,  1.,  1.],
       [ 1.,  1., -1.,  1.,  1.]])

Как заставить python возвращать матрицу 3x2, как в R? Или как мне выбрать правильные значения в матрице Python 5x5, чтобы они соответствовали результату R?

Возможно, это делает как cor(cbind(X, Y))

akrun 18.03.2022 22:19

Я предполагаю, что причина упоминается здесь. Возможный дубликат

akrun 18.03.2022 22:44
Почему в Python есть оператор "pass"?
Почему в Python есть оператор "pass"?
Оператор pass в Python - это простая концепция, которую могут быстро освоить даже новички без опыта программирования.
Некоторые методы, о которых вы не знали, что они существуют в Python
Некоторые методы, о которых вы не знали, что они существуют в Python
Python - самый известный и самый простой в изучении язык в наши дни. Имея широкий спектр применения в области машинного обучения, Data Science,...
Основы Python Часть I
Основы Python Часть I
Вы когда-нибудь задумывались, почему в программах на Python вы видите приведенный ниже код?
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
LeetCode - 1579. Удаление максимального числа ребер для сохранения полной проходимости графа
Алиса и Боб имеют неориентированный граф из n узлов и трех типов ребер:
Оптимизация кода с помощью тернарного оператора Python
Оптимизация кода с помощью тернарного оператора Python
И последнее, что мы хотели бы показать вам, прежде чем двигаться дальше, это
Советы по эффективной веб-разработке с помощью Python
Советы по эффективной веб-разработке с помощью Python
Как веб-разработчик, Python может стать мощным инструментом для создания эффективных и масштабируемых веб-приложений.
0
2
31
1
Перейти к ответу Данный вопрос помечен как решенный

Ответы 1

Ответ принят как подходящий

В R, когда x и y являются матрицами, cor(x,y) вернет корреляцию столбцов x (n = 3) со столбцами y (n = 2). В python вы можете нарезать результат np.corrcoef(), используя правильные индексы, которые в данном случае равны 3 и 2, для x (строки) и y (столбцы) соответственно результата.

np.corrcoef(X,Y)[0:3,0:2]

array([[ 1.,  1.],
       [ 1.,  1.],
       [-1., -1.]])

Другие вопросы по теме