- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
I have pilot data for paired non-parametric continuous variables, so I want to perform sample size calculation to detect differences using a wilcoxon signed rank test. Proc power has options for pairedmeans and twosamplewilcoxon (independent samples, NOT paired), but I can't find an option for paired non-parametric data. How do I go about doing this? Is there a way to manipulate the pairedmeans or twosamplewilcoxon statements to accomplish sample size calculation for paired non-parametric data?
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
The POWER procedure only supports the calculation of sample size for a fraction of statistical methods available in SAS. You may search on the web for formula(e) of accomplishing your goal and compile SAS codes on your own or resort to R or PASS, the latter is a statistical software dedicated to sample size calculation.
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Hello @reubendon,
In section 14.2 "One-Sample Location Problem" of the book
Chow, S.C., Shao, J., and Wang, H. (2008): Sample Size Calculations in Clinical Research, 2nd ed., Chapman & Hall/CRC Press/Taylor & Francis Group, New York
the authors derive an approximation formula using quantiles of the standard normal distribution for the sample size for a Wilcoxon signed rank test, given the significance level a, the desired power b and a series of quantities (probabilities) to be estimated based on data from pilot studies. I think this formula could be implemented in a moderately complicated DATA step reading the pilot study data. They note that their formula is valid only if there are no ties (i.e., paired differences with the same absolute value in the data of the planned study) and it would need to be modified otherwise (in a way not explained in the book).
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
diff=after-before;
and using ONESAMPLEMEANS statement performs T test.
If you need non-parameter method, you could change DIFF into ordinal value and do ttest .
diff ordianl
-1.2 1
0 2
2.3 3
4.1 4
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
I beg to disagree on substituting t-test for Wilcoxon sum-of-rank test unless the differences follow a normal distribution. I guess that assumption is highly unlikely to be tenable for @reubendon's data because if that is the case, @reubendon does not need to resort to Wilcoxon sum-of-rank test and therefore raise the question here.