options ls=78; title "Hotellings T2 - Women's Nutrition Data"; data nutrient; infile "D:\Statistics\STAT 505\data\nutrient.txt"; input id calcium iron protein a c; run; proc iml; start hotel; mu0={1000, 15, 60, 800, 75}; 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{calcium iron protein a c} into x; run hotel;