library(rbiom)
infile <- system.file("extdata", "hmp50.bz2", package = "rbiom")
biom <- read.biom(infile)
ex1 <- subset(biom, Age > 30)
ex2 <- subset(biom, `Body Site` %in% c("Saliva", "Stool"))
ex3 <- subset(biom, Age < 25 & BMI > 22)
Run the code above in your browser using DataLab