表11.3 MEANS语句与LMEANS语句计算的均值

从MEANS的角度看，药B比药A要好（假设值越大代表效果越好），但是从LSMEANS的角度看，二者无差异。一般情况下，若数据是均衡的，二者的计算结果就是相同；但如果数据不均衡，二者的结果也会不一样。从使用的角度看，若试验是不均衡的，则应该使用LSMEANS。此外，LSMEANS也用于某个因素下水平间差异的检验。例如，假设在方差分析后，我们找出了某一个因素具有显著作用，现欲在分析该因素的两两水平间的差异，这时候就可以用LSMEANS。

水平差异分析

语句LSMEANS的使用语法如下：

LSMEANS 因素，

在方差分析中，常见的选项有两个：PDIFF与SLICE，前者用于水平间的两两比较，后者用于多因素分析，当因素间有相互作用时，可固定一个因素的水平，进而分析另外一个因素的变化。有关SLICE选项的用法会在下一篇文章双因素试验中介绍。需要指出的是，与MEANS语句类似，一个PROC GLM中可以指定多个LSMEANS语句。

PDIFF的常见选项如表11.4所示。

表11.4 选项PDIFF和ADJUST取值

控制选项误差的方法大体上有两种：一种是控制两两比较的误差（Comparisonwise Error Rate，简称CER），另一种是控制整个试验的误差（Experimentwise Error Rate，简称EER）。如果选择的是前者（CER），那么可以使用语句“LSMEANS/PDIFF=ALL ADJUST = T”；若是后者，则可以使用语句“LSMEANS/PDIFF=ALL ADJUST = TUKEY”或者 “LSMEANS/PDIFF = CONTROL(‘控制组’) ADJUST = DUNNETT”，取值“TUKEY”适用于仅考虑水平间两两比较的情形。

例11.4：在例11.3中，我们发现因素Medicine的确是一个显著因素。该因素共有4个水平：A，B，C和D，现欲分析两两比较药的疗效。

示例代码如下：

提交上述代码，在结果查看窗口中，首先可以看到每个组的LSMEANS信息，如图11.8所示。

图11.8 例11.4输出的LSMEANS信息

紧接的是LSMEANS部分主要的输出结果（如图11.9所示），一个4*4阶矩阵，记为，这里的i，j=1，2，…， 4。下标值i, j对应的信息由上一个表格中的列LSMEAN号确定。原假设为i和j的LSMEANS值相等（无显著性差别），显著性水平仍然为默认值0.05。

图11.9 例11.4变量MEDICINE水平间差异分析报表

由于=0.9960 > 0.05，所以我们接受原假设，认为二者间无显著性差别；同理，由于=0.0040 < 0.05，拒绝原假设，认为二者有显著性差别。类似地可以对矩阵中的其他元素值进行分析。

语句LSMEANS中的选项“pdiff = All”要求ODS 选项是打开的，其作用是生成如图11.10所示的图表。

图11.10 例11.4不同水平间均值差异的置信区间图

在图11.10中，斜向上的对角线（虚线）是组间差异为0的参考线。斜向下的实线代表了CLASS变量下两个不同水平间的均值。本例中，CLASS变量下共有4个不同水平，可以组成6对不同水平组合。每一对水平组合下的两个水平之间的均值差异的置信区间由图中斜向下的实线表示。以图中直线L为例，直线L通过竖直线D和水平线C的交点，因此，直线L对应的是水平D和水平C均值差异的置信区间。直线L与斜向上的虚线没有交点，即0不在该置信区间内。据此，我们可以推断水平D和水平C的均值存在显著差异。简而言之，如果斜向下的实线与参考线不相交，就表示对应于该实线的两个水平的均值是有显著差异的；反之，若实线与参考线相交，则该实线对应的两水平均值没有显著差异。

本文结束，下一篇文章开始介绍“双因素试验的方差分析”。

回复【数据和代码】可以下载《深入解析SAS》一书的数据和代码。

查询往期文章，请回复下列关键字：

【安全】==>【SAS智能平台安全管理】系列文章

【Base基础】==>【Base SAS基础】系列文章

【编程概念】==>【SAS编程基本概念】系列文章

【描述性统计】==>【SAS统计分析系列：描述性统计分析】系列文章

【读外部数据】==>【读取外部数据到SAS数据集】系列文章

【假设检验】==>【参数估计与假设检验】系列文章

小贴士

读者可以从以下链接获取SAS公司提供的免费版环境：

SAS大学版（SAS® University Edition）是SAS为在校大学生免费提供的基于虚拟机和网页的SAS环境。回复关键字【大学版】，可以查看详细介绍。

下载路径：

https://www.sas.com/en_us/software/university-edition/download-software.html

SAS学术版（SAS® OnDemand for Academics）是 SAS 为学术届人士免费提供的、在线的、基于SAS 私有云上的应用服务环境。

用户首先需要注册，然后按照提示信息就可登录。

注册路径：

http://odamid.oda.sas.com

本文转自《深入解析SAS — 数据处理、分析优化与商业应用 》

作者：夏坤庄、徐唯、潘红莲、林建伟

如若转载本文，请在文章顶部标注 “本文转自SAS知识 （ID: SASAdvisor），摘自《深入解析SAS — 数据处理、分析优化与商业应用 》”

作者介绍

夏坤庄

《深入解析SAS — 数据处理、分析优化与商业应用》第一作者， SAS软件研究开发（北京）有限公司客户职能部总监。在承担研发工作的同时，夏及其团队负责对SAS非英语市场提供技术支持，并且与在美国及其它地区的团队一起，服务于SAS的SaaS/RaaS业务，同时提供和验证关于SAS产品和技术在应用领域的最佳实践。在加入SAS软件研究开发（北京）有限公司之前，夏就职于SAS中国公司，历任资深咨询顾问、项目经理、首席顾问、咨询经理，拥有丰富的咨询和项目实施经验。在长期的从业经历中，不但为SAS的金融行业客户成功实施了众多深受好评的项目，而且在近年领导实施了非金融行业的多个大数据分析项目。

欢迎大家投稿，一起分享SAS的点滴

投稿邮箱: [email protected]返回搜狐，查看更多