Dear list members,
I am having trouble understanding the output of an adonis model I am
running first with then without the argument ?strata? in the formula. I
have seen this question asked before but no answer were provided.
I would like to compare microbial communities associated with healthy and
diseased reef algae. Healthy and diseased algae have been collected in
different islands. I am more interested in the difference between the
Healthy and Diseased groups, I am not really interested in comparing the
islands but still would like to take the islands into account. From what I
read it looks like a way to do it is to include Island as a random factor
in my model.
Variable: Bray-Curtis dissimilarity matrix (dist)
Fixed factor: Disease (Dis)
Random factor: Island
adonis1<- adonis(dist ~ Dis, strata = Island, data = meta)
adonis2<- adonis(dist ~ Dis, data = meta)
However, I get the exact same output for both models:
Df SumsOfSqs MeanSqs F.Model R2 Pr(>F)
Dis_status 1 0.8934 0.89345 1.9845 0.0397 0.001 ***
Residuals 48 21.6103 0.45021 0.9603
Total 49 22.5038 1.0000
---
I tried with a subset and I had the same output besides a different
p-value.
Could someone please explain to me why the outputs are identical?
I am sorry if that sounds trivial but but I don?t really understand why
R2, SumsOfSqs and MeanSqs are exactly identical with or without the strata
argument. I?d like to understand better what?s behind it.
Thanks a lot,
Gaelle
[[alternative HTML version deleted]]