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