0
$\begingroup$

Let's say I have a model (pseudolanguage): survival = My_Factor + sex + age.

I do it to adjust the hazard ratio for sex and categorized age={<30yr; >=30yrs} (not my choice). I can get two p-values:

The one for the My_factor coefficient. That's the Wald's one. I guess it's now adjusted for the others, as the HR changed.

The overall Wald and Log-rank (called "score") test. But I believe the two are just global tests about all coefficients (b1=b2=b3=0).

Is there any way to get adjusted logrank? Or should I report just the beta's one (Wald)?

enter image description here

$\endgroup$
0

1 Answer 1

0
$\begingroup$

You are correct that the 3 tests reported at the bottom of your display are for the overall model.

Usual practice is to report the Wald test values for individual coefficients. Those are easy to get, as all you need is the covariance matrix of the coefficient estimates from the final model. They also can be combined easily into a test on multiple coefficients, for example combining all levels of a multi-level categorical predictor. The Wald test assumes, however, a symmetric shape of the likelihood profile with respect to changes in each predictor around its estimate.

With small sample sizes, a likelihood-ratio test is the gold standard. To evaluate a specific predictor of interest, you can do such a test on models with and without that predictor. The related profile-likelihood confidence intervals (CI) can also be more reliable than the Wald CI, as they make no assumptions on profile shape. Instead, they directly evaluate the likelihood as a function of predictor value (effectively refitting the model for each of a set of predictor values). Getting those CI can, however, be computationally intensive so they aren't so often calculated or reported.

Score tests on individual predictors are possible, as explained on this page that covers all 3 tests. They are done starting with a model that omits the predictor of interest, however, so there is no practical advantage over a likelihood ratio test, and it provides results that are less reliable. Can't say that I've ever seen individual-coefficient score tests reported for a multiple-regression Cox model.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.