sample function in shuffle_within_groups behaves differently if length of input is 1

[ctanes]

sample(8) for example would return sample(1:8) and not 8.

Can be fixed with
shuffle_within_groups <- function(x,g) {
ave(x, g, FUN = function(a) if(length(a)>1) sample(a) else a)
}