ř

2016-05-14 66 views
0

我正尝试创建两个不同的矩阵之间的相关性(具有p值)(操作分类单元对环境paramenters)中的Rř

第一张表是本

biotic1 biotic2 T1 1.540184 3.080025 T2 1.354927 5.012977 T3 1.449712 4.715981 T4 1.146659 2.442083 X1 1.705184 3.881878 X2 1.182721 3.014836 X3 1.536956 2.636719 X4 1.808025 4.434525 A1 1.132737 2.135737 A2 1.506048 3.114281 A3 1.285308 4.363828 A4 3.008994 7.290423

和第二表

OTU1 OTU2 OTU3 OTU4 OTU5 OTU6 OTU7 OTU8 T1 109 80 175 14 71 46 61 39 T2 102 48 26 8 23 5 35 10 T3 26 19 61 3 68 13 10 29 T4 143 56 9 11 16 13 49 24 X1 70 36 20 15 39 9 26 12 X2 39 33 12 32 15 2 11 3 X3 43 17 2 14 8 2 7 2 X4 160 60 8 26 25 7 9 15 A1 90 73 41 15 22 23 33 7 A2 344 109 18 28 22 13 93 16 A3 65 16 15 9 5 10 18 6 A4 141 140 6 86 18 3 43 4

我已经尝试COR()及其Corr .test()但它只关注关联来自第一个表的值

任何建议?

非常感谢您

˚F

回答

1

目前还不清楚我是什么导致你期望。但是,如果你想执行一个简单的相关性检验,您必须在基质中矢量格式。您可以尝试类似:

cor(c(as.matrix(your_matrix1)), c(as.matrix(your_matrix2))) 

cor.test(c(as.matrix(your_matrix1)), c(as.matrix(your_matrix2))) 

,看看其中一个选项符合您的期望。

然而,使用典型相关分析来探索数据集对我来说更有意义。使用基础R您可以使用:

cancor(matrix1, matrix2) 

你也可以使用一些程序有一系列的工具来解释结果(如library(CCA)

+0

感谢您的回答G. Cocca,我想每一个相关个别OTU与环境参数。 – FranciscoC