I also have a request for decontam: By definition, isContaminant() only works for taxa that are found in 1+ samples. Can you add a note about that in the isContaminant() documentation? When I first ran the function, my logical vector contained lots of NAs because I hadn't yet excluded taxa found in only 1 sample. This could be confusing for some users.
From Nicole: