Por sinal, aumentando o numero de niveis de CC para 1250 e usando o
codigo q te mostrei acima:
CC <- paste("H-",1:1250,sep="")
dim(a)
> system.time(a.2 <- ddply(a,c('R','CC','TF'),ma.fun,k=3,.parallel=T))
user system elapsed
65.757 30.706 91.359