[R-sig-eco] Adonis output with or without strata
Jakob Russel
ru@@e|2620 @end|ng |rom gm@||@com
Fri Apr 10 09:52:47 CEST 2020
Dear Gaëlle,
The p-values in PERMANOVA are generated by permutations. These permutations
are by default free, but they can be constrained by the strata factor.
The strata argument only affects the permutations, and hence will only
affect the p-value.
The Sum of Squares (and R2) are only based on the fixed effects.
Regards,
Jakob
Den fre. 10. apr. 2020 kl. 09.42 skrev Gaëlle Quere <GaelleQUERE using hotmail.com
>:
> 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]]
>
> _______________________________________________
> R-sig-ecology mailing list
> R-sig-ecology using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
>
[[alternative HTML version deleted]]
More information about the R-sig-ecology
mailing list