Last chance! 50% off unlimited learning
Sale ends in
locScore(x, ...)## S3 method for class 'coxph':
locScore(x, ..., all = FALSE, hypo = NULL,
ties = c("breslow", "efron", "exact"))
locLR(x, ...)
## S3 method for class 'coxph':
locLR(x, ..., all = FALSE, hypo = NULL,
ties = c("breslow", "efron", "exact"))
locWald(x, ...)
## S3 method for class 'coxph':
locWald(x, ..., all = FALSE, hypo = NULL)
coxph
TRUE
anbreslow
, efron
or exact
.locScore
a list
with the following elements,
which are data.table
s:locLR
and locWald
, a data.table
showing the hypothesis being tested and the results of the test.colSums(survival::coxph.detail(x)$score)
and x$var
respectively).
For the likelihood ratio test, the model is also refit with
the coefficients of interest excluded. The likelihood ratios
from the full model and those with coefficients excluded
are used to construct the test:
data(larynx, package="KMsurv")
c1 <- coxph(Surv(time, delta) ~ factor(stage) + age, data=larynx)
locScore(c1, all=TRUE)
locScore(c1, hypo=c(0, 0, 0, 1))
locScore(coxph(Surv(time, delta) ~ stage + age, data=larynx))
###
data(larynx, package="KMsurv")
c1 <- coxph(Surv(time, delta) ~ factor(stage) + age, data=larynx, method="breslow")
locLR(c1, all=TRUE)
locLR(c1, hypo=c(FALSE, FALSE, FALSE, TRUE))
###
data(larynx, package="KMsurv")
c1 <- coxph(Surv(time, delta) ~ factor(stage) + age, data=larynx, method="breslow")
locWald(c1, all=TRUE)
locWald(c1, hypo=c(0, 0, 0, 1))
Run the code above in your browser using DataLab