cluster_anova_sim.m

NREP = 30000;
N = 1000;
K = 5; 

p_random = nan(NREP, 1);
p_cluster = nan(NREP, 1);

parfor rep=1:NREP
    data = randn(N, 1);
    random_group = mod(randperm(N), K);
    cluster_group = kmeans(data, K);
    
    p_random(rep) = anova1(data, random_group, 'off');
    p_cluster(rep) = anova1(data, cluster_group, 'off');
end


sp = 0.025:0.05:1;

subplot(1,2,1);
hist(p_random, sp);
axis('tight');

subplot(1,2,2);
hist(p_cluster, sp);
axis('tight');