[Above Average]
x;
#x>mu;
#double mu=mean(x);
[Below Average]
x;
#x<mu;
#double mu=mean(x);
[Within 3-Sigma Band]
x;
#x>=lcl and x<=ucl;
#double mu, std, lcl, ucl;
mu=mean(x);
std=stddev(x);
lcl=mu-3*std;
ucl=mu+3*std;
[Within 95% CI]
x;
#x>=lcl and x<=ucl;
#double mu, std, tc, lcl, ucl;
int nx;
stats ix:=x mean:=mu sd:=std n:=nx;
tc=tinv(0.975,nx-1);
lcl=mu-std*tc/sqrt(nx);
ucl=mu+std*tc/sqrt(nx);
