using adonis to compare to distance matrices

Question

I have 2 two dissimilarity matrices. One with observed data comparing among 111 sites and another generated using a null model.

I would like to use the adnois function in vegan to test whether the observed dissimilarities differ significantly from those expected by the null model. However the adonis function will only take one dissimilarity matrix on the left side of the formula.

Does anyone have any idea how to model this test?

Thanks

Answer 1

The answer to this problem was:

meanjac <- function(x) mean(vegdist(x, method='jaccard', diag=TRUE))
test <- oecosimu(x, nestfun=meanjac, method="r1", nsimul = 10^3, statistic='adonis')

which passes a function to get the mean of jaccard dissimilarity matrix to oecosimu, which then uses the 'r1' method to generate null community matrices by randomly shuffling the binary community matrix but assigning the probability of species occupancies based on their observed occupancy and comparing this to the observed dissimilarity matrix.

Thanks Jari for pointing me in the right direction...