临床生信文章普遍存在的两个统计问题

临床预后模型,大多包含风险分层和独立预后分析两个验证分析,而很多文章在两个分析上存在不严谨。下面介绍一些分析的关键点。

1. 风险分层与cutoff

一般而言,预后模型会有训练集和验证集,有一些还分为了训练集、验证集和测试集。首先,在有多个数据集的情况下,对于表达谱数据,需要进行批次校正;其次,如果仅有训练集和验证集(两个独立队列),cutoff不必保证相同,而如果有训练集、验证集和测试集(其中训练集和验证集是一个队列的拆分),训练集和验证集的cutoff需要保持一致,测试集的cutoff不强求。

2. 独立预后分析

(1) PH假定

很多文章都会在构建了prognostic scoring system的基础上,将临床病理性状纳入做单因素和多因素回归,也就是独立预后分析,有的还会基于多因素Cox模型构建nomogram,但事实是大多数文章的结果不可靠。我们知道Cox回归需要满足PH假定,而对PH假定的检验主要有LML曲线法、KM曲线法和Schoenfeld's 检验法。
对于分类数据,如性别、病理分级、组织学分期、有无突变等等,用LML曲线法和KM曲线法检验,若曲线互不相交、互相分离,则一般认为满足PH假定,可以纳入Cox模型;对于连续型变量,如年龄、构建的评分系统、肿瘤大小等等,最好使用三种方法都进行检验,其中Schoenfeld's 检验的p值需要>0.05,若p<0.05,则纳入该因素的结果是不可靠的。

如果想探索的因素不满足PH假定怎么办?

网上有一些解决办法,如分层回归和时依回归,但这两种方法操作起来相对困难且结果不易解释。有一种比较简单的方法大家可以试一试,但仅针对于连续变量。可以将连续变量变为分类数据,如年龄分为>60和<=60岁,risk score分为高低风险组,此时再用LML和KM检验法进行验证。如果对于分类变量还不满足PH假定,那只能考虑剔除这个因素或者进行分层/时依回归了。

(2) 等级变量

对于分期、分级这样的等级变量,很多文章直接将数据变为0、1、2、3等,这是严谨的。应该将该变量设置成哑变量,设置其中一个哑变量为参考组 (reference),在数据科学中,称作热编码 (one-hot encoding)。

虽然以上分析比“不严谨”的方法操作起来多了一些步骤,但一旦被审稿人抓到把柄,返修起来将更加繁琐,所以还是一次做到位吧!

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,509评论 6 504
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,806评论 3 394
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 163,875评论 0 354
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,441评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,488评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,365评论 1 302
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,190评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,062评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,500评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,706评论 3 335
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,834评论 1 347
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,559评论 5 345
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,167评论 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,779评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,912评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,958评论 2 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,779评论 2 354

推荐阅读更多精彩内容

  • 1.1假设检验,估计和预测 统计学包括了对其它领域研究的设计、假设检验、估计和预测等方面。本书着重讲述估计和预测,...
    王子威PtaYoth阅读 684评论 0 0
  • 一、生存分析(survival analysis)的定义 生存分析:对一个或多个非负随机变量进行统计推断,研究...
    Seurat_Satija阅读 8,203评论 0 12
  • 统计学词汇中英文对照完整版 A Absolute deviation, 绝对离差 Absolute number,...
    生信F3阅读 5,315评论 0 4
  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 6,049评论 0 4
  • 公元:2019年11月28日19时42分农历:二零一九年 十一月 初三日 戌时干支:己亥乙亥己巳甲戌当月节气:立冬...
    石放阅读 6,879评论 0 2