Hallo List,
I'm trying to implemement a restricted permutation scheme in permutest(). More
precisely I have dependence in my data that should be allowed for in the
permutation - I simulated the problem in the example of the vegan documentation
p.24:
library(vegan)
data(varespec)
## Bray-Curtis distances between samples
dis <- vegdist(varespec)
## First 16 sites grazed, remaining 8 sites ungrazed
groups <- factor(c(rep(1,16), rep(2,8)), labels = c("grazed","ungrazed"))
## Calculate multivariate dispersions
mod <- betadisper(dis, groups)
## Perform test
anova(mod)
##simulation of dependence; blocks
blocks<-factor(c(rep(1,4),rep(2,4),rep(3,4),rep(4,4)))
##the unrestricted Permutation test for F is
permutest(mod, pairwise = TRUE)
..I wasn't able to find out the right syntax for implementation of blocks, I
tried with argument strata, which would be blocks, in permutest(mod, pairwise =
TRUE, blocks) but this doesn't work. By the way, what is 'pairwise = TRUE' for?
Thanks for any help,
Kay
---------------------------------------------
Mag. Kay Cichini
Institute of Botany
Sternwartestr. 15
A-6020 Innsbruck
Tel.: +43 (0)512 507 5938
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
