options ls=78; title "Profile Analysis - Women's Nutrition Data"; data nutrient; infile "D:\Statistics\STAT 505\data\nutrient.txt"; input id calcium iron protein a c; calcium=calcium/1000; iron=iron/15; protein=protein/60; a=a/800; c=c/75; diff1=iron-calcium; diff2=protein-iron; diff3=a-protein; diff4=c-a; run; proc iml; start hotel; mu0={0,0,0,0}; one=j(nrow(x),1,1); ident=i(nrow(x)); ybar=x`*one/nrow(x); s=x`*(ident-one*one`/nrow(x))*x/(nrow(x)-1.0); print mu0 ybar; print s; t2=nrow(x)*(ybar-mu0)`*inv(s)*(ybar-mu0); f=(nrow(x)-ncol(x))*t2/ncol(x)/(nrow(x)-1); df1=ncol(x); df2=nrow(x)-ncol(x); p=1-probf(f,df1,df2); print t2 f df1 df2 p; finish; use nutrient; read all var{diff1 diff2 diff3 diff4} into x; run hotel;