SAS-给公众号做一个秩和检验

日前，小编的公众号的关注人数终于突破了2000了！但还远远没有完成今年的年度计划，于是小编就想看一看，每日增粉数量是否与关注基数有显著的统计学意义~好吧，作为一个不会统计学的菜鸟，写这篇推送的目的是希望各位老师能够指点一下小编，并检阅一下小编给自己布置的作业是否做对了，还请各位老师不吝赐教！

嗯，于是小编从公众号上下载了自2017年11月11日-2018年03月25日的公众号每日增粉相关的数据...接着小编就开始分组了，以500人为区间，分成3个组进行对照研究（group1:

嗯，看这里的P值是小于0.0001的，得出的结论大概是有显著的统计学意义的，为啥小编这里用秩和检验呢，听说用秩和检验不用考虑数据的正态性，所以就它了。从结果可以看出，均值是逐渐增大的，那么是不是可以说明随着关注量的增大，每日增关量也会增加呢？

接着呢，小编想看一看增粉数量与星期有没有显著性差异，于是小编用了卡方检验，同样得出了P值小于0.0001，是的，有显著性差异！那么还能看出啥呢？小编就不太知道了！

在就是结果中三线表的输出，三线表的输出小编以前虽然写过相关的推送，不过吧，现在水平又提升了一丁点！也正好想找点数据来练手，于是小编就将在本文中在一次的写一写report过程，以及将以前写过的一些基础的专题与推送，简单的连贯起来进行一次实践性质的编程！

首先，我们在写程序的时候，做的一件事是建立文件夹以及逻辑库路径！

嗯，这是我建立的文件夹，数据、文档、日志、宏、输出、程序、源数据都有自己的夹子。如果有兴趣的朋友，可以在公众号对话框回复：文件夹1，既可以下载本文涉及的全部程序以及数据以及macro程序包！解压后程序可以直接运行。

回到正文上来，下载数据 -数据导入：

这段代码用到了一个系统宏变量（sysver）来获取你SAS的版本号，小编电脑上装了SAS9.2与SAS9.4，因此将宏分别执行成2个版本宏包（执行后的宏包不能夸版本），所以小编这里利用这个宏变量自动获取当前SAS软件的版本号，然后进行判断进而选择正确的宏包。嗯，下载文件夹后，解压压缩包，不改变压缩包内部文件以及文件夹，在目录下的PGM\Table\Table.sas程序打开是可以直接运行的。

数据导入后呢，数据集的转换以及建立分组等：

接着就到了描述性统计量的计算，这里用mean过程步：

计算后的结果：这里为啥小编需要给数据转换成字符型的变量呢，主要原因是这样的，后面小编还要进行数据集的追加，变成字符变量比较好操作！

看到上面的结果：是不是和前面RTF中的排版结构相差很大呢，那是因为没有转置，接下来就来转置一下：

结果是这样的：

现在就和前面RTF输出的结构差不多了，好了，然后在进行P值的计算：

运行完后的结果是这样的： 下面首列的一串字符是啥意思呢，这是RTF标记语言，实现一些排版缩进等操作！

这里一部分算是完成了，接下来在来看第二部分。卡方检验结果的输出！其实和上面的也是很类似，都整到数据集中，在进行转置啊等几步数据结构的操作，这样就可以实现想要输出结构的排版，这里就不细说，直接贴代码了！

最后在来看看要输出的final数据集：

做到这一步，整个数据的操作基本上是完成了，接下来就是数据集的输出，输出到RTF中。这里小编就用了直接写好的宏进行输出。

首先，我先导入输出RTF中的标题、备注等信息，前面可以看到小编输出的rtf里面是带有备注等信息的，其实就是通过这里控制。顺便看一看titles的Excel模板是啥样的。为啥要做Excel中呢，放在外部，其实也便于修改与管理，以及用Macro来实现自动添加这些信息等等！第二个sheet有一些常用的RTF标记的例子，方便忘记的时候随时查找！

紧接着，就是一些Macro的调用了，这里就不过多的说了，有兴趣的朋友可以私下找我讨论。

然后就到了proc report过程步了...

到这里，基本上整个程序是写完了！如对程序有兴趣，可以在可以在公众号对话框回复：文件夹1。

小编是初学者，各位老师如果觉得文中涉及统计存在问题，欢迎指正！以免其他初学者被小编带偏了！

今天就这么多了，后续内容，敬请期待~

如果觉得好或者公众号有价值：请打赏我把~

打赏包括但不限于：转发，分享，推荐，赞赏本公众号或者文章，谢谢！