如何用sas软件正确分析生物医学科研资料-xix．r×c列联表资料的统计分析与sas软件实现（二）

《如何用sas软件正确分析生物医学科研资料-xix．r×c列联表资料的统计分析与sas软件实现（二）》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、396±里垦垫竺技术2012年10月第7卷第5期ChinMedBiotechnol，October2012，Vo1．7，No．5DOh1O．3969／cmba．j．issn．1673—713X．2012．05．015·讲座·如何用SAS软件正确分析生物医学科研资料XIX．R×C列联表资料的统计分析与SAS软件实现(二)王琪，胡良平编者按生物统计学是生物学领域科学研究和实际X-作中必不可少的工具，在分子生物学迅速发展的今天，生物统计学更显示出了它的重要性。实验设计与数据统计分析是现代生物学的基石，是生物学研究者检验假说、寻找模式、建立生物学理论的有利X-

2、具，也是生物学研究者探索微观和宏观生物世界的必备基础知识。对于每天甚至是每时每刻涌现的大量的、以天文数字计量的分子遗传数据，必须借助统计学知识加以分析处理，才能从中获得有意义的信息。“生物多样性数据分析”是开展生物多样性研究的一个重要方面，数据分析能力的高低极大地影响着我们对各种生态学现象认识的深度和广度。现在，电子计算机的普及使得生物统计分析过程大大简化，生物统计分析软件包的普及将生物统计学从统计学家的书本里解放了出来，简化了生物统计分析过程，使之成为生物学研究者的常用工具。本刊特邀军事医学科学院生物医学统计学咨询中心主任胡良平教授，以“如何用SAS

3、软件正确分析生物医学科研资料”为题，撰写系列统计学讲座，希望该系列讲座能对生物医学科研工作者有所帮助。如上期所述，R×c列联表资料可以分为双向无序的料视为双向无序的R×C列联表资料，可根据资料具备的R×C列联表资料、结果变量为有序变量的单向有序R×C前提条件，选用一般检验或Fisher精确检验。列联表资料、双向有序且属性相同的R×C列联表资料和以下将通过实例向读者介绍如何用SAS软件处理双双向有序且属性不同的R×C列联表资料，上期已介绍了向有序且属性不同的R×C列联表资料。前两类。本文将继续介绍第三类，即双向有序且属性不同的【例1】为探讨在夏季三伏天采

4、用三伏灸治疗慢性呼吸系R×C列联表资料及其用SAS软件实现统计分析的全部统疾病的情况。对184例患者进行了分析，数据见表1，内容。患者按年龄分为3组，治疗效果分为三类。当R×C列联表中的两个定性变量，即原因变量与结果变量都是有序变量，并且它们的性质不同，这样的列联表表1患者年龄与治疗效果的情况资料被称为双向有序日_属性／1同的R×C列联表资料。对l『双向有序且属性不同的R×C列联表资料所采用的统计分析方法小能一概而论，应视具体的分析目的而定。一般来说，有4个可能的分析目的lI_：(1)只关心各组结果变量“平均”取值之间的差别是否具有统计学意义。此时，原

5、因变量的有序性就变得无关紧要了，可将此时的“双向有序R×C列联表资料”视为“结果变分析与解答：依据前述4个分析目的，分别实现相应量为有序变量的单向有序R×C列联表资料”，可以选用的的统计分析。统计分析方法有秩和检验、Ridit分析和有序变量的logistic(1)考察各年龄段患者疗效之间的差别有无统计学意义回归分析。SAS程序如下，设程序名为lil．sas。(2)希望考察原因变量与结果变量之间是否存在线性相关关系。此时，需要选用处理定性资料的相关分析方法，通dataal；odshtml；常采用Spearman秩相关分析。doa=lto3：proenpa

6、rlwaywilcoxon；(3)若两个有序变量之间存在的相关关系有统计学意义，dob=lto3：freqt研究者希望进一步了解这两个有序变量之间的变化关系是inputf@@；output；classa：end；end；varb：呈直线关系还是呈某种曲线关系，此时宜选用线性趋势检马令。(4)希望考察各行上的频数分布是否相同。此时，将此资作者单位：100850北京，军事医学科学院生物医学统计学咨询中心通讯作者：胡良平，Emaihlphu812@sina．com中国医药生物技术2012年10月第7卷第5期ChinMedBiotechn，!10b：!：!!!

7、：397cards：rim；程序说明：程序中调用COIT过程，用spearman选项40254odshtm1close；指定进行Spearman秩相关分析。2545lSAS程序运行结果：42515mn：Corr过程程序说明：数据步，建立数据集名为a1，建立数值型变量a、b、f，分别读入行号、列号、每格实际频数；过程步，调用单因素非参数检验过程npar1way，wilcoxon选项要求进行秩和检验，多组比较时进行Kruskal—Wallis检验。Spearman相关系数，N=184由于数据集中的变量是频数，此处用freq语句表示变量b当H0：Rho=0时

8、，Prob>(疗效)出现的频数。最后用class指定分类变量，var指定结果变量。SAS程序运