精密度实验变异系数计算的SAS宏程序(1)讲解

1、精密度实验变异系数计算的SAS宏程序】评价医院实验室检测结果的准确性和稳定性，需进行精密度实验，以确定检测结果是否处于所控制的范围内。通过计算精密度实验批内、批间、天间 以及总的变异系数，能够反映实验仪器精密度好坏。通过编写SAS宏程序，可以应用SAS统计软件直接输出以上变异系数的统计报表。【关键词】 精密度实验； 变异系数；统计分析报表；SAS宏程序1精密度实验评价医院实验室检测结果的准确性和稳定性，需进行精密度评价实验，以 确定检测结果是否处于所控制的范围内。精密度实验通常包括批内、批间以及 日间重复实验。对同一批次质控标本的重复测定，求每天在不同时间点测定同 一批次质控标本2次（2次测定

2、间隔不得少于2小时），为批内重复实验；每 次测定均做不同批次质控标本双份，为批间重复实验；一般求连续测定20天,为天间重复实验，这是对检测系统天间不精密度的观察。对精密度实验结果进行统计分析，反映实验仪器精密度好坏的指标是变异 系数（CV。CV越小精密度越好，反之则差，故也称其为不精密度。通常按以 下公式可以计算出批内、批间、天间和总 CV其中总CV最重，它代表整个分 析体系的可重复程度。S 批内=ni=12j=1 (Xij1-Xij2)24n式中:S批内为批内标准差；n为实验天数（n=20） ； i为第i天（120） ； j 为1天内的批数（1或2） ； xij1为第i天第j批的第1个结果；

3、xij2为第i天 第j批的第2个结果。A= ni=1 (Xi1-Xi2)22n式中：A为批间差异水平；Xi1为第i天第1批的结果均数；Xi2为第i天第2批的结果均数B= ni=1 (i-)2 n-1式中：B为天间差异水平；i为第i天的结果均数；为所有实验结果均数。S 总=2B2 A2 S2 批内 2式中：S总为总标准差。CV总=5总/式中：CV总为总变异系数；S总为总标准差；为所有实验结果均数。2编写SAS宏程序为直接得到如表1所示统计分析报表，编写以下 SAS宏程序。数据集名为 &database，统计变量为&var，其中第1批2次测定数据结果分别为&var.1和 &

4、amp;var.2，第2批2次测定数据结果分别为&var.3和&var.4，输出总变异系数的 数据集名为&dataout。表1精密度试验变异系数统计分析(略)2.1定义输出表的格式%macro tformat; /* 定义宏，输出统计报表的格式*/proc format;inv alue g1=202=40;%me nd tformat;2.2 计算变异系数调用proc uni variate过程计算批内标准差s、批间差异水平a、天间差异水平b和所有实验结果均数x,利用公式计算出总的变异系数 CV%macro cv(database,var,dataout); /* 定

5、义 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(kee

6、p=&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 &va

7、r.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 ，

8、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;row &

9、;p;put #row 4 "&cvar" 22 s 6.2 '%' 34 a 6.2 '%' 50 b 5.2 '%' 64cv 5.2 '%'#(row 1) 2 75*'-'end;if inds=1 then inds=0;end;%mend fc;定义宏nullset，产生一个输出表，集成已产生的批内标准差 s、批间差异 水平a、天间差异水平b和总的变异系数CV并按定义排列。将输出表存入 d 盘 sas 目录下，文件名为 &tab 的文本文件， &tab 为宏变

10、量。%macro nullset(data);data _null_;file "d: sas &tab.txt" print n=ps no titles header=head;set &data;inds 1;col=input(g,g.);【摘】 评价医院实验室检测结果的准确性和稳定性，需进行精密度实验， 以确定检测结果是否处于所控制的范围本篇论文是由3COM文档频道的网友为您在网络上收集整理饼投稿至本站 的，论文版权属原作者，请不用于商业用途或者抄袭，仅供参考学习之用，否 者后果自负，如果此文侵犯您的合法权益，请联系我们。%mend nullset;2.4 运行宏程序将以上SAS程序提交SAS系统运行，即可自动生成统计分析表1。其中数据集名为main，变量名分别为high、mid、low。%let tab=