Hello,

I am attempting to use the Adonis function (part of vegan package in R) to determine the significance of geographic location (Site) vs host species (Species) in determining Parasite species. While I get each explanatory variable, I would also like to see the significance of the interaction between the variables. However, I can't figure out how to get that result in the permanova table.

adonis2(formula = parasite ~ Species * Site, data = parasite.env)

Df SumOfSqs R2 F Pr(>F)

Species 2 1.1714 0.09079 2.6592 0.018 *

Site 5 3.5817 0.27759 3.2522 0.001 ***

Residual 37 8.1497 0.63162

Total 44 12.9028 1.00000

When the sample data is used with the same format, it gives both explanatory variables AND the interaction between the variables (Management:A1):

adonis2(formula = dune ~ Management * A1, data = dune.env)

Df SumOfSqs R2 F Pr(>F)

Management 3 1.4686 0.34161 3.2629 0.004 **

A1 1 0.4409 0.10256 2.9387 0.015 *

Management:A1 3 0.5892 0.13705 1.3090 0.201

Residual 12 1.8004 0.41878

Total 19 4.2990 1.00000

Maybe there is some discrepancy between my data and the sample data format that is causing this issue, but I cannot figure it out. Any help is appreciated, and I'm happy to provide more details if needed