Open rosebaesj opened 2 years ago
line 135 dim(all_species)
line 144 dim(all_species)
line 151 dim(all_species[apply(all_species, 1, function(bug) sum(bug >= 0.01) >= 0.1*ncol(all_species)), ])
possible answer: the metaphlan_taxonomic_profiles.tsv may have been updated
line 135 dim(all_species)
[1] 468 929
sunjeong: [1] 543 929
line 144 dim(all_species)
[1] 468 925
sunjeong: [1] 543 925
line 151 dim(all_species[apply(all_species, 1, function(bug) sum(bug >= 0.01) >= 0.1*ncol(all_species)), ])
[1] 139 925 = the 139 species after QC included in initial starr manuscripts
sunjeong: [1] 156 925