options ls=78; title "Paired Hotelling's T-Square"; data spouse; infile "D:\Statistics\STAT 505\data\spouse.txt"; input h1 h2 h3 h4 w1 w2 w3 w4; d1=h1-w1; d2=h2-w2; d3=h3-w3; d4=h4-w4; run; proc print; 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 spouse; read all var{d1 d2 d3 d4} into x; run hotel;