joshua@0: %N trials, M answers -> N+M-1 choose M combinations
joshua@0: t_actual=6.121;%change this to find best match for a reported t value
joshua@0: threshold=0.1;
joshua@0: answers=[-4 -3 -2 -1 1 2 3 4];
joshua@0: for count1=1:8
joshua@0: for count2=1:count1
joshua@0: for count3=1:count2
joshua@0: for count4=1:count3
joshua@0: for count5=1:count4
joshua@0: for count6=1:count5
joshua@0: if (count6)~=(count1)
joshua@0:     Answers=[answers(count1) answers(count2) answers(count3) answers(count4) answers(count5) answers(count6)];
joshua@0:     t=sqrt(6)*mean(Answers)/std(Answers);
joshua@0:     test=abs((t-t_actual)/t_actual);
joshua@0:     if ((std(r)>0)&&(test <= threshold)) %best match so far
joshua@0:     threshold=test;
joshua@0:     disp(Answers)
joshua@0:     disp(t)
joshua@0:     disp(t_actual)
joshua@0:     end
joshua@0: end
joshua@0: end
joshua@0: end
joshua@0: end
joshua@0: end
joshua@0: end
joshua@0: end