Thanks for the guidance, Gavin!
I've taken this approach, and I have the bootstrap running now.
I'm having the wrapper function pull 1 column of coefficients (for all
predictors/one response variable) at a time. I also made a hobbled version
of adonis that just generates the coefficients, and
Hi Jari,
Thanks so much for your clarification and expansion of Gavin's comments.
I was confused by the application of bootstrapping to a function that is
already
using a permutation procedure, but I see now that bootstrapping the
regression
coefficients would operate independently of adonis'
Rafter,
The permutation test goes nowhere near the coefficients. IIRC it works on
the sum of squares decomposition directly.
Jari has already passed on more info about the bootstrapping idea.
Basically the boot package has functions that will apply your function to
bootstrap samples of the
I should add that I'm not averse to digging around in adonis if necessary. I
modified adonis to call daisy (which accepts mixed data) instead of vegdist
(which doesn't) so that I could fit these models with predictor coefficients
in the first place.
--
View this message in context:
I apologize that this seems to be a half function-specific, half stats
question. My stats knowledge is basic in terms of working with formulae, but
pretty decent conceptually.
I'm conducting a set of multivariate analyses using adonis.
Predictor variables are mixed continuous and categorical.
Rafter,
I am not Gavin, but I step in and speak for him.
Gavin quite clearly suggested that you bootstrap adonis. The boot package is
intended to generate bootstrap samples and to estimate the bootstrap confidence
limits for *any* function. The boot package does not contain those functions,