clear all %Load datasets data1 = dlmread('YZSplineValues.csv',',',3,0); %7 skips rows(top down) and 2 skips columns (left to right) data2=data1; A = data2(:,[1]); %y15neg B = data2(:,[2]); %y10neg C = data2(:,[3]); %y5neg D = data2(:,[5]); %y5pos E = data2(:,[6]); %y10pos F = data2(:,[7]); %y15pos Array = [A B C D E F]; %ANOVA %[p,tbl,stats] = anova1(hogg); base code [p,tbl,stats] = anova1(Array); %adapted for my code [c, m, h, nms] = multcompare(stats,'alpha',0.05,'ctype','hsd'); %post-hoc tukey's analysis