-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathFigure 5
250 lines (155 loc) · 7.86 KB
/
Figure 5
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
###Spearman heatmap
require(ggpubr)
require(tidyverse)
require(Hmisc)
require(corrplot)
library("pheatmap")
library(ggpubr)
matrix_a <- Taxa
matrix_b <- Open
bb <- cor(matrix_b, matrix_a, method = "spearman", use = "pairwise.complete.obs")
nba_matrix <- data.matrix(bb)
pheatmap(nba_matrix, cluster_cols=FALSE)
###Spearman correlation
require(ggpubr)
require(tidyverse)
require(Hmisc)
require(corrplot)
library("pheatmap")
library(ggpubr)
matrix_a <- Taxa
matrix_b <- Open
##Aeromonadaceae
M1Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$sidt2, method=c("spearman"))
M1Aeromonadaceae
M2Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$rabep2, method=c("spearman"))
M2Aeromonadaceae
M3Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$ffar2, method=c("spearman"))
M3Aeromonadaceae
M4Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$aifm3, method=c("spearman"))
M4Aeromonadaceae
M5Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$ub, method=c("spearman"))
M5Aeromonadaceae
M6Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$dspa2b, method=c("spearman"))
M6Aeromonadaceae
M7Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$tmem38b, method=c("spearman"))
M7Aeromonadaceae
M8Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$manf, method=c("spearman"))
M8Aeromonadaceae
M9Aeromonadaceae <- cor.test(matrix_a$Aeromonadaceae, matrix_b$prmt3, method=c("spearman"))
M9Aeromonadaceae
####Bifidobacteriaceae
M1Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$sidt2, method=c("spearman"))
M1Bifidobacteriaceae
M2Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$rabep2, method=c("spearman"))
M2Bifidobacteriaceae
M3Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$ffar2, method=c("spearman"))
M3Bifidobacteriaceae
M4Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$aifm3, method=c("spearman"))
M4Bifidobacteriaceae
M5Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$ub, method=c("spearman"))
M5Bifidobacteriaceae
M6Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$dspa2b, method=c("spearman"))
M6Bifidobacteriaceae
M7Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$tmem38b, method=c("spearman"))
M7Bifidobacteriaceae
M8Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$manf, method=c("spearman"))
M8Bifidobacteriaceae
M9Bifidobacteriaceae <- cor.test(matrix_a$Bifidobacteriaceae, matrix_b$prmt3, method=c("spearman"))
M9Bifidobacteriaceae
###Desulfovibrionaceae
M1Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$sidt2, method=c("spearman"))
M1Desulfovibrionaceae
M2Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$rabep2, method=c("spearman"))
M2Desulfovibrionaceae
M3Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$ffar2, method=c("spearman"))
M3Desulfovibrionaceae
M4Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$aifm3, method=c("spearman"))
M4Desulfovibrionaceae
M5Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$ub, method=c("spearman"))
M5Desulfovibrionaceae
M6Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$dspa2b, method=c("spearman"))
M6Desulfovibrionaceae
M7Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$tmem38b, method=c("spearman"))
M7Desulfovibrionaceae
M8Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$manf, method=c("spearman"))
M8Desulfovibrionaceae
M9Desulfovibrionaceae <- cor.test(matrix_a$Desulfovibrionaceae, matrix_b$prmt3, method=c("spearman"))
M9Desulfovibrionaceae
###Lactobacillaceae
M1Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$sidt2, method=c("spearman"))
M1Lactobacillaceae
M2Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$rabep2, method=c("spearman"))
M2Lactobacillaceae
M3Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$ffar2, method=c("spearman"))
M3Lactobacillaceae
M4Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$aifm3, method=c("spearman"))
M4Lactobacillaceae
M5Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$ub, method=c("spearman"))
M5Lactobacillaceae
M6Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$dspa2b, method=c("spearman"))
M6Lactobacillaceae
M7Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$tmem38b, method=c("spearman"))
M7Lactobacillaceae
M8Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$manf, method=c("spearman"))
M8Lactobacillaceae
M9Lactobacillaceae <- cor.test(matrix_a$Lactobacillaceae, matrix_b$prmt3, method=c("spearman"))
M9Lactobacillaceae
###Mycoplasmataceae
M1Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$sidt2, method=c("spearman"))
M1Mycoplasmataceae
M2Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$rabep2, method=c("spearman"))
M2Mycoplasmataceae
M3Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$ffar2, method=c("spearman"))
M3Mycoplasmataceae
M4Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$aifm3, method=c("spearman"))
M4Mycoplasmataceae
M5Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$ub, method=c("spearman"))
M5Mycoplasmataceae
M6Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$dspa2b, method=c("spearman"))
M6Mycoplasmataceae
M7Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$tmem38b, method=c("spearman"))
M7Mycoplasmataceae
M8Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$manf, method=c("spearman"))
M8Mycoplasmataceae
M9Mycoplasmataceae <- cor.test(matrix_a$Mycoplasmataceae, matrix_b$prmt3, method=c("spearman"))
M9Mycoplasmataceae
###Pseudomonadaceae
M1Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$sidt2, method=c("spearman"))
M1Pseudomonadaceae
M2Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$rabep2, method=c("spearman"))
M2Pseudomonadaceae
M3Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$ffar2, method=c("spearman"))
M3Pseudomonadaceae
M4Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$aifm3, method=c("spearman"))
M4Pseudomonadaceae
M5Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$ub, method=c("spearman"))
M5Pseudomonadaceae
M6Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$dspa2b, method=c("spearman"))
M6Pseudomonadaceae
M7Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$tmem38b, method=c("spearman"))
M7Pseudomonadaceae
M8Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$manf, method=c("spearman"))
M8Pseudomonadaceae
M9Pseudomonadaceae <- cor.test(matrix_a$Pseudomonadaceae, matrix_b$prmt3, method=c("spearman"))
M9Pseudomonadaceae
###Streptococcaceae
M1Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$sidt2, method=c("spearman"))
M1Streptococcaceae
M2Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$rabep2, method=c("spearman"))
M2Streptococcaceae
M3Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$ffar2, method=c("spearman"))
M3Streptococcaceae
M4Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$aifm3, method=c("spearman"))
M4Streptococcaceae
M5Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$ub, method=c("spearman"))
M5Streptococcaceae
M6Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$dspa2b, method=c("spearman"))
M6Streptococcaceae
M7Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$tmem38b, method=c("spearman"))
M7Streptococcaceae
M8Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$manf, method=c("spearman"))
M8Streptococcaceae
M9Streptococcaceae <- cor.test(matrix_a$Streptococcaceae, matrix_b$prmt3, method=c("spearman"))
M9Streptococcaceae
#Figure 5