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
- /
- Analytics
- /
- Stat Procs
- /
- Help to understand reliability analysis

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

10-21-2015 09:28 AM

SAS Community, I found this code in an archive at work and it contains no comments. The idea is to calculate a reliable sample size. It is not working properly and I am seeking to understand each component.

The input table cy11rates contains the unique provider ID, risk adjusted rate of complication (rsrr_hwr) and total patient volumn. Is anyone familiar with this process that could help me understand?

data temp;

set cy11rates;

numerator=round(rsrr_hwr*volume);

do i=1 to volume;

d=1;

if i<= numerator then n=1;

else n=0;

output;

end;

run;

proc sql;

select sum(volume) into:cnt

from cy11rates;

quit;

proc freq data=temp;

tables provid /out=b;

run;

data c;

set b;

rtsq=count*count;

rtdv=rtsq/&cnt.;

run;

proc summary sum data=c;

var rtdv;

output out=new sum=nv;

run;

data new2;

set new;

t=(&cnt.-nv)/(_freq_-1);

run;

proc print data=new2; run;

proc anova data=temp;

ods output OverallANOVA =anova_result;

class provid;

model n=provid;

run;

quit;

proc sql;

create table reliability as

select distinct (Fvalue-1)/(Fvalue-1+t) as ICC,

t*(calculated ICC)/(1+(t-1)*(calculated ICC)) as reliability,

(0.7-0.7*(calculated ICC))/((calculated ICC)-(0.7*(calculated ICC))) as sample_size

from new2

join anova_result

on 1=1

where anova_result.source='Model';

quit;

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

10-21-2015 10:38 AM

Please describe how "it is not working properly". Error messages? Post the errors and log.

No output? Unexpected output? How is it unexpected.

Also dealing with sample size generation I would recommend looking into proc power.

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

10-21-2015 10:41 AM

Thank you for your reply. I received a sample size that was a negative number. My hope was to better understand the workings of the procedure so I could determine if adjustments are needed or if the code is bad. The hope was to get a general idea of each step. I will look into proc power though and appreciate that suggestion.

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

10-21-2015 10:48 AM

My understanding of proc power is to plan sample size for a comparison of means. In my analysis I do not have an intervention. I wish to know the minimum patient volume for a provider group so that a risk adjusted complication rate can be considered stable.