Logrank test

From Wikipedia, the free encyclopedia

In statistics, the log-rank test is a hypothesis test to compare the survival distributions of two samples. It is a nonparametric test and appropriate to use when the data are right skewed and censored (technically, the censoring must be non-informative). It is widely used in clinical trials to establish the efficacy of a new treatment in comparison with a control treatment when the measurement is the time to event (such as the time from initial treatment to a heart attack). The test is sometimes called the Mantel–Cox test, named after Nathan Mantel and David Cox. The log-rank test can also be viewed as a time-stratified Cochran–Mantel–Haenszel test.

The test was first proposed by Nathan Mantel and was named the log-rank test by Richard and Julian Peto.^[1]^[2]^[3]

Definition

The log-rank test statistic compares estimates of the hazard functions of the two groups at each observed event time. It is constructed by computing the observed and expected number of events in one of the groups at each observed event time and then adding these to obtain an overall summary across all time points where there is an event.

Let j = 1, ..., J be the distinct times of observed events in either group. For each time $j$ , let $N_{{1j}}$ and $N_{{2j}}$ be the number of subjects "at risk" (have not yet had an event or been censored) at the start of period $j$ in the two groups (often treatment vs. control), respectively. Let $N_{j}=N_{{1j}}+N_{{2j}}$ . Let $O_{{1j}}$ and $O_{{2j}}$ be the observed number of events in the groups respectively at time $j$ , and define $O_{j}=O_{{1j}}+O_{{2j}}$ .

Given that $O_{j}$ events happened across both groups at time $j$ , under the null hypothesis (of the two groups having identical survival and hazard functions) $O_{{1j}}$ has the hypergeometric distribution with parameters $N_{j}$ , $N_{{1j}}$ , and $O_{j}$ . This distribution has expected value $E_{{1j}}={\frac {O_{j}}{N_{j}}}N_{{1j}}$ and variance $V_{j}={\frac {O_{j}(N_{{1j}}/N_{j})(1-N_{{1j}}/N_{j})(N_{j}-O_{j})}{N_{j}-1}}$ .

The log-rank statistic compares each $O_{{1j}}$ to its expectation $E_{{1j}}$ under the null hypothesis and is defined as

$Z={\frac {\sum _{{j=1}}^{J}(O_{{1j}}-E_{{1j}})}{{\sqrt {\sum _{{j=1}}^{J}V_{j}}}}}.$

Asymptotic distribution

If the two groups have the same survival function, the log-rank statistic is approximately standard normal. A one-sided level $\alpha$ test will reject the null hypothesis if $Z>z_{\alpha }$ where $z_{\alpha }$ is the upper $\alpha$ quantile of the standard normal distribution. If the hazard ratio is $\lambda$ , there are $n$ total subjects, $d$ is the probability a subject in either group will eventually have an event (so that $nd$ is the expected number of events at the time of the analysis), and the proportion of subjects randomized to each group is 50%, then the log-rank statistic is approximately normal with mean $(\log {\lambda })\,{\sqrt {{\frac {n\,d}{4}}}}$ and variance 1.^[4] For a one-sided level $\alpha$ test with power $1-\beta$ , the sample size required is $n={\frac {4\,(z_{\alpha }+z_{\beta })^{2}}{d\log ^{2}{\lambda }}}$ where $z_{\alpha }$ and $z_{\beta }$ are the quantiles of the standard normal distribution.

Joint distribution

Suppose $Z_{1}$ and $Z_{2}$ are the log-rank statistics at two different time points in the same study ( $Z_{1}$ earlier). Again, assume the hazard functions in the two groups are proportional with hazard ratio $\lambda$ and $d_{1}$ and $d_{2}$ are the probabilities that a subject will have an event at the two time points where $d_{1}\leq d_{2}$ . $Z_{1}$ and $Z_{2}$ are approximately bivariate normal with means $\log {\lambda }\,{\sqrt {{\frac {n\,d_{1}}{4}}}}$ and $\log {\lambda }\,{\sqrt {{\frac {n\,d_{2}}{4}}}}$ and correlation ${\sqrt {{\frac {d_{1}}{d_{2}}}}}$ . Calculations involving the joint distribution are needed to correctly maintain the error rate when the data are examined multiple times within a study by a Data Monitoring Committee.

Relationship to other statistics

The log-rank test has been shown to be too permissive a test, allowing significant results for survivorship prediction models that have low accuracy. The F* test was developed in response to these observations and has been shown to be more critical and to track accuracy of the prediction models with higher fidelity.^[5]

The log-rank statistic can be derived as the score test for the Cox proportional hazards model comparing two groups. It is therefore asymptotically equivalent to the likelihood ratio test statistic based from that model.

The log-rank statistic is asymptotically equivalent to the likelihood ratio test statistic for any family of distributions with proportional hazard alternative. For example, if the data from the two samples have exponential distributions.

If $Z$ is the log-rank statistic, $D$ is the number of events observed, and ${\hat {\lambda }}$ is the estimate of the hazard ratio, then $\log {{\hat {\lambda }}}\approx Z\,{\sqrt {4/D}}$ . This relationship is useful when two of the quantities are known (e.g. from a published article), but the third one is needed.

The log-rank statistic can be used when observations are censored. If censored observations are not present in the data then the Wilcoxon rank sum test is appropriate.

The log-rank statistic gives all calculations the same weight, regardless of the time at which an event occurs. The Peto logrank statistic gives more weight to earlier events when there are a large number of observations.

References

↑ Mantel, Nathan (1966). "Evaluation of survival data and two new rank order statistics arising in its consideration.". Cancer Chemotherapy Reports 50 (3): 163–70. PMID 5910392.
↑ Peto, Richard; Peto, Julian (1972). "Asymptotically Efficient Rank Invariant Test Procedures". Journal of the Royal Statistical Society, Series A (Blackwell Publishing) 135 (2): 185–207. doi:10.2307/2344317. JSTOR 2344317. Cite uses deprecated parameters (help)
↑ Harrington, David (2005). "Linear Rank Tests in Survival Analysis". Encyclopedia of Biostatistics. Wiley Interscience. doi:10.1002/0470011815.b2a11047.
↑ Schoenfeld, D (1981). "The asymptotic properties of nonparametric tests for comparing survival distributions". Biometrika 68: 316–319. JSTOR 2335833.
↑ Berty, H. P.; Shi, H.; Lyons-Weiler, J. (2010). "Determining the statistical significance of survivorship prediction models". J Eval Clin Pract 16 (1): 155–165. doi:10.1111/j.1365-2753.2009.01199.x. PMID 20367827.

External links

Bland, J. M.; Altman, D. G. (2004). "The logrank test". BMJ 328 (7447): 1073. doi:10.1136/bmj.328.7447.1073. PMC 403858. PMID 15117797.

Statistics

Descriptive statistics

Continuous data

Location	Mean (Arithmetic, Geometric, Harmonic) Median Mode

Dispersion	Range Standard deviation Coefficient of variation Percentile Interquartile range

Shape	Variance Skewness Kurtosis Moments L-moments

Count data

Index of dispersion

Summary tables

Dependence

Statistical graphics

Data collection

Designing studies	Effect size Standard error Statistical power Sample size determination

Survey methodology	Sampling Stratified sampling Cluster sampling Opinion poll Questionnaire

Controlled experiment	Design of experiments Randomized experiment Random assignment Replication Blocking Factorial experiment Optimal design

Uncontrolled studies	Natural experiment Quasi-experiment Observational study

Statistical inference

Statistical theory	Sampling distribution Order statistic Scan statistic Record value Sufficiency Completeness Exponential family Permutation test (Randomization test) Empirical distribution Bootstrap U statistic Efficiency Asymptotics Robustness

Frequentist inference	Unbiased estimator (Mean unbiased minimum variance, Median unbiased) Biased estimators (Maximum likelihood, Method of moments, Minimum distance, Density estimation) Confidence interval Testing hypotheses Power Parametric tests (Likelihood-ratio, Wald, Score)

Specific tests	Z (normal) Student's t-test F Goodness of fit (Chi-squared, G, Sample source, sample normality, Skewness & kurtosis Normality, Model comparison, Model quality) Signed-rank (1-sample, 2-sample, 1-way anova) Shapiro–Wilk Kolmogorov–Smirnov

Bayesian inference	Bayesian probability Prior Posterior Credible interval Bayes factor Bayesian estimator Maximum posterior estimator

Correlation and regression analysis

Correlation	Pearson product–moment correlation Partial correlation Confounding variable Coefficient of determination

Regression analysis	Errors and residuals Regression model validation Mixed effects models Simultaneous equations models MARS

Linear regression	Simple linear regression Ordinary least squares General linear model Bayesian regression

Non-standard predictors	Nonlinear regression Nonparametric Semiparametric Isotonic Robust Heteroscedasticity Homoscedasticity

Generalized linear model	Exponential families Logistic (Bernoulli) Binomial Poisson

Partition of variance	Analysis of variance (ANOVA) Analysis of covariance Multivariate ANOVA Degrees of freedom

Categorical, multivariate, time-series, or survival analysis

Categorical data

Multivariate statistics

Time series analysis

General	Decomposition Trend Stationarity Seasonal adjustment Exponential smoothing Cointegration

Specific tests	Granger causality Q-Statistic Durbin–Watson

Time domain	ACF PACF XCF ARMA model ARIMA model ARCH Vector autoregression

Frequency domain	Spectral density estimation Fourier analysis

Survival analysis

Applications

Biostatistics	Bioinformatics Clinical trials & studies Epidemiology Medical statistics

Engineering statistics	Chemometrics Methods engineering Probabilistic design Process & Quality control Reliability System identification

Social statistics	Actuarial science Census Crime statistics Demography Econometrics National accounts Official statistics Population Psychometrics

Spatial statistics	Cartography Environmental statistics Geographic information system Geostatistics Kriging

Category
Portal
Outline
Index

This article is issued from Wikipedia. The text is available under the Creative Commons Attribution/Share Alike; additional terms may apply for the media files.