An animated depiction of major depression epidemiology
© Patten; licensee BioMed Central Ltd. 2007
Received: 28 September 2006
Accepted: 08 June 2007
Published: 08 June 2007
Epidemiologic estimates are now available for a variety of parameters related to major depression epidemiology (incidence, prevalence, etc.). These estimates are potentially useful for policy and planning purposes, but it is first necessary that they be synthesized into a coherent picture of the epidemiology of the condition. Several attempts to do so have been made using mathematical modeling procedures. However, this information is not easy to communicate to users of epidemiological data (clinicians, administrators, policy makers).
In this study, up-to-date data on major depression epidemiology were integrated using a discrete event simulation model. The mathematical model was animated in Virtual Reality Modeling Language (VRML) to create a visual, rather than mathematical, depiction of the epidemiology.
Consistent with existing literature, the model highlights potential advantages of population health strategies that emphasize access to effective long-term treatment. The paper contains a web-link to the animation.
Visual animation of epidemiological results may be an effective knowledge translation tool. In clinical practice, such animations could potentially assist with patient education and enhanced long-term compliance.
A comprehensive picture of the epidemiology of an episodic condition cannot be embodied in a single parameter such as an estimate of incidence or prevalence. The incidence of major depressive disorder conveys information about the risk of onset of an initial episode, but this is of limited value for clinical practice (where presentations occur after onset) or for public policy (since the rate of onset of new episodes does not adequately reflect the burden of disease). Prevalence is in some respects an even more obscure parameter. Prevalence is the number or proportion ("prevalence proportion") of population members who have a disorder at a point ("point prevalence") in time or during a period ("period prevalence") in time. This parameter is salient to the problem of quantifying disease burden but ultimately reflects a balance between entry into the prevalence pool (incidence or recurrence) and a other factors such as episode duration and mortality.
Attempts to integrate epidemiological estimates into a comprehensive "picture" of the epidemiology of major depression have used Markov models [1, 2] and micro-simulation approaches . The complexity of major depression epidemiology has required the development of fairly complex models, the mathematics of which may not be of great interest to many users of such information. A solution may be found in the use of modeling techniques that can be linked to visual depictions of the population dynamics. In this way, complex information can be processed visually, rather than by using equations, transition matrices etc. In this project, a set of plausible estimates for key parameters were extracted from available sources, supplemented by information from the international literature. These were combined using a simulation model, and a visual depiction of the epidemiology was produced in the form of an animation.
In a previous study, data from a longitudinal study called the National Population Health Survey (NPHS) was used to develop a prognosis calculator for major depression . This prior project was similar to the current one in its use of a discrete event simulation modeling approach, but the previous analysis only examined the issue of episode duration. The current project uses the best available data from the literature rather than using data from a single survey source. It also attempts to integrate the various estimates into a more comprehensive description of the epidemiology.
Selection of epidemiological calibration targets for the simulation model
Targets for model calibration
High in teenagers – approximately 7.1% in female adolescents, 4.4% in male adolescents. Declines with age.
Use of a Weibull distribution, to reflect declining incidence with age. Separate scale and shape parameters for men and women. Calibrate first incidence to result in lifetime prevalence of approximately 20% for men and 30% for women.
Overall, approximately 2%. Higher in women than in men, declines with age.
Relative risk of recurrence selected (see below) to predict the point prevalence.
Approximately 10% in first year after recovery from an initial episode. Approximately 85% will have at least one recurrence during their lifetime. Recurrence rate higher after multiple episodes than single episodes.
First recurrence simulated with a Weibull distribution, men and women combined together with a single scale and shape parameter.
Relative risk for multiple recurrence was estimated representing a higher risk among those with prior episodes.
The relative risk for recurrence was calibrated to achieve a point prevalence of 1.8% holding the Weibull parameters for first incidence and episode duration constant.
Approximately 15% report a 2 week duration, the time at which the diagnosis is technically possible. However, a similar proportion is not recovered after 2 years.
Weibull distribution, a single scale and shape parameter for men and women.
Model calibrated by matching the simulated proportion with a two week and 2 year duration at 15% and 85%, respectively.
The simulation model was programmed in software called Arena . This software provides a graphical interface for development of SIMAN models using discrete event simulation. The model was set up so that simulated subjects or "entities" enter at age 15. Following Kruijshaar , the simulation began by assignment of a death date to each simulated entity (depicting a person), based on sex and age-specific mortality data for Canada, which was downloaded from the Statistics Canada website using an on-line data access utility called E-STAT . Next, an incidence date was simulated using a cumulative Weibull distribution. This distribution was selected to reflect the expected decline in incidence with advancing age (see Additional file 1). When the incidence date occurred before the death date, an episode would occur during that subject's lifetime, and in this case the death date was recalculated using a relative risk for mortality of 1.1. A meta-analysis by Wulsin et al.  estimated the relative risk at 1.8, but this is almost certainly an overestimate as most of the studies included in the meta-analysis were clinical studies. A program called OptQuest, which is a part of the Arena package, was used to identify Weibull parameters (shape and scale) for incidence that were consistent both with estimates for the adolescent age group (which have been found to be very high in some studies, see Additional file 1), but diminishing at a sufficient rate to be consistent with the selected values for lifetime prevalence. OptQuest conducts a series of Monte Carlo simulations in order to identify parameters that closely fit modeling objectives.
A similar strategy was employed to simulate other elements in the model. Weibull distributions for episode duration, and the duration of remitted periods between episodes were identified. After resolution of an initial episode, the model entities moved into a first remission state. This approach resembles that of Kruijshaar et al. , except that Kruijshaar et al.  used lognormal distributions to describe episode duration. Previous work has determined that the Weibull distribution can also provide a good description . OptQuest was used to find Weibull parameters for episode duration consistent with the selected calibration targets: approximately 15% recovered before three weeks, but approximately 15% are still not recovered after 2 years, see Appendix I for a more detailed discussion. Once appropriate values for these Weibull parameters were found, the relative risk of recurrence was sought, again using OptQuest, and calibrated to a targeted point prevalence of 1.8% .
In the model, when the simulated death date was reached, the entity left the simulation. Otherwise, subjects with depressive episodes continue to cycle between depressed and remitted states until their death date is reached. Eventually, a balance between entry and exit (death) is reached, and the various calibration parameters assume steady state values. In order to reach this state, the model was run for 250,000 days. This steady state allowed the model to be calibrated, using OptQuest as described above.
Making a visual animation
Arena can save the results of a simulation model into a file type (which Arena calls a "playback" file) that can be imported into a Virtual Reality Modeling Language (VRML) program and animated in three dimensions using a program called the Arena 3-D player. Video can be captured from the animation. The animations presented in this paper were captured as "avi" files, edited with addition of captions in a program called "Camtasia" and then saved as MacroMedia Flash files. These files were uploaded to the University of Calgary's Dspace under a Creative Commons license for open access. The Flash movie contains three parts. In the first part, the structure of the visual animation is presented. In the second part, the fully calibrated model is depicted. In the third section the model is animated with the recurrence pattern in subjects with multiple episodes having been changed to that of subjects with single episodes. The third part is a way of depicting the population health impact of highly effective tertiary prevention, ie. long-term treatment substantially reducing the relapse rate. The simulations were run at 500,000 times "real time." This means that each second in the simulation represents about six days of time. The rate of entry into the population was set so that approximately 1000 simulated people (entities) would be present in the simulation at steady state. In order that subjects could be seen moving from state to state, transitions were set to last for two weeks in simulation time. The two week transition time seems clinically reasonable, and avoids a problem that would arise in the animations if the transitions occurred too quickly to be seen easily.
Optimized estimates: Weibull distributions*
Mean values for model parameters (%), after N = 100 replications
Mortality sensitivity analysis*
Strong tertiary prevention†
Incidence at age 15 (♂)
Incidence at age 15 (♀)
Total lifetime prevalence
Lifetime prevalence (♂)
Lifetime prevalence (♀)
Episodes lasting < 3 weeks
Episodes lasting ≥ 2 years (%)
Proportion MDD with single episode
Lifetime morbidity proportion**
Proportion relapse within 1 year
As expected, mortality has an effect on point prevalence, but the large change in associated mortality has only a small effect. On the other hand, changing the recurrence rate has a large effect on point prevalence.
A variety of analyses linking model variables (e.g. scale and shape parameters from the Weibull distributions) to incidence, recurrence and prevalence are presented in Appendix II, see Additional file 2. This Appendix also contains a series of simulations, each with n = 25 replications showing the effect on point prevalence of alteration to various model parameters.
The simulation model presented represents a crude synthesis of available data about the epidemiology of major depression. Rather than accept lifetime prevalence estimates made directly from epidemiological studies, thereby ignoring evidence of recall bias that is likely to affect such estimates , an estimate of the probable lifetime prevalence based on a combination of incidence and prevalence data from Kruijshaaar  was used. The mortality estimate is also somewhat arbitrary. Prior studies pointed towards a higher relative risk for mortality in subjects with major depression, but these are almost certainly over-estimates, as their samples consisted of hospitalized clinical subjects, see Additional file 1 (Appendix I). However, the estimate adopted here is somewhat arbitrary. A sensitivity analysis does not indicate that this makes a large difference to the simulation model, see Additional file 2 (Appendix II). Similarly, values for episode duration and recurrence rates that were used to calibrate the module represent reasonable values given published estimates. The rates of recovery from episodes as depicted here are lower than those reported from a Canadian data source, the National Population Health Survey . However, this study used a brief diagnostic instrument which may have picked up some sub-clinical cases .
Despite these limitations, an understanding of certain key elements of major depression epidemiology has solidified during the past decade: recurrence is greater than first incidence, incidence is high in young people, and declines with age, the recovery and remission rates are time dependent. The model reflects all of these elements of the condition. The model presents a visual impression of the interaction between a reasonable set of assumptions about the epidemiology of the condition.
The epidemiology of major depression is complex, and models that attempt to reflect this complexity may have a diminished impact because of their own inherent complexity. Depiction of a simulation model as a visual animation may assist with communicating its implications to a broad audience. Animations of this sort may also be helpful for patient education, and enhanced compliance as they illustrate some of the evidence for therapeutic recommendations. The results presented here are consistent with a conception of major depression as a chronic, recurrent condition  and suggest that relevant long-term strategies can have a large impact on population health status.
Dr. Patten is a Health Scholar with the Alberta Heritage Foundation for Medical Research and a Research Fellow with the Institute of Health Economics.
- Patten SB, Lee RC: Refining estimates of major depression incidence and episode duration in Canada using a Monte Carlo Markov model. Med Decis Making. 2004, 24: 351-358. 10.1177/0272989X04267008.View ArticlePubMedGoogle Scholar
- Patten SB, Lee RC: Describing the longitudinal course of major depression using Markov models: Data integration across three national surveys. Popul Hlth Metr. 2005, 3: 11-10.1186/1478-7954-3-11.View ArticleGoogle Scholar
- Kruijshaar ME, Barendregt J, Vos T, de Graaf R, Spijker J, Andrews G: Lifetime prevalence estimates of major depression: An indirect estimation method and a quantification of recall bias. Eur J Epidemiol. 2005, 20: 103-111. 10.1007/s10654-004-1009-0.View ArticlePubMedGoogle Scholar
- Patten SB: A major depression prognosis calculator based on episode duration. Clin Pract Epidemiol Mental Hlth. 2006, 2: 13-10.1186/1745-0179-2-13.View ArticleGoogle Scholar
- Waraich PS, Goldner EM, Somers JM, Hsu L: Prevalence and incidence studies of mood disorders: a systematic review of the literature. Can J Psychiatry. 2004, 49: 124-138.PubMedGoogle Scholar
- Andrews G, Poulton R, Skoog I: Lifetime risk of depression: restricted to a minority or waiting for most?. Br J Psychiatry. 2005, 187: 495-496. 10.1192/bjp.187.6.495.View ArticlePubMedGoogle Scholar
- Arena. 2006, Sewickley, PA, Rockwell Software Inc., [http://www.arenasimulation.com]
- Canada S: E-Stat. 2006, [http://www.statcan.ca/english/Estat/licence.htm]Google Scholar
- Wulsin LR, Vaillant GE, Wells VE: A systematic review of the mortality of depression. Psychosom Med. 1999, 61: 6-17.View ArticlePubMedGoogle Scholar
- Patten SB, Wang JL, Williams JV, Currie S, Beck CA, Maxwell CJ, El-Guebaly N: Descriptive epidemiology of major depression in Canada. Can J Psychiatry. 2006, 51: 84-90.PubMedGoogle Scholar
- Patten SB: Visual animation of major depression epidemiology. 2007, University of Calgary Dspace, [http://hdl.handle.net/1880/44649]Google Scholar
- Patten SB: An Animation of Major Depression Epidemiology (with Arena dashboards). 2007, University of Calgary DSpace, [http://hdl.handle.net/1880/44475]Google Scholar
- Andrews G, Anstey K, Brodaty H, Issakidis C, Luscombe G: Recall of depressive episode 25 years previously. Psychol Med. 1999, 29: 787-791. 10.1017/S0033291799008648.View ArticlePubMedGoogle Scholar
- Patten SB, Brandon-Christie J, Devji J, Sedmak B: Performance of the Composite International Diagnostic Interview Short Form for Major Depression in a community sample. Chron Dis Can. 2000, 21 (2): 68-72.Google Scholar
- Andrews G: Should depression be managed as a chronic disease?. Br Med J. 2001, 332: 419-421. 10.1136/bmj.322.7283.419.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-244X/7/23/prepub