精密度实验变异系数计算的SAS宏程序论文.docVIP

　　精密度实验变异系数计算的SAS宏程序论文 【摘要】 评价医院实验室检测结果的准确性和稳定性，需要进行精密度实验，以确定检测结果是否处于所控制的范围内。通过计算精密度实验批内、批间、天间以及总的变异系数，能够反映实验仪器精密度好坏。通过编写SAS宏程序，可以应用SAS统计软件直接输出以上变异系数的统计报表。 【关键词】 精密度实验； 变异系数； 统计分析报表； SAS宏程序 1 精密度实验 评价医院实验室检测结果的准确性和稳定性，需要进行精密度评价实验，以确定检测结果是否处于所控制的范围内。精密度实验通常包括批内、批间以及日间重复实验。对同一批次质控标本的重复测定，要求每天在不同时间点测定同一批次质控标本2次（2次测定间隔不得少于2小时）.freelacro tformat; /*定义宏，输出统计报表的格式*/ proc format; invalue g 1=20 2=40; %mend tformat; 2.2 计算变异系数 调用proc univariate过程计算批内标准差s、批间差异水平a、天间差异水平b和所有实验结果均数x,利用公式计算出总的变异系数CV。 %macro cv(database,var,dataout); /*定义sas宏程序cv*/ data data1; set database(keep= var.1 var.2 var.3 var.4); d1=( var.1- var.2)**2+( var.3- var.4)**2; proc univariate normal noprint; var d1; output out=d1 sum=sum n=n; data s; set d1; s=sqrt(sum/(4*n)); /*取平方根值，得到批内标准差s*/ run; data data2; set database(keep= var.1 var.2 var.3 var.4); d2=(( var.1+ var.2)/2-( var.3+ var.4)/2)**2; proc univariate normal noprint; var d2; output out=d2 sum=sum n=n; data a; set d2; a=sqrt(sum/(2*n)); /*取平方根值，得到批间差异水平a*/ run; data data3; set database(keep= var.1 var.2 var.3 var.4); d3=( var.1+ var.2+ var.3+ var.4)/4; proc univariate normal noprint; var d3; output out=b std=b mean=x; /*生成天间差异水平b和所有实验结果均数x*/ run; data cv; merge s(keep=s) a(keep=a) b(keep=b x); cv=sqrt((2*(b**2)+a**2+s**2)/2)/x; run; data dataout; set cv(keep=s a b cv) nobs=nobs; /*合并数据集，其中包含s,a，b*/ nu=nobs; s=100*s; /*由于变异系数通常用百分数表示，因此，将所得到的值乘以100*/ a=100*a; b=100*b; cv=100*cv; run; proc datasets; delete data1 data2 data3 d1 d2 s a b cv; /*删除程序中生成的临时数据集*/ quit; %mend cv; 2.3 定义输出结果的位置 定义宏FC，输出批内标准差s、批间差异水平a、天间差异水平b和总的变异系数CV的位置。 %macro fc(invar,cvar,p); if invar then do; if inds=1 then do; roend fc; 定义宏nullset，产生一个输出表，集成已产生的批内标准差s、批间差异水平a、天间差异水平b和总的变异系数CV，并按定义排列。将输出表存入d盘sas目录下，文件名为 tab的文本文件， tab为宏变量。 %macro nullset(data); data _null_; file d:＼sas＼ tab..txt print n=ps notitles header=head; set data; inds+1; col=input(g,g.); %mend nullset; 2.4 运行宏程序