Re: [R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-10 Thread Rafter
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

Re: [R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-09 Thread Rafter
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'

Re: [R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-09 Thread Gavin Simpson
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

Re: [R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-08 Thread Rafter
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:

[R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-08 Thread Rafter
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.

Re: [R-sig-eco] calculating standard error of coefficients from adonis model

2014-05-08 Thread Jari Oksanen
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,