Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Prediction of postnatal outcome in fetuses with pulmonary hypoplasia by fetal-MRI: Comparison between lung volumetry and pulmonary signal intensity values

Purpose: Pulmonary hypoplasia is a determinant parameter for extra-uterine life. In the last years, MRI appears as a complement to US in order to evaluate the degree of pulmonary hypoplasia in foetuses with congenital anomalies, by using different methods - fetal lung volumetry (FLV), lung-to-liver signal intensity ratio (LLSIR)-. But until now, information about the correlation between the MRI prediction and the real postnatal outcome is limited. Methods and materials: We retrospectively reviewed the fetal MRI performed at our Institution in the last 8 years and selected the cases with suspicion of fetal pulmonary hypoplasia (n = 30). The pulmonary volumetry data of these foetuses were collected and the lung-to-liver signal intensity ratio (LLSIR) measures performed. These data were compared with those obtained from a control group of 25 foetuses considered as normal at MRI. The data of the study group were also correlated with the autopsy records or the post-natal clinical information of the patients. Results: As expected, the control group showed higher FLV and LLSIR values than the problem group at all gestational ages. Higher values of FLV and LLSIR were associated with a better post-natal outcome. Sensitivity, specificity, positive and negative predictive values and accuracy for the relative LLSIR and the relative FLV showed no significant differences. Conclusion: Our data show that not only the FLV but also the relative LLSIR inform about the degree of fetal lung development. This information may help to predict the fetal outcome and to evaluate the need for neonatal intensive care.

A scaled difference chi-square test statistic for moment structure analysis

A family of scaling corrections aimed to improve the chi-square approximation of goodness-of-fit test statistics in small samples, large models, and nonnormal data was proposed in Satorra and Bentler (1994). For structural equations models, Satorra-Bentler's (SB) scaling corrections are available in standard computer software. Often, however, the interest is not on the overall fit of a model, but on a test of the restrictions that a null model say ${\cal M}_0$ implies on a less restricted one ${\cal M}_1$. If $T_0$ and $T_1$ denote the goodness-of-fit test statistics associated to ${\cal M}_0$ and ${\cal M}_1$, respectively, then typically the difference $T_d = T_0 - T_1$ is used as a chi-square test statistic with degrees of freedom equal to the difference on the number of independent parameters estimated under the models ${\cal M}_0$ and ${\cal M}_1$. As in the case of the goodness-of-fit test, it is of interest to scale the statistic $T_d$ in order to improve its chi-square approximation in realistic, i.e., nonasymptotic and nonnormal, applications. In a recent paper, Satorra (1999) shows that the difference between two Satorra-Bentler scaled test statistics for overall model fit does not yield the correct SB scaled difference test statistic. Satorra developed an expression that permits scaling the difference test statistic, but his formula has some practical limitations, since it requires heavy computations that are notavailable in standard computer software. The purpose of the present paper is to provide an easy way to compute the scaled difference chi-square statistic from the scaled goodness-of-fit test statistics of models ${\cal M}_0$ and ${\cal M}_1$. A Monte Carlo study is provided to illustrate the performance of the competing statistics.

Genome-Wide Prediction Methods in Highly Diverse and Heterozygous Species: Proof-of-Concept through Simulation in Grapevine.

Nowadays, genome-wide association studies (GWAS) and genomic selection (GS) methods which use genome-wide marker data for phenotype prediction are of much potential interest in plant breeding. However, to our knowledge, no studies have been performed yet on the predictive ability of these methods for structured traits when using training populations with high levels of genetic diversity. Such an example of a highly heterozygous, perennial species is grapevine. The present study compares the accuracy of models based on GWAS or GS alone, or in combination, for predicting simple or complex traits, linked or not with population structure. In order to explore the relevance of these methods in this context, we performed simulations using approx 90,000 SNPs on a population of 3,000 individuals structured into three groups and corresponding to published diversity grapevine data. To estimate the parameters of the prediction models, we defined four training populations of 1,000 individuals, corresponding to these three groups and a core collection. Finally, to estimate the accuracy of the models, we also simulated four breeding populations of 200 individuals. Although prediction accuracy was low when breeding populations were too distant from the training populations, high accuracy levels were obtained using the sole core-collection as training population. The highest prediction accuracy was obtained (up to 0.9) using the combined GWAS-GS model. We thus recommend using the combined prediction model and a core-collection as training population for grapevine breeding or for other important economic crops with the same characteristics.

A simple randomized algorithm for consistent sequential prediction of ergodic time series

We present a simple randomized procedure for the prediction of a binary sequence. The algorithm uses ideas from recent developments of the theory of the prediction of individual sequences. We show that if thesequence is a realization of a stationary and ergodic random process then the average number of mistakes converges, almost surely, to that of the optimum, given by the Bayes predictor.

On the structural difference between the evolutionary approach of Young and that of Kandori, Mailath and Rob

We provide robust examples of symmetric two-player coordination games in normal form that reveal that equilibrium selection bythe evolutionary model of Young (1993) is essentially different from equilibrium selection by the evolutionary model of Kandori, Mailath and Rob (1993).

Prediction of functional outcome 18 months after a first psychotic episode: a proton magnetic resonance spectroscopy study.

CONTEXT: Recent magnetic resonance imaging studies have attempted to relate volumetric brain measurements in early schizophrenia to clinical and functional outcome some years later. These studies have generally been negative, perhaps because gray and white matter volumes inaccurately assess the underlying dysfunction that might be predictive of outcome. OBJECTIVE: To investigate the predictive value of frontal and temporal spectroscopy measures for outcome in patients with first-episode psychoses. DESIGN: Left prefrontal cortex and left mediotemporal lobe voxels were assessed using proton magnetic resonance spectroscopy to provide the ratio of N-acetylaspartate (NAA) and choline-containing compounds to creatine and phosphocreatine (Cr) (NAA/Cr ratio). These data were used to predict outcome at 18 months after admission, as assessed by a systematic medical record audit. SETTING: Early psychosis clinic. PARTICIPANTS: Forty-six patients with first-episode psychosis. MAIN OUTCOME MEASURES: We used regression models that included age at imaging and duration of untreated psychosis to predict outcome scores on the Global Assessment of Functioning Scale, Clinical Global Impression scales, and Social and Occupational Functional Assessment Scale, as well as the number of admissions during the treatment period. We then further considered the contributions of premorbid function and baseline level of negative symptoms. RESULTS: The only spectroscopic predictor of outcome was the NAA/Cr ratio in the prefrontal cortex. Low scores on this variable were related to poorer outcome on all measures. In addition, the frontal NAA/Cr ratio explained 17% to 30% of the variance in outcome. CONCLUSIONS: Prefrontal neuronal dysfunction is an inconsistent feature of early psychosis; rather, it is an early marker of poor prognosis across the first years of illness. The extent to which this can be used to guide treatment and whether it predicts outcome some years after first presentation are questions for further research.

Effects of hatching asynchrony on sibling negotiation, begging, jostling for position and within-brood food allocation in the barn owl Tyto alba

When siblings differ markedly in their need for food, they may benefit from signalling to each other their willingness to contest the next indivisible food item delivered by the parents. This sib-sib communication system, referred to as 'sibling negotiation', may allow them to adjust optimally to investment in begging. Using barn owl (Two alba) broods. I assessed the role of within-brood age hierarchy on sibling negotiation, and in turn on jostling for position where parents predictably deliver food (i.e. nest-box entrance), begging and within-brood food allocation. More specifically, I examined three predictions derived from a game-theoretical model of sibling negotiation where a senior and a junior sibling compete for food resources (Roulin, 2002a, Johnstone and Roulin, 2003): (1) begging effort invested by the senior sibling should be less sensitive to the junior sibling's negotiation than vice versa; (2) the junior should invest less effort in sibling negotiation than its senior sibling but a similar amount of effort in begging; and (3) within-brood food allocation should be directly related to begging but only indirectly to sibling negotiation. Two-chick broods were created and vocalization in the absence (negotiation signals directed to siblings) and presence (begging signals directed to parents) of parents was recorded. In support of the first prediction, juniors begged at a low cadence after their senior sibling negotiated intensely, probably because negotiation reflects prospective investment in begging and hence willingness to compete. In contrast, the begging of senior siblings was not sensitive to their junior sibling's negotiation. In contrast to the second prediction, juniors negotiated and begged more intensely than their senior sibling apparently because they were hungrier rather than younger. In line with the third prediction, juniors monopolized food delivered by their parents when their senior sibling begged at a low level. The begging cadence of both the junior and senior sibling, the junior's negotiation cadence, the difference in age between the two nest-mates and jostling for position were not associated with the likelihood of monopolizing food. In conclusion, sibling negotiation appears to influence begging behaviour, which, in turn, affects within-brood food allocation. Juniors may negotiate to challenge their senior siblings, and thereby determine whether seniors are less hungry before deciding to beg for food. In contrast, seniors may negotiate to deter juniors from begging.

Prediction of subcutaneous fatty acid composition from ham of CLA fed boars and gilts using a NIRS fiber optic probe

The aim of this work was the use of NIR technology by direct application of a fiber optic probe on back fat to analyze the fatty acid composition of CLA fed boars and gilts. 265 animals were fed 3 different diets and the fatty acid profile of back fat from Gluteus medius was analyzed using gas chromatography and FT-NIR. Spectra were acquired using a Bruker Optics Matrix-F duplex spectrometer equipped with a fiber optic probe (IN-268-2). Oleic and stearic fatty acids were predicted accurately; myristic, vaccenic and linoleic fatty acids were predicted with lower accuracy, while palmitic and α-linolenic fatty acids were poorly predicted. The relative percentage of fatty acids and NIR spectra showed differences in fatty acid composition of back fat from pigs fed CLA which increased the relative percentage of SFA and PUFA while MUFA decreased. Results suggest that a NIR fiber optic probe can be used to predict total saturated and unsaturated fatty acid composition, as well as the percentage of stearic and oleic. NIR showed potential as a rapid and easily implemented method to discriminate carcasses from animals fed different diets.

Ruling out coronary heart disease in primary care patients with chest pain: a clinical prediction score.

ABSTRACT: BACKGROUND: Chest pain raises concern for the possibility of coronary heart disease. Scoring methods have been developed to identify coronary heart disease in emergency settings, but not in primary care. METHODS: Data were collected from a multicenter Swiss clinical cohort study including 672 consecutive patients with chest pain, who had visited one of 59 family practitioners' offices. Using delayed diagnosis we derived a prediction rule to rule out coronary heart disease by means of a logistic regression model. Known cardiovascular risk factors, pain characteristics, and physical signs associated with coronary heart disease were explored to develop a clinical score. Patients diagnosed with angina or acute myocardial infarction within the year following their initial visit comprised the coronary heart disease group. RESULTS: The coronary heart disease score was derived from eight variables: age, gender, duration of chest pain from 1 to 60 minutes, substernal chest pain location, pain increases with exertion, absence of tenderness point at palpation, cardiovascular risks factors, and personal history of cardiovascular disease. Area under the receiver operating characteristics curve was of 0.95 with a 95% confidence interval of 0.92; 0.97. From this score, 413 patients were considered as low risk for values of percentile 5 of the coronary heart disease patients. Internal validity was confirmed by bootstrapping. External validation using data from a German cohort (Marburg, n = 774) revealed a receiver operating characteristics curve of 0.75 (95% confidence interval, 0.72; 0.81) with a sensitivity of 85.6% and a specificity of 47.2%. CONCLUSIONS: This score, based only on history and physical examination, is a complementary tool for ruling out coronary heart disease in primary care patients complaining of chest pain.

Runoff measurement and prediction for a watershed under natural vegetation in central Brazil

This work aimed to measure and analyze total rainfall (P), rainfall intensity and five-day antecedent rainfall effects on runoff (R); to compare measured and simulated R values using the Soil Conservation Service Curve Number method (CN) for each rainfall event; and to establish average R/P ratios for observed R values. A one-year (07/01/96 to 06/30/97) rainfall-runoff data study was carried out in the Capetinga watershed (962.4 ha), located at the Federal District of Brazil, 47° 52' longitude West and 15° 52' latitude South. Soils of the watershed were predominantly covered by natural vegetation. Total rainfall and runoff for the period were 1,744 and 52.5 mm, respectively, providing R/P of 3% and suggesting that watershed physical characteristics favored water infiltration into the soil. A multivariate regression analysis for 31 main rainfall-runoff events totaling 781.9 and 51.0 mm, respectively, indicated that the amount of runoff was only dependent upon rainfall volume. Simulated values of total runoff were underestimated about 15% when using CN method and an area-weighted average of the CN based on published values. On the other hand, when average values of CN were calculated for the watershed, total runoff was overestimated about 39%, suggesting that CN method shoud be used with care in areas under natural vegetation.

Early multimodal outcome prediction after cardiac arrest in patients treated with hypothermia.

OBJECTIVES: Therapeutic hypothermia and pharmacological sedation may influence outcome prediction after cardiac arrest. The use of a multimodal approach, including clinical examination, electroencephalography, somatosensory-evoked potentials, and serum neuron-specific enolase, is recommended; however, no study examined the comparative performance of these predictors or addressed their optimal combination. DESIGN: Prospective cohort study. SETTING: Adult ICU of an academic hospital. PATIENTS: One hundred thirty-four consecutive adults treated with therapeutic hypothermia after cardiac arrest. MEASUREMENTS AND MAIN RESULTS: Variables related to the cardiac arrest (cardiac rhythm, time to return of spontaneous circulation), clinical examination (brainstem reflexes and myoclonus), electroencephalography reactivity during therapeutic hypothermia, somatosensory-evoked potentials, and serum neuron-specific enolase. Models to predict clinical outcome at 3 months (assessed using the Cerebral Performance Categories: 5 = death; 3-5 = poor recovery) were evaluated using ordinal logistic regressions and receiving operator characteristic curves. Seventy-two patients (54%) had a poor outcome (of whom, 62 died), and 62 had a good outcome. Multivariable ordinal logistic regression identified absence of electroencephalography reactivity (p < 0.001), incomplete recovery of brainstem reflexes in normothermia (p = 0.013), and neuron-specific enolase higher than 33 μg/L (p = 0.029), but not somatosensory-evoked potentials, as independent predictors of poor outcome. The combination of clinical examination, electroencephalography reactivity, and neuron-specific enolase yielded the best predictive performance (receiving operator characteristic areas: 0.89 for mortality and 0.88 for poor outcome), with 100% positive predictive value. Addition of somatosensory-evoked potentials to this model did not improve prognostic accuracy. CONCLUSIONS: Combination of clinical examination, electroencephalography reactivity, and serum neuron-specific enolase offers the best outcome predictive performance for prognostication of early postanoxic coma, whereas somatosensory-evoked potentials do not add any complementary information. Although prognostication of poor outcome seems excellent, future studies are needed to further improve prediction of good prognosis, which still remains inaccurate.

TRENDS IN HOSPITAL MANAGEMENT OF ACUTE MYOCARDIAL INFARCTION IN SWITZERLAND, 1998 TO 2008 THE REGION MAKES THE DIFFERENCE: DISPARITIES IN MANAGEMENT OF ACUTE MYOCARDIAL INFARCTION WITHIN SWITZERLAND

The aim of this master's thesis was to assess the ten- year trends and regional differences in management and outcome of acute myocardial infarction (AMI) within Switzerland. The thesis is composed of two articles. First, in the article "Trends in hospital management of acute myocardial infarction in Switzerland, 1998 to 2008" over 102,700 cases of AMI with corresponding management and revascularization procedures were assessed. The results showed a considerable increase in the numbers of hospital discharges for AMI, namely due to the increase of between- hospital transfers. Rates of intensive care unit admissions remained stable. All types of revascularization procedures showed an increase. In particular, overall stenting rates increased with drug-eluting stents partly replacing bare stents. Second, in the article "The region makes the difference: disparities in management of acute myocardial infarction within Switzerland" around 25,600 cases of AMI with corresponding management were assessed for the period of 2007-2008 and according to seven Swiss regions. As reported by our results, considerable regional differences in AMI management were stated within Switzerland. Although each region showed different trends regarding revascularization interventions, Leman and Ticino contrast significantly by presenting the minimum and maximum rates in almost all assessed parameters. As a consequence these two regions differ the most from the Swiss average. The impact of the changes in trends and the regional differences in AMI management on Swiss patient's outcome and economics remains to be assessed. Purpose: To assess ten-year trends in management and outcome of acute myocardial infarction (AMI) in Switzerland. Methods: Swiss hospital discharge database for the 1998 to 2008 period. AMI was defined as a primary discharge diagnosis code I21 according to the CIM-10 classification of the World Health Organization. Management and revascularization procedures were assessed. Results: Overall, 102,729 hospital discharges with a diagnosis of AMI were analyzed. The number of hospital discharges increased almost three-fold from 5530 in 1998 to 13,834 in 2008, namely due to a considerable increase in between-hospital transfers (1352 in 1998, 6494 in 2008). Relative to all hospital discharges, Intensive Care Unit admission rate was 38.0% in 1998 and remained stable (36.2%) in 2008 (p for trend=0.25). Percutaneous revascularization rates increased from 6.0% to 39.9% (p for trend<0.001). Non-drug-eluting stent use increased from 1.3% to 16.6% (p for trend<0.05). Drug eluting stents appeared in 2004 and increased to 23.5% of hospital discharges in 2008 (p for trend=0.07). Coronary artery bypass graft increased from 1.0% to 3.0% (p for trend<0.001). Circulatory assistance increased from 0.2% to 1.7% (p for trend<0.001). Thrombolysis showed no significant changes, from 0.5% to 1.9% (p for trend=0.64). Most of these trends were confirmed after multivariate adjustment. Conclusion: Between 1998 and 2008 the number of hospital discharges for AMI increased considerably in Switzerland, namely due to between-hospital transfers. Overall stenting rates increased, drug-eluting stents partly replacing bare stents. The impact of these changes on outcome and economics remains to be assessed.

Correlates of frailty, prediction of functional decline and preventative approaches - selected results from the lausanne cohort 65+ (lc65+) study (Switzerland) (Lausanne) and the longitudinal urban cohort ageing study (Lucas) (Germany)

Overall introduction.- Longitudinal studies have been designed to investigate prospectively, from their beginning, the pathway leading from health to frailty and to disability. Knowledge about determinants of healthy ageing and health behaviour (resources) as well as risks of functional decline is required to propose appropriate preventative interventions. The functional status in older people is important considering clinical outcome in general, healthcare need and mortality. Part I.- Results and interventions from lucas (longitudinal urban cohort ageing study). Authors.- J. Anders, U. Dapp, L. Neumann, F. Pröfener, C. Minder, S. Golgert, A. Daubmann, K. Wegscheider,. W. von Renteln-Kruse Methods.- The LUCAS core project is a longitudinal cohort of urban community-dwelling people 60 years and older, recruited in 2000/2001. Further LUCAS projects are cross-sectional comparative and interventional studies (RCT). Results.- The emphasis will be on geriatric medical care in a population-based approach, discussing different forms of access, too. (Dapp et al. BMC Geriatrics 2012, 12:35; http://www.biomedcentral.com/1471-2318/12/35): - longitudinal data from the LUCAS urban cohort (n = 3.326) will be presented covering 10 years of observation, including the prediction of functional decline, need of nursing care, and mortality by using a self-filling screening tool; - interventions to prevent functional decline do focus on first (pre-clinical) signs of pre-frailty before entering the frailty-cascade ("Active Health Promotion in Old Age", "geriatric mobility centre") or disability ("home visits"). Conclusions.- The LUCAS research consortium was established to study particular aspects of functional competence, its changes with ageing, to detect pre-clinical signs of functional decline, and to address questions on how to maintain functional competence and to prevent adverse outcome in different settings. The multidimensional data base allows the exploration of several further questions. Gait performance was exmined by GAITRite®-System. Supported by the Federal Ministry for Education and Research (BMBF Funding No. 01ET1002A). Part II.- Selected results from the lausanne cohort 65+ (Lc65 + ) Study (Switzerland). Authors.- Prof Santos-Eggimann Brigitte, Dr Seematter-Bagnoud Laurence, Prof Büla Christophe, Dr Rochat Stéphane. Methods.- The Lc65+ cohort was launched in 2004 with the random selection of 3054 eligible individuals aged 65 to 70 (birth year 1934-1938) in the non-institutionalized population of Lausanne (Switzerland). Results.- Information is collected about life course social and health-related events, socio-economics, medical and psychosocial dimensions, lifestyle habits, limitations in activities of daily living, mobility impairments, and falls. Gait performance are objectively measured using body-fixed sensors. Frailty is assessed using Fried's frailty phenotype. Follow-up consists in annual self-completed questionnaires, as well as physical examination and physical and mental performance tests every three years. - Lausanne cohort 65+ (Lc65 + ): design and longitudinal outcomes. The baseline data collection was completed among 1422 participants in 2004-2005 through self-completed questionnaires, face-to-face interviews, physical examination and tests of mental and physical performances. Information about institutionalization, self-reported health services utilization, and death is also assessed. An additional random sample (n = 1525) of 65-70 years old subjects was recruited in 2009 (birth year 1939-1943). - lecture no 4: alcohol intake and gait parameters: prevalent and longitudinal association in the Lc65+ study. The association between alcohol intake and gait performance was investigated.