Bootstrapping pairs in Distance-Based Regression

La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.

Individual prediction of automobile bodily injury claims liabilities

[spa] La mayoría de siniestros con daños corporales se liquidan mediante negociación, llegando a juicio menos del 5% de los casos. Una estrategia de negociación bien definida es, por tanto, fundamental para las compañías aseguradoras. En este artículo asumimos que la compensación monetaria concedida en juicio es la máxima cuantía que debería ser ofrecida por el asegurador en el proceso de negociación. Usando una base de datos real, implementamos un modelo log-lineal para estimar la máxima oferta de negociación. Perturbaciones no-esféricas son detectadas. Correlación ocurre cuando más de una siniestro se liquida en la misma sentencia judicial. Heterocedasticidad por grupos se debe a la influencia de la valoración del forense en la indemnización final.

The role of alliance in the relationship between therapist competence and outcome in brief psychodynamic psychotherapy.

Therapist competence is a key variable for psychotherapy research. Empirically, the relationship between competence and therapeutic outcome has shown contradictory results and needs to be clarified, especially with regard to possible variables influencing this relationship. A total of 78 outpatients were treated by 15 therapists in a very brief 4-session format, based on psychoanalytic theory. Data were analyzed by means of a nested design using hierarchical linear modeling. No direct link between therapist competence and outcome has been found, however, results corroborated the importance of alliance patterns as moderator in the relationship between therapist competence and outcome. Only in dyads with alliance change over the course of treatment was it clear that competence is positively related to outcome. These findings are discussed with regard to the importance for outcome of therapist competence and alliance construction processes.

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

MRI monitoring of focal cerebral ischemia in peroxisome proliferator-activated receptor (PPAR)-deficient mice.

Peroxisome proliferator-activated receptors (PPARs) are a potential target for neuroprotection in focal ischemic stroke. These nuclear receptors have major effects in lipid metabolism, but they are also involved in inflammatory processes. Three PPAR isotypes have been identified: alpha, beta (or delta) and gamma. The development of PPAR transgenic mice offers a promising tool for prospective therapeutic studies. This study used MRI to assess the role of PPARalpha and PPARbeta in the development of stroke. Permanent middle cerebral artery occlusion induced focal ischemia in wild-type, PPARalpha-null mice and PPARbeta-null mice. T(2)-weighted MRI was performed with a 7 T MRI scan on day 0, 1, 3, 7 and 14 to monitor lesion growth in the various genotypes. General Linear Model statistical analysis found a significant difference in lesion volume between wild-type and PPAR-null mice for both alpha and beta isotypes. These data validate high-resolution MRI for monitoring cerebral ischemic lesions, and confirm the neuroprotective role of PPARalpha and PPARbeta in the brain.

Modelling the influence of change in fire regime on the local distribution of a Mediterranean pyrophytic plant species (Cistus salviifolius) at its northern range limit

Aims: To assess the potential distribution of an obligate seeder and active pyrophyte, Cistus salviifolius, a vulnerable species in the Swiss Red List; to derive scenarios by changing the fire return interval; and to discuss the results from a conservation perspective. A more general aim is to assess the impact of fire as a natural factor influencing the vegetation of the southern slopes of the Alps. Locations: Alps, southern Switzerland. Methods: Presence-absence data to fit the model were obtained from the most recent field mapping of C. salviifolius. The quantitative environmental predictors used in this study include topographic, climatic and disturbance (fire) predictors. Models were fitted by logistic regression and evaluated by jackknife and bootstrap approaches. Changes in fire regime were simulated by increasing the time-return interval of fire (simulating longer periods without fire). Two scenarios were considered: no fire in the past 15 years; or in the past 35 years. Results: Rock cover, slope, topographic position, potential evapotranspiration and time elapsed since the last fire were selected in the final model. The Nagelkerke R-2 of the model for C. salviifolius was 0.57 and the Jackknife area under the curve evaluation was 0.89. The bootstrap evaluation revealed model robustness. By increasing the return interval of fire by either up to 15 years, or 35 years, the modelled C. salviifolius population declined by 30-40%, respectively. Main conclusions: Although fire plays a significant role, topography and rock cover appear to be the most important predictors, suggesting that the distribution of C. salviifolius in the southern Swiss Alps is closely related to the availability of supposedly competition-free sites, such as emerging bedrock, ridge locations or steep slopes. Fire is more likely to play a secondary role in allowing C. salviifolius to extend its occurrence temporarily, by increasing germination rates and reducing the competition from surrounding vegetation. To maintain a viable dormant seed bank for C. salviifolius, conservation managers should consider carrying out vegetation clearing and managing wild fire propagation to reduce competition and ensure sufficient recruitment for this species.

Bootstrapping pairs in Distance-Based Regression

La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.

Can Valproic Acid be an Inducer of Clozapine Metabolism?

Introduction: Prior clozapine studies indicated no effects, mild inhibition or induction of valproic acid (VPA) on clozapine metabolism. The hypotheses that (i) VPA is a net inducer of clozapine metabolism, and (ii) smoking modifies this inductive effect were tested in a therapeutic drug monitoring study. Methods: After excluding strong inhibitors and inducers, 353 steady-state total clozapine (clozapine plus norclozapine) concentrations provided by 151 patients were analyzed using a random intercept linear model. Results: VPA appeared to be an inducer of clozapine metabolism since total plasma clozapine concentrations in subjects taking VPA were significantly lower (27% lower; 95% confidence interval, 14-39%) after controlling for confounding variables including smoking (35% lower, 28-56%). Discussion: Prospective studies are needed to definitively establish that VPA may (i) be an inducer of clozapine metabolism when induction prevails over competitive inhibition, and (ii) be an inducer even in smokers who are under the influence of smoking inductive effects on clozapine metabolism.

Predicting future distributions of mountain plants under climate change: does dispersal capacity matter?

Many studies have investigated the impacts that climate change could potentially have on the distribution of plant species, but few have attempted to constrain projections through plant dispersal limitations. Instead, most studies published so far have been using the simplification of considering dispersal as either unlimited or null. However, depending on a species' dispersal capacity, landscape fragmentation, and the rate of climatic change, these assumptions can lead to serious over- or underestimation of a species' future distribution. To quantify the discrepancies between unlimited, realistic, and no dispersal scenarios, we carried out projections of future distribution over the 21st century for 287 mountain plant species in a study area of the Western Swiss Alps. For each species, simulations were run for four dispersal scenarios (unlimited dispersal, no dispersal, realistic dispersal and realistic dispersal with long-distance dispersal events) and under four climate change scenarios. Although simulations accounting for realistic dispersal limitations did significantly differ from those considering dispersal as unlimited or null in terms of projected future distribution, using the unlimited dispersal simplification nevertheless provided good approximations for species extinctions under more moderate climate change scenarios. Overall, simulations accounting for dispersal limitations produced, for our mountainous study area, results that were significantly closer to unlimited dispersal than to no dispersal. Finally, analyzing the temporal pattern of species extinctions over the entire 21st century showed that, due to the possibility of a large number of species shifting their distribution to higher elevation, important species extinctions for our study area might not occur before the 2080-2100 time periods.

A comparison of three methods of Mendelian randomization when the genetic instrument, the risk factor and the outcome are all binary.

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio.

The influence of drinking pattern, at individual and aggregate levels, on alcohol-related negative consequences.

AIM: To determine the extent drinking patterns (at the individual and country level) are associated with alcohol-related consequences over and above the total alcohol the person consumes. METHODS: Hierarchical linear models were estimated based on general population surveys conducted in 18 countries participating in the GENACIS project. RESULTS: In general, the positive association between drinking pattern scores and alcohol-related consequences was found at both the individual and country levels, independent of volume of drinking. In addition, a significant interaction effect indicated that the more detrimental the country's drinking pattern, the less steep the association between the volume of drinking and its consequences. CONCLUSION: Drinking patterns have an independent impact on consequences over and above the relationship between volume and consequences.

Influence of weather, rank, and home advantage on football outcomes in the gulf region.

PURPOSE: The objective of this study was to investigate the effects of weather, rank, and home advantage on international football match results and scores in the Gulf Cooperation Council (GCC) region. METHODS: Football matches (n = 2008) in six GCC countries were analyzed. To determine the weather influence on the likelihood of favorable outcome and goal difference, generalized linear model with a logit link function and multiple regression analysis were performed. RESULTS: In the GCC region, home teams tend to have greater likelihood of a favorable outcome (P < 0.001) and higher goal difference (P < 0.001). Temperature difference was identified as a significant explanatory variable when used independently (P < 0.001) or after adjustment for home advantage and team ranking (P < 0.001). The likelihood of favorable outcome for GCC teams increases by 3% for every 1-unit increase in temperature difference. After inclusion of interaction with opposition, this advantage remains significant only when playing against non-GCC opponents. While home advantage increased the odds of favorable outcome (P < 0.001) and goal difference (P < 0.001) after inclusion of interaction term, the likelihood of favorable outcome for a GCC team decreased (P < 0.001) when playing against a stronger opponent. Finally, the temperature and wet bulb globe temperature approximation were found as better indicators of the effect of environmental conditions than absolute and relative humidity or heat index on match outcomes. CONCLUSIONS: In GCC region, higher temperature increased the likelihood of a favorable outcome when playing against non-GCC teams. However, international ranking should be considered because an opponent with a higher rank reduced, but did not eliminate, the likelihood of a favorable outcome.

Scale-specific sex-biased dispersal in the Valais shrew unveiled by genetic variation on the Y chromosome, autosomes, and mitochondrial DNA.

We investigated sex specificities in the evolutionary processes shaping Y chromosome, autosomes, and mitochondrial DNA patterns of genetic structure in the Valais shrew (Sorex antinorii), a mountain dwelling species with a hierarchical distribution. Both hierarchical analyses of variance and isolation-by-distance analyses revealed patterns of population structure that were not consistent across maternal, paternal, and biparentally inherited markers. Differentiation on a Y microsatellite was lower than expected from the comparison with autosomal microsatellites and mtDNA, and it was mostly due to genetic variance among populations within valleys, whereas the opposite was observed on other markers. In addition, there was no pattern of isolation by distance for the Y, whereas there was strong isolation by distance on mtDNA and autosomes. We use a hierarchical island model of coancestry dynamics to discuss the relative roles of the microevolutionary forces that may induce such patterns. We conclude that sex-biased dispersal is the most important driver of the observed genetic structure, but with an intriguing twist: it seems that dispersal is strongly male biased at large spatial scale, whereas it is mildly biased in favor of females at local scale. These results add to recent reports of scale-specific sex-biased dispersal patterns, and emphasize the usefulness of the Y chromosome in conjunction with mtDNA and autosomes to infer sex specificities.

Constraining ACT-R models of decision strategies: An experimental paradigm

It has been repeatedly debated which strategies people rely on in inference. These debates have been difficult to resolve, partially because hypotheses about the decision processes assumed by these strategies have typically been formulated qualitatively, making it hard to test precise quantitative predictions about response times and other behavioral data. One way to increase the precision of strategies is to implement them in cognitive architectures such as ACT-R. Often, however, a given strategy can be implemented in several ways, with each implementation yielding different behavioral predictions. We present and report a study with an experimental paradigm that can help to identify the correct implementations of classic compensatory and non-compensatory strategies such as the take-the-best and tallying heuristics, and the weighted-linear model.

Advanced MRI unravels the nature of tissue alterations in early multiple sclerosis.

INTRODUCTION: In patients with multiple sclerosis (MS), conventional magnetic resonance imaging (MRI) provides only limited insights into the nature of brain damage with modest clinic-radiological correlation. In this study, we applied recent advances in MRI techniques to study brain microstructural alterations in early relapsing-remitting MS (RRMS) patients with minor deficits. Further, we investigated the potential use of advanced MRI to predict functional performances in these patients. METHODS: Brain relaxometry (T1, T2, T2*) and magnetization transfer MRI were performed at 3T in 36 RRMS patients and 18 healthy controls (HC). Multicontrast analysis was used to assess for microstructural alterations in normal-appearing (NA) tissue and lesions. A generalized linear model was computed to predict clinical performance in patients using multicontrast MRI data, conventional MRI measures as well as demographic and behavioral data as covariates. RESULTS: Quantitative T2 and T2* relaxometry were significantly increased in temporal normal-appearing white matter (NAWM) of patients compared to HC, indicating subtle microedema (P = 0.03 and 0.004). Furthermore, significant T1 and magnetization transfer ratio (MTR) variations in lesions (mean T1 z-score: 4.42 and mean MTR z-score: -4.09) suggested substantial tissue loss. Combinations of multicontrast and conventional MRI data significantly predicted cognitive fatigue (P = 0.01, Adj-R (2) = 0.4), attention (P = 0.0005, Adj-R (2) = 0.6), and disability (P = 0.03, Adj-R (2) = 0.4). CONCLUSION: Advanced MRI techniques at 3T, unraveled the nature of brain tissue damage in early MS and substantially improved clinical-radiological correlations in patients with minor deficits, as compared to conventional measures of disease.