Is there a more efficient way to do many ANOVAs

I want to know if there is a vectorised implementation for doing many ANOVAs.

At the moment I have time series data, with 600 time points. I want to use the same grouping matrix for all 600 of these Y values. Is there a more efficient way of doing this than using a for loop?

e.g. NTimePoints = 600; NExamples = 250; Y = rand(NExamples ,NTimePoints ); Groups = round(rand(NExamples,1)*5)+1;

for t=1:NTimePoints

anova1(Y(:,t),Groups,'off'); end

Thanks!