Incompletely resolved phylogenetic trees inflate estimates of phylogenetic conservatism.

The tendency for more closely related species to share similar traits and ecological strategies can be explained by their longer shared evolutionary histories and represents phylogenetic conservatism. How strongly species traits co-vary with phylogeny can significantly impact how we analyze cross-species data and can influence our interpretation of assembly rules in the rapidly expanding field of community phylogenetics. Phylogenetic conservatism is typically quantified by analyzing the distribution of species values on the phylogenetic tree that connects them. Many phylogenetic approaches, however, assume a completely sampled phylogeny: while we have good estimates of deeper phylogenetic relationships for many species-rich groups, such as birds and flowering plants, we often lack information on more recent interspecific relationships (i.e., within a genus). A common solution has been to represent these relationships as polytomies on trees using taxonomy as a guide. Here we show that such trees can dramatically inflate estimates of phylogenetic conservatism quantified using S. P. Blomberg et al.'s K statistic. Using simulations, we show that even randomly generated traits can appear to be phylogenetically conserved on poorly resolved trees. We provide a simple rarefaction-based solution that can reliably retrieve unbiased estimates of K, and we illustrate our method using data on first flowering times from Thoreau's woods (Concord, Massachusetts, USA).

Missing value imputation in longitudinal measures of alcohol consumption

Attrition in longitudinal studies can lead to biased results. The study is motivated by the unexpected observation that alcohol consumption decreased despite increased availability, which may be due to sample attrition of heavy drinkers. Several imputation methods have been proposed, but rarely compared in longitudinal studies of alcohol consumption. The imputation of consumption level measurements is computationally particularly challenging due to alcohol consumption being a semi-continuous variable (dichotomous drinking status and continuous volume among drinkers), and the non-normality of data in the continuous part. Data come from a longitudinal study in Denmark with four waves (2003-2006) and 1771 individuals at baseline. Five techniques for missing data are compared: Last value carried forward (LVCF) was used as a single, and Hotdeck, Heckman modelling, multivariate imputation by chained equations (MICE), and a Bayesian approach as multiple imputation methods. Predictive mean matching was used to account for non-normality, where instead of imputing regression estimates, "real" observed values from similar cases are imputed. Methods were also compared by means of a simulated dataset. The simulation showed that the Bayesian approach yielded the most unbiased estimates for imputation. The finding of no increase in consumption levels despite a higher availability remained unaltered. Copyright (C) 2011 John Wiley & Sons, Ltd.

Bayesian estimation of speciation and extinction from incomplete fossil occurrence data.

The temporal dynamics of species diversity are shaped by variations in the rates of speciation and extinction, and there is a long history of inferring these rates using first and last appearances of taxa in the fossil record. Understanding diversity dynamics critically depends on unbiased estimates of the unobserved times of speciation and extinction for all lineages, but the inference of these parameters is challenging due to the complex nature of the available data. Here, we present a new probabilistic framework to jointly estimate species-specific times of speciation and extinction and the rates of the underlying birth-death process based on the fossil record. The rates are allowed to vary through time independently of each other, and the probability of preservation and sampling is explicitly incorporated in the model to estimate the true lifespan of each lineage. We implement a Bayesian algorithm to assess the presence of rate shifts by exploring alternative diversification models. Tests on a range of simulated data sets reveal the accuracy and robustness of our approach against violations of the underlying assumptions and various degrees of data incompleteness. Finally, we demonstrate the application of our method with the diversification of the mammal family Rhinocerotidae and reveal a complex history of repeated and independent temporal shifts of both speciation and extinction rates, leading to the expansion and subsequent decline of the group. The estimated parameters of the birth-death process implemented here are directly comparable with those obtained from dated molecular phylogenies. Thus, our model represents a step towards integrating phylogenetic and fossil information to infer macroevolutionary processes.

Estimates of radiation-absorbed dose to kidneys in patients treated with 90Y-ibritumomab tiuxetan.

The aim of the present study was to retrospectively estimate the absorbed dose to kidneys in 17 patients treated in clinical practice with 90Y-ibritumomab tiuxetan for non-Hodgkin's lymphoma, using appropriate dosimetric approaches available. METHODS: The single-view effective point source method, including background subtraction, is used for planar quantification of renal activity. Since the high uptake in the liver affects the activity estimate in the right kidney, the dose to the left kidney serves as a surrogate for the dose to both kidneys. Calculation of absorbed dose is based on the Medical Internal Radiation Dose methodology with adjustment for patient kidney mass. RESULTS: The median dose to kidneys, based on the left kidney only, is 2.1 mGy/MBq (range, 0.92-4.4), whereas a value of 2.5 mGy/MBq (range, 1.5-4.7) is obtained, considering the activity in both kidneys. CONCLUSIONS: Irrespective of the method, doses to kidneys obtained in the present study were about 10 times higher than the median dose of 0.22 mGy/MBq (range, 0.00-0.95) were originally reported from the study leading to Food and Drug Administration approval. Our results are in good agreement with kidney-dose estimates recently reported from high-dose myeloablative therapy with 90Y-ibritumomab tiuxetan.

Permeability estimates from poro-elastic dispersion analyses of sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Inter-informant agreement and prevalence estimates for substance use disorders: direct interview versus family history method.

OBJECTIVES: Family studies typically use multiple sources of information on each individual including direct interviews and family history information. The aims of the present study were to: (1) assess agreement for diagnoses of specific substance use disorders between direct interviews and the family history method; (2) compare prevalence estimates according to the two methods; (3) test strategies to approximate prevalence estimates according to family history reports to those based on direct interviews; (4) determine covariates of inter-informant agreement; and (5) identify covariates that affect the likelihood of reporting disorders by informants. METHODS: Analyses were based on family study data which included 1621 distinct informant (first-degree relatives and spouses) - index subject pairs. RESULTS: Our main findings were: (1) inter-informant agreement was fair to good for all substance disorders, except for alcohol abuse; (2) the family history method underestimated the prevalence of drug but not alcohol use disorders; (3) lowering diagnostic thresholds for drug disorders and combining multiple family histories increased the accuracy of prevalence estimates for these disorders according to the family history method; (4) female sex of index subjects was associated with higher agreement for nearly all disorders; and (5) informants who themselves had a history of the same substance use disorder were more likely to report this disorder in their relatives, which entails the risk of overestimation of the size of familial aggregation. CONCLUSION: Our findings have important implications for the best-estimate procedure applied in family studies.

Quantifying uncertainty: physicians' estimates of infection in critically ill neonates and children.

To determine the diagnostic accuracy of physicians' prior probability estimates of serious infection in critically ill neonates and children, we conducted a prospective cohort study in 2 intensive care units. Using available clinical, laboratory, and radiographic information, 27 physicians provided 2567 probability estimates for 347 patients (follow-up rate, 92%). The median probability estimate of infection increased from 0% (i.e., no antibiotic treatment or diagnostic work-up for sepsis), to 2% on the day preceding initiation of antibiotic therapy, to 20% at initiation of antibiotic treatment (P<.001). At initiation of treatment, predictions discriminated well between episodes subsequently classified as proven infection and episodes ultimately judged unlikely to be infection (area under the curve, 0.88). Physicians also showed a good ability to predict blood culture-positive sepsis (area under the curve, 0.77). Treatment and testing thresholds were derived from the provided predictions and treatment rates. Physicians' prognoses regarding the presence of serious infection were remarkably precise. Studies investigating the value of new tests for diagnosis of sepsis should establish that they add incremental value to physicians' judgment.

Parent-child agreement and prevalence estimates of diagnoses in childhood: direct interview versus family history method.

Diagnostic information on children is typically elicited from both children and their parents. The aims of the present paper were to: (1) compare prevalence estimates according to maternal reports, paternal reports and direct interviews of children [major depressive disorder (MDD), anxiety and attention-deficit and disruptive behavioural disorders]; (2) assess mother-child, father-child and inter-parental agreement for these disorders; (3) determine the association between several child, parent and familial characteristics and the degree of diagnostic agreement or the likelihood of parental reporting; (4) determine the predictive validity of diagnostic information provided by parents and children. Analyses were based on 235 mother-offspring, 189 father-offspring and 128 mother-father pairs. Diagnostic assessment included the Kiddie-schedule for Affective Disorders and Schizophrenia (K-SADS) (offspring) and the Diagnostic Interview for Genetic Studies (DIGS) (parents and offspring at follow-up) interviews. Parental reports were collected using the Family History - Research Diagnostic Criteria (FH-RDC). Analyses revealed: (1) prevalence estimates for internalizing disorders were generally lower according to parental information than according to the K-SADS; (2) mother-child and father-child agreement was poor and within similar ranges; (3) parents with a history of MDD or attention deficit hyperactivity disorder (ADHD) reported these disorders in their children more frequently; (4) in a sub-sample followed-up into adulthood, diagnoses of MDD, separation anxiety and conduct disorder at baseline concurred with the corresponding lifetime diagnosis at age 19 according to the child rather than according to the parents. In conclusion, our findings support large discrepancies of diagnostic information provided by parents and children with generally lower reporting of internalizing disorders by parents, and differential reporting of depression and ADHD by parental disease status. Follow-up data also supports the validity of information provided by adolescent offspring.

Physicians' estimates of the 10 year cardiovascular risk in hypertensive patients: an evaluation in primary care physicians in training.

To evaluate how young physicians in training perceive their patients' cardiovascular risk based on the medical charts and their clinical judgment. Cross sectional observational study. University outpatient clinic, Lausanne, Switzerland. Two hundred hypertensive patients and 50 non-hypertensive patients with at least one cardiovascular risk factor. Comparison of the absolute 10-year cardiovascular risk calculated by a computer program based on the Framingham score and adapted for physicians by the WHO/ISH with the perceived risk as assessed clinically by the physicians. Physicians underestimated the 10-year cardiovascular risk of their patients compared to that calculated with the Framingham score. Concordance between methods was 39% for hypertensive patients and 30% for non-hypertensive patients. Underestimation of cardiovascular risks for hypertensive patients was related to the fact they had a stabilized systolic blood pressure under 140 mm Hg (OR = 2.1 [1.1; 4.1]). These data show that young physicians in training often have an incorrect perception of the cardiovascular risk of their patients with a tendency to underestimate the risk. However, the calculated risk could also be slightly overestimated when applying the Framingham Heart Study model to a Swiss population. To implement a systematic evaluation of risk factors in primary care a greater emphasis should be placed on the teaching of cardiovascular risk evaluation and on the implementation of quality improvement programs.

Inter-informant agreement and prevalence estimates for mood syndromes: Direct interview vs. family history method.

BACKGROUND: The use of the family history method is recommended in family studies as a type of proxy interview of non-participating relatives. However, using different sources of information can result in bias as direct interviews may provide a higher likelihood of assigning diagnoses than family history reports. The aims of the present study were to: 1) compare diagnoses for threshold and subthreshold mood syndromes from interviews to those relying on information from relatives; 2) test the appropriateness of lowering the diagnostic threshold and combining multiple reports from the family history method to obtain comparable prevalence estimates to the interviews; 3) identify factors that influence the likelihood of agreement and reporting of disorders by informants. METHODS: Within a family study, 1621 informant-index subject pairs were identified. DSM-5 diagnoses from direct interviews of index subjects were compared to those derived from family history information provided by their first-degree relatives. RESULTS: 1) Inter-informant agreement was acceptable for Mania, but low for all other mood syndromes. 2) Except for Mania and subthreshold depression, the family history method provided significantly lower prevalence estimates. The gap improved for all other syndromes after lowering the threshold of the family history method. 3) Individuals who had a history of depression themselves were more likely to report depression in their relatives. LIMITATIONS: Low proportion of affected individuals for manic syndromes and lack of independence of data. CONCLUSIONS: The higher likelihood of reporting disorders by affected informants entails the risk of overestimation of the size of familial aggregation of depression.

Optimal robust estimates using the Hellinger distance

Optimal robust M-estimates of a multidimensional parameter are described using Hampel's infinitesimal approach. The optimal estimates are derived by minimizing a measure of efficiency under the model, subject to a bounded measure of infinitesimal robustness. To this purpose we define measures of efficiency and infinitesimal sensitivity based on the Hellinger distance.We show that these two measures coincide with similar ones defined by Yohai using the Kullback-Leibler divergence, and therefore the corresponding optimal estimates coincide too.We also give an example where we fit a negative binomial distribution to a real dataset of "days of stay in hospital" using the optimal robust estimates.

Linking temperature estimates and microstructures in deformed polymineralic mantle rocks

To constrain deformation temperatures of mantle shear zones, we studied a strike-slip shear zone (Hilti massif, Semail ophiolite, Oman) and focused on the interaction between microstructural mechanisms and chemical equilibration processes. Quantitative microfabric analysis on harzburgites with different deformation intensity (porphyroclastic tectonite, mylonite, and ultramylonite) was combined with orthopyroxene geothermometry. The average grain size of all phases decreases with decreasing shear zone thickness. Dynamic recrystallization of porphyroclasts in combination with dissolution-precipitation and nucleation result in small-sized, chemically equilibrated pyroxenes. The composition of orthopyroxene was used to calculate deformation temperatures. In the case of the porphyroclastic tectonites, the chemical composition of orthopyroxene has been reset by diffusion yielding temperature estimates of 880-900 degrees C. The mylonites were deformed by dislocation creep of olivine and show a broad range of calculated temperatures, which result from a combination of grain size reduction and inheritance of equilibrium compositions from earlier high-temperature events and diffusion. In mylonites, diffusion profiles combined with geothermometry and grain size analysis indicate a mylonitic deformation temperature of 800-900 degrees C possibly followed by diffusion. In ultramylonites, the smallest grains (<30 mu m) reveal equilibration at temperatures of similar to 700 degrees C during the last stages of ductile deformation, which was dominated by diffusion creep of olivine. Our results provide a crucial link between temperature and evolution of microstructures from dislocation creep to diffusion creep in mantle shear zones.

Determining the best population-level alcohol consumption model and its impact on estimates of alcohol-attributable harms.

BACKGROUND: The goals of our study are to determine the most appropriate model for alcohol consumption as an exposure for burden of disease, to analyze the effect of the chosen alcohol consumption distribution on the estimation of the alcohol Population- Attributable Fractions (PAFs), and to characterize the chosen alcohol consumption distribution by exploring if there is a global relationship within the distribution. METHODS: To identify the best model, the Log-Normal, Gamma, and Weibull prevalence distributions were examined using data from 41 surveys from Gender, Alcohol and Culture: An International Study (GENACIS) and from the European Comparative Alcohol Study. To assess the effect of these distributions on the estimated alcohol PAFs, we calculated the alcohol PAF for diabetes, breast cancer, and pancreatitis using the three above-named distributions and using the more traditional approach based on categories. The relationship between the mean and the standard deviation from the Gamma distribution was estimated using data from 851 datasets for 66 countries from GENACIS and from the STEPwise approach to Surveillance from the World Health Organization. RESULTS: The Log-Normal distribution provided a poor fit for the survey data, with Gamma and Weibull distributions providing better fits. Additionally, our analyses showed that there were no marked differences for the alcohol PAF estimates based on the Gamma or Weibull distributions compared to PAFs based on categorical alcohol consumption estimates. The standard deviation of the alcohol distribution was highly dependent on the mean, with a unit increase in alcohol consumption associated with a unit increase in the mean of 1.258 (95% CI: 1.223 to 1.293) (R2 = 0.9207) for women and 1.171 (95% CI: 1.144 to 1.197) (R2 = 0. 9474) for men. CONCLUSIONS: Although the Gamma distribution and the Weibull distribution provided similar results, the Gamma distribution is recommended to model alcohol consumption from population surveys due to its fit, flexibility, and the ease with which it can be modified. The results showed that a large degree of variance of the standard deviation of the alcohol consumption Gamma distribution was explained by the mean alcohol consumption, allowing for alcohol consumption to be modeled through a Gamma distribution using only average consumption.

Global alcohol exposure estimates by country, territory and region for 2005--a contribution to the Comparative Risk Assessment for the 2010 Global Burden of Disease Study.

AIMS: This study aimed to estimate the prevalence of life-time abstainers, former drinkers and current drinkers, adult per-capita consumption of alcohol and pattern of drinking scores, by country and Global Burden of Disease region for 2005, and to forecast these indicators for 2010. DESIGN: Statistical modelling based on survey data and routine statistics. SETTING AND PARTICIPANTS: A total of 241 countries and territories. MEASUREMENTS: Per-capita consumption data were obtained with the help of the World Health Organization's Global Information System on Alcohol and Health. Drinking status data were obtained from Gender, Alcohol and Culture: An International Study, the STEPwise approach to Surveillance study, the World Health Survey/Multi-Country Study and other surveys. Consumption and drinking status data were triangulated to estimate alcohol consumption across multiple categories. FINDINGS: In 2005 adult per-capita annual consumption of alcohol was 6.1 litres, with 1.7 litres stemming from unrecorded consumption; 17.1 litres of alcohol were consumed per drinker, 45.8% of all adults were life-time abstainers, 13.6% were former drinkers and 40.6% were current drinkers. Life-time abstention was most prevalent in North Africa/Middle East and South Asia. Eastern Europe and Southern sub-Saharan Africa had the most detrimental pattern of drinking scores, while drinkers in Europe (Eastern and Central) and sub-Saharan Africa (Southern and West) consumed the most alcohol. CONCLUSIONS: Just over 40% of the world's adult population consumes alcohol and the average consumption per drinker is 17.1 litres per year. However, the prevalence of abstention, level of alcohol consumption and patterns of drinking vary widely across regions of the world.

Low estimates of dietary acid load are positively associated with bone ultrasound in women older than 75 years of age with a lifetime fracture

Dietary acid load from Western diets may be a risk factor for osteoporosis. It can be estimated by net endogenous acid production (NEAP). No data currently exists for NEAP estimates and bone indices in the very elderly (i.e. > or = 75 y). The aim of this study was to determine the association between NEAP estimates by using the potential renal acid load (PRAL) equation and quantitative bone ultrasound (QUS) measurements at the heel [broadband ultrasound attenuation (BUA)] in Caucasian women. We assessed NEAP and QUS in 401 very elderly Swiss ambulatory women. We evaluated dietary intake and NEAP estimates with a validated FFQ. QUS was measured using Achilles (Lunar). We identified 2 subgroups: 256 women (80.6 y +/- 3; BUA, 96.8 dB/MHz) with a fracture history and the remaining 145 (79.9 y SD 2.9; BUA, 101.7 dB/MHz) without. Women who reported having suffered a fracture had lower BUA (P < 0.001) than nonfractured women but did not differ in nutrient intakes and NEAP. Lower NEAP (P = 0.023) and higher potassium intake (P = 0.033) were correlated with higher BUA, which remained significant even after adjustment for age, BMI, and osteoporosis treatment. BUA was positively correlated with calcium (P = 0.016) and BMI (P < 0.001). Women who reported no fractures had no significant correlations between nutrient intake, NEAP, and BUA. Low nutritional acid load was correlated with higher BUA in very elderly women with a fracture history. Although relatively weak compared with age and BMI, this association was significant and may be an important additional risk factor that might be particularly relevant in frail patients with an already high fracture risk.