p*******s 发帖数: 6 | 1 这个study 比较特殊,n 个case (n<10) , k 个 rater (k>100.)
根据下面这篇文章算icc(3,1). http://www.na-mic.org/Wiki/images/4/4b/Shrout_and_fleiss_ICC.pdf
结果icc(3,1)=0.52,但confidence interval=(-0.004,0.032).
我的计算ci 的code如下(参见文章424页左下角,icc(3,1)的计算参见文章423页右
下角):
ratio= bms/ems;
f= probf(ratio,(n-1),(n-1)*(k-1));
f975f8= quantile('F',.975,(n-1),(n-1)*(k-1));
f975f9= quantile('F',.975,(n-1)*(k-1),(n-1));
fl= f/f975f8;
fu= f*f975f9;
c31l= (fl-1)/(fl+(k-1));
c31u= (fu-1)/(fu+(k-1));
我的code错了吗?还是我的study 比较特殊,不适用文章里的公式?
谢谢! |
|