running adonis on a merged phyloseq object

[pauGuas]

Note: I realized this is actually because of my OTU table, so I am asking this on the phyloseq page now!
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

Hi, I added more data to an existing phyloseq object and need to run adonis on it. I can run adonis on each object separately, but after merging I cannot. There are no duplicate sample names. The problem is that my phyloseq::distance leads to a matrix of only NA.

ps_norm <- transform_sample_counts(new.rm, function(x) x / sum(x) )
this command works

ps.dist <- phyloseq::distance(ps_norm, method="bray")
I get this error: Error in vegdist(c(0.0724891494707987, 0, 0.00399093584063314, 0.00148928857830221, : missing values are not allowed with argument 'na.rm = FALSE'

ps.dist <- phyloseq::distance(ps_norm, method="bray", na.rm=TRUE)
A matrix of the resulting object shows that every distance is "NA".

Does anyone have advice? I have successfully merged phyloseq objects before, so I'm not sure what is wrong with this one.