Short-term and long-term risk factors in gastric cancer

Abstract

While in chronic diseases, such as diabetes, mortality rates slowly increases with age, in oncological series mortality usually changes dramatically during the follow-up, often in an unpredictable pattern. For instance, in gastric cancer mortality peaks in the first two years of follow-up and declines thereafter. Also several risk factors, such as TNM stage, largely affect mortality in the first years after surgery, while afterward their effect tends to fade. Temporal trends in mortality were compared between a gastric cancer series and a cohort of type 2 diabetic patients. For this purpose, 937 patients, undergoing curative gastrectomy with D1/D2/D3 lymphadenectomy for gastric cancer in three GIRCG (Gruppo Italiano Ricerca Cancro Gastrico = Italian Research Group for Gastric Cancer) centers, were compared with 7148 type 2 diabetic patients from the Verona Diabetes Study. In the early/advanced gastric cancer series, mortality from recurrence peaked to 200 deaths per 1000 person-years 1 year after gastrectomy and then declined, becoming lower than 40 deaths per 1000 person-years after 5 years and lower than 20 deaths after 8 years. Mortality peak occurred earlier in more advanced T and N tiers. At variance, in the Verona diabetic cohort overall mortality slowly increased during a 10-year follow-up, with ageing of the type 2 diabetic patients. Seasonal oscillations were also recorded, mortality being higher during winter than during summer. Also the most important prognostic factors presented a different temporal pattern in the two diseases: while the prognostic significance of T and N stage markedly decrease over time, differences in survival among patients treated with diet, oral hypoglycemic drugs or insulin were consistent throughout the follow-up. Time variations in prognostic significance of main risk factors, their impact on survival analysis and possible solutions were evaluated in another GIRCG series of 568 patients with advanced gastric cancer, undergoing curative gastrectomy with D2/D3 lymphadenectomy. Survival curves in the two different histotypes (intestinal and mixed/diffuse) were superimposed in the first three years of follow-up and diverged thereafter. Likewise, survival curves as a function of site (fundus vs body/antrum) started to diverge after the first year. On the contrary, survival curves differed among age classes from the very beginning, due to different post-operative mortality, which increased from 0.5% in patients aged 65-74 years to 9.9% in patients aged 75-91 years; this discrepancy later disappeared. Accordingly, the proportional hazards assumption of the Cox model was violated, as regards age, site and histology. To cope with this problem, multivariable survival analysis was performed by separately considering either the first two years of follow-up or subsequent years. Histology and site were significant predictors only after two years, while T and N, although significant both in the short-term and in the long-term, became less important in the second part of follow-up. Increasing age was associated with higher mortality in the first two years, but not thereafter. Splitting survival time when performing survival analysis allows to distinguish between short-term and long-term risk factors. Alternative statistical solutions could be to exclude post-operative mortality, to introduce in the model time-dependent covariates or to stratify on variables violating proportionality assumption.

Key Words: Gastric cancer, Type 2 diabetes, Mortality, Short-term risk factors, Long-term risk factors, Survival analysis, Cox model, Proportional hazards assumption

Core tip: In gastric cancer patients mortality from recurrence peaked in the first two years after curative gastrectomy and then declined. The prognostic significance of risk factors was not stable, as the effect of T and N decreased with time, while the effect of histology became apparent after two years. Advancing age was associated with increased post-operative mortality. Hence, the proportional hazards assumption of the Cox model was violated. A possible solution was to split survival analysis in the first two years of follow-up and subsequent years. By comparison, overall mortality and prognostic significance of risk factors was rather stable in a cohort of type 2 diabetic patients.

INTRODUCTION

Declining trend in mortality after gastric cancer surgery

Despite the declining incidence, gastric cancer is currently the third leading cause of cancer death worldwide[1]. Prediction of prognosis after surgical resection for primary gastric adenocarcinoma is based mainly on the American Joint Committee on Cancer (AJCC) staging system[2] which takes into account depth of tumour invasion, pathological nodal status and the presence of distant metastasis as the only prognostic factors. Certainly, disease burden at the diagnosis is the most important predictor of cancer related-death, as up to 70% of deaths from cancer recurrence occur in the first two years after gastrectomy[3,4], mainly in patients with tumours infiltrating serosa, with extensive nodal involvement or systemic metastases.

However, after the first two years mortality rates remarkably decrease. In a Dutch-American study[4] 5-year disease-specific survival from the day of surgery presented a median improvement of 7.2%, 19.1%, and 31.6% respectively, in patients surviving 1, 2, or 3 years after surgery. Therefore, according to the US Gastric Cancer Collaborative[5] survival estimates following surgical resection of gastric adenocarcinoma can be considered dynamic: “the probability of survival increases with time already survived”.

Changing prognostic significance of risk factors during follow-up

The improvement in prognosis, recorded in gastric cancer patients surviving the first years of follow-up, is particularly large in patients with advanced stages[5,6].

In the Surveillance, Epidemiology, and End Results (SEER 17) database, five-year conditional survival (the probability to survive in the next five years) increased from 61% at baseline to 85% after 5 years in patients with localized tumors, from 34% to 79% in patients with metastases to regional lymph nodes, and from 2% to 64% in patients with distant metastases[6].

In the US Gastric Cancer Collaborative[5] the prognostic significance of T and N status was of paramount importance in the first three years after surgery: indeed three-year survival was respectively 77%, 59%, 43%, 27% in T1, T2, T3, T4, and 65%, 38% in N0, N1. However after 5 years of follow-up depth of tumor invasion and nodal status at baseline were no longer related to prognosis: indeed conditional survival in the next three years was about 78%-87%, irrespective of T or N tier[5].

In a Japanese study patients with early recurrence had more advanced T and N status, and a larger tumor size than patients with late recurrence[3]. However only depth of tumor invasion was independently associated with timing of recurrence.

Thus, in advanced gastric cancer patients, who survive over the first years but still present a notable risk of recurrence compared to early gastric cancers, factors other than TNM stage should be evaluated to correctly predict prognosis.

Problems in applying survival analysis to cancer series

As stated above, the trend in mortality changes over time and is often unpredictable, therefore analyzing mortality in the oncological field is particular difficult. Parametric models are not feasible, such as the Weibull model which assumes a monotonous increase or decrease in mortality rate over time[7]. Hence, the Cox model, which does not make a priori assumption on the temporal pattern of mortality, is usually the only choice for cancer patients. However, also the Cox model is often difficult to apply to patients with gastric cancer, due to violation of the proportional hazards assumption.

On the contrary, survival analysis encounters less problems in chronic degenerative diseases, such as diabetes, where mortality usually increases during follow-up due to ageing of the population under study.

These topics were investigated in two homogeneous GIRCG series (Gruppo Italiano di Ricerca sul Cancro Gastrico; Italian Research Group for Gastric Cancer), undergoing curative gastrectomy for gastric cancer: (1) temporal trends in mortality, either in whole series or as a function of disease stage, were inspected in a cohort of patients from Padua/Siena/Verona with early or advanced gastric cancer, undergoing curative gastrectomy with limited/extended/superextended (D1/D2/D3) lymphadenectomy[8]. Temporal trend in this gastric cancer series was compared with trend recorded in a cohort of diabetic patients, enrolled in the frame of the Verona Diabetes Study[9]; and (2) problems in survival analysis and possible solutions were investigated in a cohort of patients from Siena/Verona with advanced gastric cancer, operated on with curative gastrectomy and D2/D3 lymphadenectomy[10].

The lack of pre-operative chemotherapy and the adequate staging in both gastric cancer cohorts provided a unique opportunity to investigate both the hazard function over time and the temporal pattern of prognostic significance of main risk factors.

TREND OF MORTALITY DURING FOLLOW-UP IN A SERIES WITH EARLY/ADVANCED GASTRIC CANCER

Between 1988 and 2002, 1086 patients underwent gastrectomy for gastric cancer in three GIRCG centers: Padua (n = 270), Siena (n = 371) and Verona (n = 445)[8]. After excluding 1 neuroendocrine tumor, 1 T0 tumor, 7 multicentric tumors, 40 patients dying from postoperative causes, 100 macroscopic or microscopic residual tumors, 937 patients were considered for the study.

As regards T and N status, coded according to the UICC/AJCC TNM 6^th edition (1997), 261 patients were classified as T1, 359 as T2, 285 as T3, 32 as T4, while N status was N0, N1, N2, N3 respectively in 404, 299, 142, 92 patients.

Median follow-up in surviving patients was 106 mo (range: 36-216 mo): 3 patients were lost to follow-up. 515 patients died, 432 from tumor recurrence, 83 from other causes.

An Epanechnikov kernel function was used to calculate the weighted kernel-density estimate required to produce a smoothed hazard-function over time[11].

In this GIRCG series mortality from recurrence peaked to 200 deaths per 1000 person-years after 1 year of follow-up and then declined, becoming lower than 40 deaths per 1000 person-years after 5 years and lower than 20 deaths per person-years after 8 years (Figure 1A). This pattern could be inferred also by a shrewd inspection of the corresponding survival curve, which declined very steeply between 6 and 24 mo of follow-up, becoming rather flat thereafter (Figure 1B).

Figure 1 Temporal trend of recurrence-related mortality (A) and corresponding cancer-related survival (B) in a GIRCG series of 937 gastric cancer patients. To plot the estimated hazard function, a kernel smooth was used with a bandwidth of 2.5 mo. Survival curve was estimated by the Kaplan-Meier method.

As shown in Figure 2, mortality from recurrence varied widely as a function of T and N status in the first two years of follow-up, and tended to converge thereafter. Interestingly the mortality peak occurred earlier in more advanced T and N tiers.

Figure 2 Temporal trend in mortality from recurrence, as a function of T (1997) (A) and N (1997) (B) status in a GIRCG series of 937 patients. To plot the estimated hazard function, a kernel smooth was used with a bandwidth of 3 mo.

Comparison with mortality pattern in the Verona diabetes study

The mortality pattern in the gastric cancer series was compared with the pattern observed in 7148 type 2 diabetic patients from the Verona Diabetes Study. Briefly the Verona Diabetes Study[9] is a population-based monocentric survey, which identified 7148 type 2 diabetic patients in the Verona Social Health Unit on the 31^th December 1986, using three different sources of ascertainment (Diabetes Clinic, Family Physicians and Drug Prescription Database) with an estimated 75% completeness of ascertainment[12]. At baseline, 860 patients (12.1%) were treated with diet, 5821 (81.8%) with oral hypoglycemic drugs, 437 (6.1%) with insulin. From 1987 to 1996 2896 patients died, the overall follow-up amounting to 56410 person-years[13]. It can be appreciated that mortality slowly increased during the 10-year follow-up, with ageing of the type 2 diabetic cohort. Seasonal oscillations were also recorded: mortality rates was higher by about 10 deaths/1000 person-years during winter with respect to summer (Figure 3A). Interestingly, these oscillations were too small to be detected in the corresponding survival curve.

Figure 3 Temporal trend of all-cause mortality (A) and corresponding overall survival (B) in the 7148 type 2 diabetic patients from the Verona Diabetes Study. To plot the estimated hazard function, a kernel smooth was used with a bandwidth of 1.5 mo. Survival curve was estimated by the Kaplan-Meier method.

When hazard function was analyzed as a function of treatment, an important marker of disease progression, seasonal oscillations in overall mortality could be clearly detected in all groups but the insulin-treated group, where a very large bandwidth had to be used to smooth the estimated hazard function, due to the small size and the high mortality. Mortality rates as a function of treatment remained separated throughout the follow-up (Figure 4).

Figure 4 Temporal trend of all-cause mortality in the 7148 type 2 diabetic patients from the Verona Diabetes Study as a function of treatment, a proxy of disease progression. To plot the estimated hazard function, a kernel smooth was used with a bandwidth of 1.5 mo for the groups treated with diet and oral hypoglycemic drugs, and 4 mo for the group treated with insulin.

Of course, the comparison between the two series has some limitations. First of all, gastric cancer series represents an incident cohort, while the type 2 diabetic patients from the Verona Diabetes Study represents a prevalent cohort. Second, while mortality from gastric cancer is more easy to assess and usually corrected reported in death certificates[14], mortality from diabetes is often a matter of debate, as diabetes is an important risk factor for cardiovascular diseases, and even for certain tumors[13]. As a consequence, mortality from diabetes is often under-reported in death certificates[14]. For these reasons, we chose to analyze mortality from recurrence in gastric cancer patients and overall mortality in type 2 diabetic patients.

With all these limitations in mind, a clear difference was detected between mortality pattern in patients with gastric cancer, a disease with a very high lethality but which nevertheless can be cured, and mortality in patients with diabetes, a disease with much lower lethality but which cannot be cured. Overall mortality pattern was more stable in diabetic patients, as well as the prognostic significance of risk factors over time.

ANALYSIS OF A GIRCG SERIES WITH ADVANCED GASTRIC CANCER, UNDERGOING AT LEAST EXTENDED LYMPHADENECTOMY

Patients with advanced gastric cancer

The series comprised 791 patients undergoing gastrectomy with at least D2 lymphadenectomy for advanced gastric cancer in Verona or Siena[10]. Of these 445 patients were operated in Verona from January 1992 to May 2011 and 346 in Siena from January 1994 to June 2011. After excluding 19 subjects with Bormann IV tumour, 2 subjects with neuroendocrine tumours, 75 subjects with early gastric cancer, 127 subjects with non-curative resections, 568 subjects (312 from Verona, 256 from Siena) were left for the analysis.

Tumour invasion (pT) and lymph node status (pN) followed the UICC pathological tumour node metastasis (pTNM) criteria, 7^th edition.

Three hundred and fifty-one patients (61.8%) were male, and mean age ± SD was 66.4 ± 11.7 years, ranging from 24 to 92 years. 274 patients underwent D2 lymphadenectomy and 294 D3 lymphadenectomy, and nearly all (539/565 = 95.4%) could be adequately staged, as they had at least 15 lymph nodes retrieved. None of the patients received preoperative chemotherapy.

Patients were followed-up till December 2011. Median follow-up in surviving patients was 89 mo (range, 7-240 mo).

Statistical methods

In survival analysis both post-operative deaths and death from recurrence were considered as terminal events. Deaths caused by surgical and non-surgical complications were coded as post-operative deaths, and occurred in the first 38 d after surgery. Survival analysis was repeated by considering only death from recurrence.

Both the log-rank test and the Wilcoxon (-Breslow-Gehan) test were used to evaluate significance of differences among survival curves. The latter test places more weight to the events occurring at the beginning of follow-up, when there are more subjects at risk. The Wilcoxon test is preferred to the log-rank test when the hazard functions do not vary in a proportional way over time[7].

To test whether prognostic significance of the main risk factors changed over time, the proportional-hazards assumption of the Cox model was tested on the basis of Schoenfeld residuals. In addition, the proportionality assumption was checked by graphic methods: it was verified whether -ln[-ln[survival]] curves for each category of risk factors were parallel, when plotted vs ln[analysis time].

Univariable survival analysis

During the follow-up 315 (55.5%) patients died, 18 (3.2%) from post-operative causes, 242 (42.6%) from cancer recurrence and 53 (9.3%) from other causes. Cause of death was missing in 2 patients and timing of death in 3 patients.

One hundred and sixty-five deaths from recurrence (68.75%) occurred in the first two years of follow-up, 220 (91.7%) in the first four years, and 234 (97.5%) in the first six years. The highest hazard of death from recurrence was observed in the first 2 years of follow-up. No death from recurrence was observed after 12 years.

Survival curves were spread as a function of age in the first years, but not afterwards (Figure 5A). This was mainly due to different post-operative mortality, which was low in patients aged 25-64 years (1.3% = 3/231) or 65-74 years (0.5% = 1/196) but substantial in patients aged 75-92 years (9.9% = 14/141) (P < 0.001). Indeed when post-operative mortality was not considered, survival curves were nearly superimposed (Figure 5B).

Figure 5 Survival curves as a function of age (25-64, 65-74, 75-91 years), estimated by Kaplan-Meier method as on a GIRCG series of 568 patients undergoing R0 gastrectomy for gastric cancer. Post-operative mortality is taken into account in (A) but not in (B).

At variance survival curves of the two different histotypes (intestinal and mixed/diffuse) were superimposed in the first three years, but substantially diverged thereafter. Likewise, survival curves of the two different sites (fundus vs body/antrum) diverged after the first year (Figure 6).

Figure 6 Survival curves, estimated by Kaplan-Meier method on a GIRCG series of 568 patients undergoing curative gastrectomy for gastric cancer. The A, B, C and D display respectively survival curves as a function of site (fundus vs antrum/body), Lauren histotype (intestinal vs mixed/diffuse), depth of tumor invasion and nodal status.

Hence, survival curves as a function of age were significantly different according to the Wilcoxon test, which places more weight to first part of follow-up, but not according to the log-rank test. On the contrary, differences among survival curves of different histotypes were nearly significant when evaluated by the log-rank test, but not when reassessed by the Wilcoxon test.

At variance, survival curves of different T and N tiers were well separated throughout the follow-up (Figure 6).

Multivariable survival analysis; checking proportional hazards assumptions

First of all, the proportional hazards assumption was checked by graphic methods: it was verified whether -ln[-ln(survival)] curves for each category of risk factors were parallel, when plotted vs ln(analysis time). It can be appreciated that the proportional-hazards assumption was violated for age, site, histology as the curves were not parallel (Figure 7). As regards age, the curves were initially well separated and converged at the end of the follow-up. The reverse pattern was observed when considering either site or histology: the curves, initially superimposed, diverged at the end of follow-up. Parallel curves were observed as regards T and N status, except for the very beginning of the follow-up.

Figure 7 Log-log plots for age (A), site (B) and histology (C).

Proportional hazards assumption was also verified by applying a Cox regression model and then the test based on Schoenfeld residuals, which turned out to be significant both as a whole (P = 0.001) and as regards centre (P = 0.027), age (P = 0.005 for age 75-91 years with respect to age < 65 years), site (P = 0.007) and histology (P = 0.011). Hence, the proportional hazards assumptions of the Cox model were violated, as regards centre, age, site and histology. Hence the Cox model, which assumes the hazard ratio to be constant over time, was not feasible over the whole follow-up.

Interestingly, when considering as terminal event only death from recurrence but not post-operative death, the test based on Schoenfeld residuals was still fully significant as regards centre (P = 0.014), site (P = 0.023) and histology (P = 0.025), but not as regards age (P = 0.060 for age 75-91 years with respect to age < 65 years).

Multivariable survival analysis in the first two years of follow-up and thereafter

To distinguish between short-term and long-term risk factors, multivariable survival analysis was performed by considering either only the first two years of follow-up or only the subsequent years (Table 1). If one considers only the first 24 mo, site and histology were not significant predictors of mortality, while T and N exerted a strong effect. Histology and site were significant predictors only after two years, while T and N, although significant also in the long-term, lost a great deal of their prognostic significance. Increasing age was associated with higher mortality in the first two years, but not thereafter. Female sex was protecting both in the first two years and thereafter, although significant only in the former period.

Table 1 Prognostic significance of different risk factors in the first two years of follow-up and thereafter.

	First 2 years of follow-up		Subsequent years
	HR (95%CI)	P value	HR (95%CI)	P value
Centre (Siena vs Verona)	0.62 (0.44-0.88)	0.006	1.82 (1.06-3.14)	0.029
Sex (women vs men)	0.69 (0.49-0.97)	0.028	0.79 (0.48-1.30)	0.345
Age		< 0.001		0.586
65-74 yr vs < 65 yr	1.67 (1.14-2.44)		1.18 (0.68-2.04)
≥ 75 yr vs < 65 yr	2.28 (1.50-3.46)		0.81 (0.38-1.72)
Site (body/antrum vs fundus)	0.84 (0.60-1.18)	0.318	0.30 (0.18-0.51)	< 0.001
Histology (mix/diff vs intest)	0.80 (0.57-1.12)	0.198	1.91 (1.09-3.36)	0.023
T stage		< 0.001		0.043
T3 vs T2	1.69 (0.82-3.50)		1.25 (0.59-2.68)
T4a vs T2	3.13 (1.60-6.13)		2.21 (1.10-4.44)
T4b vs T2	5.14 (2.34-11.27)		3.26 (0.87-12.17)
N stage		< 0.001		0.007
N1 vs N0	1.49 (0.70-3.14)		0.92 (0.35-2.40)
N2 vs N0	3.47 (1.94-6.22)		2.55 (1.29-5.03)
N3 vs N0	5.28 (3.09-9.04)		2.36 (1.22-4.57)
Lymphadenectomy: D3 vs D2	0.94 (0.66-1.33)	0.717	0.92 (0.54-1.57)	0.753

P values were computed by the LR test. Hazard ratios (HR) were computed by a Cox regression model, controlling for all the other variables in the model.

DISCUSSION

The main results of the present study are: (1) mortality from recurrence peaked 1 year after curative gastrectomy in gastric cancer patients and then declined. Mortality peak occurred earlier in more advanced T and N tiers. At variance, in the Verona diabetic cohort overall mortality slowly increased during a 10-year follow-up, with ageing of the type 2 diabetic patients. Seasonal oscillations were also recorded, mortality being higher during winter; (2) depth of tumor invasion and nodal metastases gave a baseline snapshot of tumor progression at the time of diagnosis or tumor excision. In gastric cancer they are the most important predictors of mortality throughout the follow-up, but their prognostic significance tends to fade with time. On the contrary, biological intrinsic tumor characteristics, such as Lauren histotype, exerts a delayed effect on survival. The effect of age on mortality is even more limited in time, being mainly restricted to the post-operative period; and (3) accordingly, in gastric cancer series the proportional hazards assumption of the Cox model was violated, as regards age, site and histology. To cope with this problem, multivariable survival analysis was separately performed by considering only the first two years of follow-up or only the subsequent years.

Clinical implications

Hazard function over time can give additional information with respect to survival function. For instance, it allows to detect seasonal oscillations in overall mortality in type 2 diabetic patients. An excess mortality during winter has been recorded also in the general population, especially in elderly men[15]. Recently a GIRCG study recommended to focus oncological follow-up after radical surgery for gastric cancer in the first 3 years, although only 3.2% of patients with recurrence could be treated with potentially radical intent[16]. The pattern of mortality over time seems to support this statement, as mortality from recurrence peaked at 200 deaths per 1000 person-years after 1 year of follow-up and decreased below 90 deaths after 3 years.

Depth of tumor invasion and nodal metastases give a baseline snapshot of tumor progression at the time of diagnosis or tumor excision. They are the most important predictors of mortality throughout the follow-up, but their prognostic significance tends to fade with time. The effect of age on mortality is even more limited in time, being mainly restricted to the post-operative period.

On the contrary, biological intrinsic tumor characteristics, such as Lauren histotype, can affect prognosis in the short term by determining more advanced tumor stage, but they also display a delayed effect on survival. Indeed as regards tumor histology, in the first years of follow-up its prognostic significance is usually masked by N status, as Lauren diffuse cancers are more prone to give lymph node metastases. For instance, in the Verona series Lauren histology was an independent prognostic factor when N classification was based on site (TNM 1987) but not when based on number of positive nodes (TNM 1997)[17].

A similar pattern was recorded in adenocarcinoma arising from the ampulla of Vater, where “depth of invasion was the most important prognostic variable in the first 2 years of follow-up”, while allelic losses (LOH) at chromosome 17p, containing the TP53 gene, became the most important variable thereafter[18] (Figure 8). Indeed, “when considering only the first 20 mo of follow-up, statistical significance was retained by T stage (P < 0.001), but not by chromosome 17p LOH (P = 0.060); the reverse pattern was observed when excluding the first 20 mo of follow-up (P = 0.542 for T stage and P < 0.001 for chromosome 17p LOH)”[18].

Figure 8 Kaplan-Meier survival curves, according to both T-stage and chromosome 17p allelic status in ampullary adenocarcinoma. Modified from Iacono et al[18]. 17p-: Chromosome 17p allelic loss; 17p+: Chromosome 17p allelic retention.

Statistical implications

In the present series of gastric cancer patients the proportional hazards assumption of the Cox model was violated, as regards centre, age, site and histology.

Several solutions are available[7] but none seems fully satisfactory: (1) the Cox model could be stratified by those variables which do not satisfy the proportional hazards assumption. As a drawback, these variables are virtually excluded from survival analysis; indeed their impact on survival cannot be estimated nor significance of interactions with other variables. In the present series, it would be practical to stratify by centre, but not by age, site or histology, as the latter variables could be useful to tailor treatment for patient subgroups; (2) post-operative mortality could be excluded from the analysis. This would be particularly suited for age, as post-operative mortality preferentially affects older people, and this hampers the proportional hazards assumption. However, post-operative mortality is necessary to evaluate different treatments, such as survival after different surgical procedures[19,20] or chemoradiotherapy[21], or the quality of surgical procedures[8]; and (3) time-dependent covariates could be used, but they are somewhat difficult to interpret. Alternatively follow-up can be split into an early and a late period, and two separate Cox models can be employed for each period.

Survival analysis is simpler in patients affected by chronic degenerative diseases, such as diabetes, where the trend in mortality is rather stable over time. Hence parametric models could be suited for these patients, such as the Weibull model, assuming a monotonous increase or decrease in mortality rate over time[7]. Indeed in the Verona Diabetes Study observed deaths in different time intervals were not significantly different from expected deaths according to the Weibull model[22].

CONCLUSION

It can be concluded that in gastrointestinal cancer biological characteristics, such as histology or genetic features, exerts a delayed effect on prognosis. On the contrary, the actual spread of the tumor has an immediate effect on survival.

According to the Authors’ opinion methodological research in the surgical field is necessary in order to develop new statistical methods to compare groups with non-proportional hazards.