Multi-level zero-inflated Poisson regression modelling of correlated count data with excess zeros

Count data with excess zeros relative to a Poisson distribution are common in many biomedical applications. A popular approach to the analysis of such data is to use a zero-inflated Poisson (ZIP) regression model. Often, because of the hierarchical Study design or the data collection procedure, zero-inflation and lack of independence may occur simultaneously, which tender the standard ZIP model inadequate. To account for the preponderance of zero counts and the inherent correlation of observations, a class of multi-level ZIP regression model with random effects is presented. Model fitting is facilitated using an expectation-maximization algorithm, whereas variance components are estimated via residual maximum likelihood estimating equations. A score test for zero-inflation is also presented. The multi-level ZIP model is then generalized to cope with a more complex correlation structure. Application to the analysis of correlated count data from a longitudinal infant feeding study illustrates the usefulness of the approach.

Suicide in urban New South Wales, Australia 1985-1994: Socio-economic and migrant interactions

Variation of suicide with socio-economic status (SES) in urban NSW (Australia) during 1985-1994, by sex and country or region of birth, was examined using Poisson regression analysis of vital statistics and population data (age greater than or similar to 15 yr). Quintiles of SES were defined by municipality of residence and comparisons of suicide by SES were adjusted for age and country (or region) of birth (COB), and examined by COB. Risk of suicide in females was 28% that of males for all adults and 21% for youth (age 15-24 yr). Suicide risk was lower in males from southern Europe, Middle East and Asia, and higher in northern and eastern European males, compared to the Australian-born. Risks for suicide increased significantly with decreasing SES in males, but not in females. The relationship of male suicide and SES was stronger when controlled for COB. For males, the relative risk of suicide, adjusted for age and COB, was 66% higher in the lowest SES quintile compared to the highest quintile, and 39% higher for youth (age 15-24 yr). For male suicide, the population attributable fraction for SES (less than the highest quintile) was 27%. Analysis of SES differentials in male suicide according to COB indicated a significant inverse suicide gradient in relation to SES for the Australian-born and those burn in New Zealand and the United Kingdom or fire. but not in non-English speaking COB groups, except for Asia. For Australian-born males, suicide risk was 71% higher in the lowest SES group (compared to the highest), adjusted for age. These findings indicate that SES plays an important role in male suicide rates among the Australian-born and migrants from English-speaking countries and Asia, and among youth; but not in female suicide, nor suicide in most non-English speaking migrant groups. Reduction in SES differentials through economic and social policies may reduce male suicide in lower SES groups and should be seen to be at least as important as individual level interventions. (C) 1998 Elsevier Science Ltd. All rights reserved.

Hysterectomy fractions in New South Wales, 1971-2006

Hysterectomy fractions by age group for particular periods are of interest for: estimating proper population denominators for calculation of disease and procedure rates affecting the cervix and uterus; estimating the target population for Pap test programs, and response rates; and as a way of displaying the cumulative consequences of hysterectomies in a population. Hysterectomy fractions for populations can be determined by direct inquiry via a representative sample survey, or, as in this study, from prior hysterectomy rates of the cohorts of women which compose each age bracket. Hysterectomy data 1979-93 were obtained from the hospital In-patients Statistics Collection (ISC) which covers both public and private hospitals in NSW. Annual population denominators of women were obtained from Census data. Data were modelled by Poisson regression, using five.-year age group (15-greater than or equal to 85 years), annual period, and five-year birth cohort (APC model). Forward- and back-projection of the period effects were undertaken. The resultant NSW hysterectomy fractions by age and period are consistent with fractions obtained from modelled hysterectomy rates for Western Australia (1980-84), and fractions from national representative sample surveys (1989/90 and 1995) for younger women, but not for women aged greater than or equal to 70 years in 1995, which revealed higher hysterectomy fractions than modelled hysterectomy data would suggest. Hysterectomy fractions for NSW women by five-year age group for quinquennia centred on 1971 to 2006 are provided.

Socio-economic, migrant and geographic differentials in coronary heart disease occurrence in New South Wales

Objective: This study examines the variation in coronary heart disease (CHD) mortality and acute myocardial infarction (AMI) by socio-economic status (SES), country of birth (COB) and geography (urban/rural) in the total population of New South Wales (Australia) in 1991-95. Method: CHD deaths and AMI are from complete enumerations of deaths and hospital admissions, respectively; and population denominators are from census information. Data are examined separately by sex, and comparisons of SES groups (based on municipalities), COB and region are analysed using Poisson regression, after adjustment for age. Results: The study identified higher risk for AMI admissions and CHD mortality in lower SES populations with significant linear trends, for both sexes, adjusted for age, region and COB. According to the population attributable fractions (PAF), 23-41% of the risk of CHD occurrence is due to SES lower than the highest quartile. The higher age-adjusted risk for CHD occurrence in rural and remote populations for both sexes, compared with urban communities, was lessened by adjustment for COB, and all but abolished when also adjusted for SES. COB analysis indicated significantly lower age-adjusted AMI admissions and CHD mortality compared with the Australian-born, Conclusions: Higher risks for CHD in rural populations compared with the capital city (Sydney) are due, in part, to lower SES, lesser migrant composition. Implications: Strategies for reducing CHD differentials should consider demographic factors and the fundamental need to reduce socio-economic inequalities, as well as targeting appropriate prevention measures.

Urban and rural suicide differentials in migrants and the Australian-born, New South Wales, Australia 1985-1994

We estimated risk of suicide in adults in New South Wales (NSW) by sex, country of birth and rural/urban residence, after adjusting for age; we also examined youth suicide (age 15-24 years). The study population was the entire population of NSW, Australia, aged greater than or equal to 15 years during the period 1985-1994. Poisson regression was used to examine the relationship between predictor variables and the risk of suicide, with the focus on migrant status and area of residence. A significantly higher risk of suicide was found in male migrants from Northern Europe and Eastern Europe/former USSR, compared to Australian-born males; a significantly lower suicide risk occurred in males from Southern Europe, the Middle East and Asia. In female migrants, those from UK/Eire, Northern Europe, Eastern Europe/former USSR and New Zealand exhibited a significantly higher risk of suicide compared to Australian-born females. A significantly lower risk of suicide occurred in females from the Middle East. Male migrants overall were at significantly lower risk of suicide than the Australian-born, while female migrants overall had a significantly higher risk of suicide than Australian-born females. Among migrant males overall, the rural-urban suicide risk differential was significantly higher for those living in non-metropolitan areas (RR = 1.9; 95% CI: 1.7-2.1). Suicide risk was significantly higher in non-metropolitan male immigrants from the UK/Eire (RR = 1.4; 95% CI: 1.1-1.7), Southern Europe (RR = 1.7; 95% CI: 1.2-2.4), Northern/Western Europe (1.5; 95% CI: 1.2-1.9), the Middle East (RR = 3.8; 95% CI: 1.9-7.8), New :Zealand (RR = 1.4; 95% CI: 1.0-1.8) and 'other' (RR = 2.6; 95% CI: 1.9-3.5), when compared to their urban counterparts. There was no statistically significant difference in suicide risk between rural and urban Australian-born males. For female suicide, significantly lower risk was found in female immigrants living in non-metropolitan areas who were from Northern/Western Europe (RR = 0.7; 95% CI: 0.4-0.96), as well as the Australian-born (RR = 0.7; 95% CI: 0.6-0.8), when compared to their urban counterparts. The non-metropolitan/metropolitan relative risk for suicide in female migrants overall was not significantly different from one. Among male youth there was a significantly higher suicide risk in non-metropolitan areas, with a relative risk estimate of 1.4 for Australian-born youth (95% CI: 1.2-1.5) and 1.7 for migrant youth (95% CI: 1.2-2.4), when compared with metropolitan counterparts. We conclude that suicide among migrant males living in non-metropolitan areas accounts for most of the excess of male suicide in rural NSW, and the significantly lower risk of suicide for non-metropolitan Australian-born women does not apply to migrant women. (C) 1999 Published by Elsevier Science Ltd. All rights reserved.

Trends in the incidence, severity, and short-term outcome of stroke in Perth, Western Australia

Background and Purpose-This report describes trends in the key indices of cerebrovascular disease over 6 years from the end of the 1980s in a geographically defined segment of the city of Perth, Western Australia. Methods-Identical methods were used to find and assess all cases of suspected stroke in a population of approximately 134 000 residents in a triangular area of the northern suburbs of Perth. Case fatality was measured as vital status at 28 days after the onset of symptoms. Data for first-ever strokes and for all strokes for equivalent periods of 12 months in 1989-1990 and 1995-1996 were compared by age-standardized rates and proportions and Poisson regression. Results-There were 355 strokes in 328 patients and 251 first-ever strokes (71%) for 1989-1990 and 290 events in 281 patients and 213 first-ever strokes (73%) for 1995-1996. In Poisson models including age and period, overall trends in the incidence of both first-ever strokes (rate ratio = 0.75; 95% confidence limits, 0.63, 0.90) and all strokes (rate ratio = 0.73; 95% confidence limits, 0.62, 0.85) were obviously significant, but only the changes in men were independently significant. Case fatality did not change, and the balance between hemorrhagic and occlusive strokes in 1995-1996 was almost indistinguishable from that observed in 1989-1990. Conclusions-Our results, which are the only longitudinal population-based data available for Australia for key indices of stroke, suggest that it is a change in the frequency of stroke, rather than its outcome, that is chiefly responsible nationally for the fall in mortality from cerebrovascular disease.

Dental set-vice provision in private general practice for an insured population in New South Wales, Australia

Objective To determine patterns of dental set-vices provided to a cohort of the insured population 18 years and over, in private general practice in New South Wales, Australia. Basic research design A cohort study using the person-years method and Poisson regression for analysis. Setting Data were derived from claims records submitted by members of a health insurance fund (Government Employees Health Fund-GEHF) for rebates during the study period 1 January 1992-31 December 1995. Participants There were 133,467 members aged 18 years and over from New South Wales. Main outcome measures To determine, by age group, for those members who used private general practice and made a claim (referred to as 'patients') the annual number of visits, total number of services received per year and number of services received at a visit, Results The mean number of visits per patient was 2.4 per year with patients under 45 years making fewer visits than the 45-54 age group reference category. Mean number of services utilised per patient-year was 5.9, with services provided increasing from 3.5 for the 18-24-year-old group, reaching a plateau of approximately 6.2 for those aged 45 years or more. The number of services received per visit was 2.4 and there were no differences by age. Service mix was dominated by restorative (35%), diagnostic (27%), and preventive services (18%); with age specific variations observed. Conclusions Age was found to be an important determinant in the use of dental services, independent of membership duration and gender.

Increased rates of ENT surgery among young children: Have clinical guidelines made a difference?

Objectives: To examine the association between introduction of paediatric ear, nose and throat (ENT) surgery guidelines and population procedure rates. To determine changes in children's risk of undergoing ENT surgery. Methods: Trend analysis of incidence of myringotomy, tonsillectomy and adenoidectomy among New South Wales (NSW) children aged 0-14 between 1981 and mid 1999. Poisson regression models were used to estimate annual rates of change pre and postguidelines introduction and age/gender specific rates, and lifetable methods to determine risk of undergoing an ENT procedure by age 15. Results: ENT surgery rates increased by 21% over the study period. Children's risk of surgery increased from 17.9% in 1981 to 20.2% in 1998/99. Guideline introduction was associated with moderate short-term decreases in rates. For tonsillectomy, rates decreased between 1981 and 1983, but then rose continually until the introduction of myringotomy guidelines in 1993, when they fell, only to recommence rising until the end of the study period. For myringotomy, rates rose annually from 1981 to 1992/93 and fell in the 3 years following guideline introduction, after which they rose again. Increases were almost exclusively restricted to children aged 0-4 and correspond with increased use of formal childcare. The prevalence of myringotomy by the age of 5 years rose from 5.6% of children born in 1988/89 to 6.4% of those born in 1994/95, and the prevalence of tonsillectomy from 2.4% to 2.7%. Conclusions: The risk of young Australian children undergoing ENT surgery increased significantly over the last two decades despite the introduction of guidelines and no evidence of an increase in otitis media, one condition prompting surgery. Surgery increased most among the very young. We hypothesize this is related to increasing use of childcare.

Mental health and socio-economic variations in Australian suicide

This paper investigates the relationship between suicide rates and prevalence of mental disorder and suicide attempts, across socio-economic status (SES) groups based on area of residence. Australian suicide data (1996-1998) were analysed in conjunction with area-based prevalences of mental disorder derived from the National Survey of Mental Health and Well-Being (1997). Poisson regression models of suicide risk included age, quintile of area-based SES, urban-rural residence, and country of birth (COB), with males and females analysed separately. Analysis focussed on the association between suicide and prevalences of (ICD-10) affective disorders, anxiety disorders, substance use disorders and suicide attempts by SES group. Prevalences of other psychiatric symptomatology, substance use problems, health service utilisation, stressful life-events and personality were also investigated. Significant increasing gradients were evident from high to low SES groups for prevalences of affective disorders, anxiety disorders (females only), and substance use disorders (males only); sub-threshold drug and alcohol problems and depression; and suicide attempts and suicide (males only). Prevalences of mental disorder, other sub-threshold mental health items and suicide attempts were significantly associated with suicide, but in most cases associations were reduced in magnitude and became statistically non-significant after adjustment for COB, urban-rural residence, and SES. For male suicide the relative risk (RR) in the lowest SES group compared to the highest was 1.40 (95% CI 1.29-1.52, p < 0.001) for all ages, and 1.46 (95% CI 1.27-1.67, p < 0.001) for male youth (20-34 years). This relationship was not substantially modified in males when regression models included prevalences of affective disorders, and other selected mental health variables and demographic factors. From a population perspective, SES remained significantly associated with suicide after controlling for the prevalence of mental disorders and other psychiatric symptomatology. Mental conditions and previous suicidal behaviour may play an intermediary role between SES and suicide, but this study suggests that an independent relationship between suicide and SES also exists. (c) 2005 Elsevier Ltd. All rights reserved.

Cancer detection and mammogram volume of radiologists in a population-based screening programme

This study investigates the relationship between the number of screening mammograms read by radiologists and the screening breast cancer detection rate. Cancer detection rates for incident screens (all women aged >= 40 years) were compared by increasing categories of reader volume using Poisson regression. Data from New South Wales (NSW) for a 2 year period (2000-2001) were obtained from the BreastScreen NSW programme. Cancer detection rates increased with the number of mammograms read in the programme, reaching a plateau of approximately 40 per 10,000 after 1375 mammograms per year. No significant differences in cancer detection were evident above 875 mammograms (compared to below 875 mammograms) per year (RR = 0.79, 95% CI 0.63-0.99). (c) 2005 Elsevier Ltd. All rights reserved.

Anthropometric measures in relation to Basal Cell Carcinoma: a longitudinal study

Background: The relationship between anthropometric indices and risk of basal cell carcinoma ( BCC) is largely unknown. We aimed to examine the association between anthropometric measures and development of BCC and to demonstrate whether adherence to World Health Organisation guidelines for body mass index, waist circumference, and waist/ hip ratio was associated with risk of BCC, independent of sun exposure. Methods: Study participants were participants in a community- based skin cancer prevention trial in Nambour, a town in southeast Queensland ( latitude 26 degrees S). In 1992, height, weight, and waist and hip circumferences were measured for all 1621 participants and weight was remeasured at the end of the trial in 1996. Prevalence proportion ratios were calculated using a log- binomial model to estimate the risk of BCC prior to or prevalent in 1992, while Poisson regression with robust error variances was used to estimate the relative risk of BCC during the follow- up period. Results: At baseline, 94 participants had a current BCC, and 202 had a history of BCC. During the 5- year follow- up period, 179 participants developed one or more new BCCs. We found no significant association between any of the anthropometric measures or indices and risk of BCC after controlling for potential confounding factors including sun exposure. There was a suggestion that short- term weight gain may increase the risk of developing BCC for women only. Conclusion: Adherence to World Health Organisation guidelines for body mass index, waist circumference and waist/ hip ratio is not significantly associated with occurrence of basal cell carcinomas of the skin.

Asthma mortality in Australia 1920-94: Age, period, and cohort effects

Study objective-To investigate asthma mortality during 1920-94 in Australia in order to assess the relative role of period and birth cohort effects. Design-Asthma mortality (both sexes) was age standardised and examined for changes over time. The data were also examined for age, period, and cohort (APC) effects using Poisson regression modelling. Setting-National Australian mortality data. Participants-Population (both sexes) aged 15-34 years, 1920-94. Main results-Age adjusted period rates indicate an increase in asthma mortality during the 1950s, and increases and subsequent falls (epidemics) during the mid 1960s and late 1980s. APC modelling suggested an increasing cohort effect (adjusted for both age and period) from the birth cohort 1950-54 onwards. Period effects (adjusted for age and cohort) are characterised by an increase in the 1950s (possibly due to changes in diagnostic labelling), minimal or no increases in the mid 1960s and late 1980s (where period peaks had been noted when data were adjusted for age only), and declines in mortality risk subsequent to the periods where age-period analysis had noted increases. Thus, in Australia, some of the mid 1960s epidemic in asthma deaths, and all of the late 1980s mortality increase, seem to be attributable to cohort effects. Conclusions-The increase in asthma mortality cohort effect is consistent with empirical evidence of recent increases in prevalence (and presumably incidence) of asthma in Australia, and suggests the need for more research into the underlying environmental aetiology of this condition.

Assessing the impacts of grazing levels on bird density in woodland habitat: a Bayesian approach using expert opinion

Many studies on birds focus on the collection of data through an experimental design, suitable for investigation in a classical analysis of variance (ANOVA) framework. Although many findings are confirmed by one or more experts, expert information is rarely used in conjunction with the survey data to enhance the explanatory and predictive power of the model. We explore this neglected aspect of ecological modelling through a study on Australian woodland birds, focusing on the potential impact of different intensities of commercial cattle grazing on bird density in woodland habitat. We examine a number of Bayesian hierarchical random effects models, which cater for overdispersion and a high frequency of zeros in the data using WinBUGS and explore the variation between and within different grazing regimes and species. The impact and value of expert information is investigated through the inclusion of priors that reflect the experience of 20 experts in the field of bird responses to disturbance. Results indicate that expert information moderates the survey data, especially in situations where there are little or no data. When experts agreed, credible intervals for predictions were tightened considerably. When experts failed to agree, results were similar to those evaluated in the absence of expert information. Overall, we found that without expert opinion our knowledge was quite weak. The fact that the survey data is quite consistent, in general, with expert opinion shows that we do know something about birds and grazing and we could learn a lot faster if we used this approach more in ecology, where data are scarce. Copyright (c) 2005 John Wiley & Sons, Ltd.

Caregiver burden, time spent caring and health status in the first 12 months following stroke

Objective: To quantify time caring, burden and health status in carers of stroke patients after discharge from rehabilitation; to identify the potentially modifiable sociodemographic and clinical characteristics associated with these outcomes. Methods: Patients and carers prospectively interviewed 6 (n = 71) and 12 (n = 57) months after discharge. Relationships of carer and patient variables with burden, health status and time analysed by Gaussian and Poisson regression. Results: Carers showed considerable burden at 6 and 12 months. Carers spent 4.6 and 3.6 hours per day assisting patients with daily activities at 6 and 12 months, respectively. Improved patient motor and cognitive function were associated with reductions of up to 20 minutes per day in time spent in daily activities. Better patient mental health and cognitive function were associated with better carer mental health. Conclusions: Potentially modifiable factors such as these may be able to be targeted by caregiver training, support and education programmes and outpatient therapy for patients.