Data organization_OTU

Data input

library(readxl)
OTU_xulie_zhenhe <- read_excel("D:/(R) plots collection/R markdown/Data organization_OTU_gene/OTU xulie zhenhe.xlsx")
OTU_xulie_zhenhe

## # A tibble: 3,250 x 1
##    Sequence                                                                     
##    <chr>                                                                        
##  1 >OTU1 MJ12_100225                                                            
##  2 AGCTCCAACAGCGTATATAAAAGTTGTTGCGATTAAACAGCTCGTAGTCGGACTTCTGCTGAGGATGACCGGTCCG~
##  3 >OTU10 MJ12_7529                                                             
##  4 AGCTCCAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTCTGACAGGAGCGACCGGTCCG~
##  5 >OTU100 MJ12_41755                                                           
##  6 AGCTCCAATAGCGTATATTAAAGTTGTTGCGGTTAAAAAGCTCGTAGTCGAATTTCAATTATTGTCAGTGATTTAG~
##  7 >OTU1000 MJ4_9516                                                            
##  8 AGCTCCAATAGTGTATCTTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGAATTTCTGGCAATACAGCTGATCTTT~
##  9 >OTU1001 MJ4_10787                                                           
## 10 AGCTCTAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTTAGAGCTTGACAGACGGTCCG~
## # ... with 3,240 more rows

otu_table_5_12_17_29_top100 <- read_excel("D:/(R) plots collection/R markdown/Data organization_OTU_gene/otu_table_5,12,17,29_top100.xlsx")
otu_table_5_12_17_29_top100

## # A tibble: 100 x 5
##    `OTU ID`   MJ5  MJ12  MJ17  MJ29
##    <chr>    <dbl> <dbl> <dbl> <dbl>
##  1 OTU805     126     0   236 32281
##  2 OTU646     392    77  1293 18152
##  3 OTU160       0     0 13520    96
##  4 OTU151       0 11984    28     1
##  5 OTU472      11 11814    14    27
##  6 OTU2         3 10379    35   449
##  7 OTU289       0   101  9982    58
##  8 OTU286      71  3372  4406   147
##  9 OTU439       0  7811    11    46
## 10 OTU1267      0  5897    37  1321
## # ... with 90 more rows

Data organization

a = data.frame(aa = c("a", "b"))
b = data.frame(bb = c("a", "b"))
z = nrow(OTU_xulie_zhenhe)
z = z/2
for (i in 1:z) {
  a[i,1] = OTU_xulie_zhenhe[2*i-1,1]
  b[i,1] = OTU_xulie_zhenhe[2*i,1]
}
nrow(a)

## [1] 1625

nrow(b)

## [1] 1625

head(a)

##                   aa
## 1  >OTU1 MJ12_100225
## 2   >OTU10 MJ12_7529
## 3 >OTU100 MJ12_41755
## 4  >OTU1000 MJ4_9516
## 5 >OTU1001 MJ4_10787
## 6  >OTU1002 MJ4_1279

head(b)

##                                                                                                                                                                                                                                                                                                                                                                                                          bb
## 1           AGCTCCAACAGCGTATATAAAAGTTGTTGCGATTAAACAGCTCGTAGTCGGACTTCTGCTGAGGATGACCGGTCCGCCCTCTGCGTGAGCATCTGGCTTAGCCTTAGCATTTTCTTGGGCAGCGATATTGCACTTGACTGTGTGGTGTGGTACACAAGGCTTTTAGCTTGTGGAAATCAGAGCGCTTCCAGCAGGCATATGTCTTGAATACGTTACCATTGAATAATAAGATAAGACTTTATGTTTTATTGTTGTTGGTTTTTAGAATTGGGGTAATGCTTGATAGGGATACTTAGGGGCATTTCTGTTTAACTGTCAGAGGTGAACTTCTTTGATTTGTTACAGACAAACTATTGCTAAAGCATTTGCCAATGATGTTTTCC
## 2        AGCTCCAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTCTGACAGGAGCGACCGGTCCGAGCCTTCGGGTTCGTGTACTTTGTGTTGTCTCTAGTCATCCTTATGGAGAGCGTGTCTTTCATTAATTTGTTGGGCATGGGATCCATATCGTTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGGCAATCGCTGTGAATATATTAGCATGGAATAATAGGATAGGACTTTGGTACTATTTTGTTGGTTTGCATACCAAGGTAATGATTAAGAGGGACAGTTGGGGGTATTCGTATTTAATTGTCAGAGGTGAAATTCTTGGATTTATGAAAGACGAACTACTGCGAAAGCATTTACCAAGGATGTTTTCA
## 3                                                  AGCTCCAATAGCGTATATTAAAGTTGTTGCGGTTAAAAAGCTCGTAGTCGAATTTCAATTATTGTCAGTGATTTAGTCTATTAGGCTCAAGTTGCTGGCATAATTTATCTCATCACGAGATAATATTACTTTGAGGAAATTAGAGTGTTTAAAGCAAGCGATTGTTTTGAATAATTTAGCATGGAATAATGATATAGGACCTTATTTTATATTCATTGGTTTTACGAAATGAGGTAATGATTAATGGAAACAGACAGGGGCATTCGTATTTAACTGTCAGAGGTGGAATTCTTGGATTTGTTAAAGACGATCAACTGCGAAAGCGTTTGCCAGGAATGTCTTCA
## 4                                   AGCTCCAATAGTGTATCTTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGAATTTCTGGCAATACAGCTGATCTTTTCGAGAAGCACGTATATGTCATCCGCATTGGAAGCACATTCGATAGCAATATCAGTGTGTGGAAGGTGCATTTTACTTTGAAGAAATTAGAGTGTTCCAGGCAGGCTATGCTTGAATATTACAGCATGGAATAATGGAATAGGACAAGTTCATTTTATTGGTTATTGAATTTGTAATGATTAACAGGGACAGTTGGGGGCATTCGTATTGGGTAGTCAGAGGTGAAATTCTTAGATTTACTCAAGACGCACTAATGCGAAAGCATTTGCCAAGGATGTTTTCA
## 5 AGCTCTAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTTAGAGCTTGACAGACGGTCCGCTTGTTTTGCAGGTGAGTACTGTTCCTGTTGGCTTTATTCTTTTCTGGGGGCCGGCGTTCCCTTAATCAGGAGCGTCATTGGCCGTACAGATCGTTTACCTTGAAAAAATTAGAGTGTTCAAAGCAGGCTTCCGCAATTTTTGAATACATTAGCATGGAATAATGGAATAGGACTTTGGTCTTATTTTGTTGGTTTAAAGGCCGAGGTAATGATTAATAGGGATAGTTGGGGGCGTTCATATTCCGTTGTCAGAGGTGAAATTCTTGGATTTACGGAAGATGAACTACTGCGAAAGCATTCGCCAAGGATGTTTTCA
## 6                                AGCTCCAGTAGCGTATATTAAAGTTGTTGCAGTTAAAAAGCTCGTAGTTGAACTTCTGTGTTGCGATGCTGAGTCCTCCGGGACAAACTTTGCGACACATCCATCTGAATAGCATTTCGGCAATTCAAACGGCTGTTGTGTGATCAGATATTTTACCTTGAGAAAAATAGAGTGTTCCAGGCAGGCACCCGCCTGAATAATAAAGCATGGAATAATACAATAGGACTGCGTTCTATTGGTGATAGAACATAGTAATGATTAATAGGGACAGTTGGGGGCATTCGTATTTAATTGTCAGAGGTGAAATTCTTTGATTTATTAAAGACGAACCTATGCGAAAGCATTTGCCAAGGATGTTTTCA

c = cbind(a,b)
head(c)

##                   aa
## 1  >OTU1 MJ12_100225
## 2   >OTU10 MJ12_7529
## 3 >OTU100 MJ12_41755
## 4  >OTU1000 MJ4_9516
## 5 >OTU1001 MJ4_10787
## 6  >OTU1002 MJ4_1279
##                                                                                                                                                                                                                                                                                                                                                                                                          bb
## 1           AGCTCCAACAGCGTATATAAAAGTTGTTGCGATTAAACAGCTCGTAGTCGGACTTCTGCTGAGGATGACCGGTCCGCCCTCTGCGTGAGCATCTGGCTTAGCCTTAGCATTTTCTTGGGCAGCGATATTGCACTTGACTGTGTGGTGTGGTACACAAGGCTTTTAGCTTGTGGAAATCAGAGCGCTTCCAGCAGGCATATGTCTTGAATACGTTACCATTGAATAATAAGATAAGACTTTATGTTTTATTGTTGTTGGTTTTTAGAATTGGGGTAATGCTTGATAGGGATACTTAGGGGCATTTCTGTTTAACTGTCAGAGGTGAACTTCTTTGATTTGTTACAGACAAACTATTGCTAAAGCATTTGCCAATGATGTTTTCC
## 2        AGCTCCAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTCTGACAGGAGCGACCGGTCCGAGCCTTCGGGTTCGTGTACTTTGTGTTGTCTCTAGTCATCCTTATGGAGAGCGTGTCTTTCATTAATTTGTTGGGCATGGGATCCATATCGTTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGGCAATCGCTGTGAATATATTAGCATGGAATAATAGGATAGGACTTTGGTACTATTTTGTTGGTTTGCATACCAAGGTAATGATTAAGAGGGACAGTTGGGGGTATTCGTATTTAATTGTCAGAGGTGAAATTCTTGGATTTATGAAAGACGAACTACTGCGAAAGCATTTACCAAGGATGTTTTCA
## 3                                                  AGCTCCAATAGCGTATATTAAAGTTGTTGCGGTTAAAAAGCTCGTAGTCGAATTTCAATTATTGTCAGTGATTTAGTCTATTAGGCTCAAGTTGCTGGCATAATTTATCTCATCACGAGATAATATTACTTTGAGGAAATTAGAGTGTTTAAAGCAAGCGATTGTTTTGAATAATTTAGCATGGAATAATGATATAGGACCTTATTTTATATTCATTGGTTTTACGAAATGAGGTAATGATTAATGGAAACAGACAGGGGCATTCGTATTTAACTGTCAGAGGTGGAATTCTTGGATTTGTTAAAGACGATCAACTGCGAAAGCGTTTGCCAGGAATGTCTTCA
## 4                                   AGCTCCAATAGTGTATCTTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGAATTTCTGGCAATACAGCTGATCTTTTCGAGAAGCACGTATATGTCATCCGCATTGGAAGCACATTCGATAGCAATATCAGTGTGTGGAAGGTGCATTTTACTTTGAAGAAATTAGAGTGTTCCAGGCAGGCTATGCTTGAATATTACAGCATGGAATAATGGAATAGGACAAGTTCATTTTATTGGTTATTGAATTTGTAATGATTAACAGGGACAGTTGGGGGCATTCGTATTGGGTAGTCAGAGGTGAAATTCTTAGATTTACTCAAGACGCACTAATGCGAAAGCATTTGCCAAGGATGTTTTCA
## 5 AGCTCTAATAGCGTATATTAATGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATTTTAGAGCTTGACAGACGGTCCGCTTGTTTTGCAGGTGAGTACTGTTCCTGTTGGCTTTATTCTTTTCTGGGGGCCGGCGTTCCCTTAATCAGGAGCGTCATTGGCCGTACAGATCGTTTACCTTGAAAAAATTAGAGTGTTCAAAGCAGGCTTCCGCAATTTTTGAATACATTAGCATGGAATAATGGAATAGGACTTTGGTCTTATTTTGTTGGTTTAAAGGCCGAGGTAATGATTAATAGGGATAGTTGGGGGCGTTCATATTCCGTTGTCAGAGGTGAAATTCTTGGATTTACGGAAGATGAACTACTGCGAAAGCATTCGCCAAGGATGTTTTCA
## 6                                AGCTCCAGTAGCGTATATTAAAGTTGTTGCAGTTAAAAAGCTCGTAGTTGAACTTCTGTGTTGCGATGCTGAGTCCTCCGGGACAAACTTTGCGACACATCCATCTGAATAGCATTTCGGCAATTCAAACGGCTGTTGTGTGATCAGATATTTTACCTTGAGAAAAATAGAGTGTTCCAGGCAGGCACCCGCCTGAATAATAAAGCATGGAATAATACAATAGGACTGCGTTCTATTGGTGATAGAACATAGTAATGATTAATAGGGACAGTTGGGGGCATTCGTATTTAATTGTCAGAGGTGAAATTCTTTGATTTATTAAAGACGAACCTATGCGAAAGCATTTGCCAAGGATGTTTTCA

c$aa = gsub(" .*$", "", c$aa)
head(c$aa)

## [1] ">OTU1"    ">OTU10"   ">OTU100"  ">OTU1000" ">OTU1001" ">OTU1002"

selected = data.frame(
  a1 = c("a", "b"), a2 = c("a", "b")
)
for (i in 1:100) {
  selected[i,] = subset(c, grepl(otu_table_5_12_17_29_top100$`OTU ID`[i], c$aa))
}

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU160", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU160", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU151", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU151", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU2", : replacement element 1 has 111 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU2", : replacement element 2 has 111 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU46", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU46", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU159", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU159", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU19", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU19", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU5", : replacement element 1 has 111 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU5", : replacement element 2 has 111 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU75", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU75", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU69", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU69", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU84", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU84", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU117", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU117", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU123", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU123", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU61", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU61", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU40", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU40", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU30", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU30", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU86", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU86", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU80", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU80", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU137", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU137", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU98", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU98", : replacement element 2 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU58", : replacement element 1 has 11 rows to replace 1 rows

## Warning in `[<-.data.frame`(`*tmp*`, i, , value = structure(list(aa =
## c(">OTU58", : replacement element 2 has 11 rows to replace 1 rows

head(selected)

##           a1
## 1    >OTU805
## 2    >OTU646
## 669  >OTU160
## 569  >OTU151
## 1041 >OTU472
## 738    >OTU2
##                                                                                                                                                                                                                                                                                                                                                                                                      a2
## 1    AGCTCCAATAGCGTATATTAAAGTTGCTGCAGTTAAAAAGCTCGTAGTTGGATCTCGGGGTACCTGGTCTAGCCGGTATAGCTCACGCTAAAACTGGTCTGACCCCCCAACGTGAGGTATTCGTCTCCGGGTGCACTTCATTGGGTGCTCGTGAGGGGCCTCGTAGTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGGCTCATATGTTCGGATAGTCGTGCATGGAATAACTGAATACGTGACTGGCGCTTTTCAGTTGGTCTCGGCGACTGGTCAAATGATTAAGAGGGACTGCCGGGGGCATTGGTATGGCGGGGTTAGAGGTGAAATTCTTGGATCCTCGCCAGACCAACAAAAGCGAAAGCATTTGCCAAGAATGTCTTCA
## 2         AGCTCCAATAGCGTATATTAAAGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATCTCGGGTTCGGACTGGCGGTCCTACTCACGTAGGTACTGCTTGTTCCGACCTCCCTGCCGGTTTGTCCCTTGGTGCTCTTCGCTGAGTGCCTTGGGTGGCCGGAACGTTTACTTTGAAAAAATTAGAGTGTTCAAGGCAGGCGTCTCGCCTGAATATTGTTGCATGGAATAATGGAATAGGACCTCGGTTCTATTTTGTTGGTTTTCGGAACTTGAGGTAATGATTAAGAGGGACAGACGGGGGCATTAGTATTACGGTGTTAGAGGTGAAATTCATGGATCGCCGTAAGACTAACTACTGCGAAAGCATTTGCCAAGAATGTTTTCA
## 669                     AGCTCCAAAAGCGTATGCTAAAATTGTTGCGGTTAAAAAGTTCGTAGTTGAAGTTTAACTCATAGGAATCATATCCTCCAGCTGGAGAAGGTGGTTTCATGAGTATGGTGATCTTATATGTTACTCTTGGTAACATAGAGGGAAAGCCCCTTTACTTTGAAGAAATTAGAGTGTTAAAAGCAAGGTGAAGCCTGAATATGCTAGCATGGAATAATGAAATACAGTAATAGGCTTGTTTTATCGGTTCGAATCTATTATGATAATTAACAGGAGTAGGTAGGGACATCCGTATTCGCTTGTGAGAGGTGGAATTCTTGGATCTTGCGGGGACGATCAACTGCGAAAGCATTTGTCAGATATGCTTTCT
## 569           AGCTCCAATAGCGTATATTTTAGTTGTTGCAGTTAAAAAGCTCGTAGTTGTATTTCTGGTGGGAGAATGCGGTCCGACTTATGTGGGTACTTGCACTCTCCTGCCATCTTTGGGTGGAGCCAGTTTGGCATTAAGTTGTCGAGCTGGGGAGACCCATCGTTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGACTTATGTCGTTGAATATATTAGCATGGAATAATTAGATATAACAGAAATATTTCTAGTTGGTTGTGTATTTTTGTAATGATAAATAGGAACAGTTGTGGTTATTCGTATTCACTTGTAAGAGGTGAAATTCTTAGATTAAGGGAAGACGAACTACTGCGAAAGCATTTAACAAGGATGTTTTCA
## 1041      AGCTCCAATAGCATATATTAAAGTTGTTGCGCTTAAAAAGCTCGTAGTTGGACTTCTGCCAAGGACGACAGGTCCGCCCTCTGGGTGTGCATTTGCTTCGTCCCGGGTATCTTCTTGAAGGGTGTGACTGCTCTTTGCTGGGTGGTTCATTCTTCTGGACTTTTACTTTGAGGAAATTAGAGTGCTTCAAGCAGGCATATGCTTTGAATACTTTGGCATGGAATAATGCATAAGGATTTCGGTTCGAGTTTGTTGGCTTCAAGGGCTGAAATAATGATTAATAGGGATGGTTGGGGGCATTCGTATTTAACTGTCAGAGGTGAAATTCTTGGATTTGTTAAAGACGGACTACTGCGAAAGCATTTGCCAAGGATGTTTTCA
## 738       AGCTCCAATAGCGTATATTAAAGTTGTTGCGGTTAAAAAGCTCGTAGTTGGATTTCTGCTGAGGACGACCGGTCCGCCCTCTGGGTGAGTATCTGGCTCGGCCTGGGCATCTTCTTGGAGAACGTAGCTGCACTTGACTGTGTGGTGCGGTATCCAAGACTTTTACTTTGAGGAAATTAGAGTGTTTCAAGCAGGCACACGCCTTGAATACATTAGCATGGAATAATAAGATAGGACCTTGGTTCTATTTTGTTGGTTTCTAGAGCTGAGGTAATGATTAATAGGGATAGTTGGGGGCATTCGTATTTAACTGTCAGAGGTGAAATTCTTGGATTTGTTAAAGACGGACTACTGCGAAAGCATTTGCCAAGGATGTTTTCA

c_copy = c
c_copy$aa = gsub(">", "", c_copy$aa)
head(c_copy$aa)

## [1] "OTU1"    "OTU10"   "OTU100"  "OTU1000" "OTU1001" "OTU1002"

selected_2 = c_copy[match(otu_table_5_12_17_29_top100$`OTU ID`, c_copy$aa), ]
head(selected_2)

##          aa
## 1411 OTU805
## 1234 OTU646
## 669  OTU160
## 569  OTU151
## 1041 OTU472
## 738    OTU2
##                                                                                                                                                                                                                                                                                                                                                                                                      bb
## 1411 AGCTCCAATAGCGTATATTAAAGTTGCTGCAGTTAAAAAGCTCGTAGTTGGATCTCGGGGTACCTGGTCTAGCCGGTATAGCTCACGCTAAAACTGGTCTGACCCCCCAACGTGAGGTATTCGTCTCCGGGTGCACTTCATTGGGTGCTCGTGAGGGGCCTCGTAGTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGGCTCATATGTTCGGATAGTCGTGCATGGAATAACTGAATACGTGACTGGCGCTTTTCAGTTGGTCTCGGCGACTGGTCAAATGATTAAGAGGGACTGCCGGGGGCATTGGTATGGCGGGGTTAGAGGTGAAATTCTTGGATCCTCGCCAGACCAACAAAAGCGAAAGCATTTGCCAAGAATGTCTTCA
## 1234      AGCTCCAATAGCGTATATTAAAGTTGTTGCAGTTAAAAAGCTCGTAGTTGGATCTCGGGTTCGGACTGGCGGTCCTACTCACGTAGGTACTGCTTGTTCCGACCTCCCTGCCGGTTTGTCCCTTGGTGCTCTTCGCTGAGTGCCTTGGGTGGCCGGAACGTTTACTTTGAAAAAATTAGAGTGTTCAAGGCAGGCGTCTCGCCTGAATATTGTTGCATGGAATAATGGAATAGGACCTCGGTTCTATTTTGTTGGTTTTCGGAACTTGAGGTAATGATTAAGAGGGACAGACGGGGGCATTAGTATTACGGTGTTAGAGGTGAAATTCATGGATCGCCGTAAGACTAACTACTGCGAAAGCATTTGCCAAGAATGTTTTCA
## 669                     AGCTCCAAAAGCGTATGCTAAAATTGTTGCGGTTAAAAAGTTCGTAGTTGAAGTTTAACTCATAGGAATCATATCCTCCAGCTGGAGAAGGTGGTTTCATGAGTATGGTGATCTTATATGTTACTCTTGGTAACATAGAGGGAAAGCCCCTTTACTTTGAAGAAATTAGAGTGTTAAAAGCAAGGTGAAGCCTGAATATGCTAGCATGGAATAATGAAATACAGTAATAGGCTTGTTTTATCGGTTCGAATCTATTATGATAATTAACAGGAGTAGGTAGGGACATCCGTATTCGCTTGTGAGAGGTGGAATTCTTGGATCTTGCGGGGACGATCAACTGCGAAAGCATTTGTCAGATATGCTTTCT
## 569           AGCTCCAATAGCGTATATTTTAGTTGTTGCAGTTAAAAAGCTCGTAGTTGTATTTCTGGTGGGAGAATGCGGTCCGACTTATGTGGGTACTTGCACTCTCCTGCCATCTTTGGGTGGAGCCAGTTTGGCATTAAGTTGTCGAGCTGGGGAGACCCATCGTTTACTGTGAAAAAATTAGAGTGTTCAAAGCAGACTTATGTCGTTGAATATATTAGCATGGAATAATTAGATATAACAGAAATATTTCTAGTTGGTTGTGTATTTTTGTAATGATAAATAGGAACAGTTGTGGTTATTCGTATTCACTTGTAAGAGGTGAAATTCTTAGATTAAGGGAAGACGAACTACTGCGAAAGCATTTAACAAGGATGTTTTCA
## 1041      AGCTCCAATAGCATATATTAAAGTTGTTGCGCTTAAAAAGCTCGTAGTTGGACTTCTGCCAAGGACGACAGGTCCGCCCTCTGGGTGTGCATTTGCTTCGTCCCGGGTATCTTCTTGAAGGGTGTGACTGCTCTTTGCTGGGTGGTTCATTCTTCTGGACTTTTACTTTGAGGAAATTAGAGTGCTTCAAGCAGGCATATGCTTTGAATACTTTGGCATGGAATAATGCATAAGGATTTCGGTTCGAGTTTGTTGGCTTCAAGGGCTGAAATAATGATTAATAGGGATGGTTGGGGGCATTCGTATTTAACTGTCAGAGGTGAAATTCTTGGATTTGTTAAAGACGGACTACTGCGAAAGCATTTGCCAAGGATGTTTTCA
## 738       AGCTCCAATAGCGTATATTAAAGTTGTTGCGGTTAAAAAGCTCGTAGTTGGATTTCTGCTGAGGACGACCGGTCCGCCCTCTGGGTGAGTATCTGGCTCGGCCTGGGCATCTTCTTGGAGAACGTAGCTGCACTTGACTGTGTGGTGCGGTATCCAAGACTTTTACTTTGAGGAAATTAGAGTGTTTCAAGCAGGCACACGCCTTGAATACATTAGCATGGAATAATAAGATAGGACCTTGGTTCTATTTTGTTGGTTTCTAGAGCTGAGGTAATGATTAATAGGGATAGTTGGGGGCATTCGTATTTAACTGTCAGAGGTGAAATTCTTGGATTTGTTAAAGACGGACTACTGCGAAAGCATTTGCCAAGGATGTTTTCA

nrow(selected_2)

## [1] 100

selected_2$aa = paste0(">", selected_2$aa)
head(selected_2$aa)

## [1] ">OTU805" ">OTU646" ">OTU160" ">OTU151" ">OTU472" ">OTU2"

selected_2$aa == paste0(">", otu_table_5_12_17_29_top100$`OTU ID`)

##   [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [31] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [46] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [61] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [76] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [91] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE

selected$a1 == selected_2$aa

##   [1] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [16] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [31] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [46] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [61] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [76] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE
##  [91] TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE TRUE

Data organization_OTU_gene

CFY-the Great Outlaw (grady1108@163.com)

2021/5/7

Data input

Data organization

It is better to use the “match” method. Thank you for watching~~~