How can I calculate the 95% confidence intervals of the hazard ratio from coxphfit?

22 次查看（过去 30 天）

Alexander MacFarlane IV 2018-8-20

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/415574-how-can-i-calculate-the-95-confidence-intervals-of-the-hazard-ratio-from-coxphfit

评论： Victor Fernandes 2023-4-5

Hi everybody.

I am using coxphfit to compute P-values and hazard ratios for data related to clinical trials and need to be able to compute a the 2-sided 95% confidence interval of the hazard ratio.

Any help will be appreciated

5 个评论
显示 3更早的评论隐藏 3更早的评论

Alexander MacFarlane IV 2018-8-21

Wow, thanks!

I see now why you have MVP next to your name.

Star Strider 2018-8-21

My pleasure!

It's an interesting problem, since MATLAB usually provides these, or an easy way to calculate them.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Mau Coen Moraga 2018-8-22

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/415574-how-can-i-calculate-the-95-confidence-intervals-of-the-hazard-ratio-from-coxphfit#answer_333746

编辑：Mau Coen Moraga 2018-8-22

在 MATLAB Online 中打开

Hi Alexander,

My name is Mau and I can confirm that there is no documentation provided to calculate the confidence intervals for the Cox Proportional Hazards Model. Upon further deliberation with my colleagues, we will consider adding the confidence interval functionality to “coxphfit” in a future release of MATLAB.

In the meantime, in order to calculate a confidence interval from the fitted Cox Proportional Hazards Model, you can follow the commands below:

% Generate Weibull data depending on predictor |X|. 
rng('default') % for reproducibility
X = 4*rand(100,1);
A = 50*exp(-0.5*X); 
B = 2;
y = wblrnd(A,B);
% Fit a Cox proportional hazards model. 
[b,logL,H,stats] = coxphfit(X,y);
% Hazard Ratio exp(b)
exp(b)
ans = 
   2.5623
% Calculate 95% confidence interval of Hazard Ratio
CI = exp(b + [-1 1]*1.96*stats.se)
CI =
  2.0030    3.2778