采用均数和统一标准差进行方差分析的方法及SAS实现.docVIP

采取均数和标准差进行方差分析方法及SAS实现 广东省深圳市中医院郭志武 【摘要】目标探讨采取均数和标准差进行方差分析方法及SAS程序。方法经过推导，得到基于均数和标准差方差分析计算公式，并用简练SAS程序实现。结果实现了相关算法及SAS程序。结论采取均数和标准差进行方差分析是完全可行，编制SAS程序能够有用实现相关计算，得到可靠结果。 【关键词】均数标准差方差分析SAS 在实际工作中，通常会碰到某些研究资料只有均数和标准差而缺失原始数据情况，比如年代久远研究资料、期刊文件资料等等。对于这些资料，若需进行t检验是很轻易，有现成公式可用；但若需对多个样本资料做方差分析就不易了，现有计算公式和统计软件算法全部是基于原始数据，所以有必需探讨基于均数和标准差方差分析算法，以满足实际工作需要。 1.算法推导 对于单原因方差分析，传统算法是先计算总变异，然后分别算出组间变异和组内变异，最后求出F值和对应P值。这全部是基于原始数据来计算。若只有均数和标准差，则需要经过推导得出计算组内变异和组间变异公式。 组间变异理论公式为： （1） 而，所以： （2） 组内变异理论公式为： （3） 因为组内标准差，所以，从而： （4） 这样，经过（2）式和（4）式，能够很方便地计算出和，从而深入计算出F值。 至于双原因方差分析和多原因方差分析，可作类似推导，这里不再赘述。 多个样本两两比较能够遵照现有公式完成。 2.SAS程序实现 下面以《卫生统计学》第四版（倪宗瓒主编）第51页样例数据为例，介绍采取均数和标准差进行方差分析SAS程序实现。 /*创建数据集，输入各组组号、样本例数、均数、标准差，组号需按ABCD...次序输入*/ dataa; inputgp$nXS; cards; A11105.4510.87 B9102.3914.55 C10122.8017.07 ; /*采取sql过程进行方差分析计算*/ procsql; createtablebas selectsum(n)asN,/*计算总例数N*/ count(n)ask,/*计算组数k*/ sum(n*X**2)-sum(n*X)**2/sum(n)asSSZJ,/*计算组间变异*/ sum(S**2*(n-1))asSSZN/*计算组内变异*/ froma; quit; datac; setb; MSZJ=SSZJ/(k-1);/*计算组间均方*/ MSZN=SSZN/(N-k);/*计算组内均方*/ F=MSZJ/MSZN;/*计算F值*/ P=1-probf(F,k-1,N-k);/*计算P值*/ procprintdata=c; varFP;/*显示F值，P值*/ run; /*以下进行两两比较计算，采取LSD-t检验*/ datad; setc; keepNkMSZN; procsql; createtableeas selecta1.gpasgp1,a2.gpasgp2,/*交叉匹配生成两两比较数据集*/ a1.nasn1,a2.nASn2, a1.XASX1,a2.XasX2 fromaasa1,aasa2 wherea1.gpa2.gp; createtablefas selecte.*,d.*frome,d;/*生成两两比较计算表*/ quit; datag; setf; dmean=x1-x2; sd=sqrt(MSZN*(1/n1+1/n2)); t=dmean/sd; p=(1-probt(abs(t),n-k))*2; procprintdata=g; vargp1gp2n1n2x1x2dmeantp;/*显示两两比较结果*/ run; 3.讨论 参考文件 倪宗瓒主编.卫生统计学[M].第四版.北京：人民卫生出版社， .