实操：通过SEER数据库发一篇SCI系列（十）

本周在吴军的《态度》一书里，读到他给他小女儿的一封信，里面这样说：人在年轻的时候，需要做成几件事，并且通过这些成功的过程，学会取得成功的方法。一个人如果做一件事失败了，虽然可以总结经验，吸取教训，但是第二次哪怕他离成功再近，都有可能在最后时刻功亏一篑。

人只有成功过一次，才更容易成功第二次、第三次。因此，失败不是成功之母，成功才是成功之母。

快到三十而立之际，回想过往，竟找不到一件能让自己值得骄傲的事，内心不禁一阵伤感。从小到大，我的人生都是在随大流，偶尔落后，从没领先过。

直至我有幸加入《007-不写就出局》社群，目前已持续践行677天，每周写一篇，输出90+篇文章。作为一名临床医生，虽然可能被周围人认为有些“不务正业”，但我认为这是我最值得骄傲的一件事，因为我和一群陌生而熟悉的人一起在持续做一件事——写作。

大家来自不同地区、干着不同职业，但我们有着共同的目标：7天1篇写7年，7年后一起去南极。

上周把所有肺癌的资料从STAT上按年份下载下来，然后将每年的肺癌TXT数据复制到EXCEL里面，再利用STATA一键筛选出肺类癌病例，最后合并，整理出1996-2015年之间的10043例肺类癌数据。

本周继续上周内容，例如，我想研究下针对I期肺类癌的亚肺叶与肺叶切除的临床预后的比较，可以分三个步骤。

1

第一步，选出原发病，不要合并其他的疾病；

输入命令：

gen ksequencenumber=9

replace ksequencenumber =0 if strmatch(Sequencenumber, '*One*')

drop if ksequencenumber ==9

剩下的“One primary only”就是单发的肺类癌，无其他合并症。

2

第二步，筛选出I期的病例；

输入命令：

gen kDerivedAJCCStageGroup7thed=9

replace kDerivedAJCCStageGroup7thed =0 if strmatch(DerivedAJCCStageGroup7thed, '*IA*')|strmatch(DerivedAJCCStageGroup7thed, '*IB*')

drop if kDerivedAJCCStageGroup7thed ==9

这里抓字段需要注意的地方就是，如果抓 '*IA*'，程序就会连'*IIA*'、'*IIIA*'一起抓了，因此，抓完后还得调整，这里建议先抓'*IIIA*'、'*IIIB*'，再drop，再抓'*IIA*'、'*IIB*'，再drop，这样就剩下'*IA*'、'*IB*'了。

3

第三步，手术方式筛选，亚肺叶VS肺叶，通过查阅STAT手术编码手册，确定“21”、“22”为亚肺叶切除，“30”、“33”为肺叶切除；

输入命令：

gen kRXSummSurgPrimSite1998=9

replace kRXSummSurgPrimSite1998 =1 if strmatch(RXSummSurgPrimSite1998, '*21*')|strmatch(RXSummSurgPrimSite1998, '*22')

replace kRXSummSurgPrimSite1998 =2 if strmatch(RXSummSurgPrimSite1998, '*30*')

replace kRXSummSurgPrimSite1998 =2 if strmatch(RXSummSurgPrimSite1998, '*33*')

drop if kRXSummSurgPrimSite1998 ==9

这样两组数据就出来了，“1”表示亚肺叶切除的病例，“2”表示肺叶切除的病例，这样就可以将两组数据导入SPSS，分析比较它们的五年生存率了，还可以将其他变量做一个多因素回归分析，并做一个列线图，一篇文章是不是就成功了。。

不知不觉，文章都写到系列10了，道理、操作我都会，是该到自己完成自己的第一篇SCI的时刻了。下周我将重点捣弄数据，找出有差异性的研究，直至写一篇SCI投稿。

周末毒鸡汤：学习，从来都是创造惊喜、创造好运的最优路径。

END