У меня есть два фрейма данных, выглядящих так: (1)
GeneID OTU B_F
Zm00001d002422 OTU_3525 B
Zm00001d008299 OTU_3525 B
Zm00001d008956 OTU_3525 B
Zm00001d021674 OTU_8672 B
Zm00001d025402 OTU_8672 B
Zm00001d003553 OTU_211 B
Zm00001d004592 OTU_211 B
Zm00001d004721 OTU_211 B
и (2)
GeneID OTU B_F
Zm00001d002422 OTU_29 F
Zm00001d003793 OTU_29 F
Zm00001d004162 OTU_29 F
Zm00001d019227 OTU_867 F
Zm00001d023804 OTU_867 F
Zm00001d024027 OTU_867 F
Zm00001d008305 OTU_87 F
Zm00001d011642 OTU_87 F
Zm00001d013177 OTU_87 F
Zm00001d003553 OTU_87 F
Единственное отличие состоит в том, что (1) - это B (2)F и имена OTU разные.Фактические файлы намного длиннее.
Некоторые из идентификаторов генов перекрываются между таблицами.
Я хочу получить значения корреляции между OTU таблицы (1) и таблицы (2), поэтому матрицы корреляции результатов должны быть примерно такими:что:
P-val OTU_3525 OTU_8672 OTU_211
OTU_29 0.03 1 1
OTU_867 1 1 1
OTU_87 1 1 0.04
Corr. coeff OTU_3525 OTU_8672 OTU_211
OTU_29 0.3 0 0
OTU_867 1 0 0
OTU_87 0 0 0.2
Есть ли способ в R выполнить это вычисление?