The statistical analysis of truncated data: application to the Sverdlovsk anthrax outbreak

An outbreak of anthrax occurred in the city of Sverdlovsk in Russia in the spring of 1979. The outbreak was due to the inhalation of spores that were accidentally released from a military microbiology facility. In response to the outbreak a public health intervention was mounted that included distribution of antibiotics and vaccine. The objective of this paper is to develop and apply statistical methodology to analyse the Sverdlovsk outbreak, and in particular to estimate the incubation period of inhalational anthrax and the number of deaths that may have been prevented by the public health intervention. The data available for analysis from this common source epidemic are the incubation periods of reported deaths. The statistical problem is that incubation periods are truncated because some individuals may have had their deaths prevented by the public health interventions and thus are not included in the data. However, it is not known how many persons received the intervention or how efficacious was the intervention. A likelihood function is formulated that accounts for the effects of truncation. The likelihood is decomposed into a binomial likelihood with unknown sample size and a conditional likelihood for the incubation periods. The methods are extended to allow for a phase-in of the intervention over time. Assuming a lognormal model for the incubation period distribution, the median and mean incubation periods were estimated to be 11.0 and 14.2 days respectively. These estimates are longer than have been previously reported in the literature. The death toll from the Sverdlovsk anthrax outbreak could have been about 14% larger had there not been a public health intervention; however, the confidence intervals are wide (95% CI 0-61%). The sensitivity of the results to model assumptions and the parametric model for the incubation period distribution are investigated. The results are useful for determining how long antibiotic therapy should be continued in suspected anthrax cases and also for estimating the ultimate number of deaths in a new outbreak in the absence of any public health interventions.     

The mitigating effect of masks on the spread of Covid-19

We use US state-level data from early in the pandemic —March 15, 2020 to November 15, 2020— to estimate the effects of mask mandates and compliance with mandates on Covid-19 cases and deaths, conditional on mobility. A one-standard-deviation increase in mobility is associated with a 6 to 20 percent increase in the cases growth rate; a mask mandate can offset about one third of this increase with our most conservative estimates. Also, mask mandates are more effective in states with higher compliance. Given realized mobility, our estimates imply that total infections in the US on November 15, 2020 would have been 23.7 to 30.4 percent lower if a national mask mandate had been enacted on May 15, 2020. This reduction in cases translates to a 25 to 35 percent smaller decline in aggregate hours worked over the same period relative to a 2019 baseline.     

Studies of the mortality of atomic bomb survivors. report 12, part I. Cancer: 1950-1990

This continues the series of periodic general reports on cancer mortality in the cohort of A-bomb survivors followed by the Radiation Effects Research Foundation. The follow-up is extended by the 5 years 1986-1990, and analysis includes an additional 10,500 survivors with recently estimated radiation doses. Together these extensions add about 550,000 person-years of follow-up. The cohort analyzed consists of 86,572 subjects, of which about 60% have dose estimates of at least 0.005 Sv. During 1950-1990 there have been 3086 and 4741 cancer deaths for the less than and greater than 0.005 Sv groups, respectively. It is estimated that among these there have been approximately 420 excess cancer deaths during 1950-1990, of which about 85 were due to leukemia. For cancers other than leukemia (solid cancers), about 25% of the excess deaths in 1950-1990 occurred during the last 5 years; for those exposed as children this figure is nearly 50%. For leukemia only about 3% of the excess deaths in 1950-1990 occurred in the last 5 years. Whereas most of the excess for leukemia occurred in the first 15 years after exposure, for solid cancers the pattern of excess risk is apparently more like a life-long elevation of the natural age-specific cancer risk. Taking advantage of the lengthening follow-up, increased attention is given to clarifying temporal patterns of the excess cancer risk. Emphasis is placed on describing these patterns in terms of absolute excess risk, as well as relative risk. For example: (a) although it is becoming clearer that the excess relative risk for those exposed as children has declined over the follow-up, the excess absolute risk has increased rapidly with time; and (b) although the excess relative risk at a given age depends substantially on sex and age at exposure, the age-specific excess absolute risk depends little on these factors. The primary estimates of excess risk are now given as specific to sex and age at exposure, and these include projections of dose-specific lifetime risks for this cohort. The excess lifetime risk per sievert for solid cancers for those exposed at age 30 is estimated at 0.10 and 0.14 for males and females, respectively. Those exposed at age 50 have about one-third these risks. Projection of lifetime risks for those exposed at age 10 is more uncertain. Under a reasonable set of assumptions, estimates for this group range from about 1.0-1.8 times the estimates for those exposed at age 30. The excess life-time risk for leukemia at 1 Sv for those exposed at either 10 or 30 years is estimated as about 0.015 and 0.008 for males and females, respectively. Those exposed at age 50 have about two-thirds that risk. Excess risks for solid cancer appear quite linear up to about 3 Sv, but for leukemia apparent nonlinearity in dose results in risks at 0.1 Sv estimated at about 1/20 of those for 1.0 Sv. Site-specific risk estimates are given, but it is urged that great care be taken in interpreting these, because most of their variation can be explained simply by imprecision in the estimates.     

Nonparametric estimation of the incubation period of AIDS based on a prevalent cohort with unknown infection times.

Estimation of the incubation period distribution of human immunodeficiency virus based on prevalent cohorts of subjects, already infected at the time of recruitment, is complicated by the absence of information on the original times of infection. Here, we overcome this difficulty by using a prior distribution for the infection times, based on external data. Our estimate is nonparametric, but uses smoothness assumptions to avoid instability. The method is illustrated on two prevalent cohorts from San Francisco, separately and combined. The estimates produced agree with other published estimates of the incubation period distribution.     

Reductions in Cardiovascular,Cerebrovascular, and Respiratory Mortality following the National Irish Smoking Ban: Interrupted Time-Series Analysis

Background

Previous studies have shown decreases in cardiovascular mortality following the implementation of comprehensive smoking bans. It is not known whether cerebrovascular or respiratory mortality decreases post-ban. On March 29, 2004, the Republic of Ireland became the first country in the world to implement a national workplace smoking ban. The aim of this study was to assess the effect of this policy on all-cause and cause-specific, non-trauma mortality.

Methods

A time-series epidemiologic assessment was conducted, utilizing Poisson regression to examine weekly age and gender-standardized rates for 215,878 non-trauma deaths in the Irish population, ages ≥35 years. The study period was from January 1, 2000, to December 31, 2007, with a post-ban follow-up of 3.75 years. All models were adjusted for time trend, season, influenza, and smoking prevalence.

Results

Following ban implementation, an immediate 13% decrease in all-cause mortality (RR: 0.87; 95% CI: 0.76–0.99), a 26% reduction in ischemic heart disease (IHD) (RR: 0.74; 95% CI: 0.63–0.88), a 32% reduction in stroke (RR: 0.68; 95% CI: 0.54–0.85), and a 38% reduction in chronic obstructive pulmonary disease (COPD) (RR: 0.62; 95% CI: 0.46–0.83) mortality was observed. Post-ban reductions in IHD, stroke, and COPD mortalities were seen in ages ≥65 years, but not in ages 35–64 years. COPD mortality reductions were found only in females (RR: 0.47; 95% CI: 0.32–0.70). Post-ban annual trend reductions were not detected for any smoking-related causes of death. Unadjusted estimates indicate that 3,726 (95% CI: 2,305–4,629) smoking-related deaths were likely prevented post-ban. Mortality decreases were primarily due to reductions in passive smoking.

Conclusions

The national Irish smoking ban was associated with immediate reductions in early mortality. Importantly, post-ban risk differences did not change with a longer follow-up period. This study corroborates previous evidence for cardiovascular causes, and is the first to demonstrate reductions in cerebrovascular and respiratory causes.     

Epidemiological determinants of the pattern and magnitude of the vCJD epidemic in Great Britain.

Understanding the epidemiology and aetiology of new-variant Creutzfeldt-Jakob (vCJD) disease in humans has become increasingly important given the scientific evidence linking it to bovine spongiform encephalopathy (BSE) in cattle and hence the wide exposure of the population of Great Britain (GB) to potentially infectious tissue. The recent analysis undertaken to determine the risk to the population from dorsal route ganglia illustrated the danger in presenting point estimates rather than ranges of scenarios in the face of uncertainty. We present a mathematical template that relates the past pattern of the BSE epidemic in cattle to the future course of any vCJD epidemic in humans, and use extensive scenario analysis to explore the wide range of possible outcomes given the uncertainty in epidemiological determinants. We demonstrate that the average number of humans infected by one infectious bovine and the incubation period distribution are the two epidemiological factors that have the greatest impact on epidemic size and duration. Using the time-series of the BSE epidemic and the cases seen to date, we show that the minimum length of the incubation period is approximately nine years, and that at least 20% of the cases diagnosed to date were exposed prior to 1986. We also demonstrate that the current age distribution of vCJD cases can only arise if younger people were either exposed to a greater extent, more susceptible to infection, or have shorter incubation periods. Extensive scenario analyses show that given the information currently available, the very high degree of uncertainty in the future size of the epidemic will remain for the next 3-5 years. Furthermore, we demonstrate that this uncertainty is unlikely to be reduced by mass screening for late-stage infection.     

Estimation of cancer mortality rates: a Bayesian analysis of small frequencies

A Bayesian method is presented for estimating mortality rates of specific diseases when the frequency of deaths over a specified time period is assumed to have a Poisson distribution with mean proportional to the population size. The estimators use information from related populations, each having its own rate which is assumed distributed according to a common prior distribution about which some information is available. The study was motivated by an epidemiological study on the geographic variation of cancer mortality in the state of Missouri. Data from this study are used to illustrate the method and to compare it to a somewhat simpler empirical Bayes method.     

Estimating the size of the HIV epidemic by using mortality data

Evidence that more people are dying as a result of HIV infection than is reflected by the number of deaths among reported cases meeting the WHO definition of AIDS is derived from mortality data. Ninety-five causes of death likely to be associated with HIV infection were selected. Standardized mortality ratios due to these causes increased for single men aged 15-54 years from 100 in 1984 to 118 in 1987. The age, sex, marital status, temporal and geographic distribution of these excess deaths suggest that they are HIV-associated. It is estimated that 58% of excess deaths due to HIV-related causes were among cases reported to the CDSC AIDS Surveillance Programme in 1987. Some of these deaths may have been among HIV-positive people who did not meet the WHO definition at the time of death. There is a need for surveillance to be extended to include HIV-positive people who die before meeting the WHO definition if the full extent of the HIV epidemic is to be identified.     

Time trends of thyroid cancer incidence in Belarus after the Chernobyl accident.

The rates of childhood thyroid cancer incidence observed in Belarus during the period 1986 to 1995 are described as a function of time after exposure, age at exposure, and sex. Conclusions are drawn for the excess absolute risk function. After a minimum latent period of about 3 years after exposure, this risk function has a linear increase with time for at least 6 years. After correction for the dependence of average doses on age, the radiation-induced absolute thyroid risk in Gomel is about a factor of 3 higher for children up to age 10 at exposure compared to older ones; this may be due in part to different case-collection quality. In addition, in the group up to 10 years at exposure, the thyroid of girls is more sensitive to radiation by a factor of about 1.5 than the thyroid of boys on an absolute scale. Risk estimates from external exposure are consistent with risk estimates from Gomel assuming that the increase in excess cases reaches a plateau soon.     

Determination of the spread of HIV from the AIDS incidence history

The relation between the incidence of HIV in the general population, the number of AIDS cases, and the incubation period for the disease is examined. The number of AIDS cases can be expressed in terms of a convolution integral over the incubation period distribution and the temporal history of HIV incidence. In order to determine the level of HIV incidence it is necessary to invert the convolution. In this manner, it is possible to determine the spread of HIV up to the present time from knowledge of the AIDS incidence history and the incubation period. We describe the inversion of the convolution in terms of a Laplace transform technique that is applicable for any given incubation period distribution. Substantial simplifications in the technique are found in the case of an Erlang distribution for the probability density. The spread of HIV infections in the United States is charted through 1988 using AIDS incidence data that are corrected for both the revised AIDS case definition and reporting time delays. The results are consistent with current estimates of the HIV incidence in the United States and show no evidence of saturation in the rate of new infections. Indeed, the rate of new infections still appears to be climbing as of that date. While the technique is unable to predict the future course of the epidemic, it may provide a useful benchmark for comparison with mathematical models of the epidemic. The techniques are conceptually applicable to diseases other than AIDS.     

Subjective duration of every three-year period for 3 to 18 years of age, estimated by students

Results indicated that both subjective time durations, one day and three years, are well described by the normal distribution if the time scale is logarithmic rather than linear. A remarkable finding obtained here is that in reference to the mean subjective time duration of the junior high school (three years), the mean subjective time durations for each of the three-year periods before elementary school were significantly longer than unity (about 1.2), and the subjective duration of the three year period of the senior high school is significantly shorter (about 0.8). When subjects transferred to other elementary schools at 6 to 11 years of age, the mean subjective time duration of the periods below 11 years of age was much longer, more than 1.3.     

Improving estimates of the basic reproductive ratio: using both the mean and the dispersal of transition times

In both within-host and epidemiological models of pathogen dynamics, the basic reproductive ratio, R(0), is a powerful tool for gauging the risk associated with an emerging pathogen, or for estimating the magnitude of required control measures. Techniques for estimating R(0), either from incidence data or in-host clinical measures, often rely on estimates of mean transition times, that is, the mean time before recovery, death or quarantine occurs. In many cases, however, either data or intuition may provide additional information about the dispersal of these transition times about the mean, even if the precise form of the underlying probability distribution remains unknown. For example, we may know that recovery typically occurs within a few days of the mean recovery time. In this paper we elucidate common situations in which R(0) is sensitive to the dispersal of transition times about their respective means. We then provide simple correction factors that may be applied to improve estimates of R(0) when not only the mean but also the standard deviation of transition times out of the infectious state can be estimated.     

Assessment of the prevalence of vCJD through testing tonsils and appendices for abnormal prion protein

The objective of this study was to determine the age group or groups which will provide the most information on the potential size of the vCJD epidemic in Great Britain via the sampling of tonsil and appendix material to detect the presence of abnormal prion protein (PrP(Sc)). A subsidiary aim was to determine the degree to which such an anonymous age-stratified testing programme will reduce current uncertainties in the size of the epidemic in future years. A cohort- and time-stratified model was used to generate epidemic scenarios consistent with the observed vCJD case incidence. These scenarios, together with data on the age distribution of tonsillectomies and appendectomies, were used to evaluate the optimal age group and calendar time for undertaking testing and to calculate the range of epidemic sizes consistent with different outcomes. The analyses suggested that the optimal five-year age group to test is 25-29 years, although a random sample of appendix tissue from all age groups is nearly as informative. A random sample of tonsil tissue from all age groups is less informative, but the information content is improved if sampling is restricted to tissues removed from those over ten years of age. Based on the assumption that the test is able to detect infection in the last 75% of the incubation period, zero detected infections in an initial random sample of 1000 tissues would suggest that the epidemic will be less than 870,000 cases. If infections are detected, then the model prediction suggests that both relatively small epidemics (800+ cases if one is detected or 8300+ if two are detected) and larger epidemics (21,000+ cases if three or more are detected) are possible. It was concluded that testing will be most informative if undertaken using appendix tissues or tonsil tissues removed from those over ten years of age. Large epidemics can only be excluded if a small number of infections are detected and the test is able to detect infection early in the incubation period.     

Mortality rates among wild chimpanzees

In order to compare evolved human and chimpanzees' life histories we present a synthetic life table for free-living chimpanzees, derived from data collected in five study populations (Gombe, Ta?, Kibale, Mahale, Bossou). The combined data from all populations represent 3711 chimpanzee years at risk and 278 deaths. Males show higher mortality than females and data suggest some inter-site variation in mortality. Despite this variation, however, wild chimpanzees generally have a life expectancy at birth of less than 15 years and mean adult lifespan (after sexual maturity) is only about 15 years. This is considerably lower survival than that reported for chimpanzees in zoos or captive breeding colonies, or that measured among modern human hunter-gatherers. The low mortality rate of human foragers relative to chimpanzees in the early adult years may partially explain why humans have evolved to senesce later than chimpanzees, and have a longer juvenile period.     

The impact of future summer temperature on public health in Barcelona and Catalonia,Spain

Several epidemiological studies have reported associations between increases in summer temperatures and risks of premature mortality. The quantitative implications of predicted future increases in summer temperature, however, have not been extensively characterized. We have quantified these effects for the four main cities in Catalonia, Spain (Barcelona, Tarragona, Lleida, Girona). We first used case-crossover analysis to estimate the association between temperature and mortality for each of these cities for the period 1983 to 2006. These exposure–response (ER) functions were then combined with local measures of current and projected changes in population, mortality and temperature for the years 2025 and 2050. Predicted daily mean temperatures were based on the A1B greenhouse gas emission, “business-as-usual” scenario simulations derived from the ENSEMBLES project. Several different ER functions were examined and significant associations between temperature and mortality were observed for all four cities. For these four cities, the age-specific piecewise linear model predicts 520 (95%CI  340, 720) additional annual deaths attributable to the change in temperature in 2025 relative to the average from the baseline period of 1960–1990. For 2050, the estimate increases to 1,610 deaths per year during the warm season. For Catalonia as a whole, the point estimates for those two years are 720 and 2,330 deaths per year, respectively, or about 2 and 3% of the warm season. In comparing these predicted impacts with current causes of mortality, they clearly represent significant burdens to public health in Catalonia.     

The age of fine-root carbon in three forests of the eastern United States measured by radiocarbon

Using a new approach involving one-time measurements of radiocarbon (¹⁴C) in fine (<2 mm diameter) root tissues we have directly measured the mean age of fine-root carbon. We find that the carbon making up the standing stock of fine roots in deciduous and coniferous forests of the eastern United States has a mean age of 3-18 years for live fine roots, 10-18 years for dead fine roots, and 3-18 years for mixed live+dead fine roots. These ¹⁴C-derived mean ages represent the time C was stored in the plant before being allocated for root growth, plus the average lifespan (for live roots), plus the average time for the root to decompose (for dead roots and mixtures). Comparison of the ¹⁴C content of roots known to have grown within 1 year with the ¹⁴C of atmospheric CO₂ for the same period shows that root tissues are derived from recently fixed carbon, and the storage time prior to allocation is <2 years and likely <1 year. Fine-root mean ages tend to increase with depth in the soil. Live roots in the organic horizons are made of C fixed 3-8 years ago compared with 11-18 years in the mineral B horizons. The mean age of C in roots increases with root diameter and also is related to branching order. Our results differ dramatically from previous estimates of fine-root mean ages made using mass balance approaches and root-viewing cameras, which generally report life spans (mean ages for live roots) of a few months to 1-2 years. Each method for estimating fine-root dynamics, including this new radiocarbon method, has biases. Root-viewing approaches tend to emphasize more rapidly cycling roots, while radiocarbon ages tend to reflect those components that persist longest in the soil. Our ¹⁴C-derived estimates of long mean ages can be reconciled with faster estimates only if fine-root populations have varying rates of root mortality and decomposition. Our results indicate that a standard definition of fine roots, as those with diameters of <2 mm, is inadequate to determine the most dynamic portion of the root population. Recognition of the variability in fine-root dynamics is necessary to obtain better estimates of belowground C inputs.     

The Very High Premature Mortality Rate among Active Professional Wrestlers Is Primarily Due to Cardiovascular Disease

Purpose

Recently, much media attention has been given to the premature deaths in professional wrestlers. Since no formal studies exist that have statistically examined the probability of premature mortality in professional wrestlers, we determined survival estimates for active wresters over the past quarter century to establish the factors contributing to the premature mortality of these individuals.

Methods

Data including cause of death was obtained from public records and wrestling publications in wrestlers who were active between January 1, 1985 and December 31, 2011. 557 males were considered consistently active wrestlers during this time period. 2007 published mortality rates from the Center for Disease Control were used to compare the general population to the wrestlers by age, BMI, time period, and cause of death. Survival estimates and Cox hazard regression models were fit to determine incident premature deaths and factors associated with lower survival. Cumulative incidence function (CIF) estimates given years wrestled was obtained using a competing risks model for cause of death.

Results

The mortality for all wrestlers over the 26-year study period was.007 deaths/total person-years or 708 per 100,000 per year, and 16% of deaths occurred below age 50 years. Among wrestlers, the leading cause of deaths based on CIF was cardiovascular-related (38%). For cardiovascular-related deaths, drug overdose-related deaths and cancer deaths, wrestler mortality rates were respectively 15.1, 122.7 and 6.4 times greater than those of males in the general population. Survival estimates from hazard models indicated that BMI is significantly associated with the hazard of death from total time wrestling (p<0.0001).

Conclusion

Professional wrestlers are more likely to die prematurely from cardiovascular disease compared to the general population and morbidly obese wrestlers are especially at risk. Results from this study may be useful for professional wrestlers, as well as wellness policy and medical care implementation.     

The effect of pharmaceutical innovation on longevity: Evidence from the U.S. and 26 high-income countries

This study examines the impact that pharmaceutical innovation, which accounts for most private biomedical research expenditure, has had on longevity. We perform two types of two-way fixed-effects analyses, which control for the effects of many potentially confounding variables. First, we analyze long-run (2006–2018) changes in longevity associated with different diseases in a single country: the U.S. Then, we analyze relative longevity levels associated with different diseases in 26 high-income countries during a single time period (2006–2016). The measure of longevity we analyze, mean age at time of death, is strongly positively correlated across countries with life expectancy at birth. The measure of pharmaceutical innovation we use is the mean vintage (year of initial world launch) of the drugs used to treat each disease in each country. Changes in the vintage distribution of drugs are due to both entry of new drugs and exit of old drugs. Our analysis of U.S. data indicates that the diseases for which there were larger increases in drug vintage tended to have larger increases in the longevity of Americans of all races and both sexes. In other words, the lower the mean age of the drugs, the higher the mean age at death. We test, and are unable to reject, the “parallel trends” hypothesis. We estimate that the 2006–2018 increase in drug vintage increased the mean age at death of Americans by about 6 months (66% of the observed increase). Controlling for sex, race, and education has only a small effect on the estimate of the vintage coefficient. The estimates indicate that drug vintage did not have a significant effect on the mean age at death of decedents with less than 9 years of education. Drug vintage had a positive and significant effect on the mean age at death of decedents with at least 9 years of education, and a larger effect on the mean age at death of decedents with at least 13 years of education. The finding that pharmaceutical innovation has a larger effect on the longevity of people with more education is consistent with previous evidence that more educated people are more likely to use newer drugs. Our analysis of data on 26 high-income countries indicates that the higher the vintage of drugs available to treat a disease in a country, the higher mean age at death was, controlling for fixed disease and country effects. The increase in drug vintage is estimated to have increased mean age at death in the 26 countries by 1.23 years between 2006 and 2016—73% of the observed increase. We obtain estimates of the cost of pharmaceutical innovation—its impact on drug expenditure—as well as estimates of an important benefit of pharmaceutical innovation—the number of life-years gained from it—and of their ratio, i.e., the incremental cost-effectiveness ratio. Estimates of the cost per life-year gained for the U.S. and the 26 countries are $35,817 and $13,904, respectively. Both figures are well below per capita GDP in the respective regions, suggesting that, overall, pharmaceutical innovation was highly cost-effective.     

Freshwater sediment pesticide biodegradation potential as an ecological indicator of microbial recovery following a decrease in chronic pesticide exposure: A case study with the herbicide diuron

The aim of this study was to evaluate the use of freshwater sediment biodegradation potential as an ecological indicator for monitoring microbial recovery following a decrease in chronic pesticide exposure. For this purpose, a four-year case study (2008–2011) was conducted in a small stream (Morcille river) long exposed to high diuron concentrations, increasing from upstream to downstream. Our results show that the ban on diuron in December 2008 resulted in a progressive decrease in its concentrations in the Morcille river over the survey period. However, diuron remained present in the water three years after the ban. The spatio-temporal variations in the sediment biodegradation potential were assessed by radiorespirometry using [ring-U-¹⁴C] diuron to estimate diuron mineralization potentials. Between autumn 2008 and autumn 2011, mean diuron mineralization percentage after 15 weeks of incubation decreased by 65% downstream and by 82% in the intermediate sector, and mean 10% diuron dissipation time values increased between 143% (downstream) and 210% (intermediate). Thus the decrease in the level of chronic diuron exposure in the river also caused a fall in sediment diuron-mineralizing capacities, revealing a corresponding recovery of microbial communities. Our results show that the use of freshwater sediment biodegradation potential may be useful for assessing microbial recovery after a decrease in chronic exposure to pollutants, opening prospects for developing a new class of ecological indicator to monitor the recovery of biological quality of water resources. In this way, the use of molecular approaches based on direct extraction of nucleic acids from environmental matrices and their subsequent analysis by PCR-based approaches to quantify the abundance of pesticide-degrading communities could represent a promising alternative.     

Long-term survival of patients with breast cancer: a study of the curability of the disease.

A retrospective analysis was made of 3878 cases of breast carcinoma first seen in Edinburgh from 1954 to 1964. During this time there was a policy to treat breast cancer by simple mastectomy and x-ray therapy, and over 90% of cases classified as international stages I and II were so treated. The mortality in these women was compared with that in an equivalent normal population using Scottish national age-specific death rates. For every year of follow-up within 20 years of initial treatment there was an excess mortality from all causes. There was an overall excess mortality of 58% among patients with breast cancer 15-20 years after initial treatment, and 20 times more deaths occurred in this period from breast cancer than in a normal population. For patients disease-free after 15 years there was still a 28% excess mortality from all causes. Factors known to be of major prognostic significance for five-year survivorship had less influence than might have been expected when the ratio of observed to expected deaths was considered for longer periods of follow-up. The effect of clinical staging (I, II, or III), though initially marked, largely disappeared by the 10th year of follow-up, and after allowing for age there was no evidence beyond 10 years of an effect on survival of the original stage of the disease. Similarly, the effect of tumour size on survival disappeared after 10 years. Women who were premenopausal at presentation still had a significant excess of deaths in the fourth quinquennium of follow-up. In the menopausal and postmenopausal groups combined there was still a small non-significant excess of deaths from all causes after 15 years but this almost disappeared when patients who had already relapsed were excluded. In terms of overall mortality only patients who have undergone the menopause before presentation and who are disease-free 15 years after primary treatment may prove to be cured by conventional techniques such as simple mastectomy and postoperative radiotherapy.