%karma
A=[ 31 34 29 26 32 35 38 34 31 29 32 31];
B=[26 24 28 29 30 29 31 29 32 26 28 32];
C=[30 30 32 31 29 27 25 30 31 32 34 33];
D=[32 25 31 26 32 27 28 29 29 28 23 25];

figure(1)
subplot(221)
normplot(A)
subplot(222)
normplot(B)
subplot(223)
normplot(C)
subplot(224)
normplot(D)
M=[A' B' C' D']

[P,ANOVATAB,STATS]=kruskalwallis(M);%anova1(M)
multcompare(STATS)