A comparison of the alr and ilr transformations for kernel density estimation of compositional data

In a seminal paper, Aitchison and Lauder (1985) introduced classical kernel density estimation techniques in the context of compositional data analysis. Indeed, they gave two options for the choice of the kernel to be used in the kernel estimator. One of these kernels is based on the use the alr transformation on the simplex SD jointly with the normal distribution on RD-1. However, these authors themselves recognized that this method has some deficiencies. A method for overcoming these dificulties based on recent developments for compositional data analysis and multivariate kernel estimation theory, combining the ilr transformation with the use of the normal density with a full bandwidth matrix, was recently proposed in Martín-Fernández, Chacón and Mateu- Figueras (2006). Here we present an extensive simulation study that compares both methods in practice, thus exploring the finite-sample behaviour of both estimators

Mental health of internally displaced preschool children: a cross-sectional study conducted in Bogotá, Colombia

A causa de los conflictos armados, como el de Colombia, se han desplazado por la fuerza a millones de personas, entre ellas una importante parte de la población infantil. Este estudio tuvo como objetivo evaluar la salud mental de los niños desplazados internos en edad preescolar en Bogotá Colombia, e identificar los determinantes de la salud mental en estos niños. Métodos: Estudio transversal realizado entre 279 niños que asisten a cuatro jardines infantiles en un barrio marginal de Bogotá. La salud mental infantil se evaluó con el instrumento validado de Comportamiento Infantil (CBCL) 1,5-5 años, aplicados a padres y cuidadores. Se realizo un análisis univariado y multivariado de regresión logística para evaluar la asociación entre el desplazamiento y la salud mental de los niños y para identificar las relaciones con la salud mental en los niños desplazados. Resultados: los Niños desplazados (n = 90) se identificaron con más frecuencia sobre los puntos de corte límite para las escalas CBCL que los no desplazados (n = 189) (por ejemplo, problemas totales 46,7 vs 22,8%;p \ 0,001). La asociación entre el desplazamiento y la presencia de problemas CBCL totales se mantuvo después del ajuste por factores socio-demográficos (OR Ajustado 3.3 del 95%: 1,5; 6,9). Donde la salud mental del cuidador explica en parte la asociación. En los niños desplazados, la salud mental del cuidador (p \ 0,01) y el funcionamiento familiar (p \ 0,01) se asociaron independientemente con la salud mental de los niños. La exposición a eventos traumáticos y el apoyo social también se asociaron con la salud mental del niño, sin embargo, las asociaciones no fueron independientes. Conclusión: En este barrio marginal de Bogotá, los niños en edad preescolar registrados como desplazados internos presentan peor salud mental que los no desplazados. El funcionamiento familiar y la salud mental del cuidador fueron fuerte e independientemente asociados con la salud mental de los niños y niñas desplazados.

Relación de exposición a polvos de carbón, sílice cristalina y neumoconiosis, en trabajadores de minas de socavón, Cundinamarca, 2014.

Objetivo: Determinar el nivel de riesgo de la exposición por fracción respirable a polvo de carbón y sílice cristalina y la prevalencia de neumoconiosis en trabajadores de minas de socavón del departamento de Cundinamarca. Métodos: estudio de corte transversal, en grupos de exposición similar (GES) en las minas seleccionadas, el tamaño muestral fue constituido por 11 empresas y 215 trabajadores en donde se realizó un muestreo ambiental para medir los niveles de polvo de carbón y sílice cristalina. Resultados: La edad promedio del grupo fue de 46±9,5 años y género masculino (97,2%), se encontró una asociación significativa entre polvo de carbón y neumoconiosis (p =0,050) y no fue significativa con exposición a sílice cristalina (p = 0,537). El modelo de regresión logística mostró asociación significativa con la escala de nivel de riesgo de carbón medio (OR=10.4, IC 95%:1.50, 71.41, p=0,02), ajustando con variables significativas como: tamaño de la empresa mediana (OR = 2,67, IC 95%:1.07, 6.66, p=0,04), antigüedad mayor o igual a 30 años (OR = 7,186, IC 95%:2.98, 17.29, p=0,001) y habito tabáquico por más de un año (OR = 4,437, IC 95%:2.06, 9.55, p=0,001) para sílice cristalina no hubo asociación en el modelo multivariado. Conclusión: El riesgo de exposición a carbón de nivel medio está relacionado con la prevalencia de neumoconiosis y otros factores adicionales como tamaño de la empresa mediana, antigüedad mayor o igual a 30 años y habito tabáquico por más de un año para los trabajadores de minería de socavón en Cundinamarca. Para los niveles de sílice cristalina no se encontró asociación significativa.

Evaluación del riesgo cardiovascular relacionado con la actividad laboral: Experiencia en una empresa del estado

La incidencia y prevalencia de enfermedad y riesgo cardiovascular (RCV) se incrementan con los años, como consecuencia de la falta de control en los factores de riesgo modificables, por ejemplo el sedentarismo, principalmente observado en trabajadores de oficina. El objetivo del presente trabajo fue identificar los factores asociados con el incremento del RCV en trabajadores de una empresa del estado en Bogotá, Colombia en el año 2013, a través de un estudio descriptivo de corte transversal a partir de una base de datos suministrada por la empresa con información de 272 trabajadores. Se incluyeron variables sociodemográficas, perfil ocupacional, factores de riesgo, historia clínica y medidas metabólicas. Los datos fueron estudiados a través de análisis univariado, bivariado y multivariado de regresión logística binaria. El 100% de los empleados tiene un contrato a término indefinido, siendo el género femenino más predominante. Se identificó que el RCV presente en el 11.8% de la población se asocia principalmente con la presencia de diabetes mellitus tipo 2 (ORA 9.97; IC95% 2.14-14.96, p=0.019), la alteración en el índice de masa corporal (ORA 5.67; IC95% 4.48-9.19, p=0.026), la hipertensión arterial sistólica (ORA 3.44; IC95% 2.21-4.01, p=0.037. Además hubo una relación inversa respecto al puntaje de la escala Framingham, donde menores puntajes se asociaron a menor RCV (ORA 0.04; IC95% 0.02-0.71, p=0.029), una vez se ajustó el modelo por edad, género y antigüedad en la empresa. No se encontró relación estadísticamente significativa entre el RCV, el cargo y la antigüedad laboral. Se concluye que en esta población trabajadora, independientemente de la edad, tiempo de antigüedad en la empresa y el género, los factores de riesgo clásicos para RCV están presentes y por lo tanto se deben iniciar medidas de promoción y prevención en aras de disminuir la probabilidad que el RCV encontrado se traduzca en un evento cardiovascular y de ésta manera optimizar la productividad en esta empresa.

A geostatistical analysis of the spatial variation of soil mineral nitrogen and potentially available nitrogen within an arable field.

The technology for site-specific applications of nitrogen (N) fertilizer has exposed a gap in our knowledge about the spatial variation of soil mineral N, and that which will become available during the growing season within arable fields. Spring mineral N and potentially available N were measured in an arable field together with gravimetric water content, loss on ignition, crop yield, percentages of sand, silt, and clay, and elevation to describe their spatial variation geostatistically. The areas with a larger clay content had larger values of mineral N, potentially available N, loss on ignition and gravimetric water content, and the converse was true for the areas with more sandy soil. The results suggest that the spatial relations between mineral N and loss on ignition, gravimetric water content, soil texture, elevation and crop yield, and between potentially available N and loss on ignition and silt content could be used to indicate their spatial patterns. Variable-rate nitrogen fertilizer application would be feasible in this field because of the spatial structure and the magnitude of variation of mineral N and potentially available N.

The impact of fallowing and green manuring on soil conditions and the growth of sugarcane

There are currently concerns within some sugar industries that long-term monoculture has led to soil degradation and consequent yield decline. An investigation was conducted in Swaziland to assess the effects of fallowing and green manuring practices, over a seven-month period, on sugarcane yields and the physical properties of a poorly draining clay soil. In the subsequent first sugarcane crop after planting, yields were improved from 129 t ha(-1) under continuous sugarcane to 141-144 t ha(-1) after fallowing and green manuring, but there were no significant responses in the first and second ratoon crops. Also, in the first crop after planting, root length index increased from 3.5 km m(-2) under continuous sugarcane to 5.2-6.8 km m(-2) after fallowing, and improved rooting was still evident in the first ratoon crop where there had been soil drying during the fallow period. Soil bulk density, total porosity and water-holding capacity were not affected by the fallowing practices. However, air-filled porosity increased from 11% under continuous sugarcane to 16% after fallowing, and steady state ponded infiltration rates were increased from 0.61 mm h(-1) to 1.34 mm h(-1), but these improvements were no longer evident after a year back under sugarcane. Levels of soil organic matter were reduced in all cases, probably as a result of the tillage operations involved. In the plant crop, root length was well correlated with air-filled porosity, indicating the importance of improving belowground air supply for crop production on poorly draining clay soils.

Tropical and extratropical responses of the North Atlantic atmospheric circulation to a sustained weakening of the MOC

The tropospheric response to a forced shutdown of the North Atlantic Ocean’s meridional overturning circulation (MOC) is investigated in a coupled ocean–atmosphere GCM [the third climate configuration of the Met Office Unified Model (HadCM3)]. The strength of the boreal winter North Atlantic storm track is significantly increased and penetrates much farther into western Europe. The changes in the storm track are shown to be consistent with the changes in near-surface baroclinicity, which can be linked to changes in surface temperature gradients near regions of sea ice formation and in the open ocean. Changes in the SST of the tropical Atlantic are linked to a strengthening of the subtropical jet to the north, which, combined with the enhanced storm track, leads to a pronounced split in the jet structure over Europe. EOF analysis and stationary box indices methods are used to analyze changes to the North Atlantic Oscillation (NAO). There is no consistent signal of a change in the variability of the NAO, and while the changes in the mean flow project onto the positive NAO phase, they are significantly different from it. However, there is a clear eastward shift of the NAO pattern in the shutdown run, and this potentially has implications for ocean circulation and for the interpretation of proxy paleoclimate records.

Evidence of trends in daily climate extremes over southern and west Africa

[ 1] There has been a paucity of information on trends in daily climate and climate extremes, especially from developing countries. We report the results of the analysis of daily temperature ( maximum and minimum) and precipitation data from 14 south and west African countries over the period 1961 - 2000. Data were subject to quality control and processing into indices of climate extremes for release to the global community. Temperature extremes show patterns consistent with warming over most of the regions analyzed, with a large proportion of stations showing statistically significant trends for all temperature indices. Over 1961 to 2000, the regionally averaged occurrence of extreme cold ( fifth percentile) days and nights has decreased by - 3.7 and - 6.0 days/decade, respectively. Over the same period, the occurrence of extreme hot (95th percentile) days and nights has increased by 8.2 and 8.6 days/decade, respectively. The average duration of warm ( cold) has increased ( decreased) by 2.4 (0.5) days/decade and warm spells. Overall, it appears that the hot tails of the distributions of daily maximum temperature have changed more than the cold tails; for minimum temperatures, hot tails show greater changes in the NW of the region, while cold tails have changed more in the SE and east. The diurnal temperature range (DTR) does not exhibit a consistent trend across the region, with many neighboring stations showing opposite trends. However, the DTR shows consistent increases in a zone across Namibia, Botswana, Zambia, and Mozambique, coinciding with more rapid increases in maximum temperature than minimum temperature extremes. Most precipitation indices do not exhibit consistent or statistically significant trends across the region. Regionally averaged total precipitation has decreased but is not statistically significant. At the same time, there has been a statistically significant increase in regionally averaged daily rainfall intensity and dry spell duration. While the majority of stations also show increasing trends for these two indices, only a few of these are statistically significant. There are increasing trends in regionally averaged rainfall on extreme precipitation days and in maximum annual 5-day and 1-day rainfall, but only trends for the latter are statistically significant.

Global observed changes in daily climate extremes of temperature and precipitation

A suite of climate change indices derived from daily temperature and precipitation data, with a primary focus on extreme events, were computed and analyzed. By setting an exact formula for each index and using specially designed software, analyses done in different countries have been combined seamlessly. This has enabled the presentation of the most up-to-date and comprehensive global picture of trends in extreme temperature and precipitation indices using results from a number of workshops held in data-sparse regions and high-quality station data supplied by numerous scientists world wide. Seasonal and annual indices for the period 1951-2003 were gridded. Trends in the gridded fields were computed and tested for statistical significance. Results showed widespread significant changes in temperature extremes associated with warming, especially for those indices derived from daily minimum temperature. Over 70% of the global land area sampled showed a significant decrease in the annual occurrence of cold nights and a significant increase in the annual occurrence of warm nights. Some regions experienced a more than doubling of these indices. This implies a positive shift in the distribution of daily minimum temperature throughout the globe. Daily maximum temperature indices showed similar changes but with smaller magnitudes. Precipitation changes showed a widespread and significant increase, but the changes are much less spatially coherent compared with temperature change. Probability distributions of indices derived from approximately 200 temperature and 600 precipitation stations, with near-complete data for 1901-2003 and covering a very large region of the Northern Hemisphere midlatitudes (and parts of Australia for precipitation) were analyzed for the periods 1901-1950, 1951-1978 and 1979-2003. Results indicate a significant warming throughout the 20th century. Differences in temperature indices distributions are particularly pronounced between the most recent two periods and for those indices related to minimum temperature. An analysis of those indices for which seasonal time series are available shows that these changes occur for all seasons although they are generally least pronounced for September to November. Precipitation indices show a tendency toward wetter conditions throughout the 20th century.

Toward a combined seasonal weather and crop productivity forecasting system: Determination of the working spatial scale

A methodology is presented for the development of a combined seasonal weather and crop productivity forecasting system. The first stage of the methodology is the determination of the spatial scale(s) on which the system could operate; this determination has been made for the case of groundnut production in India. Rainfall is a dominant climatic determinant of groundnut yield in India. The relationship between yield and rainfall has been explored using data from 1966 to 1995. On the all-India scale, seasonal rainfall explains 52% of the variance in yield. On the subdivisional scale, correlations vary between variance r(2) = 0.62 (significance level p < 10(-4)) and a negative correlation with r(2) = 0.1 (p = 0.13). The spatial structure of the relationship between rainfall and groundnut yield has been explored using empirical orthogonal function (EOF) analysis. A coherent, large-scale pattern emerges for both rainfall and yield. On the subdivisional scale (similar to 300 km), the first principal component (PC) of rainfall is correlated well with the first PC of yield (r(2) = 0.53, p < 10(-4)), demonstrating that the large-scale patterns picked out by the EOFs are related. The physical significance of this result is demonstrated. Use of larger averaging areas for the EOF analysis resulted in lower and (over time) less robust correlations. Because of this loss of detail when using larger spatial scales, the subdivisional scale is suggested as an upper limit on the spatial scale for the proposed forecasting system. Further, district-level EOFs of the yield data demonstrate the validity of upscaling these data to the subdivisional scale. Similar patterns have been produced using data on both of these scales, and the first PCs are very highly correlated (r(2) = 0.96). Hence, a working spatial scale has been identified, typical of that used in seasonal weather forecasting, that can form the basis of crop modeling work for the case of groundnut production in India. Last, the change in correlation between yield and seasonal rainfall during the study period has been examined using seasonal totals and monthly EOFs. A further link between yield and subseasonal variability is demonstrated via analysis of dynamical data.

Probabilistic simulations of crop yield over western India using the DEMETER seasonal hindcast ensembles

Process-based integrated modelling of weather and crop yield over large areas is becoming an important research topic. The production of the DEMETER ensemble hindcasts of weather allows this work to be carried out in a probabilistic framework. In this study, ensembles of crop yield (groundnut, Arachis hypogaea L.) were produced for 10 2.5 degrees x 2.5 degrees grid cells in western India using the DEMETER ensembles and the general large-area model (GLAM) for annual crops. Four key issues are addressed by this study. First, crop model calibration methods for use with weather ensemble data are assessed. Calibration using yield ensembles was more successful than calibration using reanalysis data (the European Centre for Medium-Range Weather Forecasts 40-yr reanalysis, ERA40). Secondly, the potential for probabilistic forecasting of crop failure is examined. The hindcasts show skill in the prediction of crop failure, with more severe failures being more predictable. Thirdly, the use of yield ensemble means to predict interannual variability in crop yield is examined and their skill assessed relative to baseline simulations using ERA40. The accuracy of multi-model yield ensemble means is equal to or greater than the accuracy using ERA40. Fourthly, the impact of two key uncertainties, sowing window and spatial scale, is briefly examined. The impact of uncertainty in the sowing window is greater with ERA40 than with the multi-model yield ensemble mean. Subgrid heterogeneity affects model accuracy: where correlations are low on the grid scale, they may be significantly positive on the subgrid scale. The implications of the results of this study for yield forecasting on seasonal time-scales are as follows. (i) There is the potential for probabilistic forecasting of crop failure (defined by a threshold yield value); forecasting of yield terciles shows less potential. (ii) Any improvement in the skill of climate models has the potential to translate into improved deterministic yield prediction. (iii) Whilst model input uncertainties are important, uncertainty in the sowing window may not require specific modelling. The implications of the results of this study for yield forecasting on multidecadal (climate change) time-scales are as follows. (i) The skill in the ensemble mean suggests that the perturbation, within uncertainty bounds, of crop and climate parameters, could potentially average out some of the errors associated with mean yield prediction. (ii) For a given technology trend, decadal fluctuations in the yield-gap parameter used by GLAM may be relatively small, implying some predictability on those time-scales.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Dietary supplements of whole linseed and vitamin E to increase levels of alpha-linolenic acid and vitamin E in bovine milk

The potential to increase the concentrations of n-3 polyunsaturated fatty acids (PUFAs) in milk fat was investigated by studying the effects of feeding a xylose-treated, whole cracked linseed supplement ( rich in alpha-linolenic acid) to dairy cows. Also the effect of increasing the dietary intake of vitamin E on the vitamin E status of milk was investigated. The effect of pasteurisation on milk fatty acid composition was also examined. Using a 3 x 2 factorial design, a total of 60 Holstein dairy cows were fed a total mixed ration based on grass silage supplemented with one of three levels of whole cracked linseed (78, 142 or 209 g . kg(-1) diet dry matter (DM); designated LL, ML or HL, respectively) in combination with one of two levels of additional dietary vitamin E intake ( 6 or 12 g vitamin E . animal(-1) . day(-1); designated LE or HE, respectively). Increasing lipid supplementation reduced (P < 0.01) diet DM intake and milk yield, and increased (P < 0.001) the overall content of oleic, vaccenic, alpha-linolenic and conjugated linoleic acids, and total PUFAs and monounsaturated fatty acids (MUFA). Myristic and palmitic acids in milk fat were reduced ( P < 0.001) through increased lipid supplementation. While α-linolenic acid concentrations were substantially increased this acid only accounted for 0.02 of total fatty acids in milk at the highest level of supplementation (630 g α-linolenic acid &BULL; animal(-1) &BULL; day(-1) for HL). Conjugated linoleic acid concentrations in milk fat were almost doubled by increasing the level of lipid supplementation (8.9, 10.4 and 16.1 g &BULL; kg(-1) fatty acids for LL, ML and HL, respectively). Although milk vitamin E contents were generally increased there was no benefit (P > 0.05) of increasing vitamin E intake from 6 to 12 g . animal(-1) . day(-1). The fatty acid composition of milk was generally not affected by pasteurisation.

Genetic parameters and evaluations from single- and multiple-trait analysis of dairy cow fertility and milk production

Genetic parameters and breeding values for dairy cow fertility were estimated from 62 443 lactation records. Two-trait analysis of fertility and milk yield was investigated as a method to estimate fertility breeding values when culling or selection based on milk yield in early lactation determines presence or absence of fertility observations in later lactations. Fertility traits were calving interval, intervals from calving to first service, calving to conception and first to last service, conception success to first service and number of services per conception. Milk production traits were 305-day milk, fat and protein yield. For fertility traits, range of estimates of heritability (h(2)) was 0.012 to 0.028 and of permanent environmental variance (c(2)) was 0.016 to 0.032. Genetic correlations (r(g)) among fertility traits were generally high ( > 0.70). Genetic correlations of fertility with milk production traits were unfavourable (range -0.11 to 0.46). Single and two-trait analyses of fertility were compared using the same data set. The estimates of h(2) and c(2) were similar for two types of analyses. However, there were differences between estimated breeding values and rankings for the same trait from single versus multi-trait analyses. The range for rank correlation was 0.69-0.83 for all animals in the pedigree and 0.89-0.96 for sires with more than 25 daughters. As single-trait method is biased due to selection on milk yield, a multi-trait evaluation of fertility with milk yield is recommended. (C) 2002 Elsevier Science B.V. All rights reserved.

Effect of replacing calcium salts of palm oil distillate with extruded linseeds on milk fatty acid composition in Jersey and Holstein cows

Clinical and biomedical studies have provided evidence for the critical role of n-3 fatty acids on the reduction of chronic disease risk in humans, including cardiovascular disease. In the current experiment, the potential to enhance milk n-3 content in two breeds with inherent genetic differences in mammary lipogenesis and de novo fatty acid synthesis was examined using extruded linseeds. Six lactating cows (three Holstein and three Jersey) were used in a two-treatment switchback design with 3 × 21-day experimental periods to evaluate the effect of iso-energetic replacement of calcium salts of palm oil distillate (CPO) in the diet (34 g/kg dry matter (DM)) with 100 g/kg DM extruded linseeds (LIN). For both breeds, replacing CPO with LIN had no effect (P > 0.05) on DM intake or milk yield, but reduced (P < 0.05) milk fat and protein yield (on average, from 760 to 706 and 573 to 552 g/day, respectively). Relative to CPO, the LIN treatment reduced (P < 0.01) total saturated fatty acid content and enhanced (P < 0.001) 18:3n-3 in milk, whereas breed by diet interactions were significant for milk fat 16:0, total trans fatty acid and conjugated linoleic acid concentrations. Increases in 18:3n-3 intake derived from LIN in the diet were transferred into milk with a mean marginal transfer efficiency of 1.8%. Proportionate changes in milk fatty acid composition were greater in the Jersey, highlighting the importance of diet–genotype interactions on mammary lipogenesis. More extensive studies are required to determine the role of genotype on milk fat composition responses to oilseeds in the diet.