Question

我有两个3D（12x12x10）矩阵，它们来自CONN Software中的功能连接分析和.mat格式。每个3D矩阵由12个感兴趣区域的10个单独矩阵组成。一个是考虑休息条件而另一个是任务条件。我想比较FC中执行R中两个3D矩阵之间相关性的差异，但是我不知道如何让R理解我有一个3D矩阵！它混合在一个奇怪的2D矩阵中。使用以下代码：

# Load connectivity matrix
mat<-read.table("R/Matriz/neural", header = FALSE)
View(mat)
r<-corr.test(mat,mat)

尝试计算只有1个值的相关矩阵，我得到了一个完全不同的矩阵：

Call:corr.test(x = mat, y = mat)
Correlation matrix 
       V1    V2    V3    V4    V5    V6    V7    V8    V9   V10   V11   V12
V1   1.00  0.84  0.43 -0.14  0.02 -0.30 -0.20 -0.08 -0.04 -0.20 -0.46 -0.51
V2   0.84  1.00  0.55 -0.03  0.02 -0.23 -0.12 -0.02 -0.04 -0.13 -0.49 -0.50
V3   0.43  0.55  1.00  0.15  0.20 -0.03  0.14  0.35  0.09 -0.08 -0.31 -0.23
V4  -0.14 -0.03  0.15  1.00  0.54  0.45  0.57  0.51  0.23 -0.09  0.20  0.19
V5   0.02  0.02  0.20  0.54  1.00 -0.18  0.04  0.16  0.80  0.12  0.37  0.39
V6  -0.30 -0.23 -0.03  0.45 -0.18  1.00  0.68  0.51 -0.44 -0.31 -0.20 -0.25
V7  -0.20 -0.12  0.14  0.57  0.04  0.68  1.00  0.69 -0.20 -0.11  0.01  0.02
V8  -0.08 -0.02  0.35  0.51  0.16  0.51  0.69  1.00 -0.04 -0.11 -0.13  0.02
V9  -0.04 -0.04  0.09  0.23  0.80 -0.44 -0.20 -0.04  1.00  0.40  0.55  0.60
V10 -0.20 -0.13 -0.08 -0.09  0.12 -0.31 -0.11 -0.11  0.40  1.00  0.45  0.51
V11 -0.46 -0.49 -0.31  0.20  0.37 -0.20  0.01 -0.13  0.55  0.45  1.00  0.87
V12 -0.51 -0.50 -0.23  0.19  0.39 -0.25  0.02  0.02  0.60  0.51  0.87  1.00

Answer 1

我假设你想要两个不同矩阵之间的相关性（与代码不同），并将所有值视为独立数据。您描述数据的方式可能暗示这两者之间的相关性可能不是比较的最佳方式，但如果您想比较，一种方法是展开您的矩阵，因为它们是3D的事实与计算相关他们之间的相关性如果你把它们做成3D，那么R假设变量和结构。

试

r<-corr.test(as.vector(mat),as.vector(mat2))

使用R软件与MatLab的3D矩阵相关联

1 个答案: