Effect of nonnormal random components on level and power in group-randomized trials

The use of group-randomized trials is particularly widespread in the evaluation of health care, educational, and screening strategies. Group-randomized trials represent a subset of a larger class of designs often labeled nested, hierarchical, or multilevel and are characterized by the randomization of intact social units or groups, rather than individuals. The application of random effects models to group-randomized trials requires the specification of fixed and random components of the model. The underlying assumption is usually that these random components are normally distributed. This research is intended to determine if the Type I error rate and power are affected when the assumption of normality for the random component representing the group effect is violated. ^ In this study, simulated data are used to examine the Type I error rate, power, bias and mean squared error of the estimates of the fixed effect and the observed intraclass correlation coefficient (ICC) when the random component representing the group effect possess distributions with non-normal characteristics, such as heavy tails or severe skewness. The simulated data are generated with various characteristics (e.g. number of schools per condition, number of students per school, and several within school ICCs) observed in most small, school-based, group-randomized trials. The analysis is carried out using SAS PROC MIXED, Version 6.12, with random effects specified in a random statement and restricted maximum likelihood (REML) estimation specified. The results from the non-normally distributed data are compared to the results obtained from the analysis of data with similar design characteristics but normally distributed random effects. ^ The results suggest that the violation of the normality assumption for the group component by a skewed or heavy-tailed distribution does not appear to influence the estimation of the fixed effect, Type I error, and power. Negative biases were detected when estimating the sample ICC and dramatically increased in magnitude as the true ICC increased. These biases were not as pronounced when the true ICC was within the range observed in most group-randomized trials (i.e. 0.00 to 0.05). The normally distributed group effect also resulted in bias ICC estimates when the true ICC was greater than 0.05. However, this may be a result of higher correlation within the data. ^

Apoptosis of hippocampal neurons in organotypic slice culture models: direct effect of bacteria revisited

Neurons of the hippocampal dentate gyrus selectively undergo programmed cell death in patients suffering from bacterial meningitis and in experimental models of pneumococcal meningitis in infant rats. In the present study, a membrane-based organotypic slice culture system of rat hippocampus was used to test whether this selective vulnerability of neurons of the dentate gyrus could be reproduced in vitro. Apoptosis was assessed by nuclear morphology (condensed and fragmented nuclei), by immunochemistry for active caspase-3 and deltaC-APP, and by proteolytic caspase-3 activity. Co-incubation of the cultures with live pneumococci did not induce neuronal apoptosis unless cultures were kept in partially nutrient-deprived medium. Complete nutrient deprivation alone and staurosporine independently induced significant apoptosis, the latter in a dose-response way. In all experimental settings, apoptosis occurred preferentially in the dentate gyrus. Our data demonstrate that factors released by pneumococci per se failed to induce significant apoptosis in vitro. Thus, these factors appear to contribute to a multifactorial pathway, which ultimately leads to neuronal apoptosis in bacterial meningitis.

Effect sizes and standardization in neighbourhood models of forest stands: potential biases and misinterpretations

Effects of conspecific neighbours on survival and growth of trees have been found to be related to species abundance. Both positive and negative relationships may explain observed abundance patterns. Surprisingly, it is rarely tested whether such relationships could be biased or even spurious due to transforming neighbourhood variables or influences of spatial aggregation, distance decay of neighbour effects and standardization of effect sizes. To investigate potential biases, communities of 20 identical species were simulated with log-series abundances but without species-specific interactions. No relationship of conspecific neighbour effects on survival or growth with species abundance was expected. Survival and growth of individuals was simulated in random and aggregated spatial patterns using no, linear, or squared distance decay of neighbour effects. Regression coefficients of statistical neighbourhood models were unbiased and unrelated to species abundance. However, variation in the number of conspecific neighbours was positively or negatively related to species abundance depending on transformations of neighbourhood variables, spatial pattern and distance decay. Consequently, effect sizes and standardized regression coefficients, often used in model fitting across large numbers of species, were also positively or negatively related to species abundance depending on transformation of neighbourhood variables, spatial pattern and distance decay. Tests using randomized tree positions and identities provide the best benchmarks by which to critically evaluate relationships of effect sizes or standardized regression coefficients with tree species abundance. This will better guard against potential misinterpretations.

Effect of Health Risk Assessment and Counselling on Health Behaviour and Survival in Older People: A Pragmatic Randomised Trial.

BACKGROUND Potentially avoidable risk factors continue to cause unnecessary disability and premature death in older people. Health risk assessment (HRA), a method successfully used in working-age populations, is a promising method for cost-effective health promotion and preventive care in older individuals, but the long-term effects of this approach are unknown. The objective of this study was to evaluate the effects of an innovative approach to HRA and counselling in older individuals for health behaviours, preventive care, and long-term survival. METHODS AND FINDINGS This study was a pragmatic, single-centre randomised controlled clinical trial in community-dwelling individuals aged 65 y or older registered with one of 19 primary care physician (PCP) practices in a mixed rural and urban area in Switzerland. From November 2000 to January 2002, 874 participants were randomly allocated to the intervention and 1,410 to usual care. The intervention consisted of HRA based on self-administered questionnaires and individualised computer-generated feedback reports, combined with nurse and PCP counselling over a 2-y period. Primary outcomes were health behaviours and preventive care use at 2 y and all-cause mortality at 8 y. At baseline, participants in the intervention group had a mean ± standard deviation of 6.9 ± 3.7 risk factors (including unfavourable health behaviours, health and functional impairments, and social risk factors) and 4.3 ± 1.8 deficits in recommended preventive care. At 2 y, favourable health behaviours and use of preventive care were more frequent in the intervention than in the control group (based on z-statistics from generalised estimating equation models). For example, 70% compared to 62% were physically active (odds ratio 1.43, 95% CI 1.16-1.77, p = 0.001), and 66% compared to 59% had influenza vaccinations in the past year (odds ratio 1.35, 95% CI 1.09-1.66, p = 0.005). At 8 y, based on an intention-to-treat analysis, the estimated proportion alive was 77.9% in the intervention and 72.8% in the control group, for an absolute mortality difference of 4.9% (95% CI 1.3%-8.5%, p = 0.009; based on z-test for risk difference). The hazard ratio of death comparing intervention with control was 0.79 (95% CI 0.66-0.94, p = 0.009; based on Wald test from Cox regression model), and the number needed to receive the intervention to prevent one death was 21 (95% CI 12-79). The main limitations of the study include the single-site study design, the use of a brief self-administered questionnaire for 2-y outcome data collection, the unavailability of other long-term outcome data (e.g., functional status, nursing home admissions), and the availability of long-term follow-up data on mortality for analysis only in 2014. CONCLUSIONS This is the first trial to our knowledge demonstrating that a collaborative care model of HRA in community-dwelling older people not only results in better health behaviours and increased use of recommended preventive care interventions, but also improves survival. The intervention tested in our study may serve as a model of how to implement a relatively low-cost but effective programme of disease prevention and health promotion in older individuals. TRIAL REGISTRATION International Standard Randomized Controlled Trial Number: ISRCTN 28458424.

Effect of routine claw trimming on claw temperature in dairy cows measured by infrared thermography

Infrared thermography (IRT) was used to assess the effect of routine claw trimming on claw temperature. In total, 648 IRT observations each were collected from 81 cows housed in 6 tiestalls before and 3 wk after claw trimming. The feet were classified as either healthy (nonlesion group, n = 182) or affected with infectious foot disorders (group IFD, n = 142). The maximal surface temperatures of the coronary band and skin and the difference of the maximal temperatures (ΔT) between the lateral and medial claws of the respective foot were assessed. Linear mixed models, correcting for the hierarchical structure of the data, ambient temperature, and infectious status of the claws, were developed to evaluate the effect of time in relation to the trimming event (d 0 versus d 21) and claw (medial versus lateral). Front feet and hind feet were analyzed separately. Ambient temperature and infectious foot status were identified as external and internal factors, respectively, that significantly affected claw temperature. Before claw trimming, the lateral claws of the hind feet were significantly warmer compared with the medial claws, whereas such a difference was not evident for the claws of the front feet. At d 21, ΔT of the hind feet was reduced by ≥ 0.25 °C, whereas it was increased by ≤ 0.13 °C in the front feet compared with d 0. Therefore, trimming was associated with a remarkable decrease of ΔT of the hind claws. Equalizing the weight bearing of the hind feet by routine claw trimming is associated with a measurable reduction of ΔT between the paired hind claws.

Effect of an intervertebral disk spacer on stiffness after monocortical screw/polymethylmethacrylate fixation in simulated and cadaveric canine cervical vertebral columns.

OBJECTIVE To determine the biomechanical effect of an intervertebral spacer on construct stiffness in a PVC model and cadaveric canine cervical vertebral columns stabilized with monocortical screws/polymethylmethacrylate (PMMA). STUDY DESIGN Biomechanical study. SAMPLE POPULATION PVC pipe; cadaveric canine vertebral columns. METHODS PVC model-PVC pipe was used to create a gap model mimicking vertebral endplate orientation and disk space width of large-breed canine cervical vertebrae; 6 models had a 4-mm gap with no spacer (PVC group 1); 6 had a PVC pipe ring spacer filling the gap (PCV group 2). Animals-large breed cadaveric canine cervical vertebral columns (C2-C7) from skeletally mature dogs without (cadaveric group 1, n = 6, historical data) and with an intervertebral disk spacer (cadaveric group 2, n = 6) were used. All PVC models and cadaver specimens were instrumented with monocortical titanium screws/PMMA. Stiffness of the 2 PVC groups was compared in extension, flexion, and lateral bending using non-destructive 4-point bend testing. Stiffness testing in all 3 directions was performed of the unaltered C4-C5 vertebral motion unit in cadaveric spines and repeated after placement of an intervertebral cortical allograft ring and instrumentation. Data were compared using a linear mixed model approach that also incorporated data from previously tested spines with the same screw/PMMA construct but without disk spacer (cadaveric group 1). RESULTS Addition of a spacer increased construct stiffness in both the PVC model (P < .001) and cadaveric vertebral columns (P < .001) compared to fixation without a spacer. CONCLUSIONS Addition of an intervertebral spacer significantly increased construct stiffness of monocortical screw/PMMA fixation.

Protective effect of a germline, IL-17-neutralizing antibody in murine models of autoimmune inflammatory disease.

Monoclonal antibodies (mAbs) inhibiting cytokines have recently emerged as new drug modalities for the treatment of chronic inflammatory diseases. Interleukin-17 (IL-17) is a T-cell-derived central mediator of autoimmunity. Immunization with Qβ-IL-17, a virus-like particle based vaccine, has been shown to produce autoantibodies in mice and was effective in ameliorating disease symptoms in animal models of autoimmunity. To characterize autoantibodies induced by vaccination at the molecular level, we generated mouse mAbs specific for IL-17 and compared them to germline Ig sequences. The variable regions of a selected hypermutated high-affinity anti-IL-17 antibody differed in only three amino acid residues compared to the likely germline progenitor. An antibody, which was backmutated to germline, maintained a surprisingly high affinity (0.5 nM). The ability of the parental hypermutated antibody and the derived germline antibody to block inflammation was subsequently tested in murine models of multiple sclerosis (experimental autoimmune encephalomyelitis), arthritis (collagen-induced arthritis), and psoriasis (imiquimod-induced skin inflammation). Both antibodies were able to delay disease onset and significantly reduced disease severity. Thus, the mouse genome unexpectedly encodes for antibodies with the ability to functionally neutralize IL-17 in vivo.

Evaluation of the Hosmer -Lemeshow global goodness -of-fit statistic for mixed variable logistic regression models

The performance of the Hosmer-Lemeshow global goodness-of-fit statistic for logistic regression models was explored in a wide variety of conditions not previously fully investigated. Computer simulations, each consisting of 500 regression models, were run to assess the statistic in 23 different situations. The items which varied among the situations included the number of observations used in each regression, the number of covariates, the degree of dependence among the covariates, the combinations of continuous and discrete variables, and the generation of the values of the dependent variable for model fit or lack of fit.^ The study found that the $\rm\ C$g* statistic was adequate in tests of significance for most situations. However, when testing data which deviate from a logistic model, the statistic has low power to detect such deviation. Although grouping of the estimated probabilities into quantiles from 8 to 30 was studied, the deciles of risk approach was generally sufficient. Subdividing the estimated probabilities into more than 10 quantiles when there are many covariates in the model is not necessary, despite theoretical reasons which suggest otherwise. Because it does not follow a X$\sp2$ distribution, the statistic is not recommended for use in models containing only categorical variables with a limited number of covariate patterns.^ The statistic performed adequately when there were at least 10 observations per quantile. Large numbers of observations per quantile did not lead to incorrect conclusions that the model did not fit the data when it actually did. However, the statistic failed to detect lack of fit when it existed and should be supplemented with further tests for the influence of individual observations. Careful examination of the parameter estimates is also essential since the statistic did not perform as desired when there was moderate to severe collinearity among covariates.^ Two methods studied for handling tied values of the estimated probabilities made only a slight difference in conclusions about model fit. Neither method split observations with identical probabilities into different quantiles. Approaches which create equal size groups by separating ties should be avoided. ^

Promoting mixed stands through conversion treatments. Effect of holm oak coppices thinning on black pine regeneration

The effects of conversion treatments, depending on ecological factors and silvicultural parameters (thinning intensity, thinning type and rotation, among others) have been studied during the last fifteen years in an experimental trial in Central Spain. The general climate is continental Mediterranean; soils are low depth and limy; vegetation is an homogeneous dense coppices of Quercus ilex with isolated Pinus nigra trees. The experimental design (three locations) includes different thinning intensities (from 0 to 100% of extracted basal area). Inventories have been carried out in 1994 and 2010; thinning treatments were done in 1995 and 2011. Analysis of the effects of the conversion treatment show the increment of diameter and height growth rates, the canopy recovery and the stand resprouting, finding differences in these effects between thinning treatments. Besides the induced changes at holm oak stand, the application of conversion treatment clearly changed the woodland dynamics. Fifteen years after the thinnings, floristic composition varied and an abundant pine regeneration was installed in the woodland. In this work we describe the changes between inventories in tree species composition and diameter distribution, specially in the case of black pine. The conversion treatment caused changes in forest dynamics in the short term, increasing biodiversity and diversifying the forest structure. The fast installation of Pinus regeneration suggests the potential of the zone for the establishment of multipurpose mixed Quercus-Pinus stands in wide areas where Quercus species were favoured by human populations for firewood production. Conversion treatment of coppices, with the creation of mixed stands, constitutes a good management alternative for extensive areas and an interesting technique to adaptation to global change.

Small-scale variation of vegetation in a mixed forest understorey is partly controlled by the effect of overstory composition on litter accumulation

We investigated how richness and composition of vascular plant species in the understory of a mixed hardwood forest stand varied with respect to the abundance and composition of the overstory. The stand is in central Spain and represents the southernmost range of distribution of several tree and herbaceous species in Europe. Understory species were identified in 46 quadrats (0.25 m2) where variables litter depth and light availability were measured. In addition, we estimated tree density, basal area, and percent basal area by tree species within 6-m-radius areas around each plot. Species richness and composition were studied using path analysis and scale-dependent geostatistical methods, respectively. We found that the relative abundance of certain trees species in the overstory was more important than total overstory abundance in explaining understory species richness. Richness decreased as soil litter depth increased, and soil litter increased as the relative proportion of Fagus sylvatica in the overstory increased, which accounted for a negative, indirect effect of Fagus sylvatica on richness. Regarding understory species composition, we found that some species distributed preferentially below certain tree species. For example, Melica uniflora was most frequent below Fagus sylvatica and Quercus petraea while the increasing proportion of Q. pyrenaica in the overstory favored the presence of Cruciata glabra, Arenaria montana, Prunus avium, Conopodium bourgaei, Holcus mollis, Stellaria media and Galium aparine in the understory. Overall, these results emphasize the importance of individual tree species in controlling the assemblage and richness of understory species in mixed stands. We conclude that soil litter accumulation is one way through which overstory composition shapes the understory community.

An empirical mixed model to quantify climate influence on the growth of Pinus halepensis Mill. stands in South-Eastern Spain

The influence of climate on forest stand composition, development and growth is undeniable. Many studies have tried to quantify the effect of climatic variables on forest growth and yield. These works become especially important because there is a need to predict the effects of climate change on the development of forest ecosystems. One of the ways of facing this problem is the inclusion of climatic variables into the classic empirical growth models. The work has a double objective: (i) to identify the indicators which best describe the effect of climate on Pinus halepensis growth and (ii) to quantify such effect in several scenarios of rainfall decrease which are likely to occur in the Mediterranean area. A growth mixed model for P. halepensis including climatic variables is presented in this work. Growth estimates are based on data from the Spanish National Forest Inventory (SNFI). The best results are obtained for the indices including rainfall, or rainfall and temperature together, with annual precipitation, precipitation effectiveness, Emberger?s index or free bioclimatic intensity standing out among them. The final model includes Emberger?s index, free bioclimatic intensity and interactions between competition and climate indices. The results obtained show that a rainfall decrease about 5% leads to a decrease in volume growth of 5.5?7.5% depending on site quality.

Mixing effect on volume growth of Fagus sylvatica and Pinus sylvestris is modulated by stand density

Despite the increasing relevance of mixed stands due to their potential benefits; little information is available with regard to the effect of mixtures on yield in forest systems. Hence, it is necessary to study inter-specific relationships, and the resulting yield in mixed stands, which may vary with stand development, site or stand density, etc. In Spain, the province of Navarra is considered one of the biodiversity reservoirs; however, mixed forests occupy only a small area, probably as a consequence of management plans, in which there is an excessive focus on the productivity aspect, favoring the presence of pure stands of the most marketable species. The aim of this paper is to study how growth efficiencies of beech (Fagus sylvatica) and pine (Pinus sylvestris) are modified by the admixture of the other species and to determine whether stand density modifies interspecific relationships and to what extent. Two models were fitted from Spanish National Forest Inventory data, for P. sylvestris and F. sylvatica respectively, which relate the growth efficiency of the species, i.e. the volume increment of the species divided by the species proportion by area, with dominant height, quadratic mean diameter, stocking degree, and the species proportions by area of each species. Growth efficiency of pine increased with the admixture of beech, decreasing this positive effect when stocking degree increased. However, the positive effect of pine admixture on beech growth was greater at higher stocking degrees. Growth efficiency of beech was also dependent on stand dominant height, resulting in a net negative mixing effect when stand dominant heights and stocking degrees were simultaneously low. There is a relatively large range of species proportions and stocking degrees which results in transgressive overyielding: higher volume increments in mixed stands than that of the most productive pure pine stands. We concluded that stocking degree is a key factor in between-species interactions, being the effects of mixing not always greater at higher stand densities, but it depends on species composition.

Mixed models for short-run forecasting of electricity prices: Application for the Spanish market

Short-run forecasting of electricity prices has become necessary for power generation unit schedule, since it is the basis of every profit maximization strategy. In this article a new and very easy method to compute accurate forecasts for electricity prices using mixed models is proposed. The main idea is to develop an efficient tool for one-step-ahead forecasting in the future, combining several prediction methods for which forecasting performance has been checked and compared for a span of several years. Also as a novelty, the 24 hourly time series has been modelled separately, instead of the complete time series of the prices. This allows one to take advantage of the homogeneity of these 24 time series. The purpose of this paper is to select the model that leads to smaller prediction errors and to obtain the appropriate length of time to use for forecasting. These results have been obtained by means of a computational experiment. A mixed model which combines the advantages of the two new models discussed is proposed. Some numerical results for the Spanish market are shown, but this new methodology can be applied to other electricity markets as well

Analyzing size-symmetric vs. size-asymmetric and intra-vs. inter-specific competition in beech (Fagus sylvatica L.) mixed stands

In mixed stands, inter-specific competition can be lower than intra-specific competition when niche complementarity and/or facilitation between species prevail. These positive interactions can take place at belowground and/or aboveground levels. Belowground competition tends to be size symmetric while the aboveground competition is usually for light and almost always size-asymmetric. Interactions between forest tree species can be explored analyzing growth at tree level by comparing intra and inter-specific competition. At the same time, possible causes of niche complementarity can be inferred relating intra and inter-specific competition with the mode of competition, i.e. size-symmetric or sizeasymmetric. The aim of this paper is to further our understanding of the interactions between species and to detect possible causes of competition reduction in mixed stands of beech (Fagus sylvatica L.) with other species: pine?beech, oak?beech and fir?beech. To test whether species growth is better explained by size-symmetric and/or size-asymmetric competition, five different competition structures where included in basal area growth models fitted using data from the Spanish National Forest Inventory for the Pyrenees. These models considered either size-symmetry only (Reineke?s stand density index, SDI), size-asymmetry only (SDI of large trees or SDI of small trees), or both combined. In order to assess the influence of the admixture, these indices were introduced in two different ways, one of which was to consider that trees of all species compete in a similar way, and the other was to split the stand density indices into intra- and inter-specific competition components. The results showed that in pine?beech mixtures, there is a slightly negative effect of beech on pine basal area growth while beech benefitted from the admixture of Scots pine; this positive effect being greater as the proportion of pine trees in larger size classes increases. In oak?beech mixtures, beech growth was also positively influenced by the presence of oaks that were larger than the beech trees. The growth of oak, however, decreased when the proportion of beech in SDI increased, although the presence of beech in larger size classes promoted oak growth. Finally, in fir?beech mixtures, neither fir nor beech basal area growth were influenced by the presence of the other species. The results indicate that size-asymmetric is stronger than size-symmetric competition in these mixtures, highlighting the importance of light in competition. Positive species interactions in size-asymmetric competition involved a reduction of asymmetry in tree size-growth relationships.

Natural regeneration in mediterranean mixed forests: effect of environmental conditions on the physiological optimum niche of the species = Regeneración natural en masas mixtas mediterráneas: efecto de las condiciones ambientales en el nicho óptimo fisiológico de las especies

Los montes Mediterráneos han experimentado múltiples cambios en las últimas décadas (tanto en clima como en usos), lo que ha conducido a variaciones en la distribución de especies. El aumento previsto de las temperaturas medias junto con la mayor variabilidad intra e inter anual en cuanto a la ocurrencia de eventos extremos o disturbios naturales (como periodos prolongados de sequía, olas de frío o calor, incendios forestales o vendavales) pueden dañar significativamente al regenerado, llevándolo hasta la muerte, y jugando un papel decisivo en la composición de especies y en la dinámica del monte. La amplitud ecológica de muchas especies forestales puede verse afectada, de forma que se esperan cambios en sus nichos actuales de regeneración. Sin embargo, la migración latitudinal de las especies en busca de mejores condiciones, podría ser una explicación demasiado simplista de un proceso mucho más complejo de interacción entre la temperatura y la precipitación, que afectaría a cada especie de un modo distinto. En este sentido tanto la capacidad de adaptación al estrés ambiental de una determinada especie, así como su habilidad para competir por los recursos limitados, podría significar variaciones dentro de una comunidad. Las características fisiológicas y morfológicas propias de cada especie se encuentran fuertemente relacionadas con el lugar donde cada una puede surgir, qué especies pueden convivir y como éstas responden a las condiciones ambientales. En este sentido, el conocimiento sobre las distintas respuestas ecofisiológicas observadas ante cambios ambientales puede ser fundamentales para la predicción de variaciones en la distribución de especies, composición de la comunidad y productividad del monte ante el cambio global. En esta tesis investigamos el grado de tolerancia y sensibilidad que cada una de las tres especies de estudio, coexistentes en el interior peninsular ibérico (Pinus pinea, Quercus ilex y Juniperus oxycedrus), muestra ante los factores abióticos de estrés típicos de la región Mediterránea. Nuestro trabajo se ha basado en la definición del nicho óptimo fisiológico para el regenerado de cada especie a través de la investigación en profundidad del efecto de la sequía, la temperatura y el ambiente lumínico. Para ello, hemos desarrollado un modelo de predicción de la tasa de asimilación de carbono que nos ha permitido identificar las condiciones óptimas ambientales donde el regenerado de cada especie podría establecerse con mayor facilidad. En apoyo a este trabajo y con la idea de estudiar el efecto de la sequía a nivel de toda la planta hemos desarrollado un experimento paralelo en invernadero. Aquí se han aplicado dos regímenes hídricos para estudiar las características fisiológicas y morfológicas de cada especie, sobre todo a nivel de raíz y crecimiento del tallo, y relacionarlas con las diferentes estrategias en el uso del agua de las especies. Por último, hemos estudiado los patrones de aclimatación y desaclimatación al frio de cada especie, identificando los periodos de sensibilidad a heladas, así como cuellos de botella donde la competencia entre especies podría surgir. A pesar de que el pino piñonero ha sido la especie objeto de la gestión de estas masas durante siglos, actualmente se encuentra en la posición más desfavorable para combatir el cambio global, presentado el nicho fisiológico más estrecho de las tres especies. La encina sin embargo, ha resultado ser la especie mejor cualificada para afrontar este cambio, seguida muy de cerca por el enebro. Nuestros resultados sugieren una posible expansión en el rango de distribución de la encina, un aumento en la presencia del enebro y una disminución progresiva del pino piñonero a medio plazo en estas masas. ABSTRACT Mediterranean forests have undergone multiple changes over the last decades (in both climate and land use), which have lead to variations in the distribution of species. The expected increase in mean annual temperature together with the greater inter and intra-annual variability in extreme events and disturbances occurrence (such as prolonged drought periods, cold or heat waves, wildfires or strong winds) can significantly damage natural regeneration, up to causing death, playing a decisive role on species composition and forest dynamics. The ecological amplitude for adaptation of many species can be affected in such a way that changes in the current regeneration niches of many species are expected. However, the forecasted poleward migration of species seeking better conditions could be an oversimplification of what is a more complex phenomenon of interactions among temperature and precipitation, that would affect different species in different ways. In this regard, either the ability to adapt to environmental stresses or to compete for limited resources of a single species in a mixed forest could lead to variations within a community. The ecophysiological and morphological traits specific to each species are strongly related to the place where each species can emerge, which species can coexist, and how they respond to environmental conditions. In this regard, the understanding of the ecophysiological responses observed against changes in environmental conditions can be essential for predicting variations in species distribution, community composition, and forest productivity in the context of global change. In this thesis we investigated the degree of tolerance and sensitivity that each of the three studied species, co-occurring in central of the Iberian Peninsula (Pinus pinea, Quercus ilex and Juniperus oxycedrus), show against the typical abiotic stress factors in the Mediterranean region. Our work is based on the optimal physiological niche for regeneration of each species through in-depth research on the effect of drought, temperature and light environment. For this purpose, we developed a model to predict the carbon assimilation rate which allows us to identify the optimal environmental conditions where regeneration from each species could establish itself more easily. To obtain a better understanding about the effect of low temperature on regeneration, we studied the acclimation and deacclimation patterns to cold of each species, identifying period of frost sensitivity, as well as bottlenecks where competition between species can arise. Finally, to support our results about the effect of water availabilty, we conducted a greenhouse experiment with a view of studying the drought effect at the whole plant level. Here, two watering regimes were applied in order to study the physiological and morphological traits of each species, mainly at the level of the root system and stem growth, and so relate them to the different water use strategies of the species. Despite the fact that stone pine has been the target species for centuries, nowadays this species is in the most unfavorable position to cope with climate change. Holm oak, however, resulted the species that is best adapted to tolerate the predicted changes, followed closely by prickly juniper. Our results suggest a feasible expansion of the distribution range in holm oak, an increase in the prickly juniper presence and a progressive decreasing of stone pine presence in the medium term in these stone pine-holm oak-prickly juniper mixed forests.