Search

In a stratified analysis, the results from different strata if homogeneity assumption is met are pooled together to obtain a single summary estimate for the common effect parameter. However, the effect can appear homogeneous across strata using one measure but heterogeneous using another. Consequently, two researchers analyzing the same data can arrive at conflicting conclusions if they use different effect measures. In this paper, the author draws on the sufficient component cause model to develop a stratified-analysis method regarding a particular effect measure, the ‘peril ratio’. When there is no mechanistic interaction between the exposure under study and the stratifying variable (i.e., when they do not work together to complete any sufficient cause), the peril ratio is constant across strata. The author presents formulas for the estimation of such a common peril ratio. Three real data are re-analyzed for illustration. When the data is consistent with peril-ratio homogeneity in a stratified analysis, researchers can use the formulas in this paper to pool the strata. [ABSTRACT FROM AUTHOR]

As a performance measure for a prediction model, the area under the receiver operating characteristic curve (AUC) is insensitive to the addition of strong markers. A number of measures sensitive to performance change have recently been proposed; however, these relative-performance measures may lead to self-contradictory conclusions. This paper examines alternative performance measures for prediction models: the Lorenz curve-based Gini and Pietra indices, and a standardized version of the Brier score, the scaled Brier. Computer simulations are performed in order to study the sensitivity of these measures to performance change when a new marker is added to a baseline model. When the discrimination power of the added marker is concentrated in the gray zone of the baseline model, the AUC and the Gini show minimal performance improvements. The Pietra and the scaled Brier show more significant improvements in the same situation, comparatively. The Pietra and the scaled Brier indices are therefore recommended for prediction model performance measurement, in light of their ease of interpretation, clinical relevance and sensitivity to gray-zone resolving markers. [ABSTRACT FROM AUTHOR]

Knowledge of the incubation period of infectious diseases (time between host infection and expression of disease symptoms) is crucial to our epidemiological understanding and the design of appropriate prevention and control policies. Plant diseases cause substantial damage to agricultural and arboricultural systems, but there is still very little information about how the incubation period varies within host populations. In this paper, we focus on the incubation period of soilborne plant pathogens, which are difficult to detect as they spread and infect the hosts underground and above-ground symptoms occur considerably later. We conducted experiments on Rhizoctonia solani in sugar beet, as an example patho-system, and used modelling approaches to estimate the incubation period distribution and demonstrate the impact of differing estimations on our epidemiological understanding of plant diseases. We present measurements of the incubation period obtained in field conditions, fit alternative probability models to the data, and show that the incubation period distribution changes with host age. By simulating spatially-explicit epidemiological models with different incubation-period distributions, we study the conditions for a significant time lag between epidemics of cryptic infection and the associated epidemics of symptomatic disease. We examine the sensitivity of this lag to differing distributional assumptions about the incubation period (i.e. exponential versus Gamma). We demonstrate that accurate information about the incubation period distribution of a pathosystem can be critical in assessing the true scale of pathogen invasion behind early disease symptoms in the field; likewise, it can be central to model-based prediction of epidemic risk and evaluation of disease management strategies. Our results highlight that reliance on observation of disease symptoms can cause significant delay in detection of soil-borne pathogen epidemics and mislead practitioners and epidemiologists about the timing, extent, and viability of disease control measures for limiting economic loss. [ABSTRACT FROM AUTHOR]

No matching vaccine is immediately available when a novel influenza strain breaks out. Several nonvaccine-related strategies must be employed to control an influenza epidemic, including antiviral treatment, patient isolation, and immigration detection. This paper presents the development and application of two regional dynamic models of influenza with Pontryagin’s Maximum Principle to determine the optimal control strategies for an epidemic and the corresponding minimum antiviral stockpiles. Antiviral treatment was found to be the most effective measure to control new influenza outbreaks. In the case of inadequate antiviral resources, the preferred approach was the centralized use of antiviral resources in the early stage of the epidemic. Immigration detection was the least cost-effective; however, when used in combination with the other measures, it may play a larger role. The reasonable mix of the three control measures could reduce the number of clinical cases substantially, to achieve the optimal control of new influenza. [ABSTRACT FROM AUTHOR]

We consider the recently introduced edge-based compartmental models (EBCM) for the spread of susceptible-infected-recovered (SIR) diseases in networks. These models differ from standard infectious disease models by focusing on the status of a random partner in the population, rather than a random individual. This change in focus leads to simple analytic models for the spread of SIR diseases in random networks with heterogeneous degree. In this paper we extend this approach to handle deviations of the disease or population from the simplistic assumptions of earlier work. We allow the population to have structure due to effects such as demographic features or multiple types of risk behavior. We allow the disease to have more complicated natural history. Although we introduce these modifications in the static network context, it is straightforward to incorporate them into dynamic network models. We also consider serosorting, which requires using dynamic network models. The basic methods we use to derive these generalizations are widely applicable, and so it is straightforward to introduce many other generalizations not considered here. Our goal is twofold: to provide a number of examples generalizing the EBCM method for various different population or disease structures and to provide insight into how to derive such a model under new sets of assumptions. [ABSTRACT FROM AUTHOR]

Methicillin resistant Staphylococcus aureus (MRSA) is currently a major cause of skin and soft tissue infections (SSTI) in the United States. Seasonal variation of MRSA infections in hospital settings has been widely observed. However, systematic time-series analysis of incidence data is desirable to understand the seasonality of community acquired (CA)-MRSA infections at the population level. In this paper, using data on monthly SSTI incidence in children aged 0–19 years and enrolled in Medicaid in Maricopa County, Arizona, from January 2005 to December 2008, we carried out time-series and nonlinear regression analysis to determine the periodicity, trend, and peak timing in SSTI incidence in children at different age: 0–4 years, 5–9 years, 10–14 years, and 15–19 years. We also assessed the temporal correlation between SSTI incidence and meteorological variables including average temperature and humidity. Our analysis revealed a strong annual seasonal pattern of SSTI incidence with peak occurring in early September. This pattern was consistent across age groups. Moreover, SSTIs followed a significantly increasing trend over the 4-year study period with annual incidence increasing from 3.36% to 5.55% in our pediatric population of approximately 290,000. We also found a significant correlation between the temporal variation in SSTI incidence and mean temperature and specific humidity. Our findings could have potential implications on prevention and control efforts against CA-MRSA. [ABSTRACT FROM AUTHOR]

Background: Recent focus on earlier detection of pathogen introduction in human and animal populations has led to the development of surveillance systems based on automated monitoring of health data. Real- or near real-time monitoring of pre-diagnostic data requires automated classification of records into syndromes–syndromic surveillance–using algorithms that incorporate medical knowledge in a reliable and efficient way, while remaining comprehensible to end users. Methods: This paper describes the application of two of machine learning (Naïve Bayes and Decision Trees) and rule-based methods to extract syndromic information from laboratory test requests submitted to a veterinary diagnostic laboratory. Results: High performance (F1-macro = 0.9995) was achieved through the use of a rule-based syndrome classifier, based on rule induction followed by manual modification during the construction phase, which also resulted in clear interpretability of the resulting classification process. An unmodified rule induction algorithm achieved an F1-micro score of 0.979 though this fell to 0.677 when performance for individual classes was averaged in an unweighted manner (F1-macro), due to the fact that the algorithm failed to learn 3 of the 16 classes from the training set. Decision Trees showed equal interpretability to the rule-based approaches, but achieved an F1-micro score of 0.923 (falling to 0.311 when classes are given equal weight). A Naïve Bayes classifier learned all classes and achieved high performance (F1-micro = 0.994 and F1-macro = .955), however the classification process is not transparent to the domain experts. Conclusion: The use of a manually customised rule set allowed for the development of a system for classification of laboratory tests into syndromic groups with very high performance, and high interpretability by the domain experts. Further research is required to develop internal validation rules in order to establish automated methods to update model rules without user input. [ABSTRACT FROM AUTHOR]

Background: Graphical representation of data is one of the most easily comprehended forms of explanation. The current study describes a simple visualization tool which may allow greater understanding of medical and epidemiological data. Method: We propose a simple tool for visualization of data, known as a “quilt plot”, that provides an alternative to presenting large volumes of data as frequency tables. Data from the Australian Needle and Syringe Program survey are used to illustrate “quilt plots”. Conclusion: Visualization of large volumes of data using “quilt plots” enhances interpretation of medical and epidemiological data. Such intuitive presentations are particularly useful for the rapid assessment of problems in the data which cannot be readily identified by manual review. We recommend that, where possible, “quilt plots” be used along with traditional quantitative assessments of the data as an explanatory data analysis tool. [ABSTRACT FROM AUTHOR]

Disease control is of paramount importance in public health, with infectious disease extinction as the ultimate goal. Although diseases may go extinct due to random loss of effective contacts where the infection is transmitted to new susceptible individuals, the time to extinction in the absence of control may be prohibitively long. Intervention controls are typically defined on a deterministic schedule. In reality, however, such policies are administered as a random process, while still possessing a mean period. Here, we consider the effect of randomly distributed intervention as disease control on large finite populations. We show explicitly how intervention control, based on mean period and treatment fraction, modulates the average extinction times as a function of population size and rate of infection spread. In particular, our results show an exponential improvement in extinction times even though the controls are implemented using a random Poisson distribution. Finally, we discover those parameter regimes where random treatment yields an exponential improvement in extinction times over the application of strictly periodic intervention. The implication of our results is discussed in light of the availability of limited resources for control. [ABSTRACT FROM AUTHOR]

The objectives of this paper are to predict life years lost associated with obesity-related diseases (ORDs) for U.S. non-smoking adults, and to examine the relationship between those ORDs and mortality. Data from the National Health Interview Survey, 1997–2000, were used. We employed mixed proportional hazard models to estimate the association between those ORDs and mortality and used simulations to project life years lost associated with the ORDs. We found that obesity-attributable comorbidities are associated with large decreases in life years and increases in mortality rates. The life years lost associated with ORDs is more marked for younger adults than older adults, for blacks than whites, for males than females, and for the more obese than the less obese. Using U.S. non-smoking adults aged 40 to 49 years as an example to illustrate percentage of the life years lost associated with ORDs, we found that the mean life years lost associated with ORDs for U.S. non-smoking black males aged 40 to 49 years with a body mass index above 40 kg/m2 was 5.43 years, which translates to a 7.5% reduction in total life years. White males of the same age range and same degree of obesity lost 5.23 life years on average – a 6.8% reduction in total life years, followed by black females (5.04 years, a 6.5% reduction in life years), and white females (4.7 years, a 5.8% reduction in life years). Overall, ORDs increased chances of dying and lessened life years by 0.2 to 11.7 years depending on gender, race, BMI classification, and age. [ABSTRACT FROM AUTHOR]

Several approaches have been proposed for near real-time detection and prediction of the spread of influenza. These include search query data for influenza-related terms, which has been explored as a tool for augmenting traditional surveillance methods. In this paper, we present a method that uses Internet search query data from Baidu to model and monitor influenza activity in China. The objectives of the study are to present a comprehensive technique for: (i) keyword selection, (ii) keyword filtering, (iii) index composition and (iv) modeling and detection of influenza activity in China. Sequential time-series for the selected composite keyword index is significantly correlated with Chinese influenza case data. In addition, one-month ahead prediction of influenza cases for the first eight months of 2012 has a mean absolute percent error less than 11%. To our knowledge, this is the first study on the use of search query data from Baidu in conjunction with this approach for estimation of influenza activity in China. [ABSTRACT FROM AUTHOR]

Background: Evidence regarding the mortality rate after administration of the pandemic influenza A (H1N1) 2009 vaccine on patients with underlying diseases is currently scarce. We conducted a case-control study in Japan to compare the mortality rates of patients with idiopathic interstitial pneumonia after the vaccines were administered and were not administered. Methods: Between October 2009 and March 2010, we collected clinical records in Japan and conducted a 1∶1 matched case-control study. Patients with idiopathic interstitial pneumonia who died during this period were considered case patients, and those who survived were considered control patients. We determined and compared the proportion of each group that received the pandemic influenza A (H1N1) 2009 vaccine and estimated the odds ratio. Finally, we conducted simulations that compensated for the shortcomings of the study associated with adjusted severity of idiopathic interstitial pneumonia. Results: The case and control groups each comprised of 75 patients with idiopathic interstitial pneumonia. The proportion of patients who received the pandemic influenza A (H1N1) 2009 vaccine was 30.7% and 38.7% for the case and control groups, respectively. During that winter, the crude conditional odds ratio of mortality was 0.63 (95% confidence interval, 0.25–1.47) and the adjusted conditional odds ratio was 1.18 (95% confidence interval, 0.33–4.49); neither was significant. The simulation study showed more accurate conditional odds ratios of 0.63–0.71. Conclusions: In our study, we detected no evidence that the influenza A (H1N1) 2009 vaccine increased the mortality rate of patients with idiopathic interstitial pneumonia. The results, however, are limited by the small sample size and low statistical power. A larger-scale study is required. [ABSTRACT FROM AUTHOR]

Avian influenza viruses (AIVs) are considered a threat for their potential to seed human influenza pandemics. Despite their acknowledged importance, there are significant unknowns regarding AIV transmission dynamics in their natural hosts, wild birds. Of particular interest is the difference in subtype dynamics between human and bird populations–in human populations, typically only two or three subtypes cocirculate, while avian populations are capable of simultaneously hosting a multitude of subtypes. One species in particular–ruddy turnstones (Arenaria interpres)–has been found to harbour a very wide range of AIV subtypes, which could make them a key player in the spread of new subtypes in wild bird populations. Very little is known about the mechanisms that drive subtype dynamics in this species, and here we address this gap in our knowledge. Taking advantage of two independent sources of data collected from ruddy turnstones in Delaware Bay, USA, we examine patterns of subtype diversity and dominance at this site. We compare these patterns to those produced by a stochastic, multi-strain transmission model to investigate possible mechanisms that are parsimonious with the observed subtype dynamics. We find, in agreement with earlier experimental work, that subtype differences are unnecessary to replicate the observed dynamics, and that neutrality alone is sufficient. We also evaluate the role of subtype cross-immunity and find that it is not necessary to generate patterns consistent with observations. This work offers new insights into the mechanisms behind subtype diversity and dominance in a species that has the potential to be a key player in AIV dynamics in wild bird populations. [ABSTRACT FROM AUTHOR]

Osteoporosis is characterized by reduced bone mineral density (BMD) and increased fracture risk. Fat mass is a determinant of bone strength and both phenotypes have a strong genetic component. In this study, we examined the association between obesity associated polymorphisms (SNPs) with body composition, BMD, Ultrasound (QUS), fracture and biomarkers (Homocysteine (Hcy), folate, Vitamin D and Vitamin B12) for obesity and osteoporosis. Five common variants: rs17782313 and rs1770633 (melanocortin 4 receptor (MC4R); rs7566605 (insulin induced gene 2 (INSIG2); rs9939609 and rs1121980 (fat mass and obesity associated (FTO) were genotyped in 2 cohorts of Swedish women: PEAK-25 (age 25, n = 1061) and OPRA (age 75, n = 1044). Body mass index (BMI), total body fat and lean mass were strongly positively correlated with QUS and BMD in both cohorts (r2 = 0.2–0.6). MC4R rs17782313 was associated with QUS in the OPRA cohort and individuals with the minor C-allele had higher values compared to T-allele homozygotes (TT vs. CT vs. CC: BUA: 100 vs. 103 vs. 103; p = 0.002); (SOS: 1521 vs. 1526 vs. 1524; p = 0.008); (Stiffness index: 69 vs. 73 vs. 74; p = 0.0006) after adjustment for confounders. They also had low folate (18 vs. 17 vs. 16; p = 0.03) and vitamin D (93 vs. 91 vs. 90; p = 0.03) and high Hcy levels (13.7 vs 14.4 vs. 14.5; p = 0.06). Fracture incidence was lower among women with the C-allele, (52% vs. 58%; p = 0.067). Variation in MC4R was not associated with BMD or body composition in either OPRA or PEAK-25. SNPs close to FTO and INSIG2 were not associated with any bone phenotypes in either cohort and FTO SNPs were only associated with body composition in PEAK-25 (p≤0.001). Our results suggest that genetic variation close to MC4R is associated with quantitative ultrasound and risk of fracture. [ABSTRACT FROM AUTHOR]

Background: Cadmium (Cd) is a heavy metal that can cause renal tubular dysfunction in humans. Women are among the high-risk group for Cd health effects. Determining the thresholds of Cd-induced renal effects is important. Thus, in this article, we aimed to identify the benchmark dose (BMD) and its low limit (BMDL) levels as the Cd thresholds for Chinese women. Methods: Epidemiologic investigation was performed in county A and county B to obtain data on Cd exposure and its renal effect on respondents. Levels of Cd (UCd), β2-microglobulin (UB2M), and N-acetyl-β-D-glucosaminidase (UNAG) were measured in morning urine samples. The BMD approach was mainly performed. Results: Results of the BMD approach were similar whether the method was conducted for the two sets of data (collected in CA and CB, respectively) separately or cooperatively. The BMD/BMDL values of UCd for all subjects were 1.07/0.44 and 2.12/0.53 µg/g cr based on UB2M and UNAG, respectively, given a predetermined BMR of 0.05. Conclusions: The presented thresholds of Cd-induced renal effects (i.e., the BMDLs of UCd) are close to the counterpart values reported in Japan, Sweden and Belguim. [ABSTRACT FROM AUTHOR]

Background and Objectives: Nicaragua is highly endemic for hepatitis A. We aimed to provide an estimate of the change in the age-specific risk of hepatitis A virus (HAV) infection based on serological data from cross-sectional and longitudinal samples collected in León, Nicaragua, in 1995/96 (n = 979) and 2003 (n = 494). Methods: The observed age-specific prevalence of anti-HAV antibodies was correlated to the age-specific risk of infection by calculating the probability of freedom from infection at a specific age. Results: The proportion of seropositive children aged 1.5 to 6 years was 42% in 2003 compared to 67% in 1995/96. Estimated annual risk of infection for a 3-year old child was 30% (95% CI: 27.0%, 33.1%) in 1995 and 15.5% (95% CI: 12.4%, 19.0%) in 2003. There was good agreement between estimates based on cross-sectional and longitudinal data. The age-specific geometric mean of the quantified anti-HAV antibody levels assessed in 2003 was highest at age 4 and decreased steadily up to age 40. Conclusions: The substantially lower risk of HAV infection in 2003 than in 1995 for young children indicates a beginning transition from high to intermediate endemicity in León, Nicaragua. Consecutive age-stratified serosurveys are useful to assess changes in risk of infection following public health interventions. The decreasing age-specific GMC of anti-HAV antibodies during adulthood in a country with endemic HAV indirectly suggests that ongoing HAV exposure in the community has marginal boosting effect on antibody levels once protective immunity has been established by natural infection. [ABSTRACT FROM AUTHOR]

The one-size-fits-all paradigm in organized screening of breast cancer is shifting towards a personalized approach. The present study has two objectives: 1) To perform an economic evaluation and to assess the harm-benefit ratios of screening strategies that vary in their intensity and interval ages based on breast cancer risk; and 2) To estimate the gain in terms of cost and harm reductions using risk-based screening with respect to the usual practice. We used a probabilistic model and input data from Spanish population registries and screening programs, as well as from clinical studies, to estimate the benefit, harm, and costs over time of 2,624 screening strategies, uniform or risk-based. We defined four risk groups, low, moderate-low, moderate-high and high, based on breast density, family history of breast cancer and personal history of breast biopsy. The risk-based strategies were obtained combining the exam periodicity (annual, biennial, triennial and quinquennial), the starting ages (40, 45 and 50 years) and the ending ages (69 and 74 years) in the four risk groups. Incremental cost-effectiveness and harm-benefit ratios were used to select the optimal strategies. Compared to risk-based strategies, the uniform ones result in a much lower benefit for a specific cost. Reductions close to 10% in costs and higher than 20% in false-positive results and overdiagnosed cases were obtained for risk-based strategies. Optimal screening is characterized by quinquennial or triennial periodicities for the low or moderate risk-groups and annual periodicity for the high-risk group. Risk-based strategies can reduce harm and costs. It is necessary to develop accurate measures of individual risk and to work on how to implement risk-based screening strategies. [ABSTRACT FROM AUTHOR]

Background: Variable selection is an important step in building a multivariate regression model for which several methods and statistical packages are available. A comprehensive approach for variable selection in complex multivariate regression analyses within HIV cohorts is explored by utilizing both epidemiological and biostatistical procedures. Methods: Three different methods for variable selection were illustrated in a study comparing survival time between subjects in the Department of Defense’s National History Study and the Atlanta Veterans Affairs Medical Center’s HIV Atlanta VA Cohort Study. The first two methods were stepwise selection procedures, based either on significance tests (Score test), or on information theory (Akaike Information Criterion), while the third method employed a Bayesian argument (Bayesian Model Averaging). Results: All three methods resulted in a similar parsimonious survival model. Three of the covariates previously used in the multivariate model were not included in the final model suggested by the three approaches. When comparing the parsimonious model to the previously published model, there was evidence of less variance in the main survival estimates. Conclusions: The variable selection approaches considered in this study allowed building a model based on significance tests, on an information criterion, and on averaging models using their posterior probabilities. A parsimonious model that balanced these three approaches was found to provide a better fit than the previously reported model. [ABSTRACT FROM AUTHOR]

Comparison of control strategies against animal infectious diseases allows determining optimal strategies according to their epidemiological and/or economic impacts. However, in real life, the choice of a control strategy does not always obey a pure economic or epidemiological rationality. The objective of this study was to analyze the choice of a foot and mouth disease (FMD) control strategy as a decision-making process in which the decision-maker is influenced by several stakeholders (government, agro-food industries, public opinion). For each of these, an indicator of epizootic impact was quantified to compare seven control strategies. We then determined how, in France, the optimal control strategy varied according to the relative weights of stakeholders and to the perception of risk by the decision-maker (risk-neutral/risk-averse). When the scope of decision was national, whatever their perception of risk and the stakeholders' weights, decision-makers chose a strategy based on vaccination. This consensus concealed marked differences between regions, which were connected with the regional breeding characteristics. Vaccination-based strategies were predominant in regions with dense cattle and swine populations, and in regions with a dense population of small ruminants, combined with a medium density of cattle and swine. These differences between regions suggested that control strategies could be usefully adapted to local breeding conditions. We then analyzed the feasibility of adaptive decision-making processes depending on the date and place where the epizootic starts, or on the evolution of the epizootic over time. The initial conditions always explained at least half of the variance of impacts, the remaining variance being attributed to the variability of epizootics evolution. However, the first weeks of this evolution explained a large part of the impacts variability. Although the predictive value of the initial conditions for determining the optimal strategy was weak, adaptive strategies changing dynamically according to the evolution of the epizootic appeared feasible. [ABSTRACT FROM AUTHOR]

Background: A combination of clinical and routine laboratory data with biomarkers reflecting different pathophysiological pathways may help to refine risk stratification in heart failure (HF). A novel calculator (BCN Bio-HF calculator) incorporating N-terminal pro B-type natriuretic peptide (NT-proBNP, a marker of myocardial stretch), high-sensitivity cardiac troponin T (hs-cTnT, a marker of myocyte injury), and high-sensitivity soluble ST2 (ST2), (reflective of myocardial fibrosis and remodeling) was developed. Methods: Model performance was evaluated using discrimination, calibration, and reclassification tools for 1-, 2-, and 3-year mortality. Ten-fold cross-validation with 1000 bootstrapping was used. Results: The BCN Bio-HF calculator was derived from 864 consecutive outpatients (72% men) with mean age 68.2±12 years (73%/27% New York Heart Association (NYHA) class I-II/III-IV, LVEF 36%, ischemic etiology 52.2%) and followed for a median of 3.4 years (305 deaths). After an initial evaluation of 23 variables, eight independent models were developed. The variables included in these models were age, sex, NYHA functional class, left ventricular ejection fraction, serum sodium, estimated glomerular filtration rate, hemoglobin, loop diuretic dose, β-blocker, Angiotensin converting enzyme inhibitor/Angiotensin-2 receptor blocker and statin treatments, and hs-cTnT, ST2, and NT-proBNP levels. The calculator may run with the availability of none, one, two, or the three biomarkers. The calculated risk of death was significantly changed by additive biomarker data. The average C-statistic in cross-validation analysis was 0.79. Conclusions: A new HF risk-calculator that incorporates available biomarkers reflecting different pathophysiological pathways better allowed individual prediction of death at 1, 2, and 3 years. [ABSTRACT FROM AUTHOR]

Vaccination is one of the most efficient ways to control the spread of infectious diseases. Simulations are now widely used to assess how vaccination can limit disease spread as well as mitigate morbidity or mortality in susceptible populations. However, field studies investigating how much vaccines decrease the velocity of epizootic wave-fronts during outbreaks are rare. This study aimed at investigating the effect of vaccination on the propagation of bluetongue, a vector-borne disease of ruminants. We used data from the 2008 bluetongue virus serotype 1 (BTV-1) epizootic of southwest France. As the virus was newly introduced in this area, natural immunity of livestock was absent. This allowed determination of the role of vaccination in changing the velocity of bluetongue spread while accounting for environmental factors that possibly influenced it. The average estimated velocity across the country despite restriction on animal movements was 5.4 km/day, which is very similar to the velocity of spread of the bluetongue virus serotype 8 epizootic in France also estimated in a context of restrictions on animal movements. Vaccination significantly reduced the propagation velocity of BTV-1. In comparison to municipalities with no vaccine coverage, the velocity of BTV-1 spread decreased by 1.7 km/day in municipalities with immunized animals. For the first time, the effect of vaccination has been quantified using data from a real epizootic whilst accounting for environmental factors known to modify the velocity of bluetongue spread. Our findings emphasize the importance of vaccination in limiting disease spread across natural landscape. Finally, environmental factors, specifically those related to vector abundance and activity, were found to be good predictors of the velocity of BTV-1 spread, indicating that these variables need to be adequately accounted for when evaluating the role of vaccination on bluetongue spread. [ABSTRACT FROM AUTHOR]

Background: Accelerometers are designed to measure plausible human activity, however extremely high count values (EHCV) have been recorded in large-scale studies. Using population data, we develop methodological principles for establishing an EHCV threshold, propose a threshold to define EHCV in the ActiGraph GT1M, determine occurrences of EHCV in a large-scale study, identify device-specific error values, and investigate the influence of varying EHCV thresholds on daily vigorous PA (VPA). Methods: We estimated quantiles to analyse the distribution of all accelerometer positive count values obtained from 9005 seven-year old children participating in the UK Millennium Cohort Study. A threshold to identify EHCV was derived by differentiating the quantile function. Data were screened for device-specific error count values and EHCV, and a sensitivity analysis conducted to compare daily VPA estimates using three approaches to accounting for EHCV. Results: Using our proposed threshold of ≥ 11,715 counts/minute to identify EHCV, we found that only 0.7% of all non-zero counts measured in MCS children were EHCV; in 99.7% of these children, EHCV comprised < 1% of total non-zero counts. Only 11 MCS children (0.12% of sample) returned accelerometers that contained negative counts; out of 237 such values, 211 counts were equal to −32,768 in one child. The medians of daily minutes spent in VPA obtained without excluding EHCV, and when using a higher threshold (≥19,442 counts/minute) were, respectively, 6.2% and 4.6% higher than when using our threshold (6.5 minutes; p<0.0001). Conclusions: Quality control processes should be undertaken during accelerometer fieldwork and prior to analysing data to identify monitors recording error values and EHCV. The proposed threshold will improve the validity of VPA estimates in children’s studies using the ActiGraph GT1M by ensuring only plausible data are analysed. These methods can be applied to define appropriate EHCV thresholds for different accelerometer models. [ABSTRACT FROM AUTHOR]

Background: In this study, we aim to identify key climatic factors that are associated with the transmission of Japanese encephalitis virus in areas located near the Three Gorges Dam, between 1997 and 2008. Methods: We identified three geographical regions of Chongqing, based on their distance from the Three Gorges Dam. Collectively, the three regions consisted of 12 districts from which study information was collected. Zero-Inflated Poisson Regression models were run to identify key climatic factors of the transmission of Japanese encephalitis virus for both the whole study area and for each individual region; linear regression models were conducted to examine the fluctuation of climatic variables over time during the construction of the Three Gorges Dam. Results: Between 1997 and 2008, the incidence of Japanese encephalitis decreased throughout the entire city of Chongqing, with noticeable variations taking place in 2000, 2001 and 2006. The eastern region, which is closest to the Three Gorges Dam, suffered the highest incidence of Japanese encephalitis, while the western region experienced the lowest incidence. Linear regression models revealed that there were seasonal fluctuations of climatic variables during this period. Zero-Inflated Poisson Regression models indicated a significant positive association between temperature (with a lag of 1 and 3 months) and Japanese encephalitis incidence, and a significant negative association between rainfall (with a lag of 0 and 4 months) and Japanese encephalitis incidence. Conclusion: The spatial and temporal trends of Japanese encephalitis incidence that occurred in the City of Chongqing were associated with temperature and rainfall. Seasonal fluctuations of climatic variables during this period were also observed. Additional studies that focus on long-term data collection are needed to validate the findings of this study and to further explore the effects of the Three Gorges Dam on Japanese encephalitis and other related diseases. [ABSTRACT FROM AUTHOR]

Background: The basic reproduction number, , is one of the many measures of the epidemic potential of an infection in a population. We estimate HIV over 18 years in a rural population in Uganda, examine method-specific differences in estimated , and estimate behavioural changes that would reduce below one. Methods: Data on HIV natural history and infectiousness were collated from literature. Data on new sexual partner count were available from a rural clinical cohort in Uganda over 1991–2008. was estimated using six methods. Behavioural changes required to reduce below one were calculated. Results: Reported number of new partners per year was 0 to 16 (women) and 0 to 80 (men). When proportionate sexual mixing was assumed, the different methods yielded comparable estimates. Assuming totally assortative mixing led to increased estimates in the high sexual activity class while all estimates in the low-activity class were below one. Using the “effective” partner change rate introduced by Anderson and colleagues resulted in estimates all above one except in the lowest sexual activity class. could be reduced below one if: (a) medium risk individuals reduce their partner acquisition rate by 70% and higher risk individuals reduce their partner acquisition rate by 93%, or (b) higher risk individuals reduce the partner acquisition rate by 95%. Conclusions: The estimated depended strongly on the method used. Ignoring variation in sexual activity leads to an underestimation of . Relying on behaviour change alone to eradicate HIV may require unrealistically large reductions in risk behaviour, even though for a small proportion of the population. To control HIV, complementary prevention strategies such as male circumcision and HIV treatment services need rapid scale up. [ABSTRACT FROM AUTHOR]

In the Susceptible–Infectious–Recovered (SIR) model of disease spreading, the time to extinction of the epidemics happens at an intermediate value of the per-contact transmission probability. Too contagious infections burn out fast in the population. Infections that are not contagious enough die out before they spread to a large fraction of people. We characterize how the maximal extinction time in SIR simulations on networks depend on the network structure. For example we find that the average distances in isolated components, weighted by the component size, is a good predictor of the maximal time to extinction. Furthermore, the transmission probability giving the longest outbreaks is larger than, but otherwise seemingly independent of, the epidemic threshold. [ABSTRACT FROM AUTHOR]

Objectives: Severe influenza can lead to Intensive Care Unit (ICU) admission. We explored whether ICU data reflect influenza like illness (ILI) activity in the general population, and whether ICU respiratory infections can predict influenza epidemics. Methods: We calculated the time lag and correlation between ILI incidence (from ILI sentinel surveillance, based on general practitioners (GP) consultations) and percentages of ICU admissions with a respiratory infection (from the Dutch National Intensive Care Registry) over the years 2003–2011. In addition, ICU data of the first three years was used to build three regression models to predict the start and end of influenza epidemics in the years thereafter, one to three weeks ahead. The predicted start and end of influenza epidemics were compared with observed start and end of such epidemics according to the incidence of ILI. Results: Peaks in respiratory ICU admissions lasted longer than peaks in ILI incidence rates. Increases in ICU admissions occurred on average two days earlier compared to ILI. Predicting influenza epidemics one, two, or three weeks ahead yielded positive predictive values ranging from 0.52 to 0.78, and sensitivities from 0.34 to 0.51. Conclusions: ICU data was associated with ILI activity, with increases in ICU data often occurring earlier and for a longer time period. However, in the Netherlands, predicting influenza epidemics in the general population using ICU data was imprecise, with low positive predictive values and sensitivities. [ABSTRACT FROM AUTHOR]

Background: Communicable disease outbreaks of novel or existing pathogens threaten human health around the globe. It would be desirable to rapidly characterize such outbreaks and develop accurate projections of their duration and cumulative size even when limited preliminary data are available. Here we develop a mathematical model to aid public health authorities in tracking the expansion and contraction of outbreaks with explicit representation of factors (other than population immunity) that may slow epidemic growth. Methodology: The Incidence Decay and Exponential Adjustment (IDEA) model is a parsimonious function that uses the basic reproduction number R0, along with a discounting factor to project the growth of outbreaks using only basic epidemiological information (e.g., daily incidence counts). Principal Findings: Compared to simulated data, IDEA provides highly accurate estimates of total size and duration for a given outbreak when R0 is low or moderate, and also identifies turning points or new waves. When tested with an outbreak of pandemic influenza A (H1N1), the model generates estimated incidence at the i+1th serial interval using data from the ith serial interval within an average of 20% of actual incidence. Conclusions and Significance: This model for communicable disease outbreaks provides rapid assessments of outbreak growth and public health interventions. Further evaluation in the context of real-world outbreaks will establish the utility of IDEA as a tool for front-line epidemiologists. [ABSTRACT FROM AUTHOR]

Understanding the biological mechanisms underlying episodic outbreaks of infectious diseases is one of mathematical epidemiology’s major goals. Historic records are an invaluable source of information in this enterprise. Pertussis (whooping cough) is a re-emerging infection whose intermittent bouts of large multiannual epidemics interspersed between periods of smaller-amplitude cycles remain an enigma. It has been suggested that recent increases in pertussis incidence and shifts in the age-distribution of cases may be due to diminished natural immune boosting. Here we show that a model that incorporates this mechanism can account for a unique set of pre-vaccine-era data from Copenhagen. Under this model, immune boosting induces transient bursts of large amplitude outbreaks. In the face of mass vaccination, the boosting model predicts larger and more frequent outbreaks than do models with permanent or passively-waning immunity. Our results emphasize the importance of understanding the mechanisms responsible for maintaining immune memory for pertussis epidemiology. [ABSTRACT FROM AUTHOR]

Background: To evaluate the pattern of co-infection of human papillomavirus (HPV) types in both sexes in Sweden. Methods: Cell samples from genital swabs, first-void urine, and genital swabs immersed in first-void urine were collected in the present cross-sectional High Throughput HPV Monitoring study. Overall, 31,717 samples from women and 9,949 from men (mean age 25) were tested for 16 HPV types using mass spectrometry. Multilevel logistic regression was used to estimate the expected number of multiple infections with specific HPV types, adjusted for age, type of sample, and accounting for correlations between HPV types due to unobserved risk factors using sample-level random effects. Bonferroni correction was used to allow for multiple comparisons (120). Results: Observed-to-expected ratio for any multiple infections was slightly above unity in both sexes, but, for most 2-type combinations, there was no evidence of significant departure from expected numbers. HPV6/18 was found more often and HPV51/68 and 6/68 less often than expected. However, HPV68 tended to be generally underrepresented in co-infections, suggesting a sub-optimal performance of our testing method for this HPV type. Conclusions: We found no evidence for positive or negative clustering between HPV types included in the current prophylactic vaccines and other untargeted oncogenic types, in either sex. [ABSTRACT FROM AUTHOR]

Background: Health inequities in developing countries are difficult to eradicate because of limited resources. The neglect of adult mortality in Sub-Saharan Africa (SSA) is a particular concern. Advances in data availability, software and analytic methods have created opportunities to address this challenge and tailor interventions to small areas. This study demonstrates how a generic framework can be applied to guide policy interventions to reduce adult mortality in high risk areas. The framework, therefore, incorporates the spatial clustering of adult mortality, estimates the impact of a range of determinants and quantifies the impact of their removal to ensure optimal returns on scarce resources. Methods: Data from a national cross-sectional survey in 2007 were used to illustrate the use of the generic framework for SSA and elsewhere. Adult mortality proportions were analyzed at four administrative levels and spatial analyses were used to identify areas with significant excess mortality. An ecological approach was then used to assess the relationship between mortality “hotspots” and various determinants. Population attributable fractions were calculated to quantify the reduction in mortality as a result of targeted removal of high-impact determinants. Results: Overall adult mortality rate was 145 per 10,000. Spatial disaggregation identified a highly non-random pattern and 67 significant high risk local municipalities were identified. The most prominent determinants of adult mortality included HIV antenatal sero-prevalence, low SES and lack of formal marital union status. The removal of the most attributable factors, based on local area prevalence, suggest that overall adult mortality could be potentially reduced by ∼90 deaths per 10,000. Conclusions: The innovative use of secondary data and advanced epidemiological techniques can be combined in a generic framework to identify and map mortality to the lowest administration level. The identification of high risk mortality determinants allows health authorities to tailor interventions at local level. This approach can be replicated elsewhere. [ABSTRACT FROM AUTHOR]

The spread of certain diseases can be promoted, in some cases substantially, by prior infection with another disease. One example is that of HIV, whose immunosuppressant effects significantly increase the chances of infection with other pathogens. Such coinfection processes, when combined with nontrivial structure in the contact networks over which diseases spread, can lead to complex patterns of epidemiological behavior. Here we consider a mathematical model of two diseases spreading through a single population, where infection with one disease is dependent on prior infection with the other. We solve exactly for the sizes of the outbreaks of both diseases in the limit of large population size, along with the complete phase diagram of the system. Among other things, we use our model to demonstrate how diseases can be controlled not only by reducing the rate of their spread, but also by reducing the spread of other infections upon which they depend. [ABSTRACT FROM AUTHOR]

Background: Advancements in knowledge of obesity aetiology and mobile phone technology have created the opportunity to develop an electronic tool to predict an infant’s risk of childhood obesity. The study aims were to develop and validate equations for the prediction of childhood obesity and integrate them into a mobile phone application (App). Methods and Findings: Anthropometry and childhood obesity risk data were obtained for 1868 UK-born White or South Asian infants in the Born in Bradford cohort. Logistic regression was used to develop prediction equations (at 6±1.5, 9±1.5 and 12±1.5 months) for risk of childhood obesity (BMI at 2 years >91st centile and weight gain from 0–2 years >1 centile band) incorporating sex, birth weight, and weight gain as predictors. The discrimination accuracy of the equations was assessed by the area under the curve (AUC); internal validity by comparing area under the curve to those obtained in bootstrapped samples; and external validity by applying the equations to an external sample. An App was built to incorporate six final equations (two at each age, one of which included maternal BMI). The equations had good discrimination (AUCs 86–91%), with the addition of maternal BMI marginally improving prediction. The AUCs in the bootstrapped and external validation samples were similar to those obtained in the development sample. The App is user-friendly, requires a minimum amount of information, and provides a risk assessment of low, medium, or high accompanied by advice and website links to government recommendations. Conclusions: Prediction equations for risk of childhood obesity have been developed and incorporated into a novel App, thereby providing proof of concept that childhood obesity prediction research can be integrated with advancements in technology. [ABSTRACT FROM AUTHOR]

Introduction: Statistical interactions are a common component of data analysis across a broad range of scientific disciplines. However, the statistical power to detect interactions is often undesirably low. One solution is to elevate the Type 1 error rate so that important interactions are not missed in a low power situation. To date, no study has quantified the effects of this practice on power in a linear regression model. Methods: A Monte Carlo simulation study was performed. A continuous dependent variable was specified, along with three types of interactions: continuous variable by continuous variable; continuous by dichotomous; and dichotomous by dichotomous. For each of the three scenarios, the interaction effect sizes, sample sizes, and Type 1 error rate were varied, resulting in a total of 240 unique simulations. Results: In general, power to detect the interaction effect was either so low or so high at α = 0.05 that raising the Type 1 error rate only served to increase the probability of including a spurious interaction in the model. A small number of scenarios were identified in which an elevated Type 1 error rate may be justified. Conclusions: Routinely elevating Type 1 error rate when testing interaction effects is not an advisable practice. Researchers are best served by positing interaction effects a priori and accounting for them when conducting sample size calculations. [ABSTRACT FROM AUTHOR]

Objectives: To examine the factor structure and to evaluate the longitudinal measurement invariance of the demand-control-support questionnaire (DCSQ), using the Swedish Longitudinal Occupational Survey of Health (SLOSH). Methods: A confirmatory factor analysis (CFA) and multi-group confirmatory factor analysis (MGCFA) models within the framework of structural equation modeling (SEM) have been used to examine the factor structure and invariance across time. Results: Four factors: psychological demand, skill discretion, decision authority and social support, were confirmed by CFA at baseline, with the best fit obtained by removing the item repetitive work of skill discretion. A measurement error correlation (0.42) between work fast and work intensively for psychological demands was also detected. Acceptable composite reliability measures were obtained except for skill discretion (0.68). The invariance of the same factor structure was established, but caution in comparing mean levels of factors over time is warranted as lack of intercept invariance was evident. However, partial intercept invariance was established for work intensively. Conclusion: Our findings indicate that skill discretion and decision authority represent two distinct constructs in the retained model. However removing the item repetitive work along with either work fast or work intensively would improve model fit. Care should also be taken while making comparisons in the constructs across time. Further research should investigate invariance across occupations or socio-economic classes. [ABSTRACT FROM AUTHOR]

Background: Accurate estimates of the incidence and prevalence of methicillin-resistant Staphylococcus aureus (MRSA) infections are needed to inform public health policies. In Norway, where both MRSA infection and carriage are notifiable conditions, the reported incidence of MRSA is slowly increasing. However, the proportion of MRSA in relation to all S. aureus isolates is unknown, making it difficult to determine if the rising incidence is real or an artifact of an increasing number of tests performed. Aim: To characterize recent trends in MRSA infections and obtain a more complete understanding of the MRSA level in Norway. Methods: All reported cases of MRSA and methicillin-sensitive S. aureus (MSSA) from Oslo County (1997–2010) and Health Region East (2008–2008), representing approximately 11% and 36% of the Norwegian population, respectively, were analyzed using a stochastic time series analysis to characterize trends. Results: In Oslo County, the proportion of methicillin-resistant cases increased from 0.73% to 3.78% during the study period and was well modeled by an exponential growth with a doubling constant of 5.7 years (95% CI 4.5–7.4 years). In Health Region East, the proportion of MRSA cases increased from 0.4% to 2.1% from 2002 to 2008, with a best-fitting linear increase of 0.26% (95% CI 0.21–0.30%) per year. In both cases, the choice of a linear or exponential model for the time trend produced only marginally different model fits. We found no significant changes due to revised national MRSA guidelines published in June 2009. Significant variations in the increasing time trend were observed in the five hospitals within the region. The yearly reported incidence of MSSA was relatively stable in both study areas although we found seasonal patterns with peaks in August. Conclusion: The level of MRSA is increasing in Norway, and the proportion of methicillin resistance in all S. aureus isolates are higher than the reported proportion of MRSA in invasive infections. [ABSTRACT FROM AUTHOR]

Background: Temperature and humidity strongly affect the physiology, longevity, fecundity and dispersal behavior of Aedes aegypti, vector of dengue fever. Contrastingly, the statistical associations measured between time series of mosquito abundance and meteorological variables are often weak and contradictory. Here, we investigated the significance of these relationships at different time scales. Methods and Findings: A time series of the adult mosquito abundance from a medium-sized city in Brazil, lasting 109 weeks was analyzed. Meteorological variables included temperature, precipitation, wind velocity and humidity. As analytical tools, generalized linear models (GLM) with time lags and interaction terms were used to identify average effects while the wavelet analysis was complementarily used to identify transient associations. The fitted GLM showed that mosquito abundance is significantly affected by the interaction between lagged temperature and humidity, and also by the mosquito abundance a week earlier. Extreme meteorological variables were the best predictors, and the mosquito population tended to increase at values above and 54% humidity. The wavelet analysis identified non-stationary local effects of these meteorological variables on abundance throughout the study period, with peaks in the spring-summer period. The wavelet detected weak but significant effects for precipitation and wind velocity. Conclusion: Our results support the presence of transient relationships between meteorological variables and mosquito abundance. Such transient association may be explained by the ability of Ae. aegypti to buffer part of its response to climate, for example, by choosing sites with proper microclimate. We also observed enough coupling between the abundance and meteorological variables to develop a model with good predictive power. Extreme values of meteorological variables with time lags, interaction terms and previous mosquito abundance are strong predictors and should be considered when understanding the climate effect on mosquito abundance and population growth. [ABSTRACT FROM AUTHOR]

Malaria and lymphatic filariasis (LF) continue to cause a considerable public health burden globally and are co-endemic in many regions of sub-Saharan Africa. These infections are transmitted by the same mosquito species which raises important questions about optimal vector control strategies in co-endemic regions, as well as the effect of the presence of each infection on endemicity of the other; there is currently little consensus on the latter. The need for comprehensive modelling studies to address such questions is therefore significant, yet very few have been undertaken to date despite the recognised explanatory power of reliable dynamic mathematical models. Here, we develop a malaria-LF co-infection modelling framework that accounts for two key interactions between these infections, namely the increase in vector mortality as LF mosquito prevalence increases and the antagonistic Th1/Th2 immune response that occurs in co-infected hosts. We consider the crucial interplay between these interactions on the resulting endemic prevalence when introducing each infection in regions where the other is already endemic (e.g. due to regional environmental change), and the associated timescale for such changes, as well as effects on the basic reproduction number R0 of each disease. We also highlight potential perverse effects of vector controls on human infection prevalence in co-endemic regions, noting that understanding such effects is critical in designing optimal integrated control programmes. Hence, as well as highlighting where better data are required to more reliably address such questions, we provide an important framework that will form the basis of future scenario analysis tools used to plan and inform policy decisions on intervention measures in different transmission settings. [ABSTRACT FROM AUTHOR]