使用parallel包并行处理phyloseq数据：过滤低丰度OTU

library(parallel) cl <- makeCluster(detectCores()) clusterEvalQ(cl, library('phyloseq')) clusterExport(cl, 'physeq1')

clusterExport(cl, 'sum_thresh') sum_thresh <- function(x, thresh=0.0005) { ifelse(x <= sum(x)*thresh, 0, x) }

physeq1 <- parApply(colnames(physeq1@otu_table), function(i) { physeq1@otu_table[,i] <- sum_thresh(physeq1@otu_table[,i]) return(physeq1) }, cl=cl)

stopCluster(cl)