data table341 ; input treatment phosphorus @@ ; cards ; 1 13.0 1 24.1 1 11.7 1 16.3 1 15.5 1 24.5 2 42.0 2 18.0 2 14.0 2 36.0 2 11.6 2 19.0 3 15.6 3 23.8 3 24.4 3 24.0 3 21.0 3 21.1 4 35.3 4 22.5 4 16.9 4 25.0 4 23.1 4 26.0 ; proc print ; proc npar1way data = table341 wilcoxon ; class treatment ; var phosphorus ; exact wilcoxon / n = 2000 ; run ; proc boxplot data = table341 ; plot phosphorus*treatment ; run ; * Examine the printed frequency table from PROC FREQ to insure that the rows are ordered correctly ; proc freq data = table341 ; * must be ordered according to expected difference ; tables treatment*phosphorus / norow nocol nopercent jt ; * classification*response order in the tables statment ; exact jt / n = 20000 ; run ;