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

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

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

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequent in practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of the sample, thus producing a transfer of mass from some components to other ones, but preserving the total mass present in the system. This evolution is traditionally modelled as a system of ordinary di erential equations of the mass of each component. However, this kind of evolution can be decomposed into a compositional change, expressed in terms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despite of some subcompositions behaving linearly. The goal is to study the characteristics of such simplicial systems of di erential equa- tions such as linearity and stability. This is performed extracting the compositional dif ferential equations from the mass equations. Then, simplicial derivatives are expressed in coordinates of the simplex, thus reducing the problem to the standard theory of systems of di erential equations, including stability. The characterisation of stability of these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and the associated behaviour of the orbits are the main tools. For a three component system, these orbits can be plotted both in coordinates of the simplex or in a ternary diagram. A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is a radioactive decay

Using problem-based learning for introducing producer theory and market structure in intermediate microeconomics

This paper shows how instructors can use the problem‐based learning method to introduce producer theory and market structure in intermediate microeconomics courses. The paper proposes a framework where different decision problems are presented to students, who are asked to imagine that they are the managers of a firm who need to solve a problem in a particular business setting. In this setting, the instructors’ role is to provide both guidance to facilitate student learning and content knowledge on a just‐in‐time basis

Revaloració de la fracció cel·lular de la sang de porc procedent d'escorxadors industrials

La sang de porc és un subproducte comestible que es genera als escorxadors industrials durant el procés d'obtenció de la canal. Aquest subproducte es caracteritza per presentar una elevada càrrega contaminant i, degut a l'elevat volum que es genera, és necessari trobar estratègies que permetin la seva revaloració i aprofitament, a la vegada que disminuïm la contaminació ambiental i les despeses que es deriven del seu processament abans de l'abocament. La fracció cel·lular (FC) constitueix el 40 % de la sang de porc i conté principalment l'hemoglobina (Hb), que representa al voltant del 90 % del contingut en proteïna d'aquesta fracció (un 35 % aproximadament). L'elevat percentatge en proteïna i en ferro, i les seves bones propietats funcionals fan que l'aprofitament d'aquest subproducte com a primera matèria o ingredient de la indústria alimentària sigui una alternativa molt útil a l'hora de reduir les despeses de la indústria càrnia, sempre que es resolguin els problemes de l'enfosquiment i dels sabors estranys que pot conferir la FC quan s'addiciona a productes alimentaris. Una altra possible utilització de la FC és aprofitar les propietats colorants de l'Hb o del grup hemo, com a colorant d'origen natural en diversos productes alimentaris. Els objectius del present treball eren, en primer lloc, determinar les millors condicions d'aplicació del procés de conservació de la FC mitjançant la deshidratació per atomització i caracteritzar físico-químicament i microbiològica el concentrat d'Hb en pols. En segon lloc, avaluar l'eficàcia de diferents additius antioxidants i/o segrestants del ferro per prevenir l'enfosquiment que pateix la FC durant la deshidratació. En tercer lloc, aplicar tractaments d'altes pressions hidrostàtiques com a procés d'higienització i avaluar els efectes d'aquest tractament sobre la microbiota contaminant, el color i les propietats funcionals de la FC. Finalment, desenvolupar un procés d'obtenció d'hidrolitzats proteics descolorats a partir de l'Hb amb la finalitat d'utilitzar-los com a ingredients nutricionals i/o funcionals. La millor temperatura de deshidratació per atomització de la FC hemolitzada era 140ºC. La FC en pols presentava un contingut en humitat del 5,3 % i un percentatge de solubilitat proteica del 96 %. La deshidratació per atomització induïa canvis en l'estructura nativa de l'Hb i, per tant, un cert grau de desnaturalització que pot conduir a una disminució de les seves propietats funcionals. L'extracte sec de la FC en pols estava composat per un 94,6 % de proteïna, un 3 % de sals minerals i un 0,7 % de greix. Els valors CIE L*a*b* del color de la FC en pols eren força constants i reflectien el color vermell marró fosc d'aquesta, a causa de l'oxidació del ferro hèmic que es produeix durant la deshidratació. La càrrega contaminant de la FC fresca de la sang de porc era força elevada i el tractament d'hemòlisi amb ultrasons i la centrifugació posterior no produïen una reducció significativa de la microbiota contaminant, obtenint un producte amb uns recomptes microbiològics de l'ordre de 106 ufc·mL-1. La deshidratació per atomització produïa una disminució d'una unitat logarítmica dels recomptes totals de la FC hemolitzada. Tanmateix, el producte en pols encara reflectia l'elevada contaminació de la primera matèria, fet que condiciona negativament la seva utilització com a ingredient alimentari, a no ser que es millorin les condicions de recollida de la sang a l'escorxador o que aquesta o la FC es sotmeti a algun tractament d'higienització prèviament a la deshidratació. Les isotermes de sorció a 20ºC de la FC en pols tenien forma sigmoïdal i una histèresi estreta i llarga. L'equació GAB és un bon model matemàtic per ajustar les dades de sorció obtingudes experimentalment i determinar la isoterma d'adsorció de la FC deshidratada per atomització. El percentatge d'humitat de la FC deshidratada a 140ºC es corresponia a un valor d'aw a 20 ºC d'aproximadament el 0,16. Tenint en compte que estava per sota dels valors d'aw corresponents a la capa monomolecular, es pot garantir la conservació a temperatura ambient del producte, sempre que s'envasi en recipients tancats que no permetin l'entrada d'humitat de l'exterior. De l'estudi de la possible estabilització del color de la FC deshidratada per atomització mitjançant l'addició d'antioxidants i/o segrestants de ferro, es va observar que només l'àcid ascòrbic, la glucosa, l'àcid nicotínic i la nicotinamida, tenien efectes positius sobre el color del producte en pols. L'ascòrbic i la glucosa no milloraven la conservació del color de l'Hb però disminuïen l'enfosquiment que es produeix durant la deshidratació, amb la qual cosa es pot obtenir un producte en pols de color marró més clar. L'addició de dextrina o L-cisteïna no disminuïa l'enfosquiment ni evitava el canvi de color de l'Hb. L'àcid nicotínic i la nicotinamida protegien el color de l'Hb durant el procés de deshidratació i l'emmagatzematge de la FC en pols. Les millors condicions d'aplicació del tractament amb altes pressions hidrostàtiques (HHP) sobre la FC eren 400 MPa, a 20ºC, durant 15 minuts, perquè produïen una millora significativa de la qualitat microbiològica, no afectaven negativament al color, no comprometien gaire la solubilitat proteica l'Hb i, malgrat que produïen un augment de la viscositat, la FC romania fluida després del tractament. Aquest tractament permetia una reducció de la microbiota contaminant de la FC d'entre 2 i 3 unitats logarítmiques. L'aplicació de l'alta pressió i la posterior deshidratació per atomització permetien obtenir un producte en pols amb recomptes totals de l'ordre de 2,8 unitats logarítmiques. El color de la FC pressuritzada en pols era igual que el de la FC control deshidratada, perquè ambdues mostres presentaven la mateixa susceptibilitat a l'oxidació del grup hemo produïda per la deshidratació. L'alta pressió incrementava la susceptibilitat de l'Hb als efectes desnaturalitzants de la deshidratació, fonamentalment a pH 7 (PIE), ja que es va observar una disminució de la solubilitat proteica a pH neutre després dels 2 processos tecnològics. La FC en pols presentava una màxima capacitat escumant al PIE de l'Hb. L'aplicació del tractament HHP produïa una disminució de la capacitat escumant de la FC en pols, però no tenia efectes negatius sobre l'estabilitat de l'escuma formada. Tampoc es van observar efectes negatius del tractament HHP sobre l'activitat emulsionant de l'Hb. La màxima activitat emulsionant de l'Hb s'aconseguia amb una concentració de FC en pols de l'1,5 % a pH 7 i de l'1 % a pH 4,5. Les pastes obtingudes per escalfament de la FC presentaven característiques molt diferenciades depenent del pH. A pH neutre es formaven unes pastes dures i consistents, mentre que a pH àcid les pastes eren poc consistents, molt adhesives i més elàstiques que les anteriors. Aquestes tenien una capacitat de retenció d'aigua molt superior que les de pH 7, en les quals l'aigua quedava retinguda per capil·laritat. La textura i capacitat de retenció d'aigua de les pastes tampoc eren afectades pel tractament HHP. El tractament HHP incrementava l'activitat de la Tripsina sobre l'Hb quan el substrat i l'enzim es tractaven conjuntament i afavoria el procés d'obtenció d'hidrolitzats descolorats a partir de la FC, la qual cosa permetia assolir el mateix grau de descoloració amb una dosi d'enzim inferior. El tractament d'hidròlisi de la FC amb la utilització combinada de Tripsina seguida d'un tractament amb Pepsina permetia l'obtenció d'un hidrolitzat proteic d'Hb descolorat i hidrolitzava completament la globina, donant lloc a 2 pèptids de 10,8 i 7,4 KDa. Val a dir que també produïa un 60 80 % de nitrogen soluble en TCA, constituït fonamentalment per pèptids petits i aminoàcids lliures. Els hidrolitzats trípsics i pèpsics d'Hb, obtinguts a partir de FC no pressuritzada i deshidratats per atomització a 180ºC, eren de color blanc i tenien un contingut en humitat del 4,7 %, un 84,2 % de proteïna i 9,7 % de sals minerals. El procés d'hidròlisi permetia una reducció considerable de la contaminació de la FC, obtenint un producte en pols amb uns recomptes totals de l'ordre de 102-103 ufc·g-1. Pel que fa a la funcionalitat dels hidrolitzats d'Hb deshidratats per atomització, aquests presentaven una elevada solubilitat proteica a pH 5 i 7 i romanien solubles després d'un escalfament a 80ºC durant 30 min. Tanmateix, aquesta hidròlisi afectava molt negativament la capacitat de mantenir escumes estables i l'activitat emulsionant.

UV index measurement and model agreement: uncertainties and limitations

En les últimes dècades, l'increment dels nivells de radiació solar ultraviolada (UVR) que arriba a la Terra (principalment degut a la disminució d'ozó estratosfèric) juntament amb l'augment detectat en malalties relacionades amb l'exposició a la UVR, ha portat a un gran volum d'investigacions sobre la radiació solar en aquesta banda i els seus efectes en els humans. L'índex ultraviolat (UVI), que ha estat adoptat internacionalment, va ser definit amb el propòsit d'informar al públic general sobre els riscos d'exposar el cos nu a la UVR i per tal d'enviar missatges preventius. L'UVI es va definir inicialment com el valor màxim diari. No obstant, el seu ús actual s'ha ampliat i té sentit referir-se a un valor instantani o a una evolució diària del valor d'UVI mesurat, modelitzat o predit. El valor concret d'UVI està afectat per la geometria Sol-Terra, els núvols, l'ozó, els aerosols, l'altitud i l'albedo superficial. Les mesures d'UVI d'alta qualitat són essencials com a referència i per estudiar tendències a llarg termini; es necessiten també tècniques acurades de modelització per tal d'entendre els factors que afecten la UVR, per predir l'UVI i com a control de qualitat de les mesures. És d'esperar que les mesures més acurades d'UVI s'obtinguin amb espectroradiòmetres. No obstant, com que els costs d'aquests dispositius són elevats, és més habitual trobar dades d'UVI de radiòmetres eritemàtics (de fet, la majoria de les xarxes d'UVI estan equipades amb aquest tipus de sensors). Els millors resultats en modelització s'obtenen amb models de transferència radiativa de dispersió múltiple quan es coneix bé la informació d'entrada. No obstant, habitualment no es coneix informació d'entrada, com per exemple les propietats òptiques dels aerosols, la qual cosa pot portar a importants incerteses en la modelització. Sovint, s'utilitzen models més simples per aplicacions com ara la predicció d'UVI o l'elaboració de mapes d'UVI, ja que aquests són més ràpids i requereixen menys paràmetres d'entrada. Tenint en compte aquest marc de treball, l'objectiu general d'aquest estudi és analitzar l'acord al qual es pot arribar entre la mesura i la modelització d'UVI per condicions de cel sense núvols. D'aquesta manera, en aquest estudi es presenten comparacions model-mesura per diferents tècniques de modelització, diferents opcions d'entrada i per mesures d'UVI tant de radiòmetres eritemàtics com d'espectroradiòmeters. Com a conclusió general, es pot afirmar que la comparació model-mesura és molt útil per detectar limitacions i estimar incerteses tant en les modelitzacions com en les mesures. Pel que fa a la modelització, les principals limitacions que s'han trobat és la falta de coneixement de la informació d'aerosols considerada com a entrada dels models. També, s'han trobat importants diferències entre l'ozó mesurat des de satèl·lit i des de la superfície terrestre, la qual cosa pot portar a diferències importants en l'UVI modelitzat. PTUV, una nova i simple parametrització pel càlcul ràpid d'UVI per condicions de cel serens, ha estat desenvolupada en base a càlculs de transferència radiativa. La parametrització mostra una bona execució tant respecte el model base com en comparació amb diverses mesures d'UVI. PTUV ha demostrat la seva utilitat per aplicacions particulars com ara l'estudi de l'evolució anual de l'UVI per un cert lloc (Girona) i la composició de mapes d'alta resolució de valors d'UVI típics per un territori concret (Catalunya). En relació a les mesures, es constata que és molt important saber la resposta espectral dels radiòmetres eritemàtics per tal d'evitar grans incerteses a la mesura d'UVI. Aquest instruments, si estan ben caracteritzats, mostren una bona comparació amb els espectroradiòmetres d'alta qualitat en la mesura d'UVI. Les qüestions més importants respecte les mesures són la calibració i estabilitat a llarg termini. També, s'ha observat un efecte de temperatura en el PTFE, un material utilitzat en els difusors en alguns instruments, cosa que potencialment podria tenir implicacions importants en el camp experimental. Finalment, i pel que fa a les comparacions model-mesura, el millor acord s'ha trobat quan es consideren mesures d'UVI d'espectroradiòmetres d'alta qualitat i s'usen models de transferència radiativa que consideren les millors dades disponibles pel que fa als paràmetres òptics d'ozó i aerosols i els seus canvis en el temps. D'aquesta manera, l'acord pot ser tan alt dins un 0.1º% en UVI, i típicament entre menys d'un 3%. Aquest acord es veu altament deteriorat si s'ignora la informació d'aerosols i depèn de manera important del valor d'albedo de dispersió simple dels aerosols. Altres dades d'entrada del model, com ara l'albedo superficial i els perfils d'ozó i temperatura introdueixen una incertesa menor en els resultats de modelització.