/W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd">
In and Lee: Survival analysis: Part II — applied clinical data analysis

DOI: https://doi.org/10.4097/kja.19183.e1
Survival analysis: Part II — applied clinical data analysis

Junyong In^{}, Dong Kyu Lee^{}

Copyright © The Korean Society of Anesthesiologists, 2023
Korean J Anesthesiol 2019 October 72(5): 441457
This corrects the article “Survival analysis: Part II — applied clinical data analysis” on pages 453456.
The article by In and Lee entitled, “Survival analysis: Part II — applied clinical data analysis”, contained a wrong explanation in “R code for timedependent coefficient Cox regression model: step function.” On page 453, the authors described the combined results from the timedependent coefficient Cox regression model. The estimated coefficients corresponding to the segmented time periods cannot be combined; they should be interpreted as separated (Table 10). The R code described as “#Combined results” also is not a resultscombining process. The function “survfit” computes estimates of a survival curve using the established survival function. That is, the following R code produces the estimates of the timedependent coefficient Cox regression model, which is fitted before and defined in the object “fit.tdc”
R code (P454) ‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒
cfit.tdc < survfit(fit.tdc, newdata = combine.tdc, id = parameter)
‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒‒
This R code can calculate the estimates of specific populations by providing a userdefined data frame using the argument “newdata = ”. For example, we can get the estimates of the fitted Cox regression model for Antiemetic = 0 and Inopioid = 0, and so on. The final example model of the timedependent coefficient Cox regression has two categorical independent variables, each having two levels, and four types of populations can be defined. The provided data frame on page 453 named “combine.tdc” has the wrong code in the argument “parameter.” We have to set the four types of population, “parameter = rep(0:3, each = 3)” is right instead of “parameter = rep(0:1, each = 6).” Consequently, the figures 10 and 11B are also wrong, we present the corrected graphs below, having four types of populations. These new graphs are created using the same R code provided in the original manuscript.
Figure 1. Corrected Figure 10. Timedependent coefficient Cox regression curves are plotted for four types of populations.
Figure 2. Corrected Figure 11B. Timedependent coefficient Cox regression curves are plotted for four types of populations. Compared to the original figure 11B, the confidence intervals are removed to enhance the discriminability.
The authors apologize for any inconvenience these mistakes may have caused.