npar1way histogram

[chopper]

I have a task to create  permutation for rank sum and fit/overlay a normal curve to the resulting histogram. My npar1way code is below. I'm not sure how to create the histogram from the output of the npar1way. I'm using the latest university edition. Any help would be appreciated.

proc npar1way data=mydata wilcoxon;
class type;
exact wilcoxon / n=5000;
run;

[PGStats]

I don't think you can get the 5000 Monte-Carlo permutation Wilcoxon statistics from npar1way.

[Rick_SAS · Accepted solution]

I don't know how comfortable you are with the SAS/IML language, but you this article about

"Monte Carlo for Contingency Tables" contains about 80% of the computations that you will need to run your own Monte Carlo simulation and graph the 5,000 statistics. You just need to implement the Wilcoxon statistic, which you can do by using the RANK function in SAS/IML.