turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- General Programming
- /
- How does sas proc power calculate the sample size?

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-29-2014 02:41 PM

In SAS proc power, there is an example for Determining Required Sample Size for a Two-Sample *t* Test.

Here is the original code from SAS

proc power;

twosamplemeans

groupmeans = (13 14) (13 14.5) (13 15)

stddev = 1.2 1.7

groupweights = 1 | 1 2 3

power = 0.9

ntotal = .;

run;

I calculate the sample size by using formula, the result is always smaller the one obtained from SAS output.

For example, the first sample size I got is 60.466176 but SAS got 62.507429. Basically, the fractional n total from SAS is always bigger then the one I calculated based on the common formula which is

n=4*stddev^2*(1.96+1.28)^2/(14-13). I want to know how SAS compute the sample size and base on what formula. I hope someone can help me to find it out. Thanks.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

Posted in reply to saszyc

01-29-2014 03:06 PM

The documentation has a good description:

SAS/STAT(R) 9.2 User's Guide, Second Edition

Proc Power>Details>Computational Methods and Formulas>Two Sample Means

The results from proc power are the total sample size for both groups that match the results from a web calculator as well:

Your formula: 60.444

SAS : 64

UBC Calc: 62

All seem within the same ballpark, though not exactly the same.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

Posted in reply to saszyc

01-30-2014 08:54 AM

One thing to note is that you used the z distribution values (1.96, 1.28) in your calculation. Use of the t distribution values will result in the more conservative number obtained by SAS. Note that this is iterative, starting with the asymptotic normals, getting an estimate of sample size, calculating degrees of freedom and then obtaining t distribution values. (At least I think that is the explanation).

Steve Denham