Next Article in Journal
Aflatoxin Contamination in Food and Body Fluids in Relation to Malnutrition and Cancer Status in Cameroon
Next Article in Special Issue
A Systematic Review of Occupational Exposure to Particulate Matter and Cardiovascular Disease
Previous Article in Journal
Maternal Cigarette Smoking during Pregnancy and Offspring Externalizing Behavioral Problems: A Propensity Score Matching Analysis
Previous Article in Special Issue
Associations between Area-Level Unemployment, Body Mass Index, and Risk Factors for Cardiovascular Disease in an Urban Area
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Prevalence of Coronary Heart Disease for Small Areas Using Collateral Indicators of Morbidity

Department of Geography and Centre for Statistics, Queen Mary University of London, Mile End Rd, London E1 4NS, UK
Int. J. Environ. Res. Public Health 2010, 7(1), 164-177; https://doi.org/10.3390/ijerph7010164
Submission received: 20 November 2009 / Accepted: 14 January 2010 / Published: 18 January 2010
(This article belongs to the Special Issue Cardiovascular Diseases and Public Health)

Abstract

:
Different indicators of morbidity for chronic disease may not necessarily be available at a disaggregated spatial scale (e.g., for small areas with populations under 10 thousand). Instead certain indicators may only be available at a more highly aggregated spatial scale; for example, deaths may be recorded for small areas, but disease prevalence only at a considerably higher spatial scale. Nevertheless prevalence estimates at small area level are important for assessing health need. An instance is provided by England where deaths and hospital admissions for coronary heart disease are available for small areas known as wards, but prevalence is only available for relatively large health authority areas. To estimate CHD prevalence at small area level in such a situation, a shared random effect method is proposed that pools information regarding spatial morbidity contrasts over different indicators (deaths, hospitalizations, prevalence). The shared random effect approach also incorporates differences between small areas in known risk factors (e.g., income, ethnic structure). A Poisson-multinomial equivalence may be used to ensure small area prevalence estimates sum to the known higher area total. An illustration is provided by data for London using hospital admissions and CHD deaths at ward level, together with CHD prevalence totals for considerably larger local health authority areas. The shared random effect involved a spatially correlated common factor, that accounts for clustering in latent risk factors, and also provides a summary measure of small area CHD morbidity.

1. Background

Profiling geographic variations in health care need is important for equitable and effective targeting of resources that reflects inequalities in morbidity [1]. Coronary heart disease (CHD) is an important part of the overall disease burden faced by government health agencies and demonstrates considerable geographic inequality. Assessing prevalence variations between populations and areas for chronic diseases such as CHD is a central aspect of defining health care need. However, prevalence is not necessarily as well recorded as other health outcomes (e.g., mortality, hospitalisations). For example, in some countries such as the US, area prevalence estimates can only be made on the basis of health survey data. Focussed studies, such as the British Regional Heart Study (BRHS) considered by Morris et al. [2] may also provide evidence of geographic prevalence variations, but generally provide only limited geographic coverage; thus the BRHS included only 24 British towns. In recent years, prevalence of major conditions (including CHD) treated in primary care in England has been administratively recorded under a system known as the Quality Outcomes Framework (QOF), but not at a disaggregated spatial scale (e.g., for small areas with populations under 10 thousand). However, small area contrasts are important in defining variations in health need within local health authorities.
While treated prevalence totals are only available for local health agencies (known as Primary Care Trusts or PCTs, with 152 such PCTs in England), deaths and hospital admissions for coronary heart disease are available for smaller areas known as wards, of which there are over 8000. As argued above, estimating CHD prevalence at small area level is important, and this paper develops a shared random effect (or common factor) method to pool information regarding spatial morbidity contrasts over multiple indicators (deaths, hospitalizations, prevalence). This provides a summary index for representing small area CHD morbidity which is applied to estimate small area CHD prevalence totals and hence relative prevalence risks (comparing actual to expected prevalence). Geographic variations in latent constructs relevant to population health are typically spatially correlated, and this is recognised in the derivation of the common morbidity factor [3]. The shared random effect approach also incorporates differences between areas in deprivation levels and other forms of population risk (e.g., ethnic structure). That is, the common factor is partly predicted on the basis of known ecological risk factors or ”multiple causes”, so providing a spatial adaptation of a multiple indicator-multiple cause approach [4].
An illustration is provided by data for London. London may be disaggregated into 625 small areas known as wards, and into 31 Primary Care Trust areas. Some observed data on CHD are at ward level, namely hospital admissions and mortality totals. However, some data (namely CHD prevalence totals) are only observed for considerably larger PCTs. We wish to obtain a summary index of CHD morbidity (as a shared random effect) using all the observations, and use this index to estimate the disaggregated prevalence totals for the 625 wards.
A fully Bayesian approach is used in specifying the model and in the London case study application. This involves ascribing prior densities to model parameters and updating those densities via the likelihood for the observed data. Iterative Monte Carlo Markov Chain (MCMC) techniques [5] are used for estimation, as implemented in the WINBUGS program [6].

2. Modelling Latent Morbidity at the Lower Spatial Scale

Let j = 1, .., NL denote the set of lower level small areas within a particular region, and let i = 1, .., NH denote the set of aggregated higher level areas (e.g., local health authorities) within which the small areas are nested. The available data contain P observed indicators yj = (yj1, .., yjP) at the small area scale (such as small area death totals), and counts Zi = (Zi1, .., ZiQ) (e.g., disease prevalence totals) observed only at the aggregated area scale. However, one aim of the modelling process is to develop estimates zj = (zj1, .., zjQ) of these indicators at a small area scale.
It is assumed that correlations between the observed indicators can be represented by underlying common latent factors f = (f1, .., fR), where R is of typically of much smaller dimension than the total number P + Q of observed indicators. For simplicity, a univariate common factor f = (f1, ...fNL) is considered (i.e. R = 1). In the parlance of factor analysis techniques, the set of observed indicators are proxies for, or ”measures of”, the underlying latent factor.
The first set of small area measurement equations describe the relationship between the observations yjp (j = 1, .., NL, p = 1, …, P) and the latent factor. In population health applications, the indicators are typically discrete counts (e.g., deaths, hospital admissions), assumed either Poisson or binomial, so that a general linear mixed model is appropriate for the measurement equations. In the application here, mortality or admission is infrequent in relation to population at risk, and Poisson sampling is relevant. Expected mortality or admission counts Ojp are obtained by applying a standard age-sex schedule (for the entire region, providing an internal standard, or for the nation, providing an external standard) to small area populations at risk. Then one has
y j p P o ( μ j p ) μ j p = O j p ρ j p ,
where ρjp is the relative risk of outcome p in small area j. In the present application, expectations Ojp are scaled to equal the total of expected counts over all small areas, namely j y j p = j O j p, so that the region wide average relative risk ρp for indicator p is 1 if an internal standard is used.
As is conventional for Poisson responses, a log link is employed [7]. So one has measurement models for small area indicators p = 1, .., P,
l o g ( μ j p ) = log ( O j p ) + λ p f j + u j p ,
where the unique errors u j p N ( 0 , σ u 2 ) may be necessary for explaining any residual overdispersion. In substantive terms, the ujp also control for structural influences unrelated to population morbidity per se (e.g., effectiveness of health care services, hospital configuration). Intercepts are not included in (2), so providing a form of location constraint on the latent variable f [8]. The coefficients λp are typically known as loadings, the specification of which is considered below.
Variations in population morbidity, whether observed or latent, typically display spatial correlation between adjacent areas-unmeasured aspects of population structure relevant to health risk typically straddle administrative boundaries [9]. However, rather than a priori assume exclusively spatial dependence, the model here determines an appropriate mix between spatial and non-spatial (”exchangeable”) dependence in the latent morbidity construct.
There may also be observed variables (i.e., known rather than latent risk factors) that are relevant to defining the common morbidity factor. For example, many indices of health need are composites of variables such as unemployment rates, poverty rates, car ownership, etc. Here a spatial adaptation of the multiple indicators-multiple causes (MIMIC) approach is used, with L measured causes xj = (xj1, .., xjL) (such as small area socio-economic or population risk variables) of the latent morbidity index. These influence the latent morbidity index fj via regression terms
η j = β x j = β 1 x j 1 + + β L x j L
where the regression excludes an intercept, with residuals denoted
r j = f j η j .
Additionally if the xjl are standardised, the absolute size of the β coefficients measures the relative importance of different population risk factors or socio-economic variables in defining the morbidity index.
To allow a mix between spatial and non-spatial dependence in the latent index, define a spatial correlation parameter k ∈ (0, 1), and assume symmetric spatial interactions wjh. Also let f[j] = (f1,fj−1, fj+1, …, fNL) denote the collection of morbidity effects for all areas but area j. Under the scheme of Leroux et al. [10], though adapted here to include regression effects, as in (3.1) (3.2), the expected value of the latent effect in area j and its variance are
E ( f j | f [ j ] ) = η j + κ h j w j h r h / [ 1 κ + κ h j w j h ] ,
V a r ( f j | f [ j ] ) = σ f 2 [ 1 κ + κ h j w j h ] ,
where σ f 2 is a variance parameter. A value of k close to 1 indicates high spatial dependence in latent morbidity, while values near zero imply lack of spatial correlation.
The wjh may incorporate factors such as distances between areas j and h. However, in many applications the wjh simply represent adjacency, namely wjh = whj = 1 if areas h and j are adjacent, zero otherwise. In this case it is relevant to define the neighbourhood j of small area j, which contains the mj areas adjacent to area j, and one then has h j w j h = m j. The expectations are then
E ( f j | f [ j ] ) = η j + κ h j r h / [ 1 κ + κ m j ] .
To uniquely determine the scale of the f scores, constraints are needed on the loadings λp, or on the variance σ f 2 in (4.2). The first kind involves standardized factors, with σ f 2 = 1, as in the spatial factor model of Wang and Wall [11], with all loadings then unknowns. An alternative constraint involves appropriately fixed loadings, such as setting one of the loadings λp to a particular fixed value, usually 1. The variance σ f 2 is then an unknown parameter.

3. Methods: Estimating Prevalence at Small Area Level based on the Morbidity Index

We wish not just to obtain a latent morbidity index, but to use this index to estimate unknown indicator totals (zj1, …, zjQ) (e.g., prevalence totals) for small areas j = 1, …, NL. Estimation of the missing lower area scale data takes account (a) of values of the small area morbidity index f = (f1, …, fNL), and (b) of the known prevalence totals (Zi1, …, ZiQ) for the i = 1, …, NH higher level areas. The small areas are nested within the higher level areas, with Hj ∈ {1, ,NH} denoting the higher level area to which small areas j belong, and the region is defined equivalently by all the higher level areas or all the lower level areas.
Also assumed known are age-sex structures for the small area populations, and from these can be obtained expected totals Ejq of the small area counts zjq. This involves using an external schedule of prevalence rates rqsk for the qth outcome by age k and sex s, and applying this schedule to small area population estimates Pjsk, so that E j q = s κ P j s κ r q s κ.
To ensure the estimates of (zj1,…, zjQ) take account of the observed prevalence counts Ziq of the higher level areas they are located in, the Poisson means Δiq in the likelihood ZiqPo(Δiq) for the higher level observed totals Ziq are defined by totals of small area means δjq located within each higher area. Thus let
Δ i q = H j = i δ j q ,
denote the total mean prevalence counts for large areas i obtained from the small area model for the z-indicators.
The small area model (i.e., the model for the δjq) can be set up to ensure that the posterior means of the Δiq equal (to a close approximation) the known higher level totals Ziq. One way to achieve thus is via a collection of NH fixed effects γq,Hj in the model for the δjq, equivalent to using dummy variables in the small area model for each higher scale area, and providing a Poisson equivalence to the multinomial [12]. Thus the zjq for Hj = i are multinomial within Ziq. We also wish the values of the latent morbidity index fj to influence the multinomial allocation of Ziq to small areas in a manner analogous to that in Equation (2). So the small area model is
z j q P o ( δ j q ) ,
l o g ( δ j q ) = log ( E j q ) + γ q , H j + λ P + q f j ,
where λP+q are additional loadings on the latent spatial morbidity index. Whether they are set to known values or taken as unknowns depends on the identification constraint adopted for the scale of the fj.
Other priors for (γ1, …,γNH), for example, as random rather than fixed effects, in practice have a very similar consequence. that the means of the Δiq equal (to a close approximation) the known higher level total Ziq. For example, one might use random effect spatial priors, comparable to (4)–(5) but at the higher area level.
In some circumstances, there may be doubts about how far the Ziq are accurate measures of morbidity, and a constraint to reproduce them may not be advantageous. For example, the prevalence counts obtained under the QOF system in England may under-record prevalence in deprived areas, since the quality of primary care is lower in such areas [13], this may result in less effective case-finding [14]. To allow unconstrained estimation of small area prevalence counts, one may use an intercept in the model for δjq that is not specific to the higher area, namely
l o g ( δ j q ) = log ( E j q ) + γ q + λ P + q f j .
This model ensures i = 1 N H Z i q = j = 1 N L z j q, but does not guarantee that H j = i z j q = Z i q, as the constrained model does.
To recap, the model is a form of spatial structural equation model (SEM) that seeks to estimate small area health outcomes z for which only large area observations Z are available. The model works in practice by using observed small area health indicators y (e.g., mortality from a particular disease) which are likely to be correlated with the missing small area outcomes z (e.g., prevalence of the same disease). The information in the correlated multiple indicators y is summarised in a latent variable f that depends on observed area risk factors x, but is also spatially correlated, reflecting spatial clustering in unobserved area risk factors. The decomposition of large area totals Z to small areas is based on the latent variable f, and the decomposition can be constrained so that total small area prevalences zj sum to the known Zi for large area i within which areas j are located. It seems reasonable to use socioeconomic variables x as causes of variability in f, but another strategy would be to use small area socioeconomic variables as additional indicators of the latent variable.

4. CHD Morbidity in London: Data

The motivating case study illustrating the above methodology involves derivation of a univariate index of CHD morbidity (i.e. R = 1) for London small areas using P = 4 observed small area health indicators, and a single health indicator (Q = 1) observed only at an aggregated area scale. The two area scales are wards and Primary Care Trusts (PCTs): there are NL = 625 wards and NH = 31 PCTs in London. The first two small area indicators (yj1, yj2) are male and female CHD deaths over 2004–2006, while (yj3, yj4) are male and female hospitalisations for CHD over three financial years 2003–2004 to 2005–2006. Expected deaths and hospitalisations Ojp in (2) are based on London wide death and hospitalisation rates specific to gender and five year age bands.
CHD prevalence totals Zi (for 2004–2005 and 2005–2006 combined) are observed only for PCTs, but one goal of the model is to estimate missing small area CHD prevalence totals zj. Expected CHD prevalence totals Ejq = Ej at ward level in (7.2) are obtained with an external schedule of CHD prevalence rates by age and sex, and applying this schedule to small area population estimates (here 2005 intercensal estimates of ward populations developed by the UK Office of National Statistics). The external schedule used is based on the 2003 Health Survey for England [15], with the expectations Ej scaled so that the London wide standard prevalence ratio is 1 (i.e., the total of observed prevalence counts Zi across all London PCTs equals the total of expected prevalence counts Ej over all London wards).
In the multiple cause sub-model (3), there are L = 3 socio-economic indicators of population CHD risk. These are x1 = average weekly household income in 2001–2002 [16], x2 = proportion of population of south Asian ethnicity, 2001 Census [17,18], and x3 = estimated ward level smoking prevalences [19]. These predictors are converted to standardised form so that their relative importance can be assessed.

5. CHD Morbidity in London: Models

Two models are compared. One assumes intercepts in the small area prevalence model that vary by PCT, as in (7.2). The other is unconstrained, as in (8). Identifiability is achieved by setting λP+1 = λ5 = 1, so that σ f 2 is an unknown, the inverse variance 1 / σ f 2 is accordingly assigned a Gamma(1,1) prior. To ensure the model produces a positive index of CHD morbidity, the remaining λp parameters also follow Gamma(1,1) priors [20]. Fixed effect parameters, namely β parameters in (3) and γ parameters in (7.2) and (8) are assigned diffuse N(0, 100) priors, while a uniform prior kU(0, 1) is assumed for the spatial correlation coefficient in (4)–(5).
Comparisons of model fit use the deviance information criterion (DIC) [21], obtained as the average deviance plus a complexity measure. The focus is on goodness of fit for the y-indicators (deaths and hospital admissions). Model 1 will automatically fit the higher level (PCT) prevalence data better as it has separate intercepts for each PCT. Model checking is based on the posterior predictive density, p(yrep|y), under a mixed predictive approach [22], where sampled replicates yrep are based on model means that include replicate samples from random effects (f and u effects). Then a mixed-predictive test for area j and outcome p has the form
p j p , m i x = P r ( y r e p , j p > y j p | y ) + 0.5 P r ( y r e p , j p = y j p | y ) ,
with extreme tail values indicating poorly fitted cases. One may compare the proportion of cases under-fitted (pjp,mix < 0:05) or over-fitted (pjp,mix > 0:95) with the expected proportions in these two tails (namely 0.05 in each).
Inferences are based on the second halves of two chain runs of 10000 iterations with convergence before iteration 5000 assessed using Gelman-Rubin scale reduction factors [23]. Table 1 presents model fit and checking criteria. It can be seen that model 1 (ward totals constrained to reproduce the QOF totals at PCT level) has a lower DIC, and satisfactory predictive performance. Table 2 summarises parameter estimates under the two models.
The estimated k coefficients in Table 2 indicate a high spatial correlation in the latent CHD morbidity index under both models. The estimated βl parameters from the multiple cause regression (3) show income differences between wards to be the most important known influence on the index, though concentrations of south Asian ethnic groups are also important. As expected, higher income levels are negatively associated with morbidity (so the 95% interval for the coefficient β1 is confined to negative values). The importance of area socioeconomic status to CHD outcomes is confirmed by other studies [24,25].
The income effect is weaker in the constrained model. This is likely to reflect discrepancies in some deprived parts of London between officially recorded prevalence (used as a constraint in model 1), and what would be expected on the basis of socioeconomic structure. Examples are the apparently low prevalence in some deprived areas in inner South East London. The consequence is that the effect of income is deflated, providing an example of measurement error affecting regression estimates. Table 3 compares prevalence (in the higher level PCT areas) based on the official QOF totals, with average income levels in such areas (weekly income in hundreds of pounds). Outliers in the negative relation between prevalence and income (there is a-0.50 correlation between PCT ranks for prevalence risk and for income, even though using official CHD prevalence data) include the deprived inner SE London area of Southwark. The latter area has the sixth lowest income, but also low measured prevalence. So while the DIC criterion prefers model 1, the geographic prevalence pattern implied by the unconstrained model 2 might be preferred on the basis of epidemiological arguments.
Figure 1 shows the spatial patterning of the CHD morbidity scores in model 1, higher values are in inner east (though not central) London and in certain parts of west London. Figure 2 maps the estimated ward level prevalences in terms of relative risks under the constrained model 1, namely the posterior means of ξj = zj/Ej. For policy purposes, the probability that a small area has significantly higher relative risk and thus possibly needs special resources is important. Therefore the marginal variance ω2 = var(ξj) is monitored during the MCMC run, and the standardized relative risks (SRRs)
S R R j = ( ξ j ξ ¯ ) / ω
are also monitored. High posterior means for these SRRs (e.g., SRRs above 1 or 2, namely 1 or 2 standard deviations above average) indicate significantly elevated prevalence, while low values (under −1, or under −2) indicate significantly low prevalence. Figure 3 maps the three categories: SRRj > 1, −1 < SRRj < 1 and SRRj < −1. Clusters of elevated risk are now clearly apparent.

6. Discussion

Estimates of prevalence at small area level are often necessary, as prevalence is less likely to be routinely reported for such areas, whereas outcomes such as mortality and hospitalization often are. Prevalence totals may, however, be reported for relatively aggregated large areas, either from health surveys, or (in the case of the UK) systems of chronic disease monitoring in primary care. The present paper has employed a common spatial factor model to disaggregate large area CHD prevalence totals to small areas. Various forms of common spatial factor model have been proposed for spatial health outcomes epidemiology. For example, Hogan & Tchernis [26] develop a small area measurement model for a socioeconomic deprivation score, while Liu et al. [27] develop a spatial structural equation model linking health outcomes to spatially correlated latent indices.
The present paper includes three main extensions on such work: first, it allows for known risk factors to influence the composite index via a multiple cause sub-model, second, it allows the data to determine the extent of spatial correlation rather than presuming a priori that latent factors are necessarily spatially dependent, and third, it applies the model to estimate missing health outcomes at a lower area scale (CHD prevalence for wards in the London case study), when observations on such outcomes are only available at a higher scale (PCTs in the London study). The essence of the method is to use all available lower scale information (both from levels of related health outcomes y and from measures of socioeconomic structure x) to provide a reasonable imputation of the missing outcomes z at the lower scale.
The case study has considered deaths and hospital admissions for CHD as the lower scale observed data (the y-variables), and a single higher scale outcome (CHD prevalence), with x-variables (causes) being income, ethnicity and smoking. It has also had a primarily urban focus, being confined to London. Under a broader geographic focus (including rural small areas), it might be relevant to consider adding an urban-rural indicator to the x-variables.
As demonstrated in Figure 3, one application of the modelling scheme is to highlight small areas with significantly elevated prevalence. This is important for prioritizing resourcing and intervention, and is based on a method that seeks to make use of all relevant information (comorbid outcomes, area social structure, and the spatial configuration of small areas). By contrast, many other health needs measures used to distribute resources are based simply on socioeconomic variables (e.g., the Jarman score) [28], or on regressions of single health indicators (e.g., hospitalizations) on socioeconomic variables [29], when multiple indicators may in fact be relevant. Existing methods also neglect spatial clustering in unobserved risk factors.
The methodology set out here has potential application to small area prevalence estimation for other chronic diseases, though the appropriate mix of y and x-variables would be different. For example, the Quality Outcomes Framework system in the UK monitors prevalence of several chronic diseases. In particular, PCT (higher scale) level counts of the prevalence of serious mental illness (SMI) are available, but one may seek ward level measures of SMI prevalence. The available y-indicators in this situation might be small area hospital admissions for conditions such as schizophrenia and bipolar disorder, while x-variables would include indicators of risk for psychiatric morbidity, such as small area income or deprivation, urban-rural status, social capital and so on.
Another potential application area is to use health survey information on disease prevalence, often obtained only for higher scale regional units. For example, the public release version of recent Health Surveys for England only contains prevalence rate estimates for chronic conditions included in the survey (e.g., obesity, diabetes) for 10 Strategic Health Authorities. However, one may wish to use this information in making estimates of such conditions for lower scale geographies such as the 354 local authorities in England. Using survey based regional estimates Zi of prevalence, one can estimate lower scale totals zj, using information on both socioeconomic structure (xj) and related outcomes (yj) at the lower spatial scale. The procedures outlined in the paper could in fact be used to disaggregate survey based estimates Zik which include relevant demographic stratifiers k (e.g. age, sex, ethnicity). Relevant spatial SEM coefficients (β and λ parameters) may well differ between demographic category. For example, one might seek to disaggregate survey-based regional estimates of diabetes by ethnicity to a lower spatial scale.

References

  1. Sundquist, K; Malmström, M; Johansson, S; Sundquist, J. Care need index, a useful tool for the distribution of primary health care resources. J. Epid. Comm. Health 2003, 57, 347–352. [Google Scholar]
  2. Morris, R; Whincup, P; Lampe, F; Walker, M; Wannamethee, S; Shaper, A. Geographic variation in incidence of coronary heart disease in Britain: the contribution of established risk factors. Heart 2001, 86, 277–283. [Google Scholar]
  3. Hogan, J; Tchernis, R. Bayesian factor analysis for spatially correlated data, with application to summarizing area-level material deprivation from census data. J. Amer. Stat. Assoc 2004, 99, 314–324. [Google Scholar]
  4. Kline, R. Principles and practice of structural equation modeling; Guilford Press: New York, NY, USA, 2004. [Google Scholar]
  5. Gelfand, A; Smith, A. Sampling based approaches to calculate marginal densities. J. Amer. Statist. Assoc 1990, 85, 398–409. [Google Scholar]
  6. Lunn, D; Spiegelhalter, D; Thomas, A; Best, N. The BUGS project: evolution, critique and future directions. Stat. Med 2009, 28, 3049–3067. [Google Scholar]
  7. McCullagh, P; Nelder, J. Generalized Linear Models; Chapman and Hall/CRC: New York, NY, USA, 1989. [Google Scholar]
  8. Best, N; Hansell, A. Geographic variations in risk: adjusting for unmeasured confounders through joint modeling of multiple diseases. Epidemiology 2009, 20, 400–410. [Google Scholar]
  9. Rezaeian, M; Dunn, G; St Leger, S; Appleby, L. Geographical epidemiology, spatial analysis and geographical information systems: a multidisciplinary glossary. J. Epidemiol. Community Health 2007, 61, 98–102. [Google Scholar]
  10. Leroux, B; Lei, X; Breslow, N. Estimation of disease rates in small areas: a new mixed model for spatial dependence. In Statistical Models in Epidemiology, the Environment and Clinical Trials; Halloran, M, Berry, D, Eds.; Springer-Verlag: New York, NY, USA, 1999; pp. 135–178. [Google Scholar]
  11. Wang, F; Wall, M. Generalized common spatial factor model. Biostatistics 2003, 4, 569–582. [Google Scholar]
  12. Forster, J. Bayesian Inference for Poisson and Multinomial Log-Linear Models; Southampton Statistical Sciences Research Institute: Southampton, UK; (S3RI Methodology Working Papers, M09/11), 2009. [Google Scholar]
  13. Wright, J; Martin, D; Cockings, S; Polack, C. Overall QOF scores lower in practices in deprived areas. Br. J. Gen. Prac 2006, 56, 277–279. [Google Scholar]
  14. Sigfrid, L; Turner, C; Crook, D; Ray, S. Using the UK primary care Quality and Outcomes Framework to audit health care equity: preliminary data on diabetes management. J. Publ. Health 2006, 28, 221–225. [Google Scholar]
  15. Strong, M; Maheswaran, R; Radford, J. Socioeconomic deprivation, coronary heart disease prevalence and quality of care: a practice-level analysis in Rotherham using data from the new UK general practitioner Quality and Outcomes Framework. J. Publ. Health 2006, 28, 39–42. [Google Scholar]
  16. Heady, P; Clarke, P; Brown, G; Ellis, K; Heasman, D; Hennell, S; Longhurst, J; Mitchell, B. Model-Based Small Area Estimation; Office for National Statistics: London, UK, 2003. [Google Scholar]
  17. Forouhi, N; Sattar, N; Tillin, T; McKeigue, P; Chaturvedi, N. Do known risk factors explain the higher coronary heart disease mortality in south Asian compared with European men? Prospective follow-up of the Southall and Brent studies, UK. Diabetologia 2006, 49, 2580–2588. [Google Scholar]
  18. Tziomalos, K; Weerasinghe, C; Mikhailidis, D; Seifalian, A. Vascular risk factors in South Asians. Int. J. Cardiol 2008, 128, 5–16. [Google Scholar]
  19. Scarborough, P; Allender, S; Rayner, M; Goldacre, M. Validation of model-based estimates (synthetic estimates) of the prevalence of risk factors for coronary heart disease for wards in England. Health and Place 2009, 15, 596–605. [Google Scholar]
  20. Sahu, S. Bayesian estimation and model choice in item response models. J. Stat. Comp. Sim 2002, 72, 217–232. [Google Scholar]
  21. Spiegelhalter, D; Best, N; Carlin, B; van der Linde, A. Bayesian measures of model complexity and fit. J. Roy. Stat. Soc. B 2002, 64, 583–639. [Google Scholar]
  22. Marshall, C; Spiegelhalter, D. Identifying outliers in Bayesian hierarchical models: a simulation-based approach. Bayesian Analysis 2007, 2, 1–33. [Google Scholar]
  23. Brooks, S; Gelman, A. Alternative methods for monitoring convergence of iterative simulations. J. Comp. Graph. Stat 1998, 7, 434–456. [Google Scholar]
  24. Winkleby, M; Sundquist, K; Cubbin, C. Inequities in CHD incidence and case fatality by neighborhood deprivation. Am. J. Prev. Med 2007, 32, 97–106. [Google Scholar]
  25. Sundquist, K; Winkleby, M; Ahlén, H; Johansson, S. Neighborhood socioeconomic environment and incidence of coronary heart disease: a follow-up study of 25,319 women and men in Sweden. Am. J. Epid 2004, 159, 655–662. [Google Scholar]
  26. Hogan, J; Tchernis, R. Bayesian factor analysis for spatially correlated data, with application to summarizing area-level material deprivation from census data. J. Amer. Stat. Assoc 2004, 99, 314–324. [Google Scholar]
  27. Liu, X; Wall, M; Hodges, J. Generalized spatial structural equation modeling. Biostatistics 2005, 6, 539–557. [Google Scholar]
  28. Jarman, B. Underprivileged areas: validation and distribution of scores. Brit. Med. J 1984, 289, 1587–1592. [Google Scholar]
  29. Glover, G; Robin, E; Emami, J; Arabscheibani, G. A needs index for mental health care. Soc Psychiatry Psychiatr Epidemiol 1998, 33, 89–96. [Google Scholar]
Figure 1. Latent Morbidity Index.
Figure 1. Latent Morbidity Index.
Ijerph 07 00164f1
Figure 2. Relative prevalence risk.
Figure 2. Relative prevalence risk.
Ijerph 07 00164f2
Figure 3. Significantly elevated and depressed risk.
Figure 3. Significantly elevated and depressed risk.
Ijerph 07 00164f3
Table 1. Fit measures.
Table 1. Fit measures.
Average DevianceComplexityDICProportion of y values with Pr(yrep > y) under 0.05 or over 0.95
Model 1 (Multinomial Constraint)2,5701,2903,8600.093
Model 2 (unconstrained)2,5291,5124,0410.085
Table 2. Summary of model parameter estimates.
Table 2. Summary of model parameter estimates.
MeanStdevnMonte Carlo SE2.5%97.5%
Model 1β1−0.1990.0200.002−0.242−0.162
β20.1800.0230.0020.1390.231
β30.0500.0200.0020.0090.087
κ0.9370.0550.0020.8000.998
λ10.4500.0390.0030.3710.525
λ20.4100.0440.0030.3250.495
λ30.7200.0400.0040.6380.783
λ40.7690.0410.0040.6820.893
Model 2β1−0.3940.0320.003−0.461−0.333
β20.1580.0250.00010.1080.208
β30.0620.0300.0020.0070.120
κ0.9350.0550.0020.7970.998
λ10.3030.0240.0010.2580.352
λ20.2410.0260.0010.1910.292
λ30.3660.0200.0010.3260.406
λ40.4130.0230.0010.3660.456
Table 3. PCT prevalence risks based on official prevalence total, compared with income patterns.
Table 3. PCT prevalence risks based on official prevalence total, compared with income patterns.
PCTObserved from QOFExpected using HSE 2003 as standardRR based on actual QOF prevalence recordsRank of RRAverage incomeIncome rank
Barking and Dagenham9,8009,1471.071255.32
Barnet20,16119,2871.045217.625
Bexley13,97314,7780.946126.714
Brent14,04013,5421.037206.610
Bromley19,46620,8830.932117.523
Camden9,4309,3891.004157.321
City and Hackney9,0308,7461.033195.54
Croydon17,51918,8680.92896.815
Ealing18,41015,2281.209297.322
Enfield14,83916,2320.91456.611
Greenwich12,41911,4641.083265.85
Hammersmith and Fulham7,0227,6090.92377.827
Haringey9,3189,3600.996146.612
Harrow13,68012,9491.056227.726
Havering16,65016,5381.007166.916
Hillingdon13,92914,4080.967137.219
Hounslow8,1277,6631.061246.59
Islington13,92914,4080.967137.219
Kensington and Chelsea6,9539,5060.73118.028
Kingston8,5738,4851.010178.129
Lambeth9,76810,4990.930106.613
Lewisham12,02711,3481.060236.18
Newham12,4959,4331.325314.81
Redbridge14,48814,2231.019186.917
Richmond and Twickenham7,80210,3120.75729.031
Southwark10,23311,1680.91666.06
Sutton and Merton19,30321,3850.90347.624
Tower Hamlets9,7247,5231.293305.43
Waltham Forest11,95510,6401.124276.07
Wandsworth10,90411,7630.92788.430
Westminster9,92111,0970.89437.118

Share and Cite

MDPI and ACS Style

Congdon, P. Estimating Prevalence of Coronary Heart Disease for Small Areas Using Collateral Indicators of Morbidity. Int. J. Environ. Res. Public Health 2010, 7, 164-177. https://doi.org/10.3390/ijerph7010164

AMA Style

Congdon P. Estimating Prevalence of Coronary Heart Disease for Small Areas Using Collateral Indicators of Morbidity. International Journal of Environmental Research and Public Health. 2010; 7(1):164-177. https://doi.org/10.3390/ijerph7010164

Chicago/Turabian Style

Congdon, Peter. 2010. "Estimating Prevalence of Coronary Heart Disease for Small Areas Using Collateral Indicators of Morbidity" International Journal of Environmental Research and Public Health 7, no. 1: 164-177. https://doi.org/10.3390/ijerph7010164

APA Style

Congdon, P. (2010). Estimating Prevalence of Coronary Heart Disease for Small Areas Using Collateral Indicators of Morbidity. International Journal of Environmental Research and Public Health, 7(1), 164-177. https://doi.org/10.3390/ijerph7010164

Article Metrics

Back to TopTop