干货|多因素做cox生存分析,怀疑数据共线性了怎么办?——多重线性回归(一)

昨天从一个坑爬出来,今天又掉进另一个坑里

问题就是这么个问题,生存回归单因素分析的95%CI和HR,为什么和cox多因素算出来的不一样?丁香园的人说,这是因为共线性的问题。

于是小张又踏上了解决共线性的问题的道路

此文贼长,看到最后有个小彩蛋

百度文库里有一个很详细全面的文档(认真一看,其实也不是很详细),丁香园子里还有个老师是用sas做的,但是那个文档看起来很头晕

多重共线性问题的几种解决方法 - 百度文库​wenku.baidu.com

还是我自己来,首先你要怎么的得知自己的几个变量它就共线性了呢

完全不知道这个是在说什么

让我们回到COX需要满足的条件

接着就是——多重线性回归看是否符合假设3-8

经过上述操作,dataview出现了五个新变量

这5个变量分别是未标化预测值(unstandardized predicted values,PRE_1)、学生化残差(studentized residuals SRE_1)、学生化删除残差(studentized deleted residuals,SDR_1)、Cook距离(Cook's Distance values,COO_1)以及杠杆值(leverage values,LEV_1)。

根据这5个新增变量和其它结果,我们将逐一对假设3-8进行检验。

假设3.DURBIN-WATSON检验可用于观测变量间是否有影响,一般处于0-4之间,越接近2,观测值相互独立的可能性越大

假设4.因变量和自变量之间存在线性关系

(1)因变量和所有自变量之间的关系

一般来说,如果未标化预测值(PRE_1)和学生化残差(SRE_1)的散点图呈水平带状,就说明多重线性回归中因所有自变量之间存在线性关系。结果提示,本研究满足因变量与所有自变量之间存在线性关系的假设。

(2).因变量与每一个自变量是否有关

假设5.残差的方差齐

这些就属于方差不齐

假设6.不存在多重共线性

一般来说,如果自变量之间的相关系数大于0.7,则可能存在多重共线性。本研究中,任意两个自变量的相关系数都小于0.7,提示不存在多重共线性。

点击Collinearity diagnosis,SPSS可以输出容忍度(Tolerance)和方差膨胀因子大于10。

实际上,方差膨胀因子是容忍度的倒数(1/容忍度),我们只需要判断其中一个指标即可。如果容忍度小于0.1(或方差膨胀因子大于10),提示数据存在多重共线性。

在本研究中,所有容忍度值都大于0.1,提示数据满足假设6。

假设7.没有显著异常值

1.离群值

方法1.Casewise Diagnostics检验标准是上下3倍标准差,并标记超出此范围的数据为离群值。

如果研究中没有标化残差超出3倍标准差的离群值,SPSS就不会输出该表格。本研究中没有标化残差超出3倍标准差的离群值,SPSS也未输出Casewise Diagnostics检验结果

方法2.将学生化删除残差降序(或升序)排列后,在数据栏的最上方和最下方检查是否存在超过3倍标准差的离群值。如果存在,就应进行相应的调整或剔除。

2.高杠杆点

自变量的异常值即高杠杆点。高杠杆点可以通过杠杆值检测。Distance框内点选Leverage values后,新增变量(LEV_1)就是杠杆值。Data View窗口中,右击LEV_1数据栏,选择"Sort Descending"(或Sort Ascending)。将杠杆值降序(或升序)排列后,检查杠杆值。

一般来说,如果杠杆值小于0.2,可以认为数据安全。但如果杠杆值位于0.2-0.5之间,则认为对应数据可能是高杠杆点;若杠杆值大于0.5,数据就很可能是高杠杆点了。

我发现自己的文章有一堆高杠杆点

3.强影响点

COO_1,小于1就认为ok

假设8.残差近似正态分布

离对角线越近,越认为是正态分布

PP图比直方图更直观

近似即可

太累了,不能连续打出一段话了,以上,我看完资料的三点总结

八个假设走一遭



文末的小彩蛋~啾

没了,哈哈哈

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

推荐阅读更多精彩内容

  • R中的线性回归函数比较简单,就是lm(),比较复杂的是对线性模型的诊断和调整。这里结合Statistical Le...
    真依然很拉风阅读 65,120评论 1 64
  • 1、OLS线性回归的基本原则 最优拟合曲线应该使各点到直线的距离的平方和(即残差平方和,简称RSS)最小。 2、O...
    wonphen阅读 1,552评论 0 2
  • 最近用R和python做了一些回归分析的项目,发现这两个工具各有特点,如果目的是分析研究用R即可,毕竟R语言有很多...
    阿童89阅读 5,948评论 1 2
  • title: "lesson 8 回归"author: "wintryheart"date: "2019/9/9"...
    冬之心阅读 9,935评论 0 5
  • 人说鸟儿生来自由 我说鸟儿生来忧愁 小小鸟想要飞得更高 可是天空之下有猎人的枪 天空海阔任尔遨游 无限带来的是离愁...
    无人问津的小姐姐阅读 198评论 0 0