Within the scope of what NPAR1WAY addresses, no problems, but NPAR1WAY does NOT work on weighted data.  I had to fool it by using FREQ as the weight and recalculating the p-value from _d_, n1 and n2 using a function I wrote myself.  To my surprise, for unweighted data my calculation and the one by NPAR1WAY agree to 8 significant digits!  I suspect I implemented the very same algorithm for the calculation that NPAR1WAY uses, which is why I'm wondering... ... View more

Yes, I did that, and realized NPAR1WAY was calculating things incorrectly because the frequency is not the same as the weight, and all it takes is the frequency (and I have a weight). E.g. the Kolmogorov-Smirnov p-value is always near zero because it uses the sum of the frequencies in the calculation, but the statistical precision goes as the number of entries (i.e. as 10,000, not as 10,000*W).  Same for Wilcoxon Rank Sum tests, using the weight as a frequency NPAR1WAY thinks there are 10,000*W people and thus the p-values are always miniscule.  For the K-S tests I can in principle recalculate p from D, n1 and n2.  But I cannot do anything about the Wilcoxon calculations. they are wrong and I can't recalculate them. Thus the need to create the new unweighted, resampled data set from the original weighted set. ... View more

Yes, using FREQ allows me to get correct distributions.  The problem is that NPAR1WAY doesn't do two-sample Kolmogorov-Smirnov or Wilcoxon Rank Sum tests for weighted data.  I'm trying to find a way to trick it by converting a small sample of weighted data into very large sample of unweighted data, then resampling back to its original size so that NPAR1WAY will give correct results.  I know of no other way of getting SAS to calculate these statistics correctly. So I believe I really need to resample from a weighted, complex survey sample down to a simple unweighted survey sample. ... View more

I tried that previously, it just printed exactly 0 for all the p-values (see column of zeros below).  Perhaps my weights are so large that it's hopeless to get a non-zero p-value that I can correct for the weights.  Seems like I'll just have to calculate the p-value myself. I sure would like to find a function in SAS that does the calculation... 1 Male (<20) ApoB .002233328 6.2103 0.22041 0 .000002242 17.335 0.22293 6.2813 0 2 Male (20-29) ApoB .002473620 7.3677 0.12183 0 .000003057 27.119 0.12511 7.5660 0 3 Male (30-39) ApoB .005193976 15.9733 0.16763 0 .000011119 105.162 0.16785 15.9943 0 4 Male (40-49) ApoB .003007115 9.7514 0.07357 0 .000002399 25.230 0.09987 13.2380 0 5 Male (50-59) ApoB .004811441 14.8999 0.09631 0 .000006126 58.747 0.12704 19.6538 0 6 Male (60-69) ApoB .007293936 19.9366 0.13297 0 .000009974 74.515 0.19752 29.6162 0 7 Male (70-79) ApoB .005602209 11.4631 0.09621 0 .000004495 18.821 0.15906 18.9521 0 8 Male (>=80) ApoB .005729162 7.4794 0.09921 0 .000011073 18.872 0.13626 10.2727 0 9 Female (<20) ApoB .002512096 6.9934 0.22836 0 .000002821 21.866 0.22836 6.9934 0 10 Female (20-29) ApoB .002397185 7.4033 0.09843 0 .000002358 22.488 0.09896 7.4432 0 11 Female (30-39) ApoB .005747304 17.9301 0.16782 0 .000014166 137.876 0.18721 20.0008 0 12 Female (40-49) ApoB .004009143 12.3046 0.08878 0 .000007887 74.293 0.10704 14.8363 0 13 Female (50-59) ApoB .005250789 17.1014 0.10402 0 .000006173 65.481 0.15556 25.5739 0 14 Female (60-69) ApoB .005870556 16.3198 0.10203 0 .000011481 88.728 0.15890 25.4182 0 15 Female (70-79) ApoB .004183428 9.3058 0.07472 0 .000003132 15.496 0.14135 17.6047 0 16 Female (>=80) ApoB .008706090 14.2755 0.16290 0 .000015396 41.394 0.19818 17.3680 0 ... View more

Yes, I saw the equations, I was hoping a function exists to do the work.  I understand in practice no one sums from zero to infinity, only a few terms are kept and special corrections are applied, so I'd much rather find a working function written by experts than have to research this on my own. ... View more

I looked in the output data, all the p-values are 0.  The problem is my weights are huge, as high as 200,000.  But the true p-value calculated properly should not be 0. Might EXACT help?  I fiddled with EXACT, MC, MAXTIME, etc. but it just sat there for long periods of time and I had to kill it.  Have to try N= as well. Anyway, if I could just access the function NPAR1WAY uses I'd be fine... ... View more

Thanks Sanjay, the new features in SAS9.4M2 may be exactly what I want.  Can you specify multiple variables in the GROUP= option? ... View more

The weights are essentially predicted frequencies, and I use them as such.  They are based on a full survey design and using the weights/frequencies when plotting a variable should give a distribution that is close to what one would get if one sampled the entire US civilian population where every sample had weight equal to one. The K-S test should work fine, I just want to find a way to get SAS to calculate the p-value correctly.  It gets the d-statistic correct, I believe. ... View more