-
Hi, Thanks for your help! |
Beta Was this translation helpful? Give feedback.
Replies: 3 comments
-
The idea behind Homogeneity of variances is harder to check within a linear model framework; one often resorts to graphical techniques like scatter plots of residuals or QQ plots. I suppose you could artificially create groups by cutting the continuous covariate in |
Beta Was this translation helpful? Give feedback.
-
There is no direct way of finding that information directly from library(vegan)
data(mite, mite.env)
mad2 <- adonis2(dist(decostand(mite, "hell")) ~ SubsDens + WatrCont, mite.env)
mdb <- dbrda(dist(decostand(mite, "hell")) ~ SubsDens + WatrCont, mite.env)
ordiresids(mdb, residuals="working")
ordiresids(mdb, "scale", residuals="working")
ordiresids(mdb, "qq", residuals="working") The internal mathematical models are identical in |
Beta Was this translation helpful? Give feedback.
-
Here a simple test to see how adequate the library(vegan)
data(mite, mite.env)
mbb <- betadisper(vegdist(sqrt(mite)), mite.env$Shrub, type="centroid") # default type = "median" is worse
mdb <- dbrda(vegdist(sqrt(mite)) ~ Shrub, mite.env)
dbb <- mbb$distances # distances to betadisper centroids
dw <- residuals(mdb, type="working")
dr <- as.matrix(residuals(mdb, type="response")) # needs as.matrix
plot(dbb, sqrt(rowSums(dw^2))) # not quite similar
plot(dbb, sqrt(rowSums(dr^2))) # not quite linear, but exact mapping So this would indicate that response residuals are more or less similar to |
Beta Was this translation helpful? Give feedback.
There is no direct way of finding that information directly from
adonis2
. However, you can build a correspongingdbrda
model and use tools to analyse its residuals:The internal mathematical models are identical in
dbrda
andadonis2
justifying this approach. However, the residuals concern internal ("working") form of dissimilarities (Gower double-centred dissimilarities), and for