hr.test: Hazard ratio test

In the-mad-statter/washu: Helpful Functions for WashU Statistical Consultations

Hazard ratio test

Description

Hazard ratio test

Usage

hr.test(beta1, se1, beta2, se2, r = 0, n = 10000)

Arguments

beta1	beta estimate for the first hazard ratio
se1	standard error of beta estimate for the first hazard ratio
beta2	standard error of beta estimate for the second hazard ratio
se2	standard error of beta estimate for the second hazard ratio
r	correlation between the first and second risk scores (i.e., predictors)
n	number of samples from which the hazard ratios were estimated

Value

A list with class "htest" containing the following components:

statistic	the value of the t-statistic.
parameter	the degrees of freedom for the t-statistic.
p.value	the p-value for the test.
conf.int	a confidence interval for the mean appropriate to the specified alternative hypothesis.
estimate	the estimated mean or difference in means depending on whether it was a one-sample test or a two-sample test.
null.value	the specified hypothesized value of the mean or mean difference depending on whether it was a one-sample test or a two-sample test.
stderr	the standard error of the mean (difference), used as denominator in the t-statistic formula.
alternative	a character string describing the alternative hypothesis.
method	a character string indicating what type of t-test was performed.
data.name	a character string giving the name(s) of the data.

Note

Currently only performs a directional test

See Also

hr.comp2

Examples

## Not run: 
set.seed(12345)
age <- as.numeric(rnorm(100, 50, 10) >= 50)
size <- as.numeric(rexp(100, 1) > 1)
stime <- rexp(100)
cens <- runif(100, .5, 2)
sevent <- as.numeric(stime <= cens)
stime <- pmin(stime, cens)
coxm1 <- survival::coxph(survival::Surv(stime, sevent) ~ age)
coxm2 <- survival::coxph(survival::Surv(stime, sevent) ~ size)
hr.test(
  coxm1$coefficients,
  drop(sqrt(coxm1$var)),
  coxm2$coefficients,
  drop(sqrt(coxm2$var)),
  cor(age, size),
  100
)

## End(Not run)

the-mad-statter/washu documentation built on May 11, 2023, 7:24 p.m.