文档章节

Coxph model Pvalue Select

o
 osc_y8yehimr
发布于 2019/03/20 11:13
字数 267
阅读 12
收藏 0

精选30+云产品,助力企业轻松上云!>>>

 

I am calculating cox propotional hazards models with the coxph 
function from the survival package.  My data relates to failure of 
various types of endovascular interventions.  I can successfully 
obtain the LR, Wald, and Score test p-values from the coxph.object, as 
well as the hazard ratio as follows: 

formula.obj = Surv(days, status) ~ type 
coxph.model = coxph(formula.obj, df) 
fit = summary(coxph.model) 
hazard.ratio = fit$conf.int[1] 
lower95 = fit$conf.int[3] 
upper95 = fit$conf.int[4] 
logrank.p.value = fit$sctest[3]  (Score (logrank) test)

wald.p.value = fit$waldtest[3]  (Wald test)

lr.p.value = fit$logtest[3]    ----    Likelihood ratio test (优先)



I had intended to report logrank P values with the hazard ratio and CI 
obtained from this function.  In one case the P was 0.04 yet the CI 
crossed one, which confused me, and certainly will raise questions by 
reviewers.  In retrospect I can see that the CI calculated by coxph is 
intimately related to the Wald p-value (which in this specific case 
was 0.06), so this would appear to me not a good strategy for 
reporting my results (mixing the logrank test with the HR and CIs from 
coxph). 

I can report the Wald p-values instead, but I have read that the Wald 
test is inferior to the score test or LR test.  My questions for 
survival analysis jockeys out there / TT: 

1. Should I just stop here and use the wald.p.value?  This appears to 
be what Stata does with the stcox function (albeit Breslow method). 

2. Should I calculate HR and CIs that "agree" with the LR or logrank 
P?  How do I do that? 

Thank you, 

o
粉丝 0
博文 500
码字总数 0
作品 0
私信 提问
加载中
请先登录后再评论。
通过sqlHelper访问数据库,查询数据

使用最简单的方式进行对数据库访问,并且查询返回需要的数据;如下步骤: 一、web.config 配置 <connectionStrings>   <add name="ABConnection" connectionString="Data Source=.;Initia......

osc_7shyb795
2019/08/24
11
0
多元逐步回归算法

python Mutiple-Stepwise-Regression-Analysis(MSRA) 先谈一下个人对多元逐步回归的理解:多元逐步回归的最本质的核心是最小二乘原理,本方法中调用smf方法。 注:输入数据的结果如下图所示:...

草木的执念
2019/04/12
72
0
第五篇,jroo自定义标签篇

1、jroo中有不少自定义标签,如数据字典标签(select、checkbox、radio、view),权限标签,分页/排序/序号标签,国际化标签,位置标签等。 2、自定义标签篇 数据字典中设置的是系统的常量数据...

yeshujun
2012/12/31
720
1
回归分析特征选择(包括Stepwise算法) python 实现

-- coding: utf-8 -- """ Created on Sat Aug 18 16:23:17 2018 @author: acadsoc """ import scipy import numpy as np import pandas as pd import matplotlib import matplotlib.pyplot a......

osc_v22siqak
2018/08/25
1
0
参数按照ASCII码从小到大排序

public static void requestPay() { Dictionary<string, string> dics = new Dictionary<string, string>(); dics.Add("amount", amount); dics.Add("callbackurl", callbackurl); dics.Add(......

osc_7c6h7z9p
2019/04/04
0
0

没有更多内容

加载失败,请刷新页面

加载更多

是否有可能从另一个git存储库中挑选一个提交? - Is it possible to cherry-pick a commit from another git repository?

问题: I'm working with a git repository that needs a commit from another git repository that knows nothing of the first. 我正在使用一个git存储库,需要从另一个不知道第一个存储库......

技术盛宴
昨天
26
0
【LeetCode】53 盛最多水的容器

题目 解题思路 双指针法: https://leetcode-cn.com/problems/container-with-most-water/solution/sheng-zui-duo-shui-de-rong-qi-by-leetcode-solution/ 代码 public class Solution { ......

JaneRoad
昨天
16
0
阿里云OSS配置CDN加速

首先购买CDN流量包 然后添加域名 添加好后 然后将域名OSS.xxxx.com 解析到 生成的CDN域名上 这样就完成了

可达鸭眉头一皱
昨天
16
0
js 整数与小数正则替换片段

说明 /(\d+)/g 整数 /(\d+\.\d+)rem/g 小数 /(\d+\.\d+|\d+)rem/g 其中 | 或 条件 例子 全局查找带 rem 单位的,替换成 px 单位 let text = text.replace(/(\d+\.\d+|\d+)rem/g, function(s......

DrChenXX
昨天
17
0
ubuntu下minicorba例子

一、开发环境安装 sudo apt install omniorb omniorb-idl omniidl libomniorb4-dev libomniorb4-2 omniorb-nameserver libomnithread4 libomnithread4-dev 二、源文件: Hi.idl module ......

wangxuwei
昨天
18
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部