spss一键单因素logistic cox分析与基线表

前几期我们讲了R、SAS一键绘制基线表和单因素多因素分析，见：

什么，你还在手动做SCI表格？

R软件一键绘制Table，解放医生双手

加上本章节，统计软件三剑客的一键基线系列就完整了。

R SAS虽然足够灵活强大，但是有些人还只会用SPSS。本章旨为解决编程能力较弱又能快速整理出论文主要基线结果，使用更加接地气临床更加常用的SPSS一键绘制基线表及单因素logistic、cox分析。

1 .基线表绘制

课本上甚至大部分SPSS教程都介绍的是explorer探索分析、freqencey描述基线。但是一个个模块一次次点是很琐碎的事情，即使用.sps语法去编辑，仍然需要重复多次并且手动录入到paper里面，后期修改起来遇到几十个上百个变量的时候稍显麻烦。演示数据采用张文彤老师的《SPSS高级教程第二版》part2部分logistic_step.sav。数据结构如下：

本章节利用“表→定制表”快速绘制基线分布及单因素p值(为t检验及X2检验)。一般来说单因素基线的定制表能解决大部分问题，但是有时候连续变量极度偏态需要考虑非参数秩和检验，单因素常常选择p<0.05 或<0.1入选多因素分析，样本量很大时也会见到<0.2或0.25。

图1 定制表界面

图2 定制表连续和分类变量勾选参数

对应.SPS代码为：

• * 定制基线表， y=low， x为其他. 
  *   连续变量格式 age [MEAN, STDDEV, MEDIAN, PTILE 25, PTILE 75] 
      分类变量格式  ht [COUNT F40.0, COLPCT.COUNT PCT40.1] .
  CTABLES 
    /VLABELS VARIABLES=age lwt race smoke ptl ht ui ftv bwt low 
      DISPLAY=LABEL 
    /TABLE age [MEAN, STDDEV, MEDIAN, PTILE 25, PTILE 75] + lwt [MEAN, STDDEV, MEDIAN, PTILE 25, 
      PTILE 75] + race [COUNT F40.0, COLPCT.COUNT PCT40.1] + smoke [COUNT F40.0, COLPCT.COUNT PCT40.1] + 
      ptl [MEAN, STDDEV, MEDIAN, PTILE 25, PTILE 75] + ht [COUNT F40.0, COLPCT.COUNT PCT40.1] + ui [COUNT 
      F40.0, COLPCT.COUNT PCT40.1] + ftv [COUNT F40.0, COLPCT.COUNT PCT40.1] + bwt [MEAN, STDDEV, MEDIAN, 
      PTILE 25, PTILE 75] BY low 
    /CATEGORIES VARIABLES=race ORDER=A KEY=VALUE EMPTY=INCLUDE TOTAL=YES POSITION=AFTER 
    /CATEGORIES VARIABLES=smoke ht ui ftv low ORDER=A KEY=VALUE EMPTY=INCLUDE 
    /CRITERIA CILEVEL=95.~~~
  

图3 一键基线分布情况与单因素p值(t、X2) 

注：t检验是双侧。卡方可选择bon、benj等多重校正。

2.一键单因素logistic表绘制

有时候研究基线展示的p值是单因素OR95% CI。如果变量少，spss菜单还算非常方便。但是变量几十上百个时候，有没有快捷方法像R、SAS一样一键绘制单因素OR%95CI呢？答案是有的。笔者利用OMS结合SPSS宏，编制了一键单因素OR95%CI的程序，仅仅需要变量视图里面赋值variable，列替换，红框位置的变量即可。同理可制作HR 95%CI。

图4 spss宏程序一键单因素logistic

OMS过程如下图。OMS是spss非常智能的一个高级功能，比如缺失值多重填补后迭代5次取均值或中位数或众数都可以用OMS进行获取计算。其他统计开发功能可根据需求自行获取中间步骤。

图5 OMS模块

本程序直接输出结果至excel表，香不香？

小结

1. spss利用“表格→定制表”可快速制作基线分布表，十分便捷。
2. spss宏程序仅需替换数据变量即可快速实现单因素logistic、cox分析。
   
   

本章节详细程序请添加笔者公众号“cardiology心思维”，发送消息：一键单因素logistic