data jain1 ; input trt plot hor HYDWT1 CARWT1 NO3WT1 NH4WT1 NITWT1 HYDWT2 CARWT2 NO3WT2 NH4WT2 NITWT2 ; dhydwt = HYDWT1 - HYDWT2 ; dcarwt = CARWT1 - CARWT2 ; dno3wt = NO3WT1 - NO3WT2 ; dnh4wt = NH4WT1 - NH4WT2 ; dnitwt = NITWT1 - NITWT2 ; dhydwt2 = sqrt(HYDWT1) - sqrt(HYDWT2) ; dcarwt2 = sqrt(CARWT1) - sqrt(CARWT2) ; dno3wt2 = sqrt(NO3WT1) - sqrt(NO3WT2) ; dnh4wt2 = sqrt(NH4WT1) - sqrt(NH4WT2) ; dnitwt2 = sqrt(NITWT1) - sqrt(NITWT2) ; shydwt2 = sqrt(HYDWT2) ; scarwt2 = sqrt(CARWT2) ; sno3wt2 = sqrt(NO3WT2) ; snh4wt2 = sqrt(NH4WT2) ; snitwt2 = sqrt(NITWT2) ; shydwt1 = sqrt(HYDWT1) ; scarwt1 = sqrt(CARWT1) ; sno3wt1 = sqrt(NO3WT1) ; snh4wt1 = sqrt(NH4WT1) ; snitwt1 = sqrt(NITWT1) ; lcarwt1 = log(CARWT1+.005) ; lno3wt1 = log(NO3WT1+.005) ; cards ; 1 901 1 3.01 23.82 0.00 0.70 1.00 0.49 5.09 0.02 0.05 0.18 1 901 2 1.99 15.49 0.00 0.00 0.44 0.97 6.95 0.00 0.00 0.25 1 901 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 901 4 10.90 31.71 0.30 1.29 1.49 25.81 226.78 0.47 0.94 7.27 1 901 5 8.43 13.49 0.34 0.67 0.01 6.91 11.93 0.38 0.44 1.13 1 959 1 3.92 31.71 0.29 1.51 0.93 0.00 0.00 0.00 0.00 0.00 1 959 2 3.19 22.99 0.23 2.41 0.82 1.07 5.34 0.06 0.53 0.27 1 959 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 959 4 7.23 23.65 0.26 6.11 1.05 7.53 12.05 0.45 1.81 0.90 1 959 5 11.73 17.60 0.35 3.64 0.01 8.97 10.76 0.45 0.45 1.08 1 972 1 2.10 16.45 0.05 1.83 0.49 1.63 13.61 0.27 1.44 0.41 1 972 2 2.77 17.94 0.07 0.86 0.51 4.31 31.85 1.18 2.42 1.01 1 972 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 972 4 6.09 11.18 0.25 5.59 0.31 8.74 32.30 0.15 0.49 1.56 1 972 5 8.21 8.21 0.21 2.19 0.00 5.94 18.90 0.32 2.65 1.13 1 976 1 2.19 17.07 0.03 0.91 0.42 2.95 1.91 0.01 0.06 0.08 1 976 2 0.71 4.63 0.00 0.00 0.46 0.61 3.85 0.00 0.00 0.17 1 976 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 976 4 6.82 16.74 5.21 3.41 0.87 9.16 19.16 1.50 2.50 1.25 1 976 5 8.44 11.25 0.28 1.83 0.01 8.20 8.20 0.37 0.45 0.67 1 984 1 1.82 14.39 0.02 0.29 0.44 0.50 3.80 0.02 0.08 0.14 1 984 2 2.94 22.90 0.27 2.10 0.65 2.28 15.68 0.09 1.58 0.54 1 984 3 5.73 11.73 0.00 0.00 0.29 0.00 0.00 0.00 0.00 0.00 1 984 4 7.21 20.96 1.11 6.88 0.98 10.05 26.39 0.71 0.91 1.32 1 984 5 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 800 1 2.69 21.50 0.00 0.00 0.50 1.91 18.70 0.00 0.00 0.49 1 800 2 2.76 14.85 0.10 1.12 0.71 0.83 6.55 0.01 0.22 0.21 1 800 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 800 4 9.19 12.53 0.25 2.25 1.09 8.54 26.84 1.59 1.53 1.28 1 800 5 8.94 8.94 0.16 0.73 0.01 8.37 9.89 0.46 0.46 0.91 1 863 1 6.54 51.48 0.42 2.52 1.35 1.51 12.67 0.02 0.20 0.40 1 863 2 1.91 13.52 0.02 0.73 0.56 1.86 12.71 0.45 4.45 0.39 1 863 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 863 4 9.17 21.15 0.00 2.75 1.41 8.92 21.55 1.11 1.04 1.26 1 863 5 8.78 15.53 0.20 3.31 0.01 10.02 14.20 0.50 0.42 1.09 1 880 1 1.09 8.03 0.01 0.51 0.28 1.47 12.58 0.03 0.36 0.47 1 880 2 1.91 13.47 0.03 1.57 0.49 1.92 11.12 0.03 0.22 0.51 1 880 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 880 4 10.56 20.30 0.32 4.95 1.54 7.72 20.58 0.39 0.64 1.22 1 880 5 9.90 9.90 0.33 2.48 0.01 8.14 12.58 0.22 1.11 1.11 1 882 1 1.00 7.52 0.09 0.58 0.28 1.04 10.66 0.11 0.23 0.35 1 882 2 2.46 15.46 0.33 1.62 0.69 1.32 8.79 0.11 2.03 0.40 1 882 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 882 4 5.58 8.11 0.30 1.52 0.71 8.80 20.00 3.20 1.28 1.52 1 882 5 7.78 7.78 0.42 2.12 0.01 7.24 11.84 0.33 0.26 0.99 1 897 1 2.58 20.96 0.14 0.34 0.56 0.58 5.84 0.01 0.12 0.21 1 897 2 3.88 29.14 0.37 2.69 0.83 1.11 8.65 0.11 2.86 0.31 1 897 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 1 897 4 8.88 25.16 3.48 6.29 1.63 10.10 18.52 0.51 2.69 1.35 1 897 5 8.15 10.66 0.13 1.38 0.01 9.85 11.49 0.25 1.23 0.99 2 384 1 1.23 9.46 0.01 0.95 0.22 0.44 6.06 0.06 0.54 0.25 2 384 2 1.07 7.02 0.05 0.68 0.26 0.00 0.00 0.00 0.00 0.00 2 384 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 384 4 5.90 24.06 0.18 1.82 1.00 21.55 88.08 0.75 6.84 6.65 2 384 5 12.85 44.06 0.37 1.93 0.02 9.48 22.60 0.51 1.60 1.60 2 385 1 0.73 5.37 0.05 0.30 0.19 0.13 2.15 0.06 0.10 0.09 2 385 2 0.54 2.74 0.02 0.15 0.10 0.63 5.33 0.21 4.85 0.23 2 385 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 385 4 13.58 46.56 0.39 6.89 2.23 9.27 42.78 0.93 5.28 1.93 2 385 5 11.67 26.04 0.36 3.23 0.01 8.40 21.84 0.76 2.18 1.26 2 389 1 0.96 6.89 0.01 0.40 0.17 0.94 10.46 0.18 0.89 0.46 2 389 2 1.16 7.90 0.01 1.25 0.25 0.00 0.00 0.00 0.00 0.00 2 389 3 0.91 7.12 0.01 0.28 0.06 0.00 0.00 0.00 0.00 0.00 2 389 4 11.13 34.87 0.30 4.08 1.11 11.31 41.76 0.61 6.18 2.00 2 389 5 11.69 32.14 0.49 2.63 0.01 9.80 27.78 0.57 1.23 1.47 2 390 1 1.14 8.81 0.01 0.20 0.46 0.21 2.64 0.02 0.22 0.10 2 390 2 0.83 5.38 0.02 0.38 0.77 0.00 0.00 0.00 0.00 0.00 2 390 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 390 4 9.61 26.60 0.30 9.61 2.96 8.82 35.29 0.40 12.51 2.09 2 390 5 13.96 32.56 0.47 1.98 0.03 7.52 17.90 0.36 1.79 1.16 2 392 1 0.52 3.77 0.00 0.38 0.12 0.69 6.36 0.06 1.59 0.23 2 392 2 1.48 10.03 0.02 0.78 0.36 0.00 0.00 0.00 0.00 0.00 2 392 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 392 4 9.07 36.29 0.30 2.19 1.51 8.45 33.79 0.46 3.76 1.69 2 392 5 9.41 25.65 0.34 2.31 0.01 7.47 25.80 0.34 1.15 1.15 2 395 1 1.56 11.89 0.05 0.87 0.28 0.37 4.10 0.09 0.87 0.15 2 395 2 2.07 16.12 0.02 0.59 0.39 0.00 0.00 0.00 0.00 0.00 2 395 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 395 4 8.08 23.56 0.27 2.22 1.01 10.05 34.79 11.36 10.36 2.09 2 395 5 10.32 25.80 0.69 2.15 0.01 9.06 22.65 3.70 1.89 1.21 2 365 1 0.56 3.81 0.00 0.28 0.11 0.40 3.26 0.02 2.06 0.07 2 365 2 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 365 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 365 4 8.61 33.10 0.26 1.19 1.19 8.90 42.07 0.89 7.77 1.94 2 365 5 10.00 25.70 0.29 1.36 0.01 9.67 27.40 0.64 3.22 1.53 2 366 1 0.57 3.78 0.01 0.30 0.10 0.48 5.13 0.20 0.62 0.20 2 366 2 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 366 3 0.97 8.12 0.01 0.32 0.07 0.93 8.89 0.34 1.56 0.07 2 366 4 8.53 35.42 0.26 2.36 0.98 5.62 23.29 2.16 3.52 1.42 2 366 5 12.11 29.26 0.40 3.13 0.01 7.91 18.98 1.98 2.69 0.79 2 369 1 0.44 3.21 0.00 0.14 0.08 0.19 1.68 0.01 0.52 0.06 2 369 2 0.94 5.59 0.01 0.46 0.35 0.00 0.00 0.00 0.00 0.00 2 369 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 369 4 9.85 27.35 0.44 2.84 2.52 6.25 21.06 0.53 3.49 1.18 2 369 5 10.76 19.37 0.43 1.18 0.02 7.40 13.22 0.53 2.03 0.79 2 372 1 0.77 5.55 0.03 0.79 0.19 0.31 6.82 0.09 0.25 0.15 2 372 2 0.62 4.14 0.03 1.64 0.20 0.67 5.30 0.06 4.31 0.23 2 372 3 0.57 5.22 0.00 0.09 0.03 0.00 0.00 0.00 0.00 0.00 2 372 4 9.05 30.45 0.33 6.09 1.81 9.98 42.63 0.82 8.89 2.27 2 372 5 10.43 28.68 0.35 3.30 0.02 8.68 20.24 0.29 1.74 1.08 2 373 1 0.85 6.40 0.01 0.15 0.13 0.12 1.55 0.01 0.02 0.07 2 373 2 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 2 373 3 3.30 31.18 0.02 2.62 0.37 0.00 0.00 0.00 0.00 0.00 2 373 4 9.92 21.82 0.40 2.08 1.19 6.42 23.89 0.60 46.28 1.49 2 373 5 9.93 18.87 0.40 1.49 0.01 10.95 46.92 1.56 3.75 2.03 3 151 1 1.03 8.16 0.01 0.16 0.17 1.10 17.59 0.14 0.21 0.52 3 151 2 0.91 4.93 0.06 0.24 0.19 0.00 0.00 0.00 0.00 0.00 3 151 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 151 4 12.00 33.75 0.98 2.63 1.35 6.72 23.52 1.46 6.50 1.46 3 151 5 20.40 33.60 3.60 2.16 0.01 8.14 18.50 2.52 1.92 1.18 3 153 1 0.35 2.72 0.01 0.06 0.05 1.73 14.32 0.89 5.47 0.47 3 153 2 2.10 16.29 3.00 0.81 0.54 0.00 0.00 0.00 0.00 0.00 3 153 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 153 4 12.80 42.40 12.16 1.36 2.16 11.55 53.90 7.55 7.08 2.85 3 153 5 12.75 16.50 10.65 1.88 0.01 15.20 29.45 25.27 1.33 2.00 3 156 1 1.42 10.89 0.05 0.15 0.22 2.26 19.56 0.06 1.04 0.48 3 156 2 2.90 22.66 3.66 2.65 0.76 0.00 0.00 0.00 0.00 0.00 3 156 3 0.17 1.36 0.03 0.38 0.02 1.93 16.02 0.30 1.59 0.37 3 156 4 11.88 37.62 2.18 1.19 1.39 14.24 77.43 1.60 5.70 3.20 3 156 5 18.81 42.57 1.98 1.49 0.02 16.20 50.40 0.63 2.61 1.71 3 160 1 2.07 16.23 0.08 0.28 0.36 1.60 16.40 0.14 2.70 0.53 3 160 2 1.28 9.09 0.78 0.69 0.36 0.00 0.00 0.00 0.00 0.00 3 160 3 0.21 1.55 0.01 0.04 0.03 3.44 31.79 0.10 4.04 0.23 3 160 4 12.74 55.86 8.04 5.10 2.45 7.98 25.25 14.54 4.95 1.72 3 160 5 12.48 33.28 0.94 1.56 0.02 14.96 46.24 15.91 3.94 2.45 3 161 1 1.71 13.76 0.24 0.59 0.35 0.00 0.00 0.00 0.00 0.00 3 161 2 0.28 1.97 0.25 0.07 0.07 0.98 6.74 0.12 1.51 0.28 3 161 3 2.97 27.36 1.15 1.37 0.28 0.00 0.00 0.00 0.00 0.00 3 161 4 13.05 45.24 9.31 2.87 1.48 12.72 33.92 0.85 7.53 2.01 3 161 5 18.45 27.06 5.41 1.11 0.01 15.99 18.45 3.32 1.35 1.35 3 171 1 0.66 5.10 0.01 0.07 0.15 1.98 17.69 0.48 6.26 0.30 3 171 2 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 171 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 171 4 7.84 20.16 0.28 1.12 0.90 4.55 20.30 0.46 8.96 1.33 3 171 5 12.60 22.50 2.16 1.17 0.01 12.61 22.31 3.49 3.10 1.94 3 172 1 1.87 13.94 0.08 0.22 0.44 1.04 8.52 0.10 2.79 0.26 3 172 2 2.53 18.50 0.94 1.14 0.68 0.49 3.27 0.19 3.11 0.19 3 172 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 172 4 12.57 43.58 13.99 2.51 2.18 11.16 47.82 1.75 9.40 2.87 3 172 5 13.40 15.47 14.85 1.34 0.01 11.63 19.38 2.13 1.94 1.45 3 173 1 3.00 22.78 0.05 0.22 0.37 0.47 6.14 0.04 0.57 0.17 3 173 2 2.76 21.27 0.30 1.01 0.63 1.82 14.58 0.59 3.59 0.36 3 173 3 3.64 35.03 0.05 5.17 0.30 6.49 53.50 0.88 2.16 0.86 3 173 4 17.32 71.33 11.82 2.96 3.36 14.99 75.90 4.03 15.93 3.84 3 173 5 19.62 76.16 5.54 12.81 0.03 7.14 20.83 2.80 2.02 1.13 3 174 1 2.68 20.34 0.64 1.40 0.58 2.02 21.55 0.17 3.42 0.67 3 174 3 4.91 41.79 3.27 3.39 0.56 8.52 77.49 0.90 14.20 0.68 3 174 4 2.26 7.33 0.53 0.86 0.34 10.47 25.76 1.69 3.14 1.45 3 174 5 17.58 31.87 9.89 1.76 0.02 19.78 32.50 1.13 2.26 2.54 3 174 8 0.62 5.41 0.03 0.85 0.04 0.00 0.00 0.00 0.00 0.00 3 182 1 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 182 2 0.00 0.00 0.00 0.00 0.00 1.12 7.17 0.19 1.59 0.32 3 182 3 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 3 182 4 7.58 29.34 5.71 2.36 1.39 6.29 19.26 11.81 6.16 1.54 3 182 5 12.67 40.32 6.22 3.92 0.02 7.72 12.55 9.94 1.93 1.16 3 182 8 6.74 57.96 0.09 1.37 0.30 2.40 21.12 0.06 4.22 0.22 ; run ; /*----------------------------------------------------------------------- The following macro is adapted from SAS code presented in "Applied Multivariate Statistics with SAS Software" by R. Khattree and D. N. Naik 1995, SAS Institute, Cary NC The argument 'data' is the name of the SAS data set containing the variables, 'var' is the list of variables on the 'data' data set, and 'p' is the number of variables in the variable list. -------------------------------------------------------------------------*/ %macro mvnplot(data, var, p) ; %let alp=%sysevalf(&p/2.); proc princomp data=&data cov std out=testb noprint ; var &var ; data qq; set testb; dsq=uss(of prin1-prin&p); run ; proc capability noprint graphics ; probplot dsq/gamma(alpha=&alp theta=0 sigma=2) square ; title2 'Multivariate normal probability plot' ; run; %mend mvnplot ; %mvnplot(jain1,sHYDWT1 sNH4WT1 sNITWT1 scarwt1 sno3wt1,5) ; /* data horiz1 ; set jain1 ; if hor = 1 ; *if plot ne ; run ; %mvnplot(horiz1, sHYDWT1 sNH4WT1 sNITWT1 scarwt1 sno3wt1,5) ; */ /* data test ; do i = 1 to 275 ; x1 = 2*rannor(0) +17 ; x2 = 3*rannor(0) +11 ; x3 = 5*rannor(0) + 7 ; x4 = 2.5*rannor(0) + 3 ; x5 = 5*rannor(0) + 9 ; x6 = 4*rannor(0) + 7 ; x7 = 3*rannor(0) + 2 ; output ; end ; %mvnplot(test,x1 x2 x3 x4 x5 x6 x7,7) ; */