 Methodology
 Open Access
 Published:
A reference relative timescale as an alternative to chronological age for cohorts with long followup
Emerging Themes in Epidemiologyvolume 12, Article number: 18 (2015)
Abstract
Background
Epidemiologists have debated the appropriate timescale for cohort survival studies; chronological age or timeonstudy being two such timescales. Importantly, assessment of risk factors may depend on the choice of timescale. Recently, chronological or attained age has gained support but a case can be made for a ‘reference relative timescale’ as an alternative which circumvents difficulties that arise with this and other scales. The reference relative time of an individual participant is the integral of a reference population hazard function between time of entry and time of exit of the individual. The objective here is to describe the reference relative timescale, illustrate its use, make comparison with attained age by simulation and explain its relationship to modern and traditional epidemiologic methods.
Results
A comparison was made between two models; a stratified Cox model with age as the timescale versus an unstratified Cox model using the reference relative timescale. The illustrative comparison used a UK cohort of cotton workers, with differing ages at entry to the study, with accrual over a time period and with long followup. Additionally, exponential and Weibull models were fitted since the reference relative timescale analysis need not be restricted to the Cox model. A simulation study showed that analysis using the reference relative timescale and analysis using chronological age had very similar power to detect a significant risk factor and both were equally unbiased. Further, the analysis using the reference relative timescale supported fullyparametric survival modelling and allowed percentile predictions and mortality curves to be constructed.
Conclusions
The reference relative timescale was a viable alternative to chronological age, led to simplification of the modelling process and possessed the defined features of a good timescale as defined in reliability theory. The reference relative timescale has several interpretations and provides a unifying concept that links contemporary approaches in survival and reliability analysis to the traditional epidemiologic methods of Poisson regression and standardised mortality ratios. The community of practitioners has not previously made this connection.
Background
In recent years, epidemiologists have debated which of the several possible timescales to use for survival analysis for longitudinal studies and it has been argued that chronological age as the timescale is preferable to the traditional timeonstudy [1, 2]. The utility of chronological age has been disputed because, it is argued, that without further covariate adjustment for age at entry to the cohort or without ‘lefttruncation’ the unadjusted age scale is inferior in performance to other models [3]. When adjustment is made for age at entry then, it has been suggested, model coefficients for risk factors of interest differ little between models using the chronological age scale and those using timeonstudy [3, 4]. However, others have reported substantial differences in the assessment of risk factors depending on the timescale used for the analysis [5–7]. The choice of correct timescale is important because the substantive findings of large cohort studies might be called into question if it were believed that an inappropriate choice leads to seriously biased estimates of hazard ratios for risk factors [8]. What is clear is that considerable care is required in the choice of timescale for any particular application; one example being when an environmental exposure is highly correlated with the calendar timescale [9]. The use of alternative timescales has been discussed widely in the context of reliability and performance and the definition of a good timescale has been proposed [10].
Although there is still no general consensus, chronological age as timescale has gained moderate acceptance for the analysis of survival data from cohort studies. With chronological age as the timescale, age is usually deemed to be truncated on the left at age of entry to the study. Cox regression is the modelling framework of choice to assess risk factors and the baseline hazard is a nonparametric function of age. This assumption is problematic for long running studies which accrue participants over a period of time since the hazard cannot be a function of age alone but must also depend on calendar time. For example, the mortality rate of a white male of 60 years in California in 1984 would not be the same as that of a white male of 60 years in California in 2004 since longevity would have improved. Using chronological age leaves unresolved the issue of calendar time when participants enter the cohort at different calendar times and this issue is not trivial for studies which span several decades. Cohorts and particularly occupational cohorts, often have these features, that participants join the cohort at different ages and at different calendar times and are followedup for long periods of time thereafter. One recommendation is to stratify on birth cohort in the Cox regression model using 5 or 10 yearly intervals but, for long running cohorts, this can result in very many strata [1, 3], although strata of varying lengths can be used to reduce the total number. Another possibility is to include age at entry as a covariate but modelling ageatentry may require a complex submodel to make adequate adjustment [5]. Stratification or covariate adjustment is not the only answer and a viable alternative, which has received little attention, is to take a relative survival approach to the analysis of cohorts with long followup. The aim of this paper is to publicise this alternative approach and to indicate how this addresses the issue of accommodating age at entry, birth cohort, differing entry times and changing longevity.
Relative survival has a long history in epidemiology [11–13]. Relative survival compares the mortality (or other event of interest) in a study population to that in an appropriate reference population. The standardised mortality ratio (SMR) is a primary traditional estimate of relative mortality and is defined as the number of events occurring in a study group relative to the number of events expected from the event rates observed in the reference population. The Poisson regression method for event counts, which uses expected number of events from the reference population in place of the person—years of observation, has been used widely in the past to model the association of risk factors with mortality or disease outcome [14–16]. Although not always recognised as such, the Poisson regression method is equivalent to modelling the SMR [17]. Poisson regression has lost popularity in recent years because it is now well recognised that modelling aggregated counts has less power to identify risk factors compared to modelling individual survival times using individual participant data (IPD). Although an extension of Poisson regression to model IPD has been described [18], Cox regression for individual survival times still comprises the regression method of first choice, since it makes no statistical distributional assumptions.
The objective of this communication is to compare, as an illustration, a standard analysis using chronological age as the timescale, lefttruncated at age of entry to the study, and with stratification for birth cohort in the Cox model with the analysis of individual survival times on a reference relative timescale. The intention is to demonstrate that this approach provides a viable alternative to the current standard analysis and is an approach which circumvents much of the controversy surrounding the choice of timescale. Also, it has appealing interpretations and is connected to the traditional methods of Poisson regression and SMR.
Methods
This section described more fully the two time scales being compared; the new reference relative timescale and the chronological age timescale. Following this the illustrative cohort data is introduced and the statistical analyses of the data are described. The section finishes with a methodological description of the simulation study which was used to establish existence of any bias and the relative power of the two timescales to detect a differential mortality risk between men and women of light smoking compared to nonsmoking.
The reference relative timescale
An individual measure of relative survival has been described previously in which, each participant’s actual survival time is transformed to a new scale which is the expected residual cumulative distribution function from a reference population [19, 20]. This measure lies in the range 0.0 to 1.0. For example, an individual participant‘s measure of 0.80 implies that this participant’s survival time is longer than 80 % of their peers in the reference population where peers are defined as people of the same gender, date of birth and other key demographics such as race and location. Modelling of the risk factors is then performed in which the response variable is the transformed survival time, Y, and the regression method is the Cox model [20]. Cox regression is valid even though the response measure, Y, is bounded above by 1.0 and does not have the fundamental attribute of time, that of being unbounded above. If the measure Y is further transformed to the measure Z where Z = −ln(1−Y) then Z has the property that the Cox model fitted to the measure Z gives identical regression estimates and fit to the Cox model fitted to the measure Y. The measure Z is also bounded below by 0.0 but Z is not bounded above and so has the attributes of a realistic timescale. The measure Z is the reference relative time.
More specifically, if an individual of gender g enters the study at age a_{0} at calendar time s_{0} and exits at age a_{1} at calendar time s_{1} then \({\text{Z}} = \int_{{s_{0} }}^{{s_{1} }} {\lambda \left( {g,a,s} \right)ds}\) where λ(g, a, s) is the hazard in the reference population for a person of gender g, at age a at calendar time s. The reference population may be defined by further demographics such as race and location, but the pertinent point is that the integral is taken over the hazard function for the individual participant’s peer group in the reference population. Thus the individual’s relative time is the cumulative hazard of the peer group in the reference population over the calendar time on study of the individual participant. The survival censoring indicator is the same as that for chronological age, which is whether the individual participant had an event at time s_{1} or was censored.
The transformation, Z, transforms both censored and uncensored survival times to a new relative timescale, called the reference relative timescale. Figure 1 shows a representation of how the reference relative timescale is obtained. The reference relative timescale has a number of interpretations, some with an intuitive appeal and some which, at first sight, may seem counterintuitive.

1.
The reference relative timescale is the expected number of events for the individual participant during their timeonstudy. A reference relative survival time of 1.0 implies that the individual reached their residual life expectancy from entry based on the reference population. A reference relative survival greater than 1.0 implies that the individual had exceeded their residual life expectancy and a value less than 1.0 that they had died or embarked before reaching residual life expectancy.

2.
The reference relative timescale is a weighted version of the real timeonstudy. Real time is ‘stretched’ when the hazard (force of mortality) in the reference population is high but time is ‘shrunk’ when the hazard is low. So a calendar time period such as 5 years lived by an older individual participant between 1990 and 1995, has a longer reference relative time value than 5 years between 1990 and 1995 lived through by a younger individual. Shrinking time is a form of accelerating time whereas stretching time, decelerates time.

3.
The reference relative timescale is the reciprocal of an individual participant’s own SMR (divided by 100 if the latter is expressed as a percentage). If the SMR is calculated for a study group of size one then the reference relative time would be the reciprocal of this value. Thus, an individual participant with an SMR of 50 % would have a reference relative survival time of 2.0. For this individual, the reference relative time on study was so long that the individual would have been ‘expected’ to have died twice, with an explanation in the following sense. On calendar time of entry the individual has an expected residual lifetime. If the individual achieves the end of this expected residual lifetime, they would then have a further (but much shorter) expected residual lifetime. If they achieve the end of this second residual lifetime then the reference relative survival time becomes 2.0. A loglinear model for the hazard on the reference relative time scale is therefore also a loglinear model for the individual participants’ SMRs.

4.
The reference relative timescale has a simple gametheoretic interpretation as the final score in a game. In each time period, the individual player joins battle against opponents that are intent on eliminating the individual. The number of opponents is proportional to the hazard in the reference population (force of opposition or mortality). If the individual eliminates all the opponents then the player banks points proportional to the number of opponents overcome. The cumulative number of banked points on exit from the game measures the individual’s total success during timeonplay in the game. Censoring applies to the banked points if on exit from the game the individual participant has not been eliminated. So, the reference relative timescale is the cumulative amount of hazard that the individual has confronted before experiencing an event and is a measure of game success.
Modelling on the reference relative timescale obtained from the transformation Z has two significant benefits. Firstly, as well as using the Cox regression, relative time can be modelled using the full range of parametric survivals models, such as exponential and Weibull models. Secondly, if risk factors and other variables are all categorical, then modelling reference relative time using an exponential distribution is equivalent to Poisson regression using observed and expected counts. If risk factors and other variables contain quantitative covariates then modelling reference relative times is an extension of Poisson regression to IPD. Thus modelling on the reference relative timescale can be understood as a generalisation of the traditional Poisson regression method to IPD and to more complex survival distributions. Reference relative survival, therefore, provides a unifying conceptual framework which encompasses both traditional and contemporary methods of analysis.
The reference relative timescale is an example of an alternative timescale or ‘operational failure timescale’ which has been discussed in the literature on reliability and performance of systems [10]. It meets the criteria for definition as an ‘ideal timescale’ which is defined by four features (1) scientific relevance (2) parsimonious and accurate description of variation in failure times under different conditions (3) a compact statistical distribution on the transformed timescale and (4) succinct and meaningful summarization of the effects of covariates of interest. An ‘ideal timescale’ arises because the hazard in the reference population provides a timevarying external continuous covariate during each individual participant’s timeonstudy.
The chronological age timescale
Currently, the conventional approach for survival analysis for cohort data with long followup, in which participants enter at different calendar times and where ages at entry are variable is to use the timescale of chronological age with left truncation for age at entry. With this approach adjustment for birth cohort is made by using a stratified Cox model with 5 or 10 year intervals [3, 5] or intervals of varying size as appropriate.
Illustrative cohort data
Analyses of a long running cohort of British cotton mill workers (1966–2007) have been previously described. The cohort was founded with the aim of understanding the longterm effects of exposure to environmental particulates on respiratory mortality and morbidity. The cohort has provided evidence that working with cotton reduced the risk of lung cancer death and that light smoking at baseline predicted higher mortality for women compared to men [21–23]. These results were obtained using the reference relative timescale as the mode of analysis and the primary inferential topic of interest here was to demonstrate that these results were robust to the choice of timescale.
Ages at entry to the cohort were in the range 15–81 years (quartiles, 32, 45 and 56 years) over the period 1966–1974 and participants’ year of birth was in the range 1885–1954. Smoking consumption was categorised at baseline as nonsmoker, light (1–14), medium (15–24), heavy (25+ cigarettes per day) smoker and former smoker. The cohort contained both men and women mostly working fulltime in all of the mill environments. Other variables measured at baseline included; the presence of the lung disease of cotton mill workers known as Byssinosis, the presence of cough and phlegm for at least 3 days per week for at least 3 months of the year, how long the worker had worked in the cotton industry and lung function expressed as forced expiratory volume in one second adjusted for age, gender and height (FEV_{1}% predicted) and the ratio of FEV_{1}% predicted to forced vital capacity (FVC). Participants alive at age 90 were censored. Table 1 illustrates the cohort data using 10 typical but hypothetical individual participants. Also, Table 1 shows the type of values produced by the transformation to reference relative survival time for the hypothetical individuals for cohort entry at one calendar time point. The participants with identifiers 1, 2 and 3 all had the same time on study but had different ages at entry. The first two participants of these three had small reference relative times because they entered at a young age, whereas the third participant entered at 40.8 years and their reference relative time on study was >1.057 showing that they had just attained their expected lifetime when they exited from the study. The participant with identifier 10 entered at age 60.7 years and their reference relative time on study was 1.756 which shows that they had well exceeded their expected residual lifetime from age at entry.
Ethics approval
Ethical approval was obtained from the University of Central Lancashire’s Faculty of Health Research Ethics Committee which accepted that the study had been granted exemption by the Department of Health’s National Information Governance Board from the need to obtain informed consent from individuals retrospectively to participate in the mortality study. In addition the Medical Research Information service at the NHS Information Centre granted permission for the study to receive vital status data.
Statistical analyses
To test the utility of the reference relative survival timescale, two regression analyses were compared using the cotton mill workers’ cohort data. The first analysis was the conventional analysis that used age as the timescale with left truncation in a stratified Cox model with 24 strata, representing 12 5year birth cohorts per gender as suggested in the literature [1, 3]. The second was the analysis using the reference relative timescale calculated using population mortality rates for England and Wales. The numbers of deaths by gender and by age in 5year bands were obtained for each year from 1966 to 2007 and the corresponding midyear population size estimates for England and Wales were also obtained [24]. The numbers of deaths were divided by the corresponding population size and then by 12 to obtain a monthly mortality hazard rate. The integral of the hazard rate was computed using a time interval of one month as an adequate approximation. Figure 1 gives an example of the estimated monthly mortality hazard for male aged 45 on 1 January 1965. The hazard is a step function because mortality was published in 5year bands. The effect of improving longevity due to better health and social care is discernible in the declining hazard over a 5year interval at the plateau on each step.
The analyses estimated the hazard ratio (HR) for light, medium, heavy and former smoker relative to nonsmoker together with HRs for other risk factors measured at baseline. The analyses also estimated the women to men relative risk ratios (RRR) for the four baseline categories of cigarette consumption; light, medium, heavy and former smoking. The relative risk ratio was defined by RRR = HR of female smokers to female nonsmokers/HR of male smokers to male nonsmokers where HR was the hazard ratio. The RRRs were the single degree of freedom components of the interaction between gender and consumption category and were estimated using indicator variables. Both analyses fitted the Cox model, but the Weibull and the exponential models were also fitted for the survival times on the reference relative timescale. The analyses were performed in the R programming language and environment [25].
Simulation study
In order to compare the bias and power of the reference relative timescale to that of chronological age at risk timescale a simulation study was carried out. Each simulation included all the participants from the illustrative cohort together with the participant’s age at entry, gender and date of entry to the study being kept fixed at their observed values. All the participant’s covariate and factor values were also kept fixed at their observed values. Each participant’s lifetime from date and age of entry was then simulated according to a model which calculated the participant’s hazard of dying in each calendar month following date of entry. For each month, in turn moving forward in time, a uniform random variable was generated. The first month, in the sequence, for which the random variable was less than the monthly hazard was selected as the simulated time of death. Once simulation of time of death was completed for all participants, a model was fitted by Cox regression using the reference relative timescale and by a stratified Cox regression using the age at risk timescale.
In the simulation, the monthly hazard for each participant was calculated from a combination of the model parameters and the England and Wales reference population monthly hazards. The model assumed that the hazard of dying was proportional to smoking status; nonsmoker, light smoker, medium smoker, heavy smoker or former smoker with relative parameters α_{0}, α_{1}, α_{2}, α_{3}, α_{4} respectively and where α_{0} = 1. These parameters applied to both males and females except in the case of light smoking females for whom the relative parameter was βα_{1}. In the simulation α_{0}, α_{1}, α_{2}, α_{3}, α_{4} were kept fixed at the values of 1.00, 1.24, 1.83, 2.08 and 0.93. Each simulation was for a fixed value of β in the range 1.0–1.4. The actual hazard used in each month was obtained from these model parameters and the numbers of participants in each gender and smoking category in the illustrative cohort. For males, if there were n_{0}, n_{1}, n_{2}, n_{3}, n_{4} in the categories; nonsmoker, light smoker, medium smoker, heavy smoker and former smoker and the monthly hazard in the reference population was h_{m} then the hazard for nonsmokers, h_{0} was computed by.
The hazards for light smokers, medium smokers, heavy smokers and former smokers were given by α_{1}h_{0}, α_{2}h_{0}, α_{3}h_{0} and α_{4}h_{0} respectively. For females the calculation of the hazard was analogous except that α_{1} was replaced by βα_{1}. This method of calculating the relative hazards had the effect that the average over a cohort with the observed mix of smoking grade would equal the reference population monthly hazard. This method generated realistic lifetimes with the appropriate relative hazards given by the model parameters. This allowed the two methods to be fairly compared in terms of their power to detect a value of β greater than 1.0 which would indicate a gender difference in the relative risk of light smoking and any bias in the estimation of β. It also allowed a check on the power of the null value of β = 1. For each value of β, 200 simulations were carried out and testing to estimate power used a 5 % significance level for the test.
The simulation also provided an opportunity to compare the methods when the reference population used to obtain the reference relative times was not an appropriate population for the cohort under study. This was achieved by increasing the monthly hazard of dying in the simulation by an additional multiplying factor which was an exponential function of age; exp (0.002 × age). This increased the hazard by about 4 % at age 20 and 13 % by age 60. This produced disproportionately shorter simulated lifetimes than in the reference population used to obtain the reference relative timescale.
Results
Illustrative cohort data
A comparison of the fitted Cox model using the reference relative timescale with the stratified Cox model using chronological age as the timescale is presented in Table 2. Gender main effect was fixed at 1.0 in the stratified model due to the choice of stratification which included genderspecific strata. Previously, the proportional hazards assumption was tested for the Cox model using the reference relative timescale by calculating the correlation between the Schoenfeld residuals and the transformed survival times and the assumption was shown to be satisfied [23]. Overall, there was a good measure of agreement between the two approaches in the parameter estimates and the 95 % confidence intervals for the risk factors of interest, suggesting that the two methods of analysis were comparable.
The fitted model coefficients for two alternative fully parametric survival models, the Weibull and the exponential, when fitted to the reference relative times, are shown in Table 3. The Weibull was previously identified as a wellfitting parametric model for the reference relative survival times in the cotton mill workers’ study [23]. It is clear from the parameter estimates in Table 3 that the RRR for light smoking, 1.35 for the Weibull model, significantly exceeded the ‘no difference’ value of 1.0 and this indicates that the relative risk of light smoking compared to never having smoked was predicted at onethird greater for women than for men. Furthermore, a lung function at 10 % below the reference standard for a person of a given gender, age and height had a parameter value 1.07, indicating a significant predicted 7 % increase in mortality hazard. This increase was after taking into account the effect of smoking by including smoking and gender in the model.
The fullyparametric models included a constant term and so the predicted percentiles of the survival distribution were computed. The percentiles were back transformed using the reference population hazards in order to build mortality curves as a function of age. This was especially easy to execute since a parametric survival model had been fitted and demonstrated the advantage of the reference relative survival time followed by parametric survival modelling. Figure 2 shows an example of this application; the effect of the higher risk to women of light smoking can be visualised by the narrowing of the gender gap compared to nonsmoking. Therefore the reference relative survival approach can satisfy the aim of obtaining survival curves as a function of age.
The exponential model shown in Table 3 is identical to the model that would have been fitted had Poisson regression for IPD data been used [17, 18] and shows that the reference relative survival timescale provides an alternative route for carrying out Poisson regression for IPD. The parameter values and their interpretation for the exponential model are very similar to those for the Weibull. Practitioners who are more comfortable with traditional methods can use the reference relative time scale and be confident that their results concur with Poisson regression and that the results will be similar to analysis using chronological age.
Simulation study
For the values of the parameter β = 1.20 and β = 1.40, the estimated power was very similar when the simulated lifetimes were analysed using chronological age at risk as timescale or the reference relative timescale (Table 4; appropriate reference population). For both methods, when β = 1.00, the estimated probability of rejecting the null hypothesis of log(β) = 0 was similar to the type 1 error rate of 5 %. In all cases the mean values of the estimates of β were close to the true values and both methods showed no evidence of bias in estimation.
When lifetimes were simulated somewhat shorter than in the reference population, the reference relative timescale method performed well and performed similarly to the chronological age at risk timescale in terms of power and bias (Table 4; inappropriate population). The simulations thus supported the assertion that the reference relative timescale was a viable alternative to using chronological age at risk.
Discussion
Transforming data, to achieve conformity to a parametric statistical distribution, is a cornerstone of much data analysis and yet surprisingly it is not often used for survival data from cohort studies. It has been used frequently in reliability and performance theory where there is motivation to determine suitable usage or exposure measures which transform real time to new scales. The term ‘operational time’ has been used for a timescale obtained by integrating the hazard function of an inhomogeneous process [26] and the expression ‘operational failure time’ has been used in reliability analysis to describe a transformed failure time. The transformation to the reference relative timescale described here is a transformation to an operational failure timescale which acts to remove a large component of the inhomogeneity in the observed survival times by using a reference population containing external knowledge of past temporal and crosssectional hazard rates as a measure of this inhomogeneity. It need not be assumed that the reference population hazard rates apply directly to the cohort participants. Indeed the same rates would not apply if the cohort were not healthy individuals but were a group with a disease diagnosis. Rather it would be enough if the hazards rates could be assumed proportional for the removal of homogeneity to succeed. The transformation to the reference relative timescale is a preprocessing of the data so that modelling can go ahead without concern for complex submodels for study group demographics, such as age, or large numbers of strata in a stratified model which has concerned data analysts [5]. The modelling effort can focus on the risk factors of interest and parsimonious models can be determined. If there was concern that the preprocessing had not been totally effective, then additional covariates such as age and age squared could be included in the regression model. This would provide a test of lackoffit since the coefficients for these terms should be negligible if the preprocessing has achieved its objective. The preprocessing also makes it more likely that a suitable parametric survival model can be identified and this may provide insight into the underlying stochastic mechanisms.
In contrast, it can be argued that a primary objective of epidemiology is to understand the complex pattern of risk over age and over calendar time and that the removal of all or a part of this pattern may hinder rather than help meet this objective. If age and calendar time are no longer in the regression model then the pattern with age and calendar time cannot be visualised in the model coefficients. However, comparative patterns with age and calendar time can be obtained using a back transformation of the reference population as has been shown here but this requires additional effort beyond interpreting the regression coefficients. Also, if age and calendar time are included in the regression model fitted after preprocessing then a good deal of care would be needed in their interpretation. It is clear that transformation to the reference relative timescale will not be appropriate in all applications and this is a limitation of the methodology.
A further limitation of the methodology as described here is that no account has been taken of the precision of estimation of the reference population hazard function. In the illustrative example actual mortality in England and Wales was likely to be reasonably well estimated but population size would have a greater degree of imprecision and this would hold true for many other potential reference populations. As described here, the methodology assumes that the reference population hazard is measured without error and, if this was in doubt, a sensitivity analysis would be needed to confirm study findings.
In the illustrative example given here, all covariates and factors were determined at baseline when participants entered the study. If covariates and factors varied during followup then the analysis would follow the usual method for handling timevarying covariates. The total calendar time interval on study for each individual would need to be subdivided into consecutive time windows during which the covariates and factors were assumed constant. Each portion before the last would be censored and the final time window would reflect whether the participant finally had the event of interest or was censored at exit from the study. Then each time window would be separately transformed to the reference relative time scale and the regression analysis proceed as usual.
To reiterate, in reliability theory a definition of a ‘good timescale’ has been proposed [10] with four defining features; (1) scientific relevance (2) parsimonious and accurate description of variation in failure times under different conditions (3) a compact statistical distribution on the transformed timescale and (4) succinct and meaningful summarization of the effects of covariates of interest. The reference relative timescale demonstrated these attributes in the illustration given here. The timescale had several relevant interpretations, values on the scale could be defined for the widely differing individual participants, the Weibull provides a compact distribution on the new timescale and effects of risk factors of interest were succinctly summarised by the model coefficients.
The reference relative timescale discussed here is a transformation of the measure of individual relative survival that has been suggested previously in the literature [19]. Subsequently, the concept of relative survival seems to have evolved and relative survival is often used now to mean only the estimation of the excess mortality due to a particular disease or condition within an additive hazards model framework [27–30]. Since the disease or condition is likely to increase mortality not only from the specific condition or disease but also from other causes, allcause mortality provides an easier quantity to measure than cause specific mortality. The reference population is used as a baseline from which to estimate the excess numbers of deaths. In this context, the models fitted are additive hazards models but the concept of relative mortality is equally applicable to proportional hazards models [10]. Where cause specific mortalities are available for the reference population, a reference relative timescale can be defined for any specific cause and so relative survival, in the context described here, can be applied to both cause specific and allcause mortality. Hence, there may be applications to the modelling of competing risks since a value on each of a number of reference relative timescales for different causes of death could be computed to provide a vector of multivariate survival data in reference relative time. Further, the metaanalysis of IPD is of increasing importance [31]. The use of the reference relative timescale may have benefit for metaanalysis of IPD since cohorts from different locations could each be transformed using different reference populations and then combined in a single regression analysis provided, of course, the same covariates were recorded in each location.
The reference relative time scale could be viewed as a composite timescale which merges an age timescale with a calendar time scale. Age–periodcohort (APC) models, which index an event count by age, calendar period and birth cohort, have an extensive literature [32] and have been of interest because of the recognised failure of identifiability of effects due to colinearity. The reference relative time scale may have applications in APC modelling since it reduces the time scales from three to two and so models might become identifiable.
The use of a reference population to gain insight into the risk factors affecting a study cohort begs the question of whether results from the study group can be generalised to the reference population. This would seem likely if it could be assumed that the study group is, in some sense, a sample from the reference population, but might the generalizability be conditional upon certain other assumptions? It would be reasonable to use a reference population to create a relative timescale so long as there was a sound belief that the hazards in the study group were proportional to those in the reference population. If this was the case, would generalisation to the reference population be valid also?
It may be worthwhile considering whether the definition of the reference relative timescale as the integral over a reference population hazard could be extended to include an additional term under the integral for the quality of life of the individual as a function of age. The transformation to the reference relative timescale stretches time when the population hazard is high and this effect is greater for older compared to younger individuals. A quality of life function could act as a penalty to reduce some of the stretching and would change conclusions regarding risk factors of interest when regression models are fitted on the reference relative timescale. It is clear, therefore, from these discussions that the reference relative timescale not only provides a viable alternative to modelling chronological age or real timeonstudy but is an interesting concept in its own right that points the way to several avenues of future research which others may wish to explore.
Conclusions
The reference relative timescale was shown to provide a viable alternative to the current standard method which uses chronological age as the time scale with left truncation for age at entry and a Cox model stratified on birth cohort. Its use led to a simplification of the modelling process and the scale possessed the defined features of a good timescale as defined in reliability theory. Simulation suggested that the two methods have similar power and are equally unbiased. The reference relative timescale has several interpretations and provides a unifying concept that links contemporary approaches in survival and reliability analysis to the traditional epidemiologic methods of Poisson regression and SMRs and can be understood as an extension of these traditional methods. The community of practitioners has previously failed to make this connection.
References
 1.
Korn EL, Graubard BI, Midhune D. Timetoevent analysis of longitudinal followup of a survey: choice of timescale. Am J Epidemiol. 1997;145:72–80.
 2.
Thiébaut ACM, Bénichou J. Choice of timescale in Cox’s model analysis of epidemiologic cohort data: a simulation study. Stat Med. 2004;23:3803–20.
 3.
Pencina MJ, Larson MG, D’Agostino RB. Choice of time scale and its effect on significance of predictors in longitudinal studies. Stat Med. 2007;26:1343–59.
 4.
Ingram DD, Makuc DM, Feldman JJ. Re: “Timetoevent analysis of longitudinal followup of a survey: choice of timescale” (Letter). Am J Epidemiol. 1997;146:528–9.
 5.
Cologne J, Hsu WL, Abbott RD, Ohishi W, Grant EJ, Fujiwara S, et al. Proportional hazards regression in epidemiologic followup studies: an intuitive consideration of primary time scale. Epidemiology. 2012;23:565–73.
 6.
Chalise P, Chicken E, McGee D. Baseline age effect on parametric estimates in Cox models. J Stat Comput Simul. 2012;82:1767–74.
 7.
Chalise P, Chicken E, McGee D. Performance and prediction for varying survival time scales. Commun Stat–Simul Comput. 2013;42:636–49.
 8.
Gail MH, Graubard B, Williamson DE, Flegal KM. Comments on “Choice of time scale and its effect on significance of predictors in longitudinal studies”. Stat Med. 2009;28:1315–8.
 9.
Griffin BA, Anderson GL, Shih RA, Whitsel EA. Use of alternative time scales in Cox proportional hazards models: implications for timevarying environmental exposures. Stat Med. 2012;31:3320–7.
 10.
Duchesne T, Lawless J. Alternative time scales and failure time models. Lifetime Data Anal. 2000;6:157–79.
 11.
Ederer F, Axtell LM, Cutler SJ. The relative survival rate: a statistical methodology. Nat Cancer Inst Monogr. 1961;6:101–21.
 12.
Andersen PK, BorchJohnsen K, Deckert T, Green A, Hougaard P, Keiding N, et al. A Cox regression model for relative mortality and its application to diabetes mellitus. Biometrics. 1985;41:921–32.
 13.
Andersen PK, Væth M. Simple parametric and nonparametric models for excess and relative mortality. Biometrics. 1989;45:523–35.
 14.
Clayton D, Hills M. Statistical models in epidemiology. Oxford: Oxford University Press; 1993.
 15.
Breslow NE, Day NE. Statistical methods in cancer research, volume II: the design and analysis of cohort studies. Lyon: IARC Scientific Publications; 1987.
 16.
Frome EL. The analysis of rates using Poisson regression models. Biometrics. 1983;39:665–74.
 17.
Tom BDM, Farewell T. Statistical methods for individuallevel data in cohort mortality studies of rheumatic diseases. Commun Stat Theor Methods. 2009;38:3472–87.
 18.
Loomis D, Richardson DB, Elliott L. Poisson regression analysis of ungrouped data. Occup Environ Med. 2005;62:325–9.
 19.
Stare J, Henderson R, Pohar M. An individual measure of relative survival. J R Stat Soc Ser CAppl Stat. 2005;54:115–26.
 20.
Pohar Perme M, Henderson R, Stare J. An approach to estimation in relative survival regression. Biostatistics. 2009;10:136–46.
 21.
Hodgson JT, Jones RD. Mortality of workers in the British cotton industry in 1968–1984. Scand J Work Environ Health. 1990;16:113–20.
 22.
McElvenny DM, Hurley MA, Lenters V, Heederik D, Wilkinson S, Coggon D. Lung cancer mortality in a cohort of UK cotton workers; an extended followup. Br J Cancer. 2011;105:1054–60.
 23.
Hurley MA. Light smoking at baseline predicts a higher mortality risk to women than to men; evidence from a cohort with long followup. BMC Public Health. 2014;14:95–103.
 24.
UK Office for National Statistics. http://www.ons.gov.uk/ons/rel/vsob1/deathregsumtables/2013/index.html. Accessed 19 Jun 2015.
 25.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for statistical computing, Vienna, Austria. ISBN: 3900051070. 2010. http://www.Rproject.org. Accessed 9 Feb 2015.
 26.
Cox DR, Isham V. Point processes. London: Chapman and Hall; 1980.
 27.
Dickman PW, Sloggett A, Hills M, Hakulinen T. Regression models for relative survival. Stat Med. 2004;23:51–64.
 28.
Nelson CP, Lambert PC, Squire IB, Jones DR. Flexible parametric models for relative survival, with application in coronary heart disease. Stat Med. 2007;26:5486–98.
 29.
Nelson CP, Lambert PC, Squire IB, Jones DR. Relative survival: what can cardiovascular disease learn from cancer. Eur Heart J. 2008;29:941–7.
 30.
Mahboubi A, Abrahamowicz M, Giorgi R, Binquet C, BonithonKopp C, Quantin C. Flexible modelling of the effects of continuous prognostic factors in relative survival. Stat Med. 2011;30:1351–65.
 31.
Thompson S, Kaptoge S, White I, Wood A, Perry P, Danesh J, et al. Statistical methods for timetoevent analysis of individual participant data from multiple epidemiological studies. Int J Epidemiol. 2010;39:1345–59.
 32.
Fu WJ, Land KC, Yang Y. On the intrinsic estimator and constrained estimators in ageperiodcohort models. Sociol Methodol. 2011;40:453–66.
Acknowledgements
Grateful thanks go to Pauline Wall and colleagues at the NHSIC for issue of death registrations, and to Sam Wilkinson, Elaine Brown and Damien McElvenny for reestablishing the UK Cotton Workers’ Cohort. The updating of the cohort data was supported in part by funding from the Colt Foundation (CF/05/07).
Competing interests
The author declare that she has no competing interests.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Operational failure time
 Life expectancy
 Time transformation
 Ageatrisk
 Cumulative hazard
 Attained age