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

精密度实验变异系数计算的SAS宏程序 ???????? 【摘要】? 　　评价医院实验室检测结果的准确性和稳定性，需要进行精密度实验，以确定检测结果是否处于所控制的范围内。通过计算精密度实验批内、批间、天间以及总的变异系数，能够反映实验仪器精密度好坏。通过编写SAS宏程序，可以应用SAS统计软件直接输出以上变异系数的统计报表。 【关键词】? 精密度实验； 变异系数； 统计分析报表； SAS宏程序 　　1? 精密度实验??? 　　评价医院实验室检测结果的准确性和稳定性，需要进行精密度评价实验，以确定检测结果是否处于所控制的范围内。精密度实验通常包括批内、批间以及日间重复实验。对同一批次质控标本的重复测定，要求每天在不同时间点测定同一批次质控标本2次（2次测定间隔不得少于2小时），为批内重复实验；每次测定均做不同批次质控标本双份，为批间重复实验；一般要求连续测定20天，为天间重复实验，这是对检测系统天间不精密度的观察。??? 　　对精密度实验结果进行统计分析，反映实验仪器精密度好坏的指标是变异系数（CV）。CV越小精密度越好，反之则差，故也称其为不精密度。通常按以下公式可以计算出批内、批间、天间和总CV，其中总CV最重要，它代表整个分析体系的可重复程度。??? 　　S批内=?ni=1 ?2j=1 (Xij1-Xij2)24n??? 　　式中:S批内为批内标准差；n为实验天数(n=20)；i为第i天(1～20)；j为1天内的批数(1或2)；xij1为第i天第j批的第1个结果；xij2为第i天第j批的第2个结果。??? 　　A=?ni=1 （Xi1-Xi2)22n??? 　　式中: A为批间差异水平；Xi1为第i天第1批的结果均数；Xi2为第i天第2批的结果均数。??? 　　B=?ni=1 （i-)2n-1??? 　　式中：B为天间差异水平；i为第i天的结果均数；为所有实验结果均数。??? 　　S总=2B2+A2+S2批内2??? 　　式中：S总为总标准差。??? 　　CV总=S总/ ??? 　　式中：CV总为总变异系数；S总为总标准差；为所有实验结果均数。 　　2? 编写SAS宏程序??? 　　为直接得到如表1所示统计分析报表，编写以下SAS宏程序。数据集名为database，统计变量为var，其中第1批2次测定数据结果分别为var.1和var.2，第2批2次测定数据结果分别为var.3和var.4，输出总变异系数的数据集名为dataout。 　　表1? 精密度试验变异系数统计分析（略） 　　2.1? 定义输出表的格式 　　%macro 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 ou