function [CIs] = makeCIs(data)
     data_size = length(data);
     df = data_size - 1;
     t_critical = tinv(0.025,df);
     std_data = std(data);
     sqrtn = sqrt(data_size);
     CIrange = abs(t_critical*std_data/sqrtn);
     CIs(1) = mean(data) - CIrange;
     CIs(2) = mean(data) + CIrange;
end