求教如何运用SAS程序进行信度分析

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

如对48个对象进行两次评分，结果如下，请教各位如何用SAS编程求同质性信度、再测信度和组内相关系数ICC值（intraclass correlation coefficients,ICC），谢谢！

评分1  评分2

8   8

8   6

11  7

10  9

8   8

12  12

12  12

14  14

7   10

14  15

13  15

13  13

10  12

14  12

14  13

10  8

9   9

9   9

9   9

9   8

10  12

9   6

9   8

12  12

12  14

9   9

9   9

11  10

11  13

13  14

13  13

8   8

8   8

9   8

9   10

15  15

15  15

10  10

10  10

8   6

8   8

9   9

9   6

10  10

10  9

9   9

9   9

10  10

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：sas程序 信度分析 coefficients correlation coefficient 程序 SAS 求教

http://sugiclub.blogspot.com/2008/06/reliability-analysis-calculate-and.html

這篇文章有附上程式碼的寫法

上面提到要直接算出組間相關係數需要寫程式碼

上面貼的網址有

自行貼到SAS上吧

谢谢bandbird的解答，我将bandbird提供的资料贴上，请大家参考，也请大家验证，谢谢！

%macro Icc_sas(ds, response, subject);

/*ds：資料名稱（含 library）response：反應變數名稱 subject：ID 名稱*/

 ods output OverallANOVA =all;
 proc glm data=&ds;
     class &subject;
     model &response=&subject;
 run;
 data Icc(keep=sb sw n R R_low R_up);
     retain sb sw n;
     set all end=last;
     if source='Model' then sb=ms;
     if source='Error' then do;sw=ms; n=df; end;
     if last then do;
     R=round((sb-sw)/(sb+sw), 0.01);
     vR1=((1-R)**2)/2;
     vR2=(((1+R)**2)/n +((1-R)*(1+3*R)+4*(R**2))/(n-1));
     VR=VR1*VR2;
     L=(0.5*log((1+R)/(1-R)))-(1.96*sqrt(VR))/((1+R)*(1-R));
     U=(0.5*log((1+R)/(1-R)))+(1.96*sqrt(VR))/((1+R)*(1-R));
     R_Low=(exp(2*L)-1)/(exp(2*L)+1);
     R_Up=(exp(2*U)-1)/(exp(2*U)+1);
     output;
     end;
 run;
 proc print data=icc noobs split='*';
     var r r_low r_up;
     label r='ICC*' r_low='Lower bound*' r_up='Upper bound*';
     title 'Reliability test: ICC and its confidence limits';
 run;
%mend;