Time-series regression models to study the short-term effects of environmental factors on health

Time series regression models are especially suitable in epidemiology for evaluating short-term effects of time-varying exposures on health. The problem is that potential for confounding in time series regression is very high. Thus, it is important that trend and seasonality are properly accounted for. Our paper reviews the statistical models commonly used in time-series regression methods, specially allowing for serial correlation, make them potentially useful for selected epidemiological purposes. In particular, we discuss the use of time-series regression for counts using a wide range Generalised Linear Models as well as Generalised Additive Models. In addition, recently critical points in using statistical software for GAM were stressed, and reanalyses of time series data on air pollution and health were performed in order to update already published. Applications are offered through an example on the relationship between asthma emergency admissions and photochemical air pollutants

Markov chain montecarlo method applied to rounding zeros of compositional data: first approach

As stated in Aitchison (1986), a proper study of relative variation in a compositional data set should be based on logratios, and dealing with logratios excludes dealing with zeros. Nevertheless, it is clear that zero observations might be present in real data sets, either because the corresponding part is completely absent –essential zeros– or because it is below detection limit –rounded zeros. Because the second kind of zeros is usually understood as “a trace too small to measure”, it seems reasonable to replace them by a suitable small value, and this has been the traditional approach. As stated, e.g. by Tauber (1999) and by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000), the principal problem in compositional data analysis is related to rounded zeros. One should be careful to use a replacement strategy that does not seriously distort the general structure of the data. In particular, the covariance structure of the involved parts –and thus the metric properties– should be preserved, as otherwise further analysis on subpopulations could be misleading. Following this point of view, a non-parametric imputation method is introduced in Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2000). This method is analyzed in depth by Martín-Fernández, Barceló-Vidal, and Pawlowsky-Glahn (2003) where it is shown that the theoretical drawbacks of the additive zero replacement method proposed in Aitchison (1986) can be overcome using a new multiplicative approach on the non-zero parts of a composition. The new approach has reasonable properties from a compositional point of view. In particular, it is “natural” in the sense that it recovers the “true” composition if replacement values are identical to the missing values, and it is coherent with the basic operations on the simplex. This coherence implies that the covariance structure of subcompositions with no zeros is preserved. As a generalization of the multiplicative replacement, in the same paper a substitution method for missing values on compositional data sets is introduced

alr approach for replacing values below the detection limit

All of the imputation techniques usually applied for replacing values below the detection limit in compositional data sets have adverse effects on the variability. In this work we propose a modification of the EM algorithm that is applied using the additive log-ratio transformation. This new strategy is applied to a compositional data set and the results are compared with the usual imputation techniques

Multivariate ARIMA Compositional Time Series Analysis

A compositional time series is obtained when a compositional data vector is observed at different points in time. Inherently, then, a compositional time series is a multivariate time series with important constraints on the variables observed at any instance in time. Although this type of data frequently occurs in situations of real practical interest, a trawl through the statistical literature reveals that research in the field is very much in its infancy and that many theoretical and empirical issues still remain to be addressed. Any appropriate statistical methodology for the analysis of compositional time series must take into account the constraints which are not allowed for by the usual statistical techniques available for analysing multivariate time series. One general approach to analyzing compositional time series consists in the application of an initial transform to break the positive and unit sum constraints, followed by the analysis of the transformed time series using multivariate ARIMA models. In this paper we discuss the use of the additive log-ratio, centred log-ratio and isometric log-ratio transforms. We also present results from an empirical study designed to explore how the selection of the initial transform affects subsequent multivariate ARIMA modelling as well as the quality of the forecasts

Combined Analysis of the Short-Term Effects of Photochemical Air Pollutants on Mortality within the EMECAM Project

In recent years, some epidemiologic studies have attributed adverse effects of air pollutants on health not only to particles and sulfur dioxide but also to photochemical air pollutants (nitrogen dioxide and ozone). The effects are usually small, leading to some inconsistencies in the results of the studies. Furthermore, the different methodologic approaches of the studies used has made it difficult to derive generic conclusions. We provide here a quantitative summary of the short-term effects of photochemical air pollutants on mortality in seven Spanish cities involved in the EMECAM project, using generalized additive models from analyses of single and multiple pollutants. Nitrogen dioxide and ozone data were provided by seven EMECAM cities (Barcelona, Gijón, Huelva, Madrid, Oviedo, Seville, and Valencia). Mortality indicators included daily total mortality from all causes excluding external causes, daily cardiovascular mortality, and daily respiratory mortality. Individual estimates, obtained from city-specific generalized additive Poisson autoregressive models, were combined by means of fixed effects models and, if significant heterogeneity among local estimates was found, also by random effects models. Significant positive associations were found between daily mortality (all causes and cardiovascular) and NO2, once the rest of air pollutants were taken into account. A 10 μg/m3 increase in the 24-hr average 1-day NO2 level was associated with an increase in the daily number of deaths of 0.43% [95% confidence interval(CI), –0.003–0.86%] for all causes excluding external. In the case of significant relationships, relative risks for cause-specific mortality were nearly twice as much as that for total mortality for all the photochemical pollutants. Ozone was independently related only to cardiovascular daily mortality. No independent statistically significant relationship between photochemical air pollutants and respiratory mortality was found. The results in this study suggest that, given the present levels of photochemical pollutants, people living in Spanish cities are exposed to health risks derived from air pollution

Comparing meta-analysis and ecological-longitudinal analysis in time-series studies: a case study of the effects of air pollution on mortality in three Spanish cities

The objective of this paper is to introduce a diVerent approach, called the ecological-longitudinal, to carrying out pooled analysis in time series ecological studies. Because it gives a larger number of data points and, hence, increases the statistical power of the analysis, this approach, unlike conventional ones, allows the complementation of aspects such as accommodation of random effect models, of lags, of interaction between pollutants and between pollutants and meteorological variables, that are hardly implemented in conventional approaches. Design—The approach is illustrated by providing quantitative estimates of the short-termeVects of air pollution on mortality in three Spanish cities, Barcelona,Valencia and Vigo, for the period 1992–1994. Because the dependent variable was a count, a Poisson generalised linear model was first specified. Several modelling issues are worth mentioning. Firstly, because the relations between mortality and explanatory variables were nonlinear, cubic splines were used for covariate control, leading to a generalised additive model, GAM. Secondly, the effects of the predictors on the response were allowed to occur with some lag. Thirdly, the residual autocorrelation, because of imperfect control, was controlled for by means of an autoregressive Poisson GAM. Finally, the longitudinal design demanded the consideration of the existence of individual heterogeneity, requiring the consideration of mixed models. Main results—The estimates of the relative risks obtained from the individual analyses varied across cities, particularly those associated with sulphur dioxide. The highest relative risks corresponded to black smoke in Valencia. These estimates were higher than those obtained from the ecological-longitudinal analysis. Relative risks estimated from this latter analysis were practically identical across cities, 1.00638 (95% confidence intervals 1.0002, 1.0011) for a black smoke increase of 10 μg/m3 and 1.00415 (95% CI 1.0001, 1.0007) for a increase of 10 μg/m3 of sulphur dioxide. Because the statistical power is higher than in the individual analysis more interactions were statistically significant,especially those among air pollutants and meteorological variables. Conclusions—Air pollutant levels were related to mortality in the three cities of the study, Barcelona, Valencia and Vigo. These results were consistent with similar studies in other cities, with other multicentric studies and coherent with both, previous individual, for each city, and multicentric studies for all three cities

Measurement Quality in Indicators of Compositions: a Compositional Multitrait-Multimethod Approach

Compositional data, also called multiplicative ipsative data, are common in survey research instruments in areas such as time use, budget expenditure and social networks. Compositional data are usually expressed as proportions of a total, whose sum can only be 1. Owing to their constrained nature, statistical analysis in general, and estimation of measurement quality with a confirmatory factor analysis model for multitrait-multimethod (MTMM) designs in particular are challenging tasks. Compositional data are highly non-normal, as they range within the 0-1 interval. One component can only increase if some other(s) decrease, which results in spurious negative correlations among components which cannot be accounted for by the MTMM model parameters. In this article we show how researchers can use the correlated uniqueness model for MTMM designs in order to evaluate measurement quality of compositional indicators. We suggest using the additive log ratio transformation of the data, discuss several approaches to deal with zero components and explain how the interpretation of MTMM designs di ers from the application to standard unconstrained data. We show an illustration of the method on data of social network composition expressed in percentages of partner, family, friends and other members in which we conclude that the faceto-face collection mode is generally superior to the telephone mode, although primacy e ects are higher in the face-to-face mode. Compositions of strong ties (such as partner) are measured with higher quality than those of weaker ties (such as other network members)

Avaluació de la toxicitat de metalls pesants i arsènic en diferents models biològics

En aquest estudi, la toxicitat de diversos metalls pesants i l'arsènic va ser analitzada utilitzant diferents models biològics. En la primera part d'aquest treball, el bioassaig de toxicitat Microtox, el qual està basat en la variació de l'emissió lumínica del bacteri luminiscent Vibrio fischeri, va ser utilitzat per establir les corbes dosi-resposta de diferents elements tòxics com el Zn(II), Pb(II), Cu(II), Hg(II), Ag(I), Co(II), Cd(II), Cr(VI), As(V) i As(III) en solucions aquoses. Els experiments es varen portar a terme a pH 6.0 i 7.0 per tal de mostrar que el pH pot influir en la toxicitat final mesurada d'alguns metalls degut als canvis relacionats amb la seva especiació química. Es varen trobar diferents tipus de corbes dosi-resposta depenent del metall analitzat i el pH del medi. En el cas de l'arsènic, l'efecte del pH en la toxicitat de l'arsenat i l'arsenit es va investigar utilitzant l'assaig Microtox en un rang de pHs comprès entre pH 5.0 i 9.0. Els valors d'EC50 determinats per l'As(V) disminueixen, reflectint un augment de la toxicitat, a mesura que el pH de la solució augmenta mentre que, en el cas de l'As(III), els valors d'EC50 quasi bé no varien entre pH 6.0 i 8.0 i només disminueixen a pH 9.0. HAsO42- i H2AsO3- es varen definir com les espècies més tòxiques. Així mateix, una anàlisi estadística va revelar un efecte antagònic entre les espècies químiques d'arsenat que es troben conjuntament a pH 6.0 i 7.0. D'altra banda, els resultats de dos mètodes estadístics per predir la toxicitat i les possibles interaccions entre el Co(II), Cd(II), Cu(II), Zn(II) i Pb(II) en mescles binàries equitòxiques es varen comparar amb la toxicitat observada sobre el bacteri Vibrio fischeri. L'efecte combinat d'aquests metalls va resultar ser antagònic per les mescles de Co(II)-Cd(II), Cd(II)-Zn(II), Cd(II)-Pb(II) i Cu(II)-Pb(II), sinèrgic per Co(II)-Cu(II) i Zn(II)-Pb(II) i additiu en els altres casos, revelant un patró complex de possibles interaccions. L'efecte sinèrgic de la combinació Co(II)-Cu(II) i la forta disminució de la toxicitat del Pb(II) quan es troba en presència de Cd(II) hauria de merèixer més atenció quan s'estableixen les normatives de seguretat ambiental. La sensibilitat de l'assaig Microtox també va ser determinada. Els valors d'EC20, els quals representen la toxicitat llindar mesurable, varen ser determinats per cada element individualment i es va veure que augmenten de la següent manera: Pb(II) < Ag(I) < Hg(II)  Cu(II) < Zn(II) < As(V) < Cd(II)  Co(II) < As(III) < Cr(VI). Aquests valors es varen comparar amb les concentracions permeses en aigues residuals industrials establertes per la normativa oficial de Catalunya (Espanya). L'assaig Microtox va resultar ser suficientment sensible per detectar els elements assajats respecte a les normes oficials referents al control de la contaminació, excepte en el cas del cadmi, mercuri, arsenat, arsenit i cromat. En la segona part d'aquest treball, com a resultats complementaris dels resultats previs obtinguts utilitzant l'assaig de toxicitat aguda Microtox, els efectes crònics del Cd(II), Cr(VI) i As(V) es varen analitzar sobre la taxa de creixement i la viabilitat en el mateix model biològic. Sorprenentment, aquests productes químics nocius varen resultar ser poc tòxics per aquest bacteri quan es mesura el seu efecte després de temps d'exposició llargs. Tot i això, en el cas del Cr(VI), l'assaig d'inhibició de la viabilitat va resultar ser més sensible que l'assaig de toxicitat aguda Microtox. Així mateix, també va ser possible observar un clar fenomen d'hormesis, especialment en el cas del Cd(II), quan s'utilitza l'assaig d'inhibició de la viabilitat. A més a més, diversos experiments es varen portar a terme per intentar explicar la manca de toxicitat de Cr(VI) mostrada pel bacteri Vibrio fischeri. La resistència mostrada per aquest bacteri podria ser atribuïda a la capacitat d'aquest bacteri de convertir el Cr(VI) a la forma menys tòxica de Cr(III). Es va trobar que aquesta capacitat de reducció depèn de la composició del medi de cultiu, de la concentració inicial de Cr(VI), del temps d'incubació i de la presència d'una font de carboni. En la tercera part d'aquest treball, la línia cel·lular humana HT29 i cultius primaris de cèl·lules sanguínies de Sparus sarba es varen utilitzar in vitro per detectar la toxicitat llindar de metalls mesurant la sobreexpressió de proteines d'estrès. Extractes de fangs precedents de diverses plantes de tractament d'aigues residuals i diferents metalls, individualment o en combinació, es varen analitzar sobre cultius cel·lulars humans per avaluar el seu efecte sobre la taxa de creixement i la capacitat d'induir la síntesi de les proteïnes Hsp72 relacionades amb l'estrès cel·lular. No es varen trobar efectes adversos significatius quan els components s'analitzen individualment. Nogensmenys, quan es troben conjuntament, es produeix un afecte advers sobre tan la taxa de creixement com en l'expressió de proteins d'estrès. D'altra banda, cèl·lules sanguínies procedents de Sparus sarba es varen exposar in vitro a diferents concentracions de cadmi, plom i crom. La proteïna d'estrès HSP70 es va sobreexpressar significativament després de l'exposició a concentracions tan febles com 0.1 M. Sota les nostres condicions de treball, no es va evidenciar una sobreexpressió de metal·lotioneïnes. Nogensmenys, les cèl·lules sanguínies de peix varen resultar ser un model biològic interessant per a ser utilitzat en anàlisis de toxicitat. Ambdós models biològics varen resultar ser molt adequats per a detectar acuradament la toxicitat produïda per metalls. En general, l'avaluació de la toxicitat basada en l'anàlisi de la sobreexpressió de proteïnes d'estrès és més sensible que l'avaluació de la toxicitat realitzada a nivell d'organisme. A partir dels resultats obtinguts, podem concloure que una bateria de bioassaigs és realment necessària per avaluar acuradament la toxicitat de metalls ja que existeixen grans variacions entre els valors de toxicitat obtinguts emprant diferents organismes i molts factors ambientals poden influir i modificar els resultats obtinguts.

Stability lobes diagram identification and surface roughness monitoring in milling processes

La millora de la productivitat i la qualitat són indubtablement dues de les principals exigències del sector productiu modern i factors clau per la competitivitat i la supervivència. Dins aquest sector,la fabricació per arrancada de material juga encara avui en dia un paper protagonista tot i l'aparició de noves tècniques de conformat per addició.Indústries com l'aeronàutica, l'automobilística,la del motlle o l'energètica, depenen en bona part de les prestacions de les màquines-eina. Aquesta Tesi aborda dos aspectes rellevants quan es tracta de millorar de la productivitat i la qualitat del sector productiu: el problema del fimbrament, més conegut per la denominació anglosaxona chatter,i la monitorització de la rugositat superficial en el mecanitzat a alta velocitat.

Desenvolupament computacional de la semblança molecular quàntica

La present Tesi Doctoral, titulada desenvolupament computacional de la semblança molecular quàntica, tracta, fonamentalment, els aspectes de càlcul de mesures de semblança basades en la comparació de funcions de densitat electrònica.El primer capítol, Semblança quàntica, és introductori. S'hi descriuen les funcions de densitat de probabilitat electrònica i llur significança en el marc de la mecànica quàntica. Se n'expliciten els aspectes essencials i les condicions matemàtiques a satisfer, cara a una millor comprensió dels models de densitat electrònica que es proposen. Hom presenta les densitats electròniques, mencionant els teoremes de Hohenberg i Kohn i esquematitzant la teoria de Bader, com magnituds fonamentals en la descripció de les molècules i en la comprensió de llurs propietats.En el capítol Models de densitats electròniques moleculars es presenten procediments computacionals originals per l'ajust de funcions densitat a models expandits en termes de gaussianes 1s centrades en els nuclis. Les restriccions físico-matemàtiques associades a les distribucions de probabilitat s'introdueixen de manera rigorosa, en el procediment anomenat Atomic Shell Approximation (ASA). Aquest procediment, implementat en el programa ASAC, parteix d'un espai funcional quasi complert, d'on se seleccionen variacionalment les funcions o capes de l'expansió, d'acord als requisits de no negativitat. La qualitat d'aquestes densitats i de les mesures de semblança derivades es verifica abastament. Aquest model ASA s'estén a representacions dinàmiques, físicament més acurades, en quant que afectades per les vibracions nuclears, cara a una exploració de l'efecte de l'esmorteïment dels pics nuclears en les mesures de semblança molecular. La comparació de les densitats dinàmiques respecte les estàtiques evidencia un reordenament en les densitats dinàmiques, d'acord al que constituiria una manifestació del Principi quàntic de Le Chatelier. El procediment ASA, explícitament consistent amb les condicions de N-representabilitat, s'aplica també a la determinació directe de densitats electròniques hidrogenoides, en un context de teoria del funcional de la densitat.El capítol Maximització global de la funció de semblança presenta algorismes originals per la determinació de la màxima sobreposició de les densitats electròniques moleculars. Les mesures de semblança molecular quàntica s'identifiquen amb el màxim solapament, de manera es mesuri la distància entre les molècules, independentment dels sistemes de referència on es defineixen les densitats electròniques. Partint de la solució global en el límit de densitats infinitament compactades en els nuclis, es proposen tres nivells de aproximació per l'exploració sistemàtica, no estocàstica, de la funció de semblança, possibilitant la identificació eficient del màxim global, així com també dels diferents màxims locals. Es proposa també una parametrització original de les integrals de recobriment a través d'ajustos a funcions lorentzianes, en quant que tècnica d'acceleració computacional. En la pràctica de les relacions estructura-activitat, aquests avenços possibiliten la implementació eficient de mesures de semblança quantitatives, i, paral·lelament, proporcionen una metodologia totalment automàtica d'alineació molecular. El capítol Semblances d'àtoms en molècules descriu un algorisme de comparació dels àtoms de Bader, o regions tridimensionals delimitades per superfícies de flux zero de la funció de densitat electrònica. El caràcter quantitatiu d'aquestes semblances possibilita la mesura rigorosa de la noció química de transferibilitat d'àtoms i grups funcionals. Les superfícies de flux zero i els algorismes d'integració usats han estat publicats recentment i constitueixen l'aproximació més acurada pel càlcul de les propietats atòmiques. Finalment, en el capítol Semblances en estructures cristal·lines hom proposa una definició original de semblança, específica per la comparació dels conceptes de suavitat o softness en la distribució de fonons associats a l'estructura cristal·lina. Aquests conceptes apareixen en estudis de superconductivitat a causa de la influència de les interaccions electró-fonó en les temperatures de transició a l'estat superconductor. En aplicar-se aquesta metodologia a l'anàlisi de sals de BEDT-TTF, s'evidencien correlacions estructurals entre sals superconductores i no superconductores, en consonància amb les hipòtesis apuntades a la literatura sobre la rellevància de determinades interaccions.Conclouen aquesta tesi un apèndix que conté el programa ASAC, implementació de l'algorisme ASA, i un capítol final amb referències bibliogràfiques.