 Research article
 Open Access
 Published:
A metaregression analysis of 41 Australian problem gambling prevalence estimates and their relationship to total spending on electronic gaming machines
BMC Public Health volume 17, Article number: 495 (2017)
Abstract
Background
Many jurisdictions regularly conduct surveys to estimate the prevalence of problem gambling in their adult populations. However, the comparison of such estimates is problematic due to methodological variations between studies. Total consumption theory suggests that an association between mean electronic gaming machine (EGM) and casino gambling losses and problem gambling prevalence estimates may exist. If this is the case, then changes in EGM losses may be used as a proxy indicator for changes in problem gambling prevalence. To test for this association this study examines the relationship between aggregated losses on electronic gaming machines (EGMs) and problem gambling prevalence estimates for Australian states and territories between 1994 and 2016.
Methods
A Bayesian metaregression analysis of 41 crosssectional problem gambling prevalence estimates was undertaken using EGM gambling losses, year of survey and methodological variations as predictor variables. General population studies of adults in Australian states and territory published before 1 July 2016 were considered in scope. 41 studies were identified, with a total of 267,367 participants. Problem gambling prevalence, moderaterisk problem gambling prevalence, problem gambling screen, administration mode and frequency threshold were extracted from surveys. Administrative data on EGM and casino gambling loss data were extracted from government reports and expressed as the proportion of household disposable income lost.
Results
Money lost on EGMs is correlated with problem gambling prevalence. An increase of 1% of household disposable income lost on EGMs and in casinos was associated with problem gambling prevalence estimates that were 1.33 times higher [95% credible interval 1.04, 1.71]. There was no clear association between EGM losses and moderaterisk problem gambling prevalence estimates. Moderaterisk problem gambling prevalence estimates were not explained by the models (I ^{2} ≥ 0.97; R ^{2} ≤ 0.01).
Conclusions
The present study adds to the weight of evidence that EGM losses are associated with the prevalence of problem gambling. No patterns were evident among moderaterisk problem gambling prevalence estimates, suggesting that this measure is either subject to pronounced measurement error or lacks construct validity. The high degree of residual heterogeneity raises questions about the validity of comparing problem gambling prevalence estimates, even after adjusting for methodological variations between studies.
Background
Introduction and rationale
Total consumption theory, or single distribution theory as it is sometimes known, predicts that the incidence of gamblingrelated harm is related to the amount of time and money spent on gambling within a given jurisdiction [1, 2]. This prediction derives from a postulate of total consumption theory that, at the population level, a fixed proportion of total gambling activity will result in harm. Any growth in gambling will accordingly produce a proportionate growth in harm. If this hypothesis is correct, it implies that gamblingrelated harm would be best prevented by reducing the gambling consumption of the entire population, not just those gambling to excess. A related implication is that changes in mean gambling consumption may be used as a proxy indicator for changes in problem gambling prevalence.
A small research literature has investigated the veracity of the propositions of total consumption theory as they relate to gambling. In a pioneering analysis of household expenditure surveys, Grun and McKeigue [3] found that mean gambling losses in British geographic regions were strongly correlated with the proportion of the population losing an ‘excessive’ sum of money on gambling, both before and after the introduction of the National Lottery. Lund [2], analysing three independent Norwegian samples, found similar correlations between average gambling frequency and the proportion of the population gambling very frequently.
These studies did not specifically examine gamblingrelated harms, a shortcoming addressed by both Hansen and Rossow [4] and Markham et al. [5]. Hansen and Rossow examined problem gambling among Norwegian adolescents grouped by school, and found a strong correlation among schools between average losses on slot machines and the reported prevalence of problemgambling symptoms. Similarly, Markham et al. found that the reported prevalence of two or more problemgambling symptoms among gamblers using electronic gaming machines (EGMs) in Australian venues was correlated with the average amount lost on EGMs in those venues.
Few studies, however, have systematically examined the relationship between gamblingrelated harm and gambling losses at the spatial scale of the regulatory jurisdiction (e.g. the country, state, territory, province, Bundesland, etc.). The jurisdictional spatial scale is important since jurisdictions compromise the territorial unit at which gambling is most frequently regulated, gambling losses are usually reported e.g. [6, 7], and problem gambling surveys are usually conducted [8]. In one notable example of a jurisdictionallevel study, the Productivity Commission [9] surveyed problem gambling prevalence in all Australian states and territories and compared these prevalence estimates to total nonlottery gambling losses in the same jurisdictions, finding a positive correlation. Unfortunately, this study was constrained by design to an examination of only eight prevalence estimates, limiting its generalizability.
Counterexamples to predictions of total consumption theory have been forwarded by Abbott [10], who describes the reduction in problem gambling prevalence estimates in New Zealand over a 9 years period during which total gambling losses increased substantially. Consequently, in a series of studies [10,11,12,13], Abbott proposed an alternative hypothesis of ‘adaptation’, in which the prevalence of problem gambling tends to fall over time. The reasons for adaptation may include a decline in gambling participation as the novelty of a new gambling activity dwindles, decreased average duration of gambling problems through destigmatisation and improved treatment, changing cultural norms, increased knowledge of gamblingrelated harms, and the introduction of regulations such as invenue smoking bans and caps on EGM numbers.
One explanation for the dearth of convincing evidence about the relationships between gambling harms and gambling losses at jurisdictional scales is the problem of interstudy heterogeneity, generally thought to result from a lack of methodological consistency between prevalence studies. Interjurisdictional comparisons may be compromised because problem gambling prevalence studies tend to use heterogeneous methods that limit comparability [8, 14,15,16,17]. As Sassen et al. [15] found in their systematic review of 39 studies, the decade between 2000 and 2010 saw little methodological convergence among prevalence studies. In practice, measurement differences between prevalence studies may be so great as to render comparisons between them invalid. Nevertheless, as an examination of almost any governmentcommissioned problem gambling prevalence study will demonstrate, comparisons between prevalence estimates are routinely drawn. Despite concerns regarding validity, almost every problem gambling prevalence study seeks to benchmark prevalence estimates against those within the same jurisdiction at a previous point in time, or within other jurisdictions at a similar point in time.
The problems inherent in comparing prevalence estimates have been recognised by some scholars, who have attempted to regularise these prevalence rates to account for methodological variations, e.g. [8, 13, 18, 19]. However, the validity of comparing regularised estimates has not yet been established because the amount of residual heterogeneity among studies after adjustment is unknown. This is important because problem gambling prevalence estimates are the primary means through which gamblingrelated harm is monitored by regulators and governments. If prevalence estimates cannot be meaningfully compared, this calls into question validity of the current, routine practice of monitoring problem gambling prevalence using surveys [16].
If the total consumption theory of gambling is correct, then monitoring total gambling losses might provide an alternative means to track the changing incidence of gamblingrelated harm. If gambling losses present an accurate and precise proxy measure for problem gambling prevalence, then the necessity to routinely conduct problem gambling prevalence estimates to monitor populationlevel rates of harm might be reduced. Instead, populationlevel gambling losses could be monitored as a proxy indicator for the incidence of harm in the population.
Objectives
This study analyses the association between problem gambling prevalence estimates and gambling losses for Australian states and territories between 1994 and 2015. It aims to answer the following specific research questions:

a)
Is there an association between EGM and casino gambling losses and problem gambling prevalence estimates in Australian states and territories?

b)
What degree of heterogeneity remains in estimates of problem gambling prevalence after regularising for methodological variations, EGM gambling losses and year of survey?
Methods
A metaregression approach was used to estimate the association between EGM and casino gambling losses and problem gambling prevalence estimates for Australian states and territories.
Setting
The units of analysis for this study were the eight states and territories of Australia. EGMs were introduced to these jurisdictions in a staggered manner, with New South Wales the first to legalise EGMs in 1956 [20]. This was followed by a wave of legalisations, mostly in the 1990s, which left Western Australia as the only jurisdiction without EGMs by 1997. Other legal gambling commodities that are available in all jurisdictions include lotteries, casino table games, instant lotteries, scratch cards, and betting on races, sports and special events.
Australia was selected as a study site because its eight federal states and territories pioneered the routine conduct of problem gambling prevalence studies, alongside Canada and the United States [8]. Consequently, there have been sufficient problem gambling prevalence studies conducted in Australia to warrant a metaanalysis of their results. The study was limited to Australian states and territories rather than including jurisdictions in multiple countries to minimise the potential differences among the populations surveyed.
Data
Two sets of data were required:

1.
problem gambling prevalence estimates and the characteristics of the studies which produced these estimates, and

2.
EGM gambling losses in the state or territory that temporally match each prevalence study.
Problem gambling prevalence studies for Australian states and territories were identified through a systematic search process. Prevalence studies are most frequently published as reports in the ‘grey literature’ rather than as peerreviewed journal articles. Consequently, the search strategy primarily involved the identification of this relatively wellknown corpus of problem gambling prevalence studies from previous inventories [8, 21]. The websites of Australian government bodies that have commissioned problem gambling prevalence studies were searched to identify further studies for examination, as were the reference lists of identified studies. The search revealed one study for which the full text was unavailable. The lead author of this study was contacted by email and conducted data extraction at the current authors’ request. To be eligible for inclusion, prevalence studies had to: 1) target the general population aged 18 years or older, 2) measure 12 or 6months problem gambling using a validated problem gambling screen, 3) report results for one or more whole states or territories in Australia, 4) have been published prior to 1 July 2016, and 5) report on independent samples rather than longitudinal studies measuring change among the same respondents over time. Data were extracted independently by two coders, MS and FM. In cases where these two coders disagreed, data were coded independently by a third coder MY and a consensus meeting held, as discussed by Orwin and Vevea [22].
EGM gambling losses were selected as the predictor variable of interest in preference to total gambling losses because expenditure on different forms of gambling produces differing levels of harm. Previous research has shown that EGM and casino losses are more closely associated with problem gambling than either total gambling losses or losses on other gambling products [23, 24]. Gambling losses on casino table games were also included in the study as Australian player loss statistics for casinos do not distinguish between EGM and nonEGM gambling losses [7]. Gambling losses on EGMs during the year of survey fieldwork were extracted from Australian Gambling Statistics, 32nd Edition, a complete and authoritative administrative dataset [7]. This dataset is compiled by Queensland Treasury, on the basis of aggregate tax records provided by each Australian state or territory government.
Measures
The following measures were extracted from each problem gambling prevalence study: a) the prevalence of ‘problem gambling’, where problem gambling was defined as Problem Gambling Severity Index (PGSI) ≥ 8 or South Oaks Gambling Screen (SOGS) ≥ 5, b) the prevalence of ‘moderaterisk problem gambling’, where moderaterisk problem gambling was defined as PGSI 3 – 7 or SOGS 3 – 4, c) the jurisdiction, d) the year during which data was collected, e) the administration mode of the survey, i.e. telephone or facetoface, f) whether the SOGS [25] or the PGSI [26] was used to assess problem gambling, g) the sample size of the survey, and h) any ‘frequency threshold’ used to select which respondents would be administered the problem gambling screen. A frequency threshold is a rule by which the problem gambling screening instrument is only administered to a subset of respondents, selected on the basis of their reported gambling frequency. For example, the screen may only be administered to those who gamble at least weekly, with less frequent gamblers being imputed a problem gambling score of zero. The prevalence of problem gambling and the prevalence of moderaterisk problem gambling were the key outcome variables of interest.
Problem gambling screen and administration mode variables were coded as dummy variables. The gambling frequency threshold variable was collapsed into four categories: weekly gambling, fortnightly gambling, monthly gambling and less than monthly gambling. The year the survey was conducted was subtracted from 2015 to calculate the age of the reported survey.
The measure of EGM gambling losses used was the sum of EGM gambling losses in hotels and clubs and all gambling losses in casinos, both expressed as a percentage of total household disposable income (HDI). This is an aggregate measure, with a single number reported for each state and territory in each year. It is calculated by expressing the total gambling expenditure on EGMs and in casinos for the jurisdiction as a percentage of HDI for the jurisdiction, where HDI is the total income accruing to the household sector less household sector taxes. HDI for each state and territory is reported annually in the Australian System of National Accounts [27]. Because gambling losses were reported for fiscal years (spanning 1 July – 30 June) while survey dates were recorded for calendar years, losses for each calendar year were estimated by calculating the mean of losses for the two overlapping fiscal years.
Statistical analysis
Random effects metaregression was used to estimate the partial correlation between problem gambling prevalence and EGM and casino losses, after adjusting for methodological variations in prevalence studies and the year of the survey. Metaregression is an extension of metaanalysis, but has different aims. In general, the goal of a metaanalysis is to pool the varying results of primary studies and thereby arrive at a more accurate and precise estimate of a quantity of interest (e.g. the population prevalence of problem gambling). In contrast, a metaregression analysis aims to understand what causes variation in the findings of primary studies, using procedures developed for regression analysis [28]. In this context, we might interpret this study as contributing to an ‘epidemiology of problem gambling prevalence studies’ [29].
A random effects model is a statistical extension to fixed effects metaregression. While a fixed effects analysis relies on the assumption that the quantity of interest (e.g. the population prevalence of problem gambling) is truly consistent across all studies even if it is imperfectly measured, this is rarely the case. For example, the populations under study are unlikely to be identical in all relevant factors, even after covariates are adjusted for. A random effects specification is more conservative because it does not make this strong assumption. When applied, random effects metaregression usually results in estimates with wider confidence intervals than fixed effects metaregression [28].
The adopted statistical approach modelled the estimated prevalence of problem gambling and moderaterisk problem gambling in each study as a function of the following predictor variables: EGM gambling losses in the jurisdiction; problem gambling screen; administration mode; frequency threshold; and survey year. Variance inflation factors for predictor variables were all less than 4.0. A binomial model specification with a logistic link function was used. Following Higgins and Thompson and Borenstein and colleagues [28, 30], the random effects metaregression model was specified as follows:
where: y _{ i } is the number of problem gamblers identified in study i; ssize _{ i } is the sample size of study i; α is a constant intercept; X _{ ik }is a matrix of k predictor variables and β _{ k } is a commensurate vector of estimated regression coefficients including EGM gambling losses, study year and methodological variations; θ _{ i } is a normally distributed random effect with a standard deviation between studies of τ (or equivalently a precision of \( \frac{1}{\tau^2} \)); n is the number of studies under analysis; and I ^{2} is the proportion of residual variation in the estimates of problem gambling prevalence that is due to heterogeneity between studies (rather than sampling variation). The uniform prior distribution ranging from 0.0 to 10.0 for τ was specified on the basis of the simulations carried out by Lambert and colleagues [31]. All models were implemented using ‘gold standard’ [30] fully Bayesian estimation with R and JAGS [32, 33]. JAGS code listings for all model specifications are available in Additional file 1 (see listings 1 – 6).
The estimates of β coefficients were modelled in three different ways in order to test the robustness of results to the provision of prior information about their values. In the first model, all β coefficients were estimated from the prevalence study data set in the usual manner of regression analysis, using ‘weakly informative’ priors distributions. This method has the advantage of minimising residual heterogeneity, but risks identifying spurious correlations [34]. In the second model, coefficients – except for the intercept, the EGM loss coefficient and the year of survey coefficient – were ‘fixed’ on the basis of estimates derived from the small literature concerned with their estimation. This is equivalent to the recent practice that has been applied to compare prevalence estimates between studies e.g. [8, 13], where prevalence estimates are normalised by multiplying them by fixed adjustment factors. Fixing coefficients is advantageous as it forces control variables to be set at plausible values and reduces the effective degrees of freedom of the models. However, it admits no variance in coefficient estimates. Consequently, a third method which constitutes a compromise between the first two was also adopted, assigning ‘informative’ prior distributions to control variable parameter coefficients, regularising them within plausible ranges while still admitting uncertainty in their estimated values.
Prior distributions for β coefficients in all three models are listed in Table 1. Weakly informative priors were placed on the β coefficients for HDI loss and year of survey in all models. The standard deviations for the informative priors were inflated by a factor of four compared to those derived from metaanalysis, on the basis that the small number of studies synthesised in the metaanalyses were likely to lead to an overestimation of the precision of these distributions. The unpublished metaanalyses that formed the basis of the informative prior distributions are attached as Additional file 2. Each metaregression model was estimated twice, first with the prevalence of problem gambling as the outcome variable, and second with the prevalence of moderaterisk problem gambling as the outcome variable.
These models address the two research objectives of this study. First, the player loss β coefficient can be interpreted as an indicator of the magnitude and direction of any association between aggregate EGM gambling losses and the prevalence of problem gambling. Second, the I ^{2} estimate describes the degree of heterogeneity that remains among problem gambling prevalence estimates after accounting for both sampling variability within individual prevalence studies and the predictor variables described above. In other words, I ^{2} measures the inconstancy in prevalence estimates between studies rather than real variation in the prevalence of problem gambling. I ^{2} has a range of 0.0 to 1.0, with Higgins and colleagues suggesting that the values of 0.25, 0.5, and 0.75 indicate low, moderate and high degrees of heterogeneity respectively [35].
Results
Selection
A total of 45 problem gambling prevalence estimates were identified, including eight statelevel estimates derived the Productivity Commission’s 1999 national survey [see Additional file 3, Fig. S1 and Table S1]. Once ineligible studies were excluded, 41 studies were identified that estimated the prevalence of problem gambling, 40 of which also estimated the prevalence of moderaterisk problem gambling. All extracted data for each study are listed in Table 2.
Just over half of studies identified used the PGSI (n = 21), with the remainder using SOGS. The vast majority were administered by telephone (n = 37). The most commonly used gambling frequency thresholds for the administration of a problem gambling screen were a weekly threshold (n = 19) and the combined category of an annual threshold, a sixmonthly threshold or no threshold at all (n = 19). Only two studies used a monthly threshold, while a single study used a fortnightly threshold. The median sample size was 4303 (Interquartile range [IQR] = 1253 – 9408), with a total of 267,367 adults responding to the 41 surveys. The median survey year was 2001 (IQR = 1999 – 2008).
Outcome data and main results
The average nonregularised prevalence of problem gambling across all studies was 0.9% of adults [95% Credible Interval (Cr.I.) 0.8%, 1.1%], with the average nonregularised prevalence of moderaterisk problem gambling estimated to be 1.8% of adults (95% Cr.I. 1.5%, 2.1%) (see Additional file 3, Fig. S2 and S3). There was an extraordinarily large degree of heterogeneity among these studies, with I ^{2} for problem gambling and moderaterisk problem gambling estimated at 0.95 (95% Cr.I. 0.94, 0.95) and 0.94 (95% Cr.I. 0.93, 0.96) respectively.
An association between the prevalence of problem gambling and EGM and casino gambling losses was apparent in the metaregression model with weakly informative priors. Parameter coefficients, expressed as ‘prevalence ratios’ are displayed in Table 3. Prevalence ratios should be interpreted analogously to incidence rate ratios, and can be multiplied with the intercept to predict the value of the outcome variable for a given set of predictor variable values. Every increase of 1 % of household disposable income lost on EGMs and at casinos was associated with problem gambling prevalence estimates that were 1.35 (95% Cr.I. 1.04, 1.74) times higher. Placing informative priors on other metaregression coefficients decreased the parameter estimate slightly to 1.33 (95% Cr.I. 1.04 1.71). Fixing coefficients related to methodological variations in prevalence studies slightly reduced the estimated association between prevalence and EGM and casino gambling losses and decreased the precision of the estimate, bringing the ‘no association’ prevalence ratio of 1.0 to within the 95% credible interval (1.29, 95% Cr.I. 0.98, 1.72). Posterior estimates of the associations between prevalence and losses are visualised in Fig. 1.
Few clear associations were found between moderaterisk problem gambling and any predictor variables using the models with weakly informative priors. With the use of informative priors, only the estimated association between moderaterisk problem gambling prevalence and a weekly gambling frequency threshold became reasonably precise (0.68, 95% Cr.I. 0.50, 0.94).
The number of years the study was conducted before 2015 was positively associated with both problem gambling prevalence and moderaterisk problem gambling prevalence, but only when metaregression coefficients were fixed at the values arrived at from metaanalyses of previous withinstudy estimates. Models with fixed priors found that a 1 year increase in the age of the study was associated with prevalence estimates that were 1.04 (95% Cr.I. 1.02, 1.07) times greater for problem gambling and 1.05 (95% Cr.I. 1.02, 1.09) times greater for moderaterisk problem gambling.
A great deal of residual heterogeneity was evident among all models. While the models of problem gambling prevalence explained up to 66% of the variation among estimates, I ^{2} for these data fell in the range between 0.78 and 0.95. This means that even after adjusting for covariates, a great deal of variation remained among problem gambling prevalence estimates, with unexplained heterogeneity between studies dominating random sampling error. This unsatisfactory situation was more extreme for estimates of moderaterisk problem gambling. No model explained even 1% of the variation among moderaterisk problem gambling prevalence estimates. The lower bound of estimates of I ^{2} for models of moderaterisk problem gambling was 0.95.
Discussion
The study had three key findings. First, problem gambling prevalence was associated with EGM and casino gambling losses in models with informative and weakly informative priors. An increase of 1% of household disposable income spend on EGMs and casino gambling associated with prevalence estimates that were approximately 1.3 times greater. In models where control parameter coefficients were fixed at values derived from metaanalyses, the point estimates of the gambling loss coefficient were similar but the 95% credible intervals widened to include a prevalence ratio of 1.0. In short, these results support the total consumption theory of gambling, but should be interpreted cautiously, given the degree of uncertainty evident in estimates.
The relatively wide uncertainty interval surrounding this finding is unsurprising given the relatively modest number of studies (n = 41) and their high degree of heterogeneity with respect to problem gambling prevalence (I ^{2} > 0.85 in all cases). Consequently, this study provides only a moderate degree of confidence that EGM gambling losses and problem gambling prevalence estimates are correlated. The bestfitting model, which used informative priors, suggested that an increase in EGM gambling losses of 1% of HDI is associated with a populationlevel increase in problem gambling prevalence of around 1.33 times (95% Cr.I. 1.04 – 1.71). The width of this uncertainty band is likely to be a consequence of measurement error overwhelming true variation in problem gambling prevalence. This high level of statistical noise in measurement is likely to derive from methodological variations that were unaccounted for in this study. For example, we did not adjust for way the survey was described to potential respondents, a factor that can impact nonresponse bias [36]. Future studies using more consistentlycollected data could seek to measure the relationship between EGM losses and problem gambling prevalence with more precision. These results support the need to phase out statebased prevalence studies and transition to national problem gambling prevalence studies that are adequately powered to investigate individual jurisdictions and that remain methodologically stable over time [16].
The moderate degree of uncertainty remaining around the association between EGM losses and harm should be interpreted in the context of parallel findings at other spatial scales. EGM gambling losses are correlated with risk of developing gambling problems for individuals [23] and for populations aggregated by county, school or gambling venue [4, 5]. The present study adds to the weight of evidence that an increase in population losses on EGMs is associated with an increase in the prevalence of problem gambling.
The second key finding of the study was that a high degree of heterogeneity exists in problem gambling prevalence estimates. Only a moderate degree of variation among prevalence estimates was explained by EGM gambling losses, methodological variations, or year of study (R ^{2} ≤ 0.66). When coefficient values for methodological variations were fixed based on prior research, variance explained fell substantially (R ^{2} = 0.47). Very little of the residual variation between problem gambling prevalence estimates was due to sampling error, with a very high degree of unexplained heterogeneity (I ^{2} ≥ 0.86). Put differently, after adjusting for methodological variations, EGM losses and year of survey, no more than 14% of the residual differences between problem gambling prevalence estimates results from sampling error. These results raise questions about the validity of comparing problem gambling prevalence estimates, even after adjusting for methodological variations between studies.
The third key finding was the absence of any apparent pattern among moderaterisk problem gambling prevalence estimates. Contrary to expectations premised on total consumption theory, no link was evident between moderaterisk problem gambling prevalence estimates and EGM and casino gambling losses. Furthermore, no model explained any meaningful amount of the variation in moderaterisk problem gambling prevalence estimates, with the point estimate of R ^{2} falling below 0.01 in all cases. In other words, none of the models explained even 1% of the variation in estimates of moderaterisk problem gambling.
The lack of an apparent relationship between EGM gambling losses and moderaterisk problem gambling prevalence has two potential interpretations. It could be that there is no real relationship between moderaterisk problem gambling and EGM spending, in contradiction to total consumption theory. Alternatively, it is also plausible that the PGSI and SOGS are mismeasuring the population at moderate risk of problem gambling. Given that problem gambling screens have been developed and validated to identify problem or pathological gamblers [26], it may be that problem gambling screens are not fit for the purpose of identifying moderaterisk problem gamblers.
Several pieces of evidence offer tentative support for the mismeasurement conjecture. First, the only published validation study of the moderaterisk classification of the PGSI found that it lacked discriminant validity [37]. In particular, this study found no practical differences between those scoring 12 on the PGSI and those scoring 37. Indeed, as McCready and Adlaf [38] note, the PGSI does not include any items designed to discriminate among gamblers with less severe problems. As the authors of the PGSI acknowledge in their original study, the screen’s division between low and moderaterisk gambling categories is only tentatively supported by the survey data from which it was derived [26]. In addition, the explanatory power of the variables included in the metaregression analysis for predicting moderaterisk problem gambling prevalence estimates was exceptionally poor, with R ^{2} < 0.01 in all three models (Table 3). Similarly, residual heterogeneity was very high, with I ^{2} ≥ 0.97 across all models. This implies that either moderaterisk problem gambling prevalence estimates are not impacted by methodological variations, or that such impacts are very small when compared to other unaccounted for factors. Finally, the metaanalysis of problem gambling screen effects presented in Additional file 2 finds a much greater degree of heterogeneity is present in estimates of screen impacts on moderaterisk problem gambling prevalence estimates (I ^{2} = 0.91) than problem gambling prevalence estimates (I ^{2} = 0.69). In other words, the ‘within study’ estimates of the impact of methodological variations on prevalence estimates vary a great deal between studies of moderaterisk problem gambling.
Taken together, this suggests that measures of moderaterisk problem gambling are extremely imprecise, to the point of possibly being the equivalent of statistical white noise. The apparent inability of current screening instruments to reliably identify this population is particularly problematic given that recent evidence suggests that this population experience the greatest mass of gamblingrelated harms when measured in terms of disabilityadjusted life years [39].
These results are subject to several limitations. First, a great deal of variation existed among prevalence estimates after adjusting for five predictor variables. Thus, estimates of any association between prevalence and EGM losses are necessarily imprecise. Second, the 41 prevalence estimates analysed in this study may be insufficient for the assessment of population trends, especially among the moderaterisk problem gambler population. Third, other methodological variations that this study was unable to adjust for may impact on problem gambling prevalence estimates. Finally, it is possible that some individuals may have been sampled in surveys in multiple years, especially in smaller jurisdictions. This is unlikely to have a substantial effect on the results.
Conclusions
This study has three key implications. First, the finding of an association between EGM and casino gambling losses and problem gambling prevalence is consistent with total consumption theory. Therefore, interventions by jurisdictional governments that reduce total EGM gambling losses among the whole population are likely to effectively reduce the prevalence of problem gambling. This result was evident despite the imprecision and heterogeneity of these estimates. Nevertheless, replication using a large crossjurisdictional surveys that are consistent over time are required to confirm this association with a greater degree of confidence.
Second, this study demonstrates that using singlejurisdiction prevalence studies for making comparisons between jurisdictions or within the same jurisdiction over time is ineffectual [17]. Even after deploying sophisticated statistical adjustments, the high degree of residual heterogeneity evident in this study suggests that the validity of comparing problem gambling prevalence estimates may be poor. The situation is far worse for moderaterisk problem gambling. It appears that although 267,367 Australian adults have responded to 41 surveys, we still are unable to confidently compare problem gambling prevalence either between jurisdictions or over time, a situation that is likely to be replicated internationally. It will be necessary to undertake adequatelypowered, multijurisdictional prevalence studies – with survey instruments and data collection protocols that remain consistent over time – if the scientific value of problem gambling prevalence studies is to be increased in future.
Third, this study suggests that the suitability of PGSI and SOGS for estimating the population prevalence of moderaterisk problem gambling needs urgent investigation. The validity and reliability of the moderaterisk problem gambling classification is unclear. No patterns were evident among the estimates of the population prevalence of this group between studies, symptomatic of either extreme measurement errors or poor construct validity. The interpretation of moderaterisk problem gambling prevalence estimates should only be undertaken with extreme caution until a greater degree of conceptual and statistical clarity is brought to the identification of this population.
Abbreviations
 Cr. I.:

Credible Interval
 EGM:

Electronic gaming machine
 HDI:

Household disposable income
 IQR:

Interquartile range
 PGSI:

Problem gambling severity index
 SOGS:

South oaks gambling screen
References
Rose G, Day S. The population mean predicts the number of deviant individuals. BMJ. 1990;301:1031–4.
Lund I. The population mean and the proportion of frequent gamblers: is the theory of Total consumption valid for gambling? J Gambl Stud. 2008;24:247–56.
Grun L, McKeigue P. Prevalence of excessive gambling before and after introduction of a national lottery in the United Kingdom: another example of the single distribution theory. Addiction. 2000;95:959–66.
Hansen M, Rossow I. Adolescent gambling and problem gambling: does the Total consumption model apply? J Gambl Stud. 2008;24:135–49.
Markham F, Young M, Doran B. Gambling expenditure predicts harm: evidence from a venuelevel study. Addiction. 2014;109:1509–16.
Canadian Partnership for Responsible Gambling. Canadian Gambling Digest, 20132014. Toronto: Responsible Gambling Council; 2015. http://www.responsiblegambling.org/docs/defaultsource/defaultdocumentlibrary/cprg_canadiangamblingdigest_201314.pdf?sfvrsn=2 archived at http://www.webcitation.org/6p1YFiWuH. Accessed 27 Jan 2016.
Queensland Government Statistician’s Office. Australian Gambling Statistics, 1989–90 to 2014–15, 32nd edition. Brisbane: Queensland Treasury; 2016. http://www.qgso.qld.gov.au/products/reports/ausgamblingstats/ausgamblingstats32ndedn.pdf archived at http://pandora.nla.gov.au/pan/135611/201611111528/www.qgso.qld.gov.au/products/reports/ausgamblingstats/ausgamblingstats32ndedn.pdf. Accessed 5 Jan 2017.
Williams RJ, Volberg RA, Stevens RMG. The population prevalence of problem gambling: methodological influences, standardized rates, jurisdictional differences, and worldwide trends. Technical Report. Ontario Problem Gambling Research Centre; 2012. http://hdl.handle.net/10133/3068 archived at http://www.webcitation.org/6OK7CeYsW. Accessed 21 Mar 2014.
Productivity Commission. Australia’s gambling industries. Report. Productivity Commission: Canberra; 1999.
Abbott M. Do EGMs and problem gambling go together like a horse and carriage? Gambl Res. 2006;18:7–38.
Abbott MW, Volberg RA, Williams MM. Seven years on: a followup study of frequent and problem gamblers living in the community. Department of Internal Affairs: Wellington, New Zealand; 1999.
Abbott M. Disabling the public interest: gambling strategies and policies for Britain: a comment on Orford 2005. Addiction. 2005;100:1233–5.
Abbott M, Stone CA, Billi R, Yeung K. Gambling and problem gambling in Victoria, Australia: changes over 5 years. J Gambl Stud. 2016;32:1–32.
Shaffer HJ, Hall MN, Bilt JV. Estimating the prevalence of disordered gambling behavior in the United States and Canada: a research synthesis. Am J Public Health. 1999;89:1369–75.
Sassen M, Kraus L, Bühringer G. Differences in pathological gambling prevalence estimates: facts or artefacts? Int J Methods Psychiatr Res. 2011;20:e83–99.
Markham F, Young M. Commentary on Dowling et al. (2016): is it time to stop conducting problem gambling prevalence studies? Addiction. 2016;111:436–7.
Doughney J. Lies, damned lies and “problem gambling” prevalence rates: the example of Victoria, Australia. J Bus Syst Gov Ethics. 2009;2:41–54.
Jackson AC, Wynne H, Dowling NA, Tomnay JE, Thomas SA. Using the CPGI to determine problem gambling prevalence in Australia: measurement issues. Int J Ment Health Addict. 2010;8:570–82.
Stone CA, Romild U, Abbott M, Yeung K, Billi R, Volberg R. Effects of different screening and scoring thresholds on PGSI gambling risk segments. Int J Ment Health Addict. 2014;13:82–102.
Australian Institute for Gambling Research. Australian gambling comparative history and analysis. Melbourne: Victorian Casino and Gaming Authority; 1999. https://assets.justice.vic.gov.au/vcglr/resources/bb81f943d85440de8babb09d8bbd610f/australiangamblingcomparativehistory.pdf archived at http://www.webcitation.org/6cNpqq7Xn.
Productivity Commission. Gambling. Report. Productivity Commission: Canberra; 2010.
Orwin RG, Vevea JL. Evaluating coding decisions. In: Cooper H, Hedges LV, Valentine JC, editors. Handbook of research synthesis and metaanalysis. 2nd ed. New York: Russell Sage Foundation; 2009. p. 177–203.
Markham F, Young M, Doran B. The relationship between player losses and gamblingrelated harm: evidence from nationally representative crosssectional surveys in four countries. Addiction. 2016;111:320–30.
Welte JW, Barnes GM, Wieczorek WF, Tidwell MCO, Parker JC. Risk factors for pathological gambling. Addict Behav. 2004;29:323–35.
Lesieur HR, Blume SB. The South Oaks gambling screen (SOGS): a new instrument for the identification of problem gamblers. Am J Psychiatry. 1987;144:1184–8.
Ferris J, Wynne H. The Canadian Problem Gambling Index: Final Report. Canadian Centre on Substance Abuse; 2001. http://www.ccgr.ca/en/projects/resources/CPGIFinalReportEnglish.pdf archived at http://www.webcitation.org/66rzvX4ZE. Accessed 10 Nov 2011.
Australian Bureau of Statistics. Australian System of National Accounts, 201415. 5204.0. Canberra: Australian Bureau of Statistics; 2016. www.abs.gov.au/ausstats/abs@.nsf/mf/5204.0. Accessed 10 Feb 2016.
Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. Introduction to metaanalysis. Chichester, U.K: Wiley; 2009.
Thompson SG, Higgins JPT. How should metaregression analyses be undertaken and interpreted? Stat Med. 2002;21:1559–73.
Higgins JPT, Thompson SG. Quantifying heterogeneity in a metaanalysis. Stat Med. 2002;21:1539–58.
Lambert PC, Sutton AJ, Burton PR, Abrams KR, Jones DR. How vague is vague? A simulation study of the impact of the use of vague prior distributions in MCMC using WinBUGS. Stat Med. 2005;24:2401–28.
R Core Team. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2015. http://www.Rproject.org.
Plummer M. JAGS Version 4.0.0 user manual. Lyon, France: International Agency for Research on Cancer, World Health Organisation; 2015. http://sourceforge.net/projects/mcmcjags/files/Manuals/4.x/archived at http://www.webcitation.org/6p1YQLhhG. Accessed 12 Dec 2016.
Higgins JPT, Thompson SG. Controlling the risk of spurious findings from metaregression. Stat Med. 2004;23:1663–82.
Higgins JPT, Thompson SG, Deeks JJ, Altman DG. Measuring inconsistency in metaanalyses. BMJ. 2003;327:557–60.
Williams RJ, Volberg RA. Impact of survey description, administration format, and exclusionary criteria on population prevalence rates of problem gambling. Int Gambl Stud. 2009;9:101–17.
Currie SR, Hodgins DC, Casey DM. Validity of the problem gambling Severity Index interpretive categories. J Gambl Stud. 2013;29:311–27.
McCready J, Adlaf E. Performance and enhancement of the Canadian Problem Gambling Index: Report and recommendations. Canadian Consortium for Gambling Research; 2006. http://www.ccgr.ca/en/projects/resources/PerformanceandEnhancementoftheCanadianProblemGamblingIndexCPGIReportandRecommendations.pdf archived at http://www.webcitation.org/6p1Ypb1Cg. Accessed 7 Feb 2016.
Browne M, Langham E, Rawat V, Greer N, Li E, Rose J, et al. Assessing gamblingrelated harm in Victoria: a public health perspective. Melbourne: Victorian Responsible Gambling Foundation; 2016. https://www.responsiblegambling.vic.gov.au/informationandresources/research/recentresearch/assessinggamblingrelatedharminvictoriaapublichealthperspective archived at http://www.webcitation.org/6mMlvRWkV. Accessed 29 Nov 2016.
Acknowledgements
The authors thank Paul Delfabbro for his assistance in extracting data from an unpublished report.
Funding
No specific funding was provided for this study.
Availability of data and materials
All data generated or analysed during this study are included in this published article and its supplementary information files.
Authors’ contributions
FM conceived of this study. All authors contributed to the design of the data collection instrument. FM and MS conducted systematic literature searches and extracted data from eligible studies. MY contributed to data extraction where FM and MS disagreed about coding. FM led the design of the study, with contributions from MY and BD. FM designed and undertook data analysis. FM led the interpretation of the results, with contributions from MY and BD. FM drafted the manuscript. All authors contributed to critically revising the manuscript. All authors approved the manuscript for publication.
Competing interests
FM was supported by an Australian Postgraduate Award from the Department of Education and Training. The first three authors have received funding from, or been employed on, projects funded by the Australian Research Council, the Community Benefit Fund (Northern Territory Government), and the Northern Territory Research and Innovation Fund (Northern Territory Government). The Community Benefit Fund is a government administered grant scheme, funded by a levy on EGMs in hotels in the Northern Territory. MY has received funding from the Department of Justice (Victorian Government) and Gambling Research Australia. FM has been employed on research projects funded by the Australian Capital Territory Gambling and Racing Commission. FM has received funding from the Alberta Gambling Research Institute (AGRI), a consortium of researchers at three universities, to present a keynote at an international conference. AGRI’s funding is provided by the Ministry of Treasury Board and Finance of Alberta. MY is a member of the National Association for Gambling Studies (Australia). FM is a member of the Public Health Association of Australia. None of the authors have ever received funding from the gambling, alcohol or tobacco industries, nor are there any constraints on the publication of this manuscript.
Consent for publication
Not applicable.
Ethics approval and consent to participate
No formal ethical approval was required for this study as it involved only secondary analysis of existing data.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Authors and Affiliations
Corresponding author
Additional files
Additional file 1:
JAGS models. JAGS code listings for all six of the models presented in this paper.c. (PDF 513 kb)
Additional file 2:
Metaanalyses of previous studies of methodological variations. An unpublished manuscript which seeks to estimate the average impact of the following methodological variations of prevalence estimates: (a) choice of problem gambling screen, (b) survey administration mode, or (c) choice of frequency threshold. (PDF 902 kb)
Additional file 3:
Supplementary figures and tables. Study inclusion flow diagram; problem gambling prevalence forest plot; moderate risk problem gambling prevalence forest plot; and full bibliographic details for each study. (PDF 1173 kb)
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
Cite this article
Markham, F., Young, M., Doran, B. et al. A metaregression analysis of 41 Australian problem gambling prevalence estimates and their relationship to total spending on electronic gaming machines. BMC Public Health 17, 495 (2017). https://doi.org/10.1186/s1288901744136
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1288901744136
Keywords
 Gambling
 Total consumption theory
 Problem gambling
 Problem gambling prevalence
 Electronic gaming machines
 Moderaterisk problem gambling
 Gambling expenditure
 Gambling losses