A framework for progressively improving small area population estimates

The paper presents a framework for small area population estimation that enables users to select a method that is fit for the purpose. The adjustments to input data that are needed before use are outlined, with emphasis on developing consistent time series of inputs. We show how geographical harmonization of small areas, which is crucial to comparisons over time, can be achieved. For two study regions, the East of England and Yorkshire and the Humber, the differences in output and consequences of adopting different methods are illustrated. The paper concludes with a discussion of how data, on stream since 1998, might be included in future small area estimates.

Population monitoring for kangaroo management

In wildlife management, the program of monitoring will depend on the management objective. If the objective is damage mitigation, then ideally it is damage that should be monitored. Alternatively, population size (N) can be used as a surrogate for damage, but the relationship between N and damage obviously needs to be known. If the management objective is a sustainable harvest, then the system of monitoring will depend on the harvesting strategy. In general, the harvest strategy in all states has been to offer a quota that is a constant proportion of population size. This strategy has a number of advantages over alternative strategies, including a low risk of over- or underharvest in a stochastic environment, simplicity, robustness to bias in population estimates and allowing harvest policy to be proactive rather than reactive. However, the strategy requires an estimate of absolute population size that needs to be made regularly for a fluctuating population. Trends in population size and in various harvest statistics, while of interest, are secondary. This explains the large research effort in further developing accurate estimation methods for kangaroo populations. Direct monitoring on a large scale is costly. Aerial surveys are conducted annually at best, and precision of population estimates declines with the area over which estimates are made. Management at a fine scale (temporal or spatial) therefore requires other monitoring tools. Indirect monitoring through harvest statistics and habitat models, that include rainfall or a greenness index from satellite imagery, may prove useful.

Estimating physical activity level: The role of domestic activities

Appropriate measures of physical activity are essential for determining the population prevalence of physical activity, for tracking trends over time, and for guiding intervention efforts. Physical activity measurement is characterised by the synthesis of information on the type, frequency, intensity, and duration of activity over a specified period. To date, emphasis in physical activity assessment has been on the measurement of leisure time physical activities. However, some domestic and transport related activities entail energy expenditures equivalent to moderate intensity of 3.0–6.0 METS1 considered to be of sufficient intensity to achieve a health benefit are yet to be included in routine population level physical activity surveillance. This leads to population estimates based only on measures of leisure time physical activities.

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. IV. Abundance and conservation status

Urban encroachment on dense, coastal koala populations has ensured that their management has received increasing government and public attention. The recently developed National Koala Conservation Strategy calls for maintenance of viable populations in the wild. Yet the success of this, and other, conservation initiatives is hampered by lack of reliable and generally accepted national and regional population estimates. In this paper we address this problem in a potentially large, but poorly studied, regional population in the State that is likely to have the largest wild populations. We draw on findings from previous reports in this series and apply the faecal standing-crop method (FSCM) to derive a regional estimate of more than 59 000 individuals. Validation trials in riverine communities showed that estimates of animal density obtained from the FSCM and direct observation were in close agreement. Bootstrapping and Monte Carlo simulations were used to obtain variance estimates for our population estimates in different vegetation associations across the region. The most favoured habitat was riverine vegetation, which covered only 0.9% of the region but supported 45% of the koalas. We also estimated that between 1969 and 1995 similar to 30% of the native vegetation associations that are considered as potential koala habitat were cleared, leading to a decline of perhaps 10% in koala numbers. Management of this large regional population has significant implications for the national conservation of the species: the continued viability of this population is critically dependent on the retention and management of riverine and residual vegetation communities, and future vegetation-management guidelines should be cognisant of the potential impacts of clearing even small areas of critical habitat. We also highlight eight management implications.

Differences in prevalence of pre-existing morbility between injured and non-injured populations

Objectives To identify and examine differences in pre-existing morbidity between injured and non-injured population-based cohorts. Methods Administrative health data from Manitoba, Canada, were used to select a population-based cohort of injured people and a sample of non-injured people matched on age, gender, aboriginal status and geographical location of residence at the date of injury. All individuals aged 18-64 years who had been hospitalized between 1988 and 1991 for injury (International Classification of Diseases, Ninth Edition, Clinical Modification (ICD-9-CM) code 800-995) (n = 21032), were identified from the Manitoba discharge database. The matched non-injured comparison group comprised individuals randomly selected 1: 1 from the Manitoba population registry. Morbidity data for the 12 months prior to the date of the injury were obtained by linking the two cohorts with all hospital discharge records and physician claims. Results Compared to the non-injured group, injured people had higher Charlson Comorbidity Index scores, 1.9 times higher rates of hospital admissions and 1.7 times higher rates of physician claims in the year prior to the injury. Injured people had a rate of admissions to hospital for a mental health disorder 9.3 times higher, and physician claims for a mental health disorder 3.5 times higher, than that of non-injured people. These differences were all statistically significant (P < 0.001). Conclusion Injured people were shown to differ from the general non-injured population in terms of pre-existing morbidity. Existing population estimates of the attributable burden of injury that are obtained by extrapolating from observed outcomes in samples of injured cases may overestimate the magnitude of the problem.

Estimates of sex- and age-class-specific survival probabilities for a southern Great Barrier Reef green sea turtle population

Sex- and age-class-specific survival probabilities of a southern Great Barrier Reef green sea turtle population were estimated using a capture - mark - recapture (CMR) study and a Cormack - Jolly - Seber (CJS) modelling approach. The CMR history profiles for 954 individual turtles tagged over a 9-year period ( 1984 - 1992) were classified into three age classes ( adult, subadult, juvenile) based on somatic growth and reproductive traits. Reduced-parameter CJS models, accounting for constant survival and time-specific recapture, fitted best for all age classes. There were no significant sex-specific differences in either survival or recapture probabilities for any age class. Mean annual adult survival was estimated at 0.9482 (95% CI: 0.92 - 0.98) and was significantly higher than survival for either subadults or juveniles. Mean annual subadult survival was 0.8474 ( 95% CI: 0.79 - 0.91), which was not significantly different from mean annual juvenile survival estimated at 0.8804 ( 95% CI: 0.84 - 0.93). The time-specific adult recapture probabilities were a function of sampling effort but this was not the case for either juveniles or subadults. The sampling effort effect was accounted for explicitly in the estimation of adult survival and recapture probabilities. These are the first comprehensive sex- and age-class-specific survival and recapture probability estimates for a green sea turtle population derived from a long-term CMR program.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Estimates of global and regional potential health gains from reducing multiple major risk factors

Background Estimates of the disease burden due to multiple risk factors can show the potential gain from combined preventive measures. But few such investigations have been attempted, and none on a global scale. Our aim was to estimate the potential health benefits from removal of multiple major risk factors. Methods We assessed the burden of disease and injury attributable to the joint effects of 20 selected leading risk factors in 14 epidemiological subregions of the world. We estimated population attributable fractions, defined as the proportional reduction in disease or mortality that would occur if exposure to a risk factor were reduced to an alternative level, from data for risk factor prevalence and hazard size. For every disease, we estimated joint population attributable fractions, for multiple risk factors, by age and sex, from the direct contributions of individual risk factors. To obtain the direct hazards, we reviewed publications and re-analysed cohort data to account for that part of hazard that is mediated through other risks. Results Globally, an estimated 47% of premature deaths and 39% of total disease burden in 2000 resulted from the joint effects of the risk factors considered. These risks caused a substantial proportion of important diseases, including diarrhoea (92%-94%), lower respiratory infections (55-62%), lung cancer (72%), chronic obstructive pulmonary disease (60%), ischaemic heart disease (83-89%), and stroke (70-76%). Removal of these risks would have increased global healthy life expectancy by 9.3 years (17%) ranging from 4.4 years (6%) in the developed countries of the western Pacific to 16.1 years (43%) in parts of sub-Saharan Africa. Interpretation Removal of major risk factors would not only increase healthy life expectancy in every region, but also reduce some of the differences between regions, The potential for disease prevention and health gain from tackling major known risks simultaneously would be substantial.

Quantitative justification for target concentration intervention - parameter variability and predictive performance using population pharmacokinetic models for aminoglycosides

Aims [1] To quantify the random and predictable components of variability for aminoglycoside clearance and volume of distribution [2] To investigate models for predicting aminoglycoside clearance in patients with low serum creatinine concentrations [3] To evaluate the predictive performance of initial dosing strategies for achieving an aminoglycoside target concentration. Methods Aminoglycoside demographic, dosing and concentration data were collected from 697 adult patients (>=20 years old) as part of standard clinical care using a target concentration intervention approach for dose individualization. It was assumed that aminoglycoside clearance had a renal and a nonrenal component, with the renal component being linearly related to predicted creatinine clearance. Results A two compartment pharmacokinetic model best described the aminoglycoside data. The addition of weight, age, sex and serum creatinine as covariates reduced the random component of between subject variability (BSVR) in clearance (CL) from 94% to 36% of population parameter variability (PPV). The final pharmacokinetic parameter estimates for the model with the best predictive performance were: CL, 4.7 l h(-1) 70 kg(-1); intercompartmental clearance (CLic), 1 l h(-1) 70 kg(-1); volume of central compartment (V-1), 19.5 l 70 kg(-1); volume of peripheral compartment (V-2) 11.2 l 70 kg(-1). Conclusions Using a fixed dose of aminoglycoside will achieve 35% of typical patients within 80-125% of a required dose. Covariate guided predictions increase this up to 61%. However, because we have shown that random within subject variability (WSVR) in clearance is less than safe and effective variability (SEV), target concentration intervention can potentially achieve safe and effective doses in 90% of patients.

Adult-population incidence of severe sepsis in Australian and New Zealand intensive care units

Objective. To determine the population incidence and outcome of severe sepsis occurring in adult patients treated in Australian and New Zealand intensive care units (ICUs), and compare with recent retrospective estimates from the USA and UK. Design. Inception cohort study. Setting. Twenty-three closed multi-disciplinary ICUs of 21 hospitals (16 tertiary and 5 university affiliated) in Australia and New Zealand. Patients. A total of 5878 consecutive ICU admission episodes. Measurements and results. Main outcome measures were population-based incidence of severe sepsis, mortality at ICU discharge, mortality at 28 days after onset of severe sepsis, and mortality at hospital discharge. A total of 691 patients, 11.8 (95% confidence intervals 10.9-12.6) per 100 ICU admissions, were diagnosed with 752 episodes of severe sepsis. Site of infection was pulmonary in 50.3% of episodes and abdominal in 19.3% of episodes. The calculated incidence of severe sepsis in adults treated in Australian and New Zealand ICUs is 0.77 (0.76-0.79) per 1000 of population. 26.5% of patients with severe sepsis died in ICU, 32.4% died within 28 days of the diagnosis of severe sepsis and 37.5% died in hospital. Conclusion. In this prospective study, 11.8 patients per 100 ICU admissions were diagnosed with severe sepsis and the calculated annual incidence of severe sepsis in adult patients treated in Australian and New Zealand ICUs is 0.77 per 1000 of population. This figure for the population incidence falls in the lower range of recent estimates from retrospective studies in the U.S. and the U.K.

The distribution of Tannerella forsythia in an adolescent and adult population

Background: The fact that Tannerella forsythia, an important periopathogen, is difficult to cultivate from mixed infections has impeded precise estimates of its distribution within a given population. In order to discern T. forsythia alone from the mixed infection of plaque, the use of sensitive 16S ribosomal RNA based polymerase chain reaction (PCR) detection is necessary. Objectives: The aim of the present study was to determine the distribution of T. forsythia in an adult and in an adolescent population. Materials and methods: Subgingival plaque samples were obtained from 498 Australian adults and from 228 adolescent subjects from Manchester, UK. Tannerella forsythia was detected using PCR and confirmed by restriction analysis. Semi-quantitation of the organisms was carried out using two specific primers of differing sensitivities. Results: In the adolescent population, 25% were found to carry T. forsythia, albeit in relatively low numbers. In the adult population, a total of 37.8% and 11% were found to carry the organism with primer 2 and primer 1, respectively, suggesting that around 27% had between 10(3) and 10(7) organisms. Although there was an apparent increased proportion of T. forsythia positive subjects in those aged >= 50 years, this was not statistical significant. However, T. forsythia positive male smokers showed increased disease severity compared with T. forsythia negative subjects. Conclusion: This study has shown that at least 25% of the adolescent population carry low numbers of T. forsythia, whereas at least 37% of adults carry the organism, with some 11% having relatively high numbers. The relationship between T. forsythia and disease progression in these populations, however, remains to be determined.

Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients

Optimal sampling times are found for a study in which one of the primary purposes is to develop a model of the pharmacokinetics of itraconazole in patients with cystic fibrosis for both capsule and solution doses. The optimal design is expected to produce reliable estimates of population parameters for two different structural PK models. Data collected at these sampling times are also expected to provide the researchers with sufficient information to reasonably discriminate between the two competing structural models.

Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

Effect on trend estimates of the difference between survey respondents and non-respondents: Results from 27 populations in the WHO MONIICA Project

Introduction: In the World Health Organization (WHO) MONICA (multinational MONItoring of trends and determinants in CArdiovascular disease) Project considerable effort was made to obtain basic data on non-respondents to community based surveys of cardiovascular risk factors. The first purpose of this paper is to examine differences in socio-economic and health profiles among respondents and non-respondents. The second purpose is to investigate the effect of non-response on estimates of trends. Methods:Socio-economic and health profile between respondents and non-respondents in the WHO MONICA Project final survey were compared. The potential effect of non-response on the trend estimates between the initial survey and final survey approximately ten years later was investigated using both MONICA data and hypothetical data. Results: In most of the populations, non-respondents were more likely to be single, less well educated, and had poorer lifestyles and health profiles than respondents. As an example of the consequences, temporal trends in prevalence of daily smokers are shown to be overestimated in most populations if they were based only on data from respondents. Conclusions: The socio-economic and health profiles of respondents and non-respondents differed fairly consistently across 27 populations. Hence, the estimators of population trends based on respondent data are likely to be biased. Declining response rates therefore pose a threat to the accuracy of estimates of risk factor trends in many countries.

Lifetime prevalence estimates of major depression: An indirect estimation method and a quantification of recall bias

The measurement of lifetime prevalence of depression in cross-sectional surveys is biased by recall problems. We estimated it indirectly for two countries using modelling, and quantified the underestimation in the empirical estimate for one. A microsimulation model was used to generate population-based epidemiological measures of depression. We fitted the model to 1-and 12-month prevalence data from the Netherlands Mental Health Survey and Incidence Study (NEMESIS) and the Australian Adult Mental Health and Wellbeing Survey. The lowest proportion of cases ever having an episode in their life is 30% of men and 40% of women, for both countries. This corresponds to a lifetime prevalence of 20 and 30%, respectively, in a cross-sectional setting (aged 15-65). The NEMESIS data were 38% lower than these estimates. We conclude that modelling enabled us to estimate lifetime prevalence of depression indirectly. This method is useful in the absence of direct measurement, but also showed that direct estimates are underestimated by recall bias and by the cross-sectional setting.