Abstract
Background
Alcohol consumption is a major risk factor in the global burden of disease, with overall volume of exposure as the principal underlying dimension. Two main sources of data on volume of alcohol exposure are available: surveys and per capita consumption derived from routine statistics such as taxation. As both sources have significant problems, this paper presents an approach that triangulates information from both sources into disaggregated estimates in line with the overall level of per capita consumption.
Methods
A modeling approach was applied to the US using data from a large and representative survey, the National Epidemiologic Survey on Alcohol and Related Conditions. Different distributions (lognormal, gamma, Weibull) were used to model consumption among drinkers in subgroups defined by sex, age, and ethnicity. The gamma distribution was used to shift the fitted distributions in line with the overall volume as derived from per capita estimates. Implications for alcoholattributable fractions were presented, using liver cirrhosis as an example.
Results
The triangulation of survey data with aggregated per capita consumption data proved feasible and allowed for modeling of alcohol exposure disaggregated by sex, age, and ethnicity. These models can be used in combination with risk relations for burden of disease calculations. Sensitivity analyses showed that the gamma distribution chosen yielded very similar results in terms of fit and alcoholattributable mortality as the other tested distributions.
Conclusions
Modeling alcohol consumption via the gamma distribution was feasible. To further refine this approach, research should focus on the main assumptions underlying the approach to explore differences between volume estimates derived from surveys and per capita consumption figures.
Introduction
The volume of alcohol consumed has been shown to be causally related to more than 230 International Classification of Disease, version 10, disease codes [13]. Most of these relationships follow the component cause model [4], in which only a fraction (the socalled attributable fraction [5]), of the incidence of a disease would disappear if the causal component, in this case alcohol use [6], could be eliminated. Liver cirrhosis may be used as an example: alcohol use has been shown to have a causal impact on liver cirrhosis [79]; however, there are also cases of liver cirrhosis where alcohol use has not been involved, such as those stemming from HCV infection in nondrinkers.
The proportion of most diseases caused by alcohol in the component cause model in a population is determined by:
• The distribution of the volume of exposure
• The relative risk associated with each level of exposure, i.e., doseresponse relationship [10]
For most disease categories, the doseresponse relationship is nonlinear and varies by sex as well as age in some cases. Thus, to calculate the alcoholattributable fractions (AAFs) for estimating a population's burden of disease attributable to alcohol, we need to characterize the volume of alcohol exposure by sex and age.
There is one main problem with the characterization of the volume of alcohol exposure in populations: the best indicator, adult per capita consumption [11], is not available by sex and age. It is derived mainly from production, sales, export, and import figures, which are almost never disaggregated [12]. As an alternative, alcohol exposure can be measured from surveys, yet this has the disadvantage that adult per capita figures, the best indicator for consumption, are often severely underestimated [11]. To give but one recent example: the national Canadian Addiction Survey [13] had a coverage rate of between 30% and 40% of the adult per capita consumption. Thus, estimating overall volume of consumption in Canada based on answers from a representative survey results in a figure that is 60% to 70% lower than the figure derived from aggregate statistics mainly based on sales and taxation. The exact magnitude of underestimation depends on the alcohol measure used in the survey [14].
This underestimation of population exposure leads to two problems: first, the absolute level of exposure is incorrectly estimated, usually underestimated. And second, inter and intrapopulation distributions based on different surveys are not adequately comparable as the degree of incorrect estimation will vary. Recent large, nationally representative surveys have shown coverage rates between 25% and more than 100% [15]. Studies have thus had to explore additional means for achieving comparability. In the case of the Comparative Risk Assessment (CRA) for alcohol within the Global Burden of Disease (GBD) 2000 study, this was achieved by triangulating survey and per capita information based on sales or production [12,15]. In short, the distribution of volume by sex and age was taken from surveys, while the overall exposure was taken from adult per capita figures [15].
There is, however, more than one way to triangulate such data. The CRA 2000, for instance, used a categorical approach based on the standard categories of volume of drinking from English and colleagues [16]. For shifting the distribution of the survey to correspond to the per capita consumption, however, an additional assumption had to be made to obtain a unique solution, and the assumption chosen at the time was that the highest drinking category was a constant fraction of the percentage of the next highest category [15]. However, no empirical evidence supports this specific assumption. Therefore, a more evidencebased approach was sought for the ongoing CRA of the GBD 2005 study. Different distributions (gamma, lognormal, Weibull) were fitted to empirical data of surveys from 66 countries with the aim of identifying associations between distribution parameters that would enable shifting survey distributions to fit the volume of drinking indicated by adult per capita consumption data [17]. The twoparameter gamma distribution proved to be best suited for triangulating survey and per capita data by shifting distributions upward because:
1) As an inherent characteristic of this distribution, the means of fitted distributions are equal to those of the empirical distributions. Thus, there is no error involved in estimating means from the fitted distribution and no need for erroneous and complicated backtransformations from fitted distributions to original scales (as in the case of the lognormal distribution).
2) There was a very high correlation between means and standard deviations of the sex and agespecific, fitted twoparameter gamma distributions (r = 0.923; N = 851; p < 0.01). We then conducted a linear regression with the standard deviation as the dependent variable, and mean and sex as the independent variables, and could predict the standard deviation with precision (r = 0.971; N = 851; p < 0.01; for the derived prediction equation, see below). Thus, the shifted distributions could easily and reliably be derived from the mean (= per capita consumption) and the standard deviation as estimated by regression methods [17].
This article discusses the modeling approach described above in its application to US data, with three main objectives:
1. To model the volume of alcohol exposure in the US with three different distributions: lognormal, gamma, and Weibull distribution for different strata by sex, age and ethnicity. The volume of alcohol use was obtained from a large representative survey, the National Epidemiologic Survey on Alcohol and Related Conditions (NESARC).
2. To shift the alcohol use distribution to the level of adult per capita consumption.
3. To show the impact of shifting the alcohol use distribution on AAFs, using liver cirrhosis as an example.
Methods
Description of underlying survey (NESARC)
This analysis is based on data from the 20012002 NESARC, which was designed and sponsored by the National Institute on Alcohol Abuse and Alcoholism. The fieldwork for the NESARC was conducted by the US Census Bureau, with data collected in facetoface, computerassisted, inhome interviews. The NESARC sample represents the civilian, noninstitutionalized adult population of the United States, including the District of Columbia, Alaska, and Hawaii, and includes people living in households and military personnel living on and off base [18]. The NESARC oversampled African Americans, Hispanics, and adults aged 18 to 24. One sample adult (age 18 or older) was selected for an interview in each household. The overall response rate was 81% (N = 43,093).
The volume of ethanol consumption reflected consumption in the 12 months preceding the interview. The volume of ethanol intake was based on data summed over a separate series of questions for coolers, beer, wine, and distilled spirits. For each beverage, volume was estimated on the basis of: overall frequency of drinking; typical and largest quantities consumed; frequency of consuming the largest quantity; frequency of consuming fiveplus drinks; typical drink size; and ethanol content by volume of the brand usually consumed. The testretest reliabilities for the various measures of alcohol consumption from the 2001 to 2002 NESARC were good to excellent, with intraclass correlation coefficients ranging from 0.68 to 0.84 [19].
Methods for fitting the distributions
To find an appropriate model for alcohol consumption, we examined three distributions that were unimodal, had a density with only one maximum, and could be used to fit rightskewed empirical data: lognormal, gamma, and Weibull. The lognormal, gamma, and Weibull densities are similar in shape and mainly differ at the tails (i.e., at high levels of consumption). Alcohol consumption has been more commonly modeled using the lognormal distribution, mostly for historical reasons related to the socalled single distribution theory [20,21]. But the lognormal distribution also has been favored because it is easy to use as a transformation and results in sufficiently accurate values that permit fitting and testing hypotheses [22]. Although doubts regarding the single parametric lognormal distribution and its justification as the best approximation for the distribution of consumption have been previously raised, [23,24] it has appeared to provide good approximations for most applications [25,26]. Later developments on modeling alcohol exposure have favored more complex distributions such as gamma [27], or revealed that mixing distributions is needed to fit separate distributions for frequency of drinking and quantity of drinking [28].
The lognormal distribution, with parameters μ and σ, describes a random variable X where log x is normally distributed with mean μ and standard deviation σ. The lognormal distribution function, with parameters μ and σ, is given by:
Although alcohol consumption is frequently modeled using the lognormal distribution, empirical distributions often deviate considerably from the lognormal model [2729].
The gamma distribution has two parameters, a scale parameter θ and a shape parameter k. The gamma distribution is more adaptable than the lognormal distribution because it has the effect of stretching or compressing its range by changing the scale parameter θ. The gamma distribution has the following probability density function:
The Weibull distribution is one of the most widely used distributions in applied statistics, especially life data analysis, because of its versatility in fitting a variety of distributions. The probability density function of a random variable X having a twoparameter Weibull distribution with shape parameter γ and scale parameter θ is given by:
The shape parameter γ in the Weibull distribution gives this distribution its flexibility.
The maximum likelihood method of estimation was used to fit all three models  lognormal, gamma, and Weibull  to the data using the R language [30]. All missing values on volume of drinking  298 in total (men 185; women 113)  were excluded from the fitted models. All numerical integration used the trapezoidal rule with many subintervals to obtain more accurate estimates. The trapezoidal rule uses trapezoids instead of rectangles for approximating the definite interval over closed bounded intervals. The NewtonRaphson algorithm, having a quadratic rate of convergence, was used to optimize the likelihood equations solving for the unknown parameters with maximum likelihood estimates.
To compare different fitted distributions, we used chisquare tests by comparing expected frequencies (derived from the fitted distributions) and observed frequencies from the empirical distributions, using a bandwidth of 10 grams of ethanol for the frequencies. The use of chisquare tests to compare two distributions is a standard method [31].
Method for shifting the distribution
As indicated above, the gamma distribution has two parameters k and θ, which can be expressed in terms of means and standard deviations using the intrinsic properties of the gamma function, namely:
hence
The gamma distribution also has the welcome property of its mean being the same as that of the empirical distribution. To shift a gamma distribution, the mean and standard deviation of the shifted distribution must be known. However, per capita consumption only indicates the mean of the (up)shifted distribution. To derive the standard deviation of this distribution, we had to find a way to predict the standard deviation based on the mean. This was achieved via regression from the large global dataset.
To shift the gamma functions to fit the adult per capita level, two crucial assumptions had to be made:
1. The proportion of abstainer categories as derived from the survey reflected the true proportion of current abstainers (lifetime abstainers plus exdrinkers) in the population.
2. The overall coverage rate for the survey (i.e., the total volume of alcohol exposure derived from the survey divided by the adult per capita consumption from sales or other statistics) applied to all subpopulations as defined by age and sex.
We found a coverage rate of 0.529% between the NESARC survey and adult per capita alcohol consumption for the US, estimated at 8.75 liters per capita for 20012002 when NESARC took place based on the Global Information System for Alcohol and Health http://apps.who.int/globalatlas/default.asp webcite. Using this value, the shifted means of the drinkers for different subpopulations as described by age, sex, and ethnicity were derived as follows with the same constant for all subpopulations:
As described above, the shifted standard deviation was derived empirically via regression analyses [17]:
where sex was coded 0 for men and 1 for women. In the regression, N = 851 sex and age subpopulations from 66 countries, individually modeled to derive the above equation [17]. The multiple linear regression based on all these surveys explained 94% of variation of the dependent variable.
Method for deriving AAFs
Knowing the percentages of abstainers and former drinkers, as well as the risk relation, the AAFs based on continuous distributions were obtained using the following formula:
while the categorical value of the AAFs were obtained using:
where P_{abs }represents the proportion of lifetime abstainers, P_{form }the proportion of former drinkers, and P(x) the probability distribution function of drinkers. RR_{form }represents the relative risk for former drinkers, and RR(x) the relative risk function for a given alcohol consumption in grams per day. The subscript i denotes the groups as characterized by different categories for volume of drinking. We conducted a sensitivity analysis with consumption capped at 150 grams of pure alcohol per day.
Confidence intervals were based on simulations using the bootstrapping method [32]. For each AAF, 10,000 simulations were run.
Results
Additional File 1 gives an overview of the volume of drinking for different groups as defined by sex, age, and ethnicity. Overall, as expected, the younger the age group among adults, the higher the volume. Men consumed more than women, and Native Americans consumed on average more than other ethnicities.
Additional file 1. Description of key parameters on volume of alcohol exposure from NESARC (weighted to be representative for the population).
Format: PDF Size: 131KB Download file
This file can be viewed with: Adobe Acrobat Reader
Table 1 and Figures 1 and 2 show the fitting of different distributions to the data for nonHispanic whites as an example. All three distributions fit the data reasonably well, with some deviations mainly in the tails of the distributions. For consumption levels up to 100 grams, the Weibull and gamma distributions were very similar and fit well, but the lognormal distribution underestimated drinking at lower levels of drinking.
Table 1. Example for fitting exposure distributions for white nonHispanic Americans
Figure 1. Histogram of alcohol exposure and fitted distributions for nonHispanic white men.
Figure 2. Histogram of alcohol exposure and fitted distributions for nonHispanic white women.
Table 2 gives an overview of all the chisquare deviations for the different subgroups as defined by sex, age, and ethnicity. Clearly, the Weibull distribution fits best, but there are no options for easily shifting this distribution [17]. However, for descriptive purposes only, the Weibull distribution seems to best fit the various distributions. For shifting distributions, the gamma distribution was chosen mainly for feasibility. However, it offered a relatively good fit for descriptive purposes as well. Overall, for several subgroups, gamma had a better fit than lognormal, whereas for others, lognormal showed the better fit (Table 2).
Table 2. Fit indices for different distributions
Table 3 lists the parameters of the original and the shifted gamma distributions. Again, the unshifted and shifted distributions for nonHispanic whites are given for illustration in Figures 3 and 4.
Table 3. Parameter estimates for original and shifted gamma distributions
Figure 3. Original and shifted gamma distributions for nonHispanic white men.
Figure 4. Original and shifted gamma distributions for nonHispanic white women.
Finally, the implications of the distributional shift for the AAFs of liver cirrhosis are displayed in Table 4. Clearly, the shift in distribution results in changes in AAFs that are public healthrelevant. On the other hand, Table 4 also shows that the different distributions did not have a marked influence on AAFs, further justifying the choice of the gamma distribution.
Table 4. AAFs (in %) for liver cirrhosis for subpopulations as defined by sex and ethnicity based on different distributions
Furthermore, the differences between the continuous and the categorical approach were not pronounced. The results for the confidence intervals are listed in Table 5.
Table 5. AAFs (in %) for liver cirrhosis and corresponding 95% confidence intervals for subpopulations as defined by sex and ethnicity based on different distributions
Discussion
A procedure for triangulating survey and per capita data for deriving population exposure based on the gamma distribution for drinkers was presented and explored. This procedure proved feasible for modeling US drinking and generating AAFs based on the per capita consumption. It also allowed for the quantitative comparability of data on alcohol exposure from surveys with different coverage rates and a necessary correction for varying coverage rates.
One crucial assumption made in the triangulation process was a constant factor of underreporting for all subpopulations as defined by sex, age, and ethnicity (see formula for the shifted mean above). There is no conclusive literature on differential underreporting by different subpopulations; even the literature on underreporting by volume of drinking is not conclusive [11]. It may be hypothesized that the more irregular the occasions of heavy drinking, the more underreporting there is as such occasions are difficult to report with our standard instruments. However, research is necessary to test this hypothesis before differential upshifting factors can be used for subpopulations. The methodological framework presented here certainly allows for such modifications.
One argument against upshifting the distribution of alcohol consumed to the level of per capita data is that the risk relations usually are also derived from selfreports on alcohol exposure, i.e., that the subject's responses in epidemiological studies also underestimate real drinking. There are, however, several counterarguments. First, alcohol exposure measurement in medical epidemiological studies and in general appears to yield valid individual consumption levels [33,34]. Second, as typically found in medical epidemiological studies [33], there are higher intercorrelations with external standards when alcohol is embedded into a series of other food items. Third, it has been shown that embedded alcohol items yield higher levels of consumption compared to questionnaires specifically targeting alcohol use [35]. Thus, there are indications that questions on alcohol in medical epidemiological studies yield more consistent and higher alcohol exposure compared to those in typical national alcohol surveys. However, the degree of difference between these approaches is not clear.
There are three principal ways to explain the undercoverage of per capita consumption in surveys:
• Measurement error in surveys due to sampling;
• Measurement error in surveys due to respondents' behavior, such as underreporting, problems in averaging, forgetting, or dropping out of the survey [11];
• Measurement error in per capita consumption.
It is beyond the scope of this article to give a full review on the types of measurement error in alcohol surveys. The literature is diverse and highly speculative, and there are few systematic studies on reasons for undercoverage. In highincome countries such as the US, there are groups left out of sampling frames who are very high alcohol consumers, such as the homeless or those living in institutions. Consumption is skewed, with a small segment of the population consuming a high proportion of alcohol. In the NESARC sample, 6.7% of the heaviest white male drinkers consume 33% of the overall consumption; in the upshifted distribution for the same group, 10.2% consume this proportion [36]. It is thus possible that a large part of the undercoverage is due to sampling schemes, and hence there might be less systematic underreporting among respondents than it appears. However, as laid out above, it is unlikely that sampling explains all of the undercoverage. There is some evidence for individual underreporting as well [11].
The third explanation is that the assumption that the sales/production figures are the gold standard might be wrong. This explanation seems implausible except in circumstances where there is large unrecorded consumption. Where there is much unrecorded consumption, including that in the sales/production estimate would of course increase the differences with estimates from surveys.
It is unlikely that there is much alcohol measured in the sales/production figures that is not in fact consumed. Why would consumers pay for goods that would be wasted, e.g., as in the case of Canada cited above, a wastage of 60% to 70%? Such behavior would be contradictory to evidence from economics. Clearly, some alcohol bought is not consumed due to spillage. However, according to industry experts, this spillage should amount to less than 10%. Other alcohol may also be stocked rather than consumed in the year of purchase, but overall across regions and years, this should cancel out. Thus, the assumption of per capita consumption derived from aggregate statistics such as sales and/or production being the best estimate for overall volume of consumption seems justified. Exploring the factors involved in undercoverage by surveys should be a research priority in coming years. It seems wasteful to conduct hundreds of annual surveys with questions about alcohol in highincome countries without knowledge of why these surveys typically cover only proportions between 30% to 60% of per capita consumption. The assumption made in the present analysis of constant undercoverage in different population segments also needs to be tested and replaced by empirical estimates of differential undercoverage.
However, irrespective of reasons for undercoverage, triangulation is necessary for the comparison of alcohol exposure. If surveys are compared to each other, they should have their proportions of coverage standardized in the same way that disease rates are standardized to correct for population distributions. It would also appear irrational to continue the comparisons of surveys in which the underlying coverage rates starkly differ as the results of such comparisons are not interpretable. This reasoning is independent of the level of upshifting chosen. Based on the uncertainty about the degree of underreporting in medical epidemiological studies, we suggest the routine application of sensitivity analyses using 100%, 90%, and 80% of per capita consumption as the target levels when the actual level of population consumption is important. This method will be used not only in the CRA 2005, but, together with new empirically determined disability weights, also in the ongoing US Burden of Disease study [37].
A final point concerns the assumptions made in the upshifting. First, we assumed the proportion of abstainers and exdrinkers to be exactly as assessed by the survey. Unfortunately, there is no better information available on which different estimates could be based. The medical epidemiological literature is of no help here, as cohorts get selected based on their potential to be followed up, and this may yield some proportions of abstainers that are not at all representative of the general population. Second, we assumed that proportions by sex and age are correctly estimated by the survey. Again, there are no better data currently available. We can only speculate if and how the inclusion of nonhousehold members shifts the proportions of alcohol consumed. More research is necessary on such populations to estimate the bias introduced by relying only on households in the sampling frame.
Conclusion
Overall, the chosen methods for estimating alcohol exposure for population health proved feasible and seemed justified based on current knowledge. Further work is needed for refinement of methods and validation of assumptions.
Conflict of interest
The authors declare that they have no competing interests.
Authors' contributions
JR supervised all aspects of the work and wrote the first draft. TK programmed most of the statistical tests and contributed to the writing. GG Jr contributed to the programming, did part of the statistical analyses, and contributed to the writing. FS and BG did part of the statistical analyses and contributed to the writing. GG Sr. contributed to the programming, helped design the study, and contributed to the writing. All authors read and approved the final manuscript.
Acknowledgements
NIAAA (contract # HHSN267200700041C "Alcohol and DrugAttributable Burden of Disease and Injury in the US" to the first author), the Global Burden of Disease the and Injury 2005 Project provided financial and/or technical support for this study. With regard to the contributions by Dr. Rehm, support to the Centre for Addiction and Mental Health for the salary of scientists and infrastructure has been provided by the Ontario Ministry of Health and Long Term Care. The views expressed in this paper do not necessarily reflect those of the Ministry of Health and Long Term Care or other funders. We would like to thank Robin Room and Fotis Kanteres for helpful comments on earlier versions of the manuscript.
References

Rehm J, Room R, Graham K, Monteiro M, Gmel G, Sempos C: The relationship of average volume of alcohol consumption and patterns of drinking to burden of disease  an overview.
Addiction 2003, 98:12091228. PubMed Abstract  Publisher Full Text

Rehm J, Mathers C, Popova S, Thavorncharoensap M, Teerawattananon Y, Patra J: Global burden of disease and injury and economic cost attributable to alcohol use and alcohol use disorders.
Lancet 2009, 373:22232233. PubMed Abstract  Publisher Full Text

Rehm J, Baliunas D, Borges GLG, Graham K, Irving HM, Kehoe T, Parry CD, Patra J, Popova L, Poznyak V, Roerecke M, Room R, Samokhvalov AV, Taylor B: The relation between different dimensions of alcohol consumption and burden of disease  an overview.

Rothman KJ, Greenland S, Lash TL: Modern epidemiology. 3rd edition. Lippincott Williams & Wilkins, PA, USA; 2008.

Hanley JA: A heuristic approach to the formulas for population attributable fraction.
J Epidemiol Community Health 2001, 55:508514. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Rehm J, Taylor B, Patra J, Gmel G: Avoidable burden of disease: conceptual and methodological issues in substance abuse epidemiology.
Int J Methods Psychiatr Res 2006, 15:181191. PubMed Abstract  Publisher Full Text

Lieber CS: Biochemical and molecular basis of alcoholinduced injury to liver and other tissues.
N Engl J Med 1988, 319:16391650. PubMed Abstract

Becker U: Epidemiology and risk factors in alcohol liver disease. In Comprehensive handbook of alcohol related pathology. Edited by Preedy VR, Watson RR. Elsevier Academic Press, London, UK; 2005:467480.

Rehm J, Taylor B, Mohapatra S, Irving H, Baliunas D, Patra J, Roerecke M: Alcohol as a risk factor for liver cirrhosis  a systematic review and metaanalysis.

Walter SD: The estimation and interpretation of attributable risk in health research.
Biometrics 1976, 32:829849. PubMed Abstract  Publisher Full Text

Rehm J, Rehn N, Room R, Monteiro M, Gmel G, Jernigan D, Frick U: The global distribution of average volume of alcohol consumption and patterns of drinking.
Eur Addict Res 2003, 9:147156. PubMed Abstract  Publisher Full Text

Canadian Centre on Substance Abuse: Canadian Addiction Survey 2004: Microdata eGuide. Canadian Centre on Substance Abuse, Ottawa; 2004.

Rehm J, Patra J, Popova S: Alcoholattributable mortality and potential years of life lost in Canada 2001: Implications for prevention and policy.
Addiction 2006, 101:373384. PubMed Abstract  Publisher Full Text

Rehm J, Klotsche J, Patra J: Comparative quantification of alcohol exposure as risk factor for global burden of disease.
Int J Methods Psychiatr Res 2007, 16:6676. PubMed Abstract  Publisher Full Text

English D, Holman C, Milne E, Winter M, Hulse G, Codde G, Bower G, Corti B, de Klerk N, Knuiman M, Kurinczuk J, Lewin G, Ryan G: The quantification of drug caused morbidity and mortality in Australia 1995. Commonwealth Department of Human Services and Health, Canberra, Australia; 1995.

Kehoe T, Gmel G, Gmel G, Rehm J: Fitting different distributions to alcohol consumption among drinkers. CAMH, Toronto, ON; 2009.

Grant B, Moore T, Kaplan K: Source and accuracy statement: Wave 1 National Epidemiologic Survey on Alcoohol and Related Conditions (NESARC). National Institute on Alcohol Abuse and Alcoholism, Bethesda, MD; 2003.

Grant BF, Dawson DA, Stinson FS, Chou PS, Kay W, Pickering R: The Alcohol Use Disorder and Associated Disabilities Interview ScheduleIV (AUDADISIV): reliability of alcohol consumption, tobacco use, family history of depression and psychiatric diagnostic modules in a general population sample.
Drug Alcohol Depend 2003, 71:716. PubMed Abstract  Publisher Full Text

Ledermann S: Alcool, Alcoolisme, Alcoolisation. Volume I. Presses Universitaires de France, Paris, France; 1956.

Ledermann S: Alcool, Alcoolisme, Alcoolisation. Volume II. Presses Universitaires de France, Paris, France; 1964.

Skog OJ: The collectivity of drinking cultures: A theory of the distribution of alcohol consumption.
Br J Addict 1985, 80:8399. PubMed Abstract  Publisher Full Text

Duffy JC: The distribution of alcohol consumption  30 years on.
Br J Addict 1986, 81:735741. PubMed Abstract  Publisher Full Text

Skog OJ: The distribution of alcohol consumption. Part I. A critical discussion of the Ledermann Model. National Institute for Alcohol Research, Oslo, Norway; 1982.

Guttorp P, Hiang H: A note on the distribution of alcohol consumption.

Skog OJ: A note on the distribution of alcohol consumption; Gamma vs Lognormal distributions. A reply to Guttorp and Song.

Skog OJ: The tail of the alcohol consumption distribution.
Addiction 1993, 88:601610. PubMed Abstract  Publisher Full Text

Gruenewald P, Nephew T: Drinking in California: Theoretical and empirical analyses of alcohol consumption patterns.
Addiction 1994, 89:707723. PubMed Abstract  Publisher Full Text

Alanko T: Statistical models for estimating the distribution function of alcohol consumption; a parametric approach.
The Finnish Foundation for Alcohol Studies, Helsinki 1997.
ISBN: 9519192611.

R Development Core Team: R: A language and environment for statistical computing. [http://www.Rproject.org] webcite
R Foundation for Statistical Computing, Vienna, Austria;
ISBN 3900051070.

Fleiss J, Levin B, Cho Paik M: Statistical Methods for Rates and Proportions. John Wiley & Sons, Hoboken, New Jersey; 2003.

Kleijnen JPC: Design and analysis of simulation experiments. Springer Verlag, New York, NY; 2008.

Rehm J: Measuring alcohol consumption: how about adopting usual epidemiological standards.

Feunekes G, van Veer P, van Staveren W, Kok F: Alcohol intake assessment: the sober facts.
Am J Epidemiol 1999, 150:105112. PubMed Abstract  Publisher Full Text

King AC: Enhancing the selfreport of alcohol consumption in the community: two questionnaire formats.
Am J Public Health 1994, 84:294296. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Greenfield TK, Rogers JD: Who drinks most of the alcohol in the US? The policy implications.
J Stud Alcohol 1999, 60:7889. PubMed Abstract  Publisher Full Text

Rehm J, Frick U: Valuation of health states in the U.S. study to establish disability weights: lessons from the literature.