Importance, R^2 and p-values for single env predictors. about s-jsdm

[MaximilianPi]

Hi Julian,

there is no central function for everything:

p-values and effect estimates

summary(model)

reports the p-value and the effect estimates of the individual predictors

R2 / analysis of variance

anova(model)

calculates the analysis of variance for the (two) three groups (env, spatial, and biotic).

res = anova(model)
print(res$result)

reports the variances and the R2 values for the three groups, but not for the individual predictors. But we can calculate the variance/R2 for each site (still experimentally):

res = anova(model, individual=TRUE)
dim(res$result)

Variation partitioning

importance(model)

calculates the variation partitioning of the three groups (env, spatial, and biotic) for each species. We could do this also for each predictor, but it is not yet implemented (we will add this feature in the next sub-version)

So we cannot yet calculate the R2 for each single environmental predictor but we will work on this.

Max

[florianhartig]

Just to add to this - it would of course be possible to calculate ANOVA for each predictor (and we could implement that), but due to the way we compute the ANOVA (see help) this would be computationally costly.

I'm not quite sure if it's worth it - calculating ANOVA for interactions / collinear variables can be a bit tricky. Depending on what you want to use this for, there can be an advantage of ANOVA over standardised effect sizes, but I think for most users, standardised effect sizes (= scale variables and look at summary), which can be interpreted as relative effect sizes, will be more in line with what they want.

[JulFrey]

Thank you very much for the very helpful replies.