set.seed(0921)
response <- c(rnorm(10),rpois(10,0.2),rnorm(10,,2))
fact <- gl(3,10,labels=LETTERS[1:3])
# Not enough permutations here but it runs faster
# permutation Bartlett test
perm.bartlett.test(response~fact,nperm=49)
# Pairwise comparisons
pairwise.perm.var.test(response,fact,nperm=49)
Run the code above in your browser using DataLab