近年のSILVAでは、Chloroplast(葉緑体)、Mitochondria(ミトコンドリア)、Eukaryota(真核生物)と同定されて出てくる配列があるので、Rでそれらをはじくコードをメモ。
dataには、行にASV(OTU)、列にサンプル+分類群(Kingdom, Phylum, ...)が入った"カウント+分類群"のデータを用いる。
data2 <- data %>%
dplyr::filter(Order != "Chloroplast" | is.na(Order)) %>%
dplyr::filter(Family != "Mitochondria" | is.na(Family)) %>%
dplyr::filter(Kingdom != "Eukaryota" | is.na(Kingdom))
参考サイト:龍見史恵さんのGitHub (https://github.com/Chikae-Tatsumi/NGS-Data-Processing/blob/master/dada2_16S_IonTorrent.r )
Comentarios